Abstract
In nature, biofilms are the most common form of bacterial growth. In biofilms, bacteria display coordinated behaviour to perform specific functions. Here, we investigated electrical signalling as a possible driver in biofilm sociobiology. Using a multi-electrode array system that enables high spatio-temporal resolution, we studied the electrical activity in two biofilm-forming strains and one non-biofilm-forming strain. The action potential rates monitored during biofilm-forming bacterial growth exhibited a one-peak maximum with a long tail, corresponding to the highest biofilm development. This peak was not observed for the non-biofilm-forming strain, demonstrating that the intensity of the electrical activity was not linearly related to the bacterial density, but was instead correlated with biofilm formation. Results obtained indicate that the analysis of the spatio-temporal electrical activity of bacteria during biofilm formation can open a new frontier in the study of the emergence of collective microbial behaviour.
1. Introduction
Bacteria exist in communities and often coordinate their behaviour to perform specific functions. Biofilm formation is an example of this behaviour. Bacteria undergo large changes involving genetic pathways, such as transcriptional reprogramming [1–3] and modulation by small non-coding RNAs [4], and phenotypic modifications during their transition from a planktonic state (free-swimming) to biofilm development.
The classical phases of a biofilm cycle consist of: (i) initial attachment to a surface; (ii) production of a matrix of extracellular polymeric substances, resulting in tighter adhesion; (iii) formation of micro-colonies; (iv) maturation of micro-colonies to mature biofilms; and (v) detachment and return to planktonic growth [5].
Genetic and molecular techniques in association with microscopic imaging have revealed that biofilm development is a well-regulated process that leads to the formation of a complex community in which bacteria integrate internal and external signals, perceive the density and type of neighbours and coordinate a number of multicellular behaviours that may be associated with morphological, physiological features and other changes in cell properties [6]. The complexity of biofilm development suggests that its formation is not a stochastic process, but rather a type of microbial development in which changes in form and function play a prominent role in the life cycle of the bacteria [6].
The involvement of cell–cell signalling known as quorum sensing and the second messenger signal cyclic di-GMP are well documented in bacterial biofilm formation [7–9]. These widely conserved and important strategies, which allow bacteria to sense changes in population density and local environmental conditions and to accomplish the transition from a motile to a sedentary state, are thought to be integrated [10].
Despite recent impressive progress, some key questions about the signals involved in biofilm development remain unanswered [9]. Electrophysiological measurements, which were recently proposed for bacterial membrane potential analysis, could be a useful tool to investigate cell–cell signalling in biofilm development. The pioneering study of Kralj et al. [11], who used voltage-sensitive fluorescent membrane proteins to obtain electrophysiological measurements of individual intact bacteria, demonstrated that bacteria generate action potentials (APs) in response to different stimuli. Classical patch-clamp recording methods are not applicable to bacteria owing to the size of the cells (corresponding to the size of a typical patch pipette) and the presence of a cell wall. However, Martinac et al. [12] established a procedure to create giant bacterial spheroplasts of Escherichia coli cells, thereby enabling the electrophysiological investigation of bacterial cell membranes and the study of ion channels in prokaryotes. The multi-electrode array (MEA) system is a well-characterized tool for investigating electrical network activity in excitable cells undergoing maturation, particularly neurons. Because random electrical activity becomes synchronized with an increasing degree of correlation as the network develops [13–15], recording the electrical activity dynamics of developing biological networks enables the characterization of their spatio-temporal organization [16–19].
In this study, the MEA system was used for the first time, to the best of the authors' knowledge, to study spatio-temporal electrical activity during the development of bacterial biofilms (Bacillus licheniformis 51.1 and Pseudomonas alcaliphila 34) and during planktonic growth (E. coli HEC30) to gain deeper insights into the complex mechanism of signalling involved in biofilm formation. The observed patterns of electrical spiking during biofilm development provide new perspectives on the investigation of cell–cell communication mechanisms in bacteria.
2. Material and methods
2.1. Bacterial strains, media and growth conditions
The bacterial strains used were B. licheniformis 51.1 [20] and P. alcaliphila 34 [21,22], which were selected because of their ability to form biofilms under anaerobic and aerobic conditions, respectively, and E. coli HEC30, which was selected as a non-biofilm-forming model. The strains were grown in lysogeny broth (LB) [23] at 30°C with shaking.
For biofilm development, overnight bacterial cultures of the strains were diluted in LB to an optical density at 590 nm (OD590) of 0.7. Then, 200 µl of the cell suspension was added to each well of a 96-well plate. The development of B. licheniformis 51.1 biofilms under anaerobic conditions was stimulated by adding 1% KNO3 to the LB medium as a terminal acceptor of metabolic electrons and covering each well of the microplates with drops of sterile vaseline oil.
The plates were incubated at 30°C without shaking for at least 72 h. At each sampling time, the planktonic cells of a plate were removed by multi-channel pipette aspiration, each well was washed with 250 µl of sterile water and the plate was vigorously shaken to remove all non-adherent bacteria. The attached biofilm was stained by the addition of 1 ml of a 0.1% crystal violet (CV) solution. After incubation for 15 min, the unbound CV was removed by aspiration, and excess stain was removed by placing the plate under running tap water. After the plate was air dried, 1 ml of 96% ethanol was added to each well, and the OD590 was measured after 30 min using a DVB9990 microplate reader (GDV, Rome, Italy).
Planktonic cultures of E. coli were grown in microplates containing LB medium and were monitored at OD590 (GDV, Italy).
2.2. Multi-electrode array recording and signal processing
The strains were grown overnight in LB medium and diluted to an OD590 of 0.7 or 0.4 in the appropriate LB medium for biofilm and planktonic culture, respectively. Then, 1.5 ml of each cell suspension was added to a sterilized planar MEA (Multichannel Systems MCS GmbH, Reutlingen, Germany) with 59 electrodes plus 1 internal reference (figure 1) and incubated at 30°C in the recording incubator. In order to promote biofilm formation in B. licheniformis 51.1, which is positively influenced by anaerobic conditions, 750 µl of sterile vaseline oil was applied on top of the MEA well, using a micropipette with sterile tips.
The MEA chip was placed on the MEA board (USB-MEA60, Multichannel Systems MCS GmbH) for long-term recording of bacterial spike activity during growth and development. Throughout, sterility was ensured using a special cover (figure 1), and the temperature was maintained using a temperature controller (TC02, Multichannel Systems MCS GmbH). Datasets from three separate replicates were obtained under the same conditions. Control experiments (without cells) were also performed in which the signals from MEA chips filled with LB medium were recorded. Pharmacological perturbation was performed using sodium azide (SA; 10 mM [24]) and carbonyl cyanide m-chlorophenylhydrazone (CCCP; 50 mg ml−1 [25]) added after 6 h from the inoculum. Each treatment was replicated three times for both biofilm-forming bacteria strains.
The extracellular activity of cultured bacteria was recorded using a USB-MEA60 system. Signals were obtained simultaneously from 59 recording electrodes (30 µm diameter with 500 µm spacing) with a sampling rate of 5 kHz. The remaining electrode was used as an internal ground. Signals were continuously recorded for at least 72 h with no perturbation. MC_RACK software (Multichannel Systems MCS GmbH) was used for data acquisition and spike detection. The values of the signal-to-noise ratio (S/N) and of the root mean-squared (RMS) of the basal noise of signal were calculated over some brief time windows (500 ms) extracted at different points along each recording session. The S/N ranged between 6 and 12, with an average value of 7.8 ± 0.7 (mean ± s.d.), and the averaged RMS of the basal noise was 4.04 ± 1.85 µV (mean ± s.d.). Both values were constant during a single recording session, allowing spike detection with a threshold-based method; because of small manufacturing differences among the several MEA arrays used for the experiments, certain variability was, indeed, found in different recording sessions. In any case, such variability was never correlated with different bacterial strains. The optimal threshold, set three times higher than the RMS of the basal noise, was chosen performing more than two attempts with different threshold values for each recording. Moreover, no filter was applied to raw data. Electrical traces were analysed in order to exclude the pure stochasticity of data; for that purpose, the probability distribution of data points was plotted and compared with the Gaussian distribution. The dedicated software NeuroExplorer (Nex Technologies, Littleton, MA) was used for the timestamp extraction, spike rate analysis and spike shape analysis. Spatial and temporal analyses of the electrical activity were performed using routines written in Matlab (The MathWorks, Inc., Natick, MA). The data are reported as the mean ± standard error (s.e.m.). Comparisons of the means were performed with one-way ANOVA and Tukey's post-test using GraphPad Prism v. 5 (GraphPad Software, San Diego, CA).
2.3. Analysis of synchronized events
We described the synchronized events by counting the electrode spikes evoked together during a fixed time interval of various selected values (Δt = 100–700 ms). The values of Δt were chosen to be much smaller than the mean interspike intervals. We labelled such events as temporal clusters (figure 2a). Then, we counted the number nt of such clusters of size C and divided it by the total number of spikes Ntot registered in a given experiment, obtaining the quantity
The next step was the identification of the spatio-temporal clusters, defined as the sites that fire together during a time interval Δt, with additional restriction for the spatial separation between sites (figure 2b).
We defined a distance d as the horizontal, vertical or diagonal separation between the electrodes, calculated using the Pythagorean formula for distance
We then calculated all distances between two given sites for each temporal cluster C ≥ 2. The total number of distances for a given C and instant of time was calculated as follows:
3. Results
Using 60-channel MEAs, we studied the development of electrical activity of biofilms (B. licheniformis 51.1 and P. alcaliphila 34) and planktonic cultures (E. coli HEC30; raw data are available upon request). The probability distribution of the data points from the raw electrical traces of the experiment showed two peaks, one corresponding to the noisy background, and the second one being a contribution from the deterministic spikes; this provided an evidence of non-Gaussianity of the data (electronic supplementary material, figure S1). Representative spikes from the three species of bacteria are illustrated in figure 3. Replicate trials revealed that the spikes had a similar shape, with an approximate duration of 1 ms. It is worth noting that in extracellular MEA recordings the time properties of signals detected are represented as their first derivative [26,27], thus they do not describe the actual membrane potential of the cell under study. This is why the reported spike duration is smaller than those already observed with other and not comparable experimental protocols [11]. The biofilm-forming bacterial strains B. licheniformis 51.1 and P. alcaliphila 34 both produced a biofilm layer covering the electrodes, as illustrated in figure 1.
The spike rate during the growth of each strain is reported in table 1. The spike rate was calculated during three different growth periods (TZ1, TZ2 and TZ3, 1 h after the beginning of the recording, at the time of maximum rate and after 72 h, respectively). The spike rate differed in biofilm- and non-biofilm-forming bacteria, particularly during TZ2, which corresponded to the time of the maximum spike rate. As expected, the spike rate in control trials with LB medium alone was close to zero (data not shown).
B. licheniformis | P. alcaliphila | E. coli | |
---|---|---|---|
AP rate (spikes s−1) | |||
average | 1.84b* | 0.98ab | 0.54a |
TZ1 | 0.11a | 0.16b | 0.08a |
TZ2 | 2.89b | 1.24b | 0.32a |
TZ3 | 1.88b | 0.88ab | 0.40a |
AP amplitude (μV) | |||
average | 19.8b | 18.9b | 13.6a |
TZ1 | 18.9b | 18.88b | 12.98a |
TZ2 | 20.55b | 19.15b | 13.95a |
TZ3 | 19.95b | 18.94b | 13.90a |
During the course of the experiments, a change in the firing rate was observed. The average number of spiking electrodes increased significantly in B. licheniformis and P. alcaliphila, whereas no significant variation was observed in E. coli at 24 h after inoculation (figure 4a). Furthermore, spike rate quantification revealed different trends for the different species. In biofilm-forming bacteria, three time zones were apparent: an initial exponential increase in spike activity, followed by a stationary phase and a phase with decreased electrical activity. The maximum peak occurred between 36 and 48 h from the beginning of the recording for B. licheniformis and between 12 and 24 h for P. alcaliphila. By contrast, the spiking activity of E. coli increased slightly at the beginning of growth, and then was stationary until the end of the observation period (figure 4b).
To determine whether an increase in the number of bacterial cells influenced the trends observed in the analysis of spike rate, we monitored biofilm and planktonic culture growth. As expected, the E. coli strain did not produce a biofilm (figure 5a), whereas the planktonic culture exhibited a typical growth curve (figure 5b). The biofilm-forming bacteria, particularly B. licheniformis, produced an abundant biofilm. Moreover, the time of the maximum biofilm formation corresponded to the maximum electrical activity recorded with the MEA system. The relationship between biofilm development and electrical activity measured by linear regression was R2 = 0.9753 (p < 0.0001) and R2 = 0.7533 (p < 0.0005) for B. licheniformis and P. alcaliphila, respectively.
The effect of metabolic uncouplers SA and CCCP on the electrical activity generated by the two biofilm-forming strains was also checked. The inhibition of the adenosintriphosphatase (ATPase) by SA induced a decrease in the spike rate quite immediately after the treatment (electronic supplementary material, table S1); such effect attenuated after a few hours. The dissipation of the proton motive force by CCCP had a more drastic effect on the electrical activity detected by the MEA system (electronic supplementary material, table S1), the rate of the electrical activity remained decreased for the rest of the recording. Two representative plots of the spike rate versus time are reported in the electronic supplementary material, figure S2. The probability distribution of the data points from the raw electrical traces after the treatments was checked and compared with that of the raw traces recorded in the absence of any kind of cells (electronic supplementary material, figure S3). No difference was found and, in both cases, data showed pure Gaussian distribution.
As already mentioned, the spikes for all strains were quite similar in terms of shape and duration (figure 3) but not amplitude (table 1). On average, we observed higher spike amplitudes in B. licheniformis and P. alcaliphila than in E. coli. Moreover, the trend in the spike amplitude of the APs measured during the course of the experiment differed during B. licheniformis and P. alcaliphila biofilm development and E. coli planktonic growth (figure 6). During the first 12 h, the spike amplitude increased for both biofilm-forming bacteria; it then continued to increase in B. licheniformis, but reached a stationary level in P. alcaliphila. The spike amplitude observed in E. coli exhibited very small variations during the duration of the experiment.
During MEA recording, the amplitude of the APs detected depends on the distance between the electrodes and the site of generation of the signal: the closer the source, the stronger and clearer the signal. Furthermore, spike amplitude is affected by the simultaneous firing of one or more cells. As a consequence, it is difficult to fully reconstruct the exact amplitude of a single spike. Therefore, to better characterize the electrical signal features produced by bacteria, we used another approach to analyse the synchronization phenomena that considers the spatio-temporal correlations of the measured spikes and not their amplitudes. We performed an analysis of the synchronized events by qualitatively comparing the electrical activity and the occurrence of synchronized events in time and space for the three different bacterial species and for random data (created in silico for each experiment and used to compare the obtained results with the experimental ones), while ignoring the spike amplitudes, to determine whether the experimental data differed from random data. We observed the appearance of temporal and spatio-temporal clusters in all cultures. Figures 7 and 8 show the distributions of temporal (pt) and spatio-temporal clusters (ps) versus cluster size C for Δt = 300 ms in bacteria (figures 7a and 8a) and random data (figures 7b and 8b). The presentation of the plot obtained using 300 ms as the time threshold is arbitrary; no different trends in the results were observed for different Δt (data not shown). The time threshold influences the absolute quantity of temporal and spatio-temporal clusters, not the relative differences between strains. In the random data, there was a critical value of cluster size Cm above which the clustering events were not observed; thus, cluster size can be considered a threshold that separates random data from data with deterministic characteristics. For bacteria, the clustering was also observed for sizes larger than Cm, thus demonstrating its purely deterministic origin. These results indicate that for cluster sizes larger than a determined value (Cm), the electrical spiking generated by bacteria shows organized synchronizations that are absent in random data. Moreover, by quantifying the temporal and spatio-temporal clusters for C > Cm, we observed that biofilm-forming bacteria had a higher average number of temporal and spatio-temporal clusters than E. coli. P. alcaliphila had larger values for pt and ps throughout all cluster sizes; B. licheniformis had a greater number of spatio-temporal clusterings than temporal ones, although temporal clusterings were observed; and E. coli showed no clustering except for a cluster size equal to 30, which may be an artefact (figures 7a and 8a). By contrast, in the case of small cluster sizes, i.e. C < Cm, the synchronized events originated both from random coincidence and deterministic sources. Nevertheless, some interesting information was evident in the pt and ps values for small cluster sizes calculated using different Δt values (electronic supplementary material, figure S4). As expected, with increasing values of Δt, which marks the temporal period for the spikes contributing to a given cluster, the values of pt and ps also increased. However, in the same range of Δt, differences between biofilm-forming and non-biofilm-forming bacteria were evident, as B. licheniformis and P. alcaliphila had more numerous temporal and spatio-temporal clusterings than E. coli (figure 9).
Finally, we performed a qualitative analysis of the selected cluster size (C = 2) at different Δt for spikes synchronized temporally (electronic supplementary material, figure S5) and spatio-temporally (electronic supplementary material, figure S6) recorded during the three different temporal zones: TZ1, TZ2 and TZ3. During each period, we observed higher numbers of temporal and spatio-temporal clusters in biofilm-forming bacteria compared with non-biofilm-forming bacteria. Such differences were observed in TZ2 and TZ3, which correspond to the periods of maximum spike rate and maximum biofilm quantification, respectively (figures 10 and 11). This result demonstrates that the development of a biofilm favours the propagation of electrical pulses.
4. Discussion
As a result of biofilm formation, cells are immobilized and kept in close proximity, forming synergistic microconsortia in which intense interactions are observed [28]. The dynamics of connectivity and functionality during bacterial biofilm development are not completely understood. In this paper, the APs produced during the development of biofilms were recoded using an MEA system equipped with 59 microelectrodes. In an MEA recording, each single signal can be considered the sum of the signals generated by a large number of cells owing to the relatively large dimensions of the electrode compared with those of a single bacterium. In fact, up to hundreds of microbial cells may be in contact with just one electrode or detected by it (a 30 μm electrode detects an area of more than 100 µm in diameter [29]). The increase in the density of the culture and the natural sinking of cells on the electrode surface may be responsible for the initial increase in the spike amplitudes. In other words, the closer position of bacterial cells to the electrodes and the increase in the number of cells, typical of biofilm development, could be responsible for the detection of higher-quality signals. For example, it is possible that the stronger signals were the result of the simultaneous firing of a group of adjacent cells. The data revealed that, similar to the events that occur during the formation of neuron networks [13–15], the number of spatio-temporal synchronized events among bacterial cells increased as the biofilm developed, allowing the bacteria to interact. The stronger electrical activity observed in biofilm-forming bacteria suggests its important role in the emergence of collective behaviours.
As single cells a few micrometres in size, bacteria can perceive stimuli from a limited area in a limited time, but as a member of a complex super-organism (such as a colony or a biofilm), each bacterium must be able to sense and communicate with other cells to cooperate, first for the self-organization of the colony/biofilm, then for distribution tasks and, eventually, facing environmental challenges. Ben-Jacob described multicellular bacterial organizations as a ‘massive brain’ capable of information processing, collective memory and self-plasticity [30]. In this perspective, the experiments described here represent a further step for approaching bacterial communities as ‘intelligent’ super-organisms that make use of electric signals to communicate.
Acknowledgements
M.C. acknowledges Regione Toscana for financial support. E. Masi, S.M. and L.S. were supported by the Future and Emerging Technologies (FET) programme within the Seventh Framework Programme for