Transcriptional profiling of Klebsiella pneumoniae defines signatures for planktonic, sessile and biofilm-dispersed cells
BMC Genomics volume 17, Article number: 237 (2016)
Surface-associated communities of bacteria, known as biofilms, play a critical role in the persistence and dissemination of bacteria in various environments. Biofilm development is a sequential dynamic process from an initial bacterial adhesion to a three-dimensional structure formation, and a subsequent bacterial dispersion. Transitions between these different modes of growth are governed by complex and partially known molecular pathways.
Using RNA-seq technology, our work provided an exhaustive overview of the transcriptomic behavior of the opportunistic pathogen Klebsiella pneumoniae derived from free-living, biofilm and biofilm-dispersed states. For each of these conditions, the combined use of Z-scores and principal component analysis provided a clear illustration of distinct expression profiles. In particular, biofilm-dispersed cells appeared as a unique stage in the bacteria lifecycle, different from both planktonic and sessile states. The K-means cluster analysis showed clusters of Coding DNA Sequences (CDS) and non-coding RNA (ncRNA) genes differentially transcribed between conditions. Most of them included dominant functional classes, emphasizing the transcriptional changes occurring in the course of K. pneumoniae lifestyle transitions. Furthermore, analysis of the whole transcriptome allowed the selection of an overall of 40 transcriptional signature genes for the five bacterial physiological states.
This transcriptional study provides additional clues to understand the key molecular mechanisms involved in the transition between biofilm and the free-living lifestyles, which represents an important challenge to control both beneficial and harmful biofilm. Moreover, this exhaustive study identified physiological state specific transcriptomic reference dataset useful for the research community.
Most bacteria can live in individual or community lifestyles. In the planktonic mode of growth, bacterial cells are free to move in suspension, whereas in the sessile state, they form surface-attached multicellular communities called biofilms. This dynamic heterogenic organization confers to its residents a powerful tolerance against stresses and facilitates symbiotic relationships between members of the communities [1, 2]. The transition between the planktonic and sessile modes of growth, as well as the biofilm development process are governed by environmental cues and the coordination of various molecular pathways linked notably to secondary messenger cyclic di-GMP and quorum sensing [3, 4]. Biofilm development progresses in three stages: i) bacterial attachment to a surface and formation of a monolayer biofilm, ii) maturation of the biofilm and emergence of a three-dimensional structure and iii) dispersion from mature biofilm.
The adhesion of planktonic cells to the surface is mostly driven by surface-exposed components like flagella, fimbriae and curli as observed in many bacteria . Subsequent biofilm maturation is concomitant with the formation of an extracellular matrix composed of exopolysaccharides, DNA, lipids and proteins . In Pseudomonas aeruginosa and Escherichia coli, exopolysaccharides and extracellular DNA also play a crucial role in the maturation process as the absence of these compounds severely impairs the formation of a three-dimensional structure .
The last step of the biofilm developmental process, dispersion from mature biofilm, constitutes an essential stage because of its crucial role in bacterial dissemination and colonization of new surfaces [8, 9]. It remains therefore unclear whether bacteria dispersed from biofilms represent or not a transition stage between biofilm and the planktonic lifestyle. Dispersion occurs either as individual cells or clumps , but the molecular mechanisms and effectors behind this process are still poorly documented . Nevertheless, secreted effectors such as glycosidases in Actinobacillus actinomycetemcomitans , proteases in Pseudomonas putida , nucleases in Haemophilus influenzae  and biosurfactants in Staphylococcus [15, 16] are able to destabilize the biofilm structure and promote dispersion. Activation of prophages in P. aeruginosa and Enterococcus faecalis was also reported as inducing cell death inside microcolonies leading to biofilm dispersion [17, 18].
Despite the accumulation of data concerning the transcriptional profile of bacteria grown in different experimental models, there has been no documented overview of all states of biofilm development and dispersion. Transcriptomic approaches by microarray or RNA sequencing have attempted to address this issue in several bacterial species like E. coli, P. aeruginosa or Acinetobacter baumannii, and showed distinct expression profiles between sessile and planktonic stages. However, cells from dispersed biofilm were not included in these analyses [19–21].
The aim of this study was to identify the transcriptional landscape of the bacteria Klebsiella pneumoniae across different experimental growth states, i.e. planktonic, sessile, and spontaneously biofilm-detached bacteria. K. pneumoniae is an ubiquitous bacterium found both in nature and in clinical environments; the molecular mechanisms leading to biofilm formation have been previously investigated, mostly by punctual mutant analysis [22, 23]. In this work, comparison of the different whole transcriptomes obtained by RNA-seq showed that each lifestyle of K. pneumoniae was associated with a unique transcriptional behavior. The comprehensive overview provided by this study allowed the identification of specific transcriptional fingerprints for each state, including the biofilm-dispersed cells.
Monitoring of biofilm development in a flow-cell model
Monitoring of biofilm development by K. pneumoniae CH1034 in a flow-cell system with confocal microscopy showed initially the formation of microcolonies leading to the development of a flat structure after 7 h of incubation (T7h) (Additional file 1: Movie S1). At T9h, a three-dimensional structure was observable and potential detachment from this mature biofilm was then assessed. Bacteria in the flow-cell effluent were harvested throughout the experiment, and CFU determination of the resulting suspensions indicated that the number of viable cells decreased in the first 3 h of the experiment, from 5.106 CFU/mL (T1h) to 1.105 (T3h), owing probably to the elimination of planktonic non-adhering cells (Fig. 1a). Observation of the harvested samples by optical microscopy revealed mainly individual bacteria (data not shown). From T3h to T6h, the number of viable bacteria in the effluent increased rapidly and then progressively in the following 10 h (T6h to T16h) (Fig. 1a). Microscopic observations revealed a progressive appearance of bacterial aggregates in the effluent, which predominated over individual cells after 12 h of incubation (Fig. 1b).
Planktonic, sessile, and biofilm-detached bacteria presented distinct transcriptional profiles
Transcriptional analysis was performed with sessile bacterial cells collected before and after the formation of a three-dimensional structure, at T7h and T13h, respectively. Detached cells isolated in the flow-cell effluent (T12h-T13h), exponential and stationary growing planktonic cells were also included. RNAseq analysis indicated that 2 052 of the 5 146 CDS of K. pneumoniae, as well as 19 of the 44 annotated ncRNA genes (excluding tRNA and rRNA genes), were differentially expressed in at least one of the ten possible pairs of conditions (∣fold-change∣ > 5 and adjusted P-value < 0.01) (Fig. 2a), with fold-changes ranging from −2 780 to 2 182 (Additional file 2: Table S1; Additional file 3: Table S2). To validate the RNA-seq efficiency, 20 genes differentially expressed between the 13 h-old biofilm bacteria and the bacteria collected in the effluent (10 genes overexpressed and 10 genes under-expressed; P-value < 0.01) were randomly selected. Their relative expression levels were determined by RT-qPCR with total RNA extracted from cells harvested in two conditions: bacteria in the effluent and 13 h-old biofilm. Results indicated a high correlation between RNAseq and RT-qPCR data (r = 0.97; P-value < 0.0001; Pearson’s correlation test) (Additional file 4: Figure S1).
PCA performed with Z-score values of the 2 052 CDS and 19 ncRNA genes indicated that the first principal component (PC1) accounted for 36.52 % and the second principal component (PC2) for 27.88 % of the total variation in the dataset (Fig. 2b). A plot of these Z-score values against a heatmap (Additional file 5: Figure S2) and the proximity of points in the PCA (Fig. 2b) demonstrated the high reproducibility of the data among the replicates. In addition, such analysis clearly indicated that all bacterial states (planktonic, sessile and bacteria in the effluent) exhibited specific transcriptional profiles (Fig. 2b and Additional file 5: Figure S2), and suggests that bacterial cells in the effluent are not pieces of biofilm mechanically detached from the biomass. Hereafter they will be referred to as biofilm-dispersed cells.
The transcriptome of the biofilm-dispersed cells presented only 224 CDS and 3 ncRNA genes differentially expressed (∣fold-change∣ > 5 and adjusted P-value < 0.01) when compared with those of the 7 h-old biofilm state. In contrast, 454 CDS and 7 ncRNA genes, 486 CDS and 2 ncRNA genes, and 1 080 CDS and 6 ncRNA genes were differentially expressed (∣fold-change∣ > 5 and adjusted P-value < 0.01) when compared with those of exponential planktonic state, 13 h-old biofilm and stationary planktonic state, respectively (Fig. 2a). Hence, biofilm-dispersed cells harbored a distinct transcriptional profile, which was closer to that of bacteria from 7 h-old biofilm than to that of 13 h-old biofilm and planktonic cells.
Gene functional classification of K. pneumoniae lifestyles through K-means clustering
K-means clustering was then used to visualize the distribution of the expression levels of the 2 052 CDS and the 19 ncRNA genes differentially expressed (∣fold-change∣ > 5 and adjusted P-value < 0.01) in the different conditions (Fig. 3a and b). Owing to the high reproducibility of data, Z-score values were able to be calculated with average values from normalized DEseq counts. This clustering indicated that the clearest representation was obtained with K = 10 for the CDS analysis and K = 5 for the ncRNA genes analysis, and showed different transcriptomic profiles between conditions. In Fig. 3a, with clusters ranging from 76 to 499 CDS for clusters 8 and 10 respectively, column clustering confirmed that dispersed cells were transcriptionally closer to 7 h-old biofilm cells than to those in all the other conditions, whereas stationary phase cells were the most different group of this dataset.
In order to highlight groups of genes highly overexpressed or under-expressed in a specific condition, the mean of the Z-scores in each cluster in the Fig. 3a was calculated for each condition. Only the Z-score groups presenting a mean value > 1 or < −1, named overexpressed boxes and under-expressed boxes, respectively (framed in Fig. 3a), were considered thereafter. All clusters presented only one overexpressed box, but clusters 5, 8 and 9 also presented one under-expressed box (Fig. 3a).
Analysis of the potential function of protein-coding genes in the under-expressed and overexpressed boxes by the Clusters of Orthologous Groups (COG) classification is represented in Fig. 3c and Additional file 6: Figure S3. A large number of genes were poorly characterized and therefore categorized in the “unknown function” class. Exponential planktonic cells exhibited two overexpressed boxes (clusters 1 and 6) (Fig. 3a), containing CDS mainly involved in inorganic ion transport and metabolism (14.9 and 15.1 % of the genes present in clusters 1 and 6, respectively) (Fig. 3c and Table 1). In parallel, two under-expressed boxes (clusters 5 and 8) were identified in the exponential planktonic condition. They contained mainly CDS involved in amino acid transport and metabolism, and energy production and conversion, as defined by the COG classification. Stationary planktonic cells exhibited three overexpressed boxes (clusters 7, 8 and 10) that contained CDS mostly implied in energy production and conversion, and in amino acid and carbohydrate transport and metabolism. The 7 h-old biofilm cells exhibited two overexpressed boxes (clusters 5 and 9) (Fig. 3a), which contained CDS chiefly involved in amino acid transport and metabolism (21.7 and 24 % of the genes present in clusters 5 and 9, respectively) (Fig. 3c and Table 1). The 13 h-old biofilm cells exhibited one overexpressed box (cluster 3), with CDS chiefly involved in carbohydrate transport and metabolism (21 % of the genes present in cluster 3). Finally, dispersed cells exhibited two overexpressed boxes (clusters 2 and 4), containing CDS chiefly involved in translation, ribosomal structure and biogenesis (21.9 and 9.3 % of the genes present in clusters 2 and 4, respectively).
Identification of a set of signature genes for each condition
Since clustering suggested the existence of specific signature genes for each condition, different stringent threshold fold-changes were applied to extract the most relevant transcriptional signature genes, up- or down-regulated, for each condition (Additional file 7: Figure S4). Forty signature CDS were identified, 11 associated with the exponential and the stationary planktonic states, 4 with the 7 h-old and the 13 h-old biofilm cells, and 10 with biofilm dispersal (Table 2). In the stationary planktonic and 13 h-old biofilm conditions, all signature CDS were upregulated, and in the 7 h-old biofilm condition, all were down-regulated, whereas exponential planktonic cells and biofilm-dispersed cells displayed both up- and down-regulated signature CDS (Table 2 and Fig. 4). The Z-score values of these 40 CDS plotted against a heatmap (Fig. 4a) and their relative expression level (Fig. 4b) confirmed their signature singularity. Putative functions of these protein encoding signatures CDS are listed in Table 2 and concern mainly transport, transcriptional regulation and metabolic pathways.
In the present study, the transcriptional changes occurring in the course of K. pneumoniae biofilm formation and biofilm-detachment were characterized by RNAseq. To date, the few data available on biofilm dispersion were obtained with artificial dispersion signals such as c-di-GMP depletion [24, 25]. In contrast, we investigated spontaneous biofilm-detached cells. Results indicated that each of the tested K. pneumoniae lifestyles, i.e. planktonic (exponential and stationary phases), sessile (7 h-old and 13 h-old biofilms) and biofilm-dispersed cells, exhibit unique and specific transcriptional profiles. The comprehensive overview presented in this study allowed the analysis of the transcriptional fate of all K. pneumoniae genes in different bacteria lifestyles.
The stationary planktonic mode of growth displayed the most particular pattern with 499 genes highly overexpressed in the K-means cluster 10. Entry in the stationary phase is the result of nutrient starvation and in consequence bacteria modulate the expression level of a considerable number of genes, many of them being under the control of the stationary-phase sigma S factor (σS) . On the basis of a study referencing the 100 most RpoS-dependent genes in stationary phase of a pathogenic E. coli strain , 54 of the 82 genes present in the K. pneumoniae genome were found in the K-means cluster 10, including 4 transcriptional signature genes of the stationary phase (ygaT (also named csiD), astA, astD and astE). Overall, the predominance of σS-dependent genes upregulated in stationary phase cells emphasized the accuracy of our data. With 1 123 differentially expressed genes, stationary planktonic cells were transcriptionally different from exponential planktonic cells (Fig. 2a), as reported elsewhere . Interestingly, three genes belonging to the same operon, cydA, cydB and ybgT (also named cydX), were under-expressed in exponential planktonic cells, and two of them, cydA and cydB, were selected as signature genes. In E. coli, the cyd operon encodes the three subunits of the cytochrome bd oxygen reductase complex, whose expression is induced under stressful growth conditions [28, 29]. The non-nutrient-limited early planktonic mode of growth explains the under-expression of this complex but also, more generally, the under-expression of pathways involved in energy production and conversion (see COG affiliation of clusters 5 and 8 in Fig. 3c and Table 1).
The response regulator CsgD, a master transcriptional regulator in biofilm formation, functions by assisting bacterial cells in transitioning from the planktonic stage to the multicellular state through the activation of expression of biofilm-linked genes [30, 31]. Accordingly, CsgD encoding gene was 25.0-fold overexpressed in 7 h-old biofilm compared to stationary planktonic growing cells, although its expression did not significantly change between the two sessile conditions. However, transcriptomic profiles of the 7 h-old and 13 h-old biofilm cells contained 290 differentially expressed CDS (∣fold-change∣ > 5 and adjusted P-value < 0.01) (Fig. 2a), which shows an evolution of the biofilm structure between these two time points and validates our experimental model. These findings are in agreement with those of previous studies showing distinct transcriptomic profiles in developing and confluent biofilm states [20, 21]. Genes of clusters 5 and 9 were specifically overexpressed in 7 h-old biofilm, showing that amino acid transport and metabolism (see COG affiliation in Table 1) is an essential process during the biofilm growth, as observed previously [32–34]. The bssS gene, encoding a biofilm regulator whose inactivation leads to an increase in both the biomass and thickness of biofilm in E. coli , was an under-expressed signature gene of the 7 h-old biofilm condition. In a more mature biofilm, 13 h-old biofilm, the overexpression of genes involved in carbohydrate transport and metabolism (cluster 3; Table 1) reflect the importance of sugar in the formation of the extracellular matrix, a crucial component for biofilm maturation . The ibpA gene was identified among the overexpressed signature genes of the 13 h-old biofilm condition, and encodes a heat shock protein whose overexpression is crucial in E. coli during biofilm growth .
The transcriptional pattern of bacteria harvested in the effluent was also specific. Surprisingly, according to K-means column clustering and the number of differentially expressed genes in the different conditions, biofilm-dispersed cells were transcriptionally closer to the 7 h-old biofilm cells than to the planktonic cells. Our results showed that dispersed cells represent a distinct stage in the bacteria lifecycle, different from both the planktonic and the biofilm states. Environmental pressure could then influence the fate of these cells converting them either into planktonic cells as suggested by Chua et al.  or into new biofilm structures.
Because spontaneously dispersed-cells were analyzed, the question of any potential input signal triggering the dispersion process was assessed. Quorum-sensing signaling is important for the proper regulation of biofilm development in several species, including K. pneumoniae [7, 37]. In our study, the operons lsrACDBFG and lsrRK encoding the regulatory network for AI-2 did present a strong up-regulation between 7 h-old biofilm and 13 h-old biofilm conditions. Interestingly, these genes were significantly under-expressed in dispersed cells compared to 13 h-old biofilm cells. Since the lsrACDBFG operon is transcriptionally regulated by both the LsrR repressor and the phosphoenolpyruvate phosphotransferase system (PTS), its expression could depend on the availability of certain substrates and the global metabolic status of the cell . In this way, our data suggested that lsr gene modulation and the subsequent down-regulation of the biofilm-linked genes trigger the dispersal process. Biofilm dispersal involving high concentrations of extracellular AI-2 was recently reported in E. faecalis and has been shown to be associated with phages release by sessile cells . A biofilm dispersal mechanism mediated by filamentous prophage-induced cell death has also been reported in P. aeruginosa [17, 39]. In our study, among the 10 transcriptional signature genes of biofilm-dispersed cells, pspA and pspB, encoding phage shock proteins A and B, were overexpressed (Fig. 4 and Table 2). Since the phage-shock protein A was overproduced in E. coli during filamentous phage infection [40, 41], it is tempting to hypothesize that the overexpression of the pspABCDE operon in K. pneumoniae dispersed cells is the consequence of bacteriophage activation, which leads to local cell death and therefore biofilm dispersal.
Since c-di-GMP depletion plays an important role in the dispersal from mature biofilms in many species [4, 42], we analyzed the expression of genes encoding proteins containing GGDEF (diguanylate cyclases) and EAL domains (phosphodiesterases), which catalyze the formation and the degradation of c-di-GMP, respectively. Two diguanylate cyclases encoding genes (CH1034_220201 and CH1034_50012) and one phosphodiesterase encoding gene (CH1034_280331 or mrkJ) were, respectively, under- and overexpressed in dispersed cells compared to 13 h-old biofilm cells. The phosphodiesterase activity of MrkJ in K. pneumoniae is an important factor in the regulation of type 3 fimbriae expression, which mediates the formation and disassembly of the biofilm . Among the other candidates potentially involved in the dispersal process, some degrading matrix enzyme-encoding genes were overexpressed in dispersed cells compared to 13 h-old biofilm, such as the protease-encoding gene ycbZ, the glucosidase-encoding gene malZ and the nucleases encoding genes endA, rnhB, nth, and yihG. Interestingly, genes involved in the SOS response (dinB, dinF, dinG, dinI, sulA, recA and recX) were also overexpressed in dispersed cells compared to 13 h-biofilm cells, suggesting a role of the stress response in biofilm dispersal. Although SOS stress response had not been directly related to biofilm dispersion, several studies reported the impact of nitrosative and nutrient stress on biofilm dispersal [13, 44]. Regarding the transcriptional status of the biofilm-dispersed cells, 21.9 and 9.3 % of the overexpressed genes in the K-means clusters 2 and 4, respectively, were categorized in the “translation, ribosomal structure and biogenesis” COG group (Fig. 3c). Dispersal probably requires high metabolic activity, even higher than that of the exponential planktonic cells. Indeed, only 4.3 and 3.5 % of the genes categorized in the K-means clusters 1 and 6, respectively (and therefore overexpressed in exponential planktonic condition), also belong to this COG group (Fig. 3c). However, ribosomal proteins could act not only in protein synthesis but also as regulators of the biofilm life cycle, as recently shown with the ribosomal proteins S11 (rpsK) and S21 (rpsU) in Bacillus subtilis . Another interesting feature of dispersed cells was the overexpression of cusA (Fig. 4 and Table 2), a member of the cusCFBA operon encoding a cation tripartite efflux pump involved in the detoxification of cooper and silver ions in the periplasm of E. coli . Two cusCFBA operons are present in the K. pneumoniae CH1034 genome and both were specifically overexpressed in dispersed cells (Additional file 2: Table S1). Because efflux systems have a major role in host colonization , we can therefore hypothesize that K. pneumoniae dispersed cells display specific phenotypes with high adaptive ability to colonize a new hostile environment. This hypothesis is reinforced by the fact that RyeE and t44, ncRNA genes, were overexpressed in dispersed cells (cluster 5, Fig. 3b); RyeE is upregulated in Yersinia pestis during lung infection  and the t44 expression level increases during initial invasion of fibroblast by Salmonella serovar Typhimurium .
Several works have already described the transcriptomic profile of biofilm cells [19–21] but none of them ever considered the overall cycle of bacterial life. The present study provides an exhaustive view of the transcriptional behavior of K. pneumoniae in the course of planktonic, biofilm formation and dispersion steps. By structuring data in clusters, we achieved a clear illustration of the specific expression profiles and functions, and identified signature genes as potential biomarkers of the different bacterial states. Further research on the genes evidenced in our work will provide a better understanding of the molecular mechanisms involved in the transition between planktonic, sessile and dispersed states.
Bacterial strains and culture conditions
K. pneumoniae CH1034 was grown in Lysogeny broth (LB) or in 0.4 % glucose M63B1 minimal medium (M63B1) at 37 °C with shaking and stored at −80 °C in LB broth containing 15 % glycerol. For subsequent RNA extraction, planktonic bacteria were cultured at 37 °C in M63B1 broth under aerobic conditions and harvested at OD620 = 0.25 (exponential phase) or after overnight growth (stationary phase).
GFP-tagged strain construction
The K. pneumoniae CH1034 GFP-tagged strain was constructed after replacement of the SHV-1 β-lactamase-encoding gene (chromosomal ampicillin resistance) by the selectable aadA7-gfpmut3 cassette. Briefly, the aadA7-gfpmut3 cassette flanked by 60-bp fragments, which correspond to the encoding upstream and downstream regions of shv, was generated using pKD4 plasmid as template, primers shv-GFP-Fw and shv-GFP-Rv and Phusion high-Fidelity DNA polymerase (Thermo Fisher Scientific, Waltham, Massachusetts, USA) according to the manufacturers’ recommendations. Primers were designed on the basis of information about the K. pneumoniae CH1034 genome sequence previously deposited in the ENA/EMBL-EBI database under the accession number: PRJEB9899 . The PCR fragment was then transformed by electroporation into the 0.4 % arabinose-induced K. pneumoniae CH1034 strain harboring the pKOBEG199, which contains the lambda-red proteins encoding genes under the control of a promoter induced by l-arabinose . The K. pneumoniae CH1034 GFP-tagged strain, named K. pneumoniae CH1034-gfp, was selected onto LB agar containing spectinomycin (70 μg/mL), and the loss of the pKOBEG199 plasmid was then checked by plating onto LB agar containing tetracycline (35 μg/mL).
Two types of flow-cell devices were used in this study, a flow-cell with three individual chambers (dimension: 35 x 1 x 5 mm; 175 mm3) to monitor biofilm development by confocal laser scanning microscopy, and a flow-cell with one chamber (dimension: 54 x 19 x 6 mm; 6156 mm3) for i) quantification and microscopic observations of the bacteria detached from biofilm, and ii) bacterial recovery for RNA-extraction. On both flow-cells, a glass cover slip ensuring a surface for biofilm development was glued with silicon glue (3 M, Saint Paul, Minnesota, USA). All components of the flow-cell system, including tubing, bubble traps, medium/waste bottles and flow-cell, were assembled as described previously . Before experiments, the system was sterilized by pumping 10 % (wt/vol) hypochlorite sodium for 1 h and then ethanol 100 % (vol/vol) for 15 min. Thereafter, the system was rinsed with M63B1 medium overnight at 37 °C. The inoculum composed of an overnight culture of K. pneumoniae CH1034 in M63B1 (4.106 and 108 cells for the three- and one-chamber flow-cells, respectively) was injected with a syringe into each compartment of the flow-cells. After 1 h of incubation at 37 °C without flow to allow bacterial adhesion, M63B1 medium was pumped at a constant rate of 0.08 mL/min (three-chamber flow-cell) or 0.9 mL/min (one-chamber flow-cell) through the devices.
Biofilm development was monitored in real time with an SP5 confocal laser microscope (Leica, Wetzlar, Germany) and a x40 oil objective. Images were processed with IMARIS software (Bitplane, Belfast, United Kingdom). Bacteria present in the effluent of the one-chamber flow-cell were observed with the Leica DM1000 optical microscope (Leica) and the Leica DFC295 camera (Leica). To quantify bacteria detached from the biofilm, viable bacteria present in the effluent were counted every hour for 16 h by serial dilution and plating on LB agar. For RNA extraction, biofilms developed on glass slide were recovered after 7 h or 13 h of incubation, and bacteria detached from the biofilm were recovered in the flow-cell effluent for 1 h after 12 h of incubation.
RNA-seq and RT-qPCR
For RNA-sequencing, total RNA was extracted from biological triplicate of planktonic, sessile or biofilm-detached bacteria prepared as described below. To avoid transcriptional changes and RNA degradation, all bacteria sampled were prepared in RNAlater® solution (Thermo Fisher Scientific) and then stored at 4 °C until RNA extraction. For exponential phase and stationary phase planktonic samples, an equivalent of 1010 CFU were pelleted by centrifugation at 6 000 g for 5 min at 4 °C, and pellets were resuspended in 2 mL of RNAlater® solution. To prepare the 7 h-old biofilm and the 13 h-old biofilm samples, biofilms developed on the glass slide of the flow-cell after the defined incubation period were scrapped in 1 mL of RNAlater® solution. In order to recover biofilm-detached bacteria, effluent of the flow-cells was directly collected in RNAlater® solution. After 1 h of collection, samples were centrifuged at 6 000 g for 5 min at 4 °C, and pellets were resuspended in 2 mL of RNAlater® solution. Before RNA extraction, bacteria were washed twice with 1X PBS. Total RNA was extracted according to the method described by Toledo-Arana et al. . Briefly, bacteria were mechanically lysed with the PreCellys 24 system (Bertin Technologies, Montigny le Bretonneux, France) at speed of 6 500 rpm for two consecutive cycles of 30 s. After acid phenol (Thermo Fisher Scientific) and TRIzol® (Thermo Fisher Scientific) extraction, total RNA was precipitated with isopropanol and treated with 10 units of TURBO DNase (Thermo Fisher Scientific). After a second phenol-chloroform extraction and ethanol precipitation, RNA pellets were suspended in DEPC-treated water. RNA concentrations were quantified with the Qubit system (Thermo Fisher Scientific) and RNA qualities were determined with Agilent RNA 6000 Pico chip (Agilent Technologies, Santa Clara, California, USA). Ribosomal RNA (rRNA) were removed from each total RNA sample with the Ribo-Zero Magnetic Kit (Bacteria) (Epicentre Biotechnologies, Madison, Wisconsin, USA), and rRNA-depleted samples were checked with Agilent RNA 6000 Pico chip. RNA-sequencing (RNA-seq) was conducted by MGX GenomiX (Montpellier, France). Libraries were produced by the Illumina TruSeq Stranded messenger RNA Sample Preparation Kit, and sequenced with the HiSeq 2000 system (Illumina, San Diego, California, USA) with a single-end protocol and read lengths of 50-bp. Short reads were mapped against the genome of K. pneumoniae CH1034 with the Burrows-Wheeler Alignment-backtrack mapper (version 0.7.12-r1039) , which allows a maximum of two mismatches within the first 32-bp. Counting was performed with the software HTSeq-count using the union mode. As data come from a strand-specific assay, the read has to be mapped to the reverse strand of the gene. Analysis of the reads mapped to intergenic regions confirmed the overall quality of the genome annotation and therefore strengthen the choice to focus on CDS and ncRNA features. Differentially expressed CDS and ncRNA genes between any pair comparisons of the five groups were determined by a negative binomial test with the DESeq package of R/Bioconductor. Transcripts were considered as differentially expressed using the following criteria: P-value < 0.01 and ∣fold-change∣ > 5. Transcriptome sequencing data were deposited in the Gene Expression Omnibus (GEO) database under the GEO accession number: GSE71754.
Reverse transcription was performed with 500 ng of total RNA prepared as described above, and the absence of DNA contamination was verified by qPCRs performed with primer pair RT-cpxR-Fw/RT-cpxR-Rv and the SsoAdvanced SYBR® Green Supermix (Bio-Rad, Hercules, California, USA) according to the manufacturer’s recommendations. cDNA were prepared with the iScript cDNA Synthesis kit (Bio-Rad) under the following conditions: 5 min at 25 °C, 30 min at 42 °C and 5 min at 85 °C. qPCRs were carried out in the CFX96 Real Time System (Bio-Rad) with the SsoAdvanced SYBR® Green Supermix (Bio-Rad) under the following conditions: initial denaturation at 95 °C for 30 s, and 40 cycles of 5 s at 95 °C and 20 s at 59 °C. qPCRs were performed in 10 μL total volume per well containing 1X SYBR® Green, 625 nM of each gene-specific primer and 2 μL of 20X diluted cDNA. Primers were designed on the basis of K. pneumoniae CH1034 genome sequence information  and are listed in Additional file 8: Table S3. Melting curve analysis was used to verify the specific single-product amplification. The gene expression levels were normalized relative to the expression levels of the cpxR housekeeping gene and relative quantifications were determined with CFX Manager software (Bio-Rad) by the E(−Delta Delta C(T)). The amplification efficiency (E) of each primer pair used for the quantification was calculated from a standard amplification curve obtained by four dilution series of genomic DNA. All assays were performed in technical triplicates with three independently isolated RNA samples.
Correlation between RNAseq and RT-qPCR was analyzed using Pearson’s correlation test in GraphPad Prism. Z-scores were calculated from the normalized DESeq expression data by the following formula: (X-Y)/Z (X: normalized DESeq counts of the sample; Y: average normalized DESeq counts of all the considered samples; Z: standard error of the counts mean for all the considered samples). Z-score values were used as a matrix to perform a principal component analysis and heatmaps with packages of R/Bioconductor: FactoMineR and Heatmap.2 (gplots), respectively. Column clustering was hierarchical, and two methods were used to cluster lines: hierarchical clustering and K-means clustering methods . K-means clustering was applied with different values of K (i.e. the number of clusters): 1 to 13. The clearest representation for each condition of the dataset was obtained with K = 10 for CDS clustering and K = 5 for ncRNA genes clustering. To highlight groups of CDS highly overexpressed or under-expressed in a specific condition, the mean of the Z-scores in each cluster was calculated for each condition, and the Z-score groups presenting a mean value > 1 or < −1 were named overexpressed boxes and under-expressed boxes, respectively.
The most relevant signature genes in the dataset were extracted using two fold-change thresholds, the Identity Threshold Fold-Change and the Differential Threshold Fold-Change. These thresholds were modulated as described in Figure S4 (Additional file 7) to obtain the most stringent signature genes for each condition.
Availability of supporting data
The RNA-seq data sets supporting the results of this article have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE71754 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE71754). All the supporting data are included as Additional files.
bis-(3'-5')-cyclic dimeric Guanosine Monophosphate
Coding DNA Sequences
colony forming unit
Clusters of Orthologous Groups
Gene Expression Omnibus
green fluorescent prrotein
- OD620 :
optical density at 620 nm
phosphate buffer saline
principal component analysis
polymerase chain reaction
quantitative Polymerase Chain Reaction
high-throughput sequencing of RNA
reverse transcription-quantitative polymerase chain reaction
- σS :
sigma S factor
Costerton JW, Stewart PS, Greenberg EP. Bacterial Biofilms: A Common Cause of Persistent Infections. Science. 1999;284(5418):1318–22.
Bogino PC, de las Mercedes Oliva M, Sorroche FG, Giordano W. The Role of Bacterial Biofilms and Surface Components in Plant-Bacterial Associations. Int J Mol Sci. 2013;14(8):15838–59.
Karatan E, Watnick P. Signals, Regulatory Networks, and Materials That Build and Break Bacterial Biofilms. Microbiol Mol Biol Rev. 2009;73(2):310–47.
Petrova OE, Cherny KE, Sauer K. The Diguanylate Cyclase GcbA Facilitates Pseudomonas aeruginosa Biofilm Dispersion by Activating BdlA. J Bacteriol. 2015;197(1):174–87.
Beloin C, Roux A, Ghigo J-M. Escherichia coli biofilms. Curr Top Microbiol Immunol. 2008;322:249–89.
Flemming H-C, Wingender J. The biofilm matrix. Nat Rev Microbiol. 2010;8(9):623–33.
Laverty G, Gorman SP, Gilmore BF. Biomolecular Mechanisms of Pseudomonas aeruginosa and Escherichia coli Biofilm Formation. Pathogens. 2014;3(3):596–632.
Hall-Stoodley L, Costerton JW, Stoodley P. Bacterial biofilms: from the Natural environment to infectious diseases. Nat Rev Microbiol. 2004;2(2):95–108.
Otto M. Staphylococcal infections: mechanisms of biofilm maturation and detachment as critical determinants of pathogenicity. Annu Rev Med. 2013;64:175–88.
Kaplan JB. Biofilm Dispersal. J Dent Res. 2010;89(3):205–18.
McDougald D, Rice SA, Barraud N, Steinberg PD, Kjelleberg S. Should we stay or should we go: mechanisms and ecological consequences for biofilm dispersal. Nat Rev Microbiol. 2012;10(1):39–50.
Kaplan JB, Ragunath C, Ramasubbu N, Fine DH. Detachment of Actinobacillus actinomycetemcomitans Biofilm Cells by an Endogenous β-Hexosaminidase Activity. J Bacteriol. 2003;185(16):4693–8.
Gjermansen M, Nilsson M, Yang L, Tolker-Nielsen T. Characterization of starvation-induced dispersion in Pseudomonas putida biofilms: genetic elements and molecular mechanisms. Mol Microbiol. 2010;75(4):815–26.
Cho C, Chande A, Gakhar L, Bakaletz LO, Jurcisek JA, Ketterer M, et al. Role of the Nuclease of Nontypeable Haemophilus influenzae in Dispersal of Organisms from Biofilms. Infect Immun. 2014. doi: 10.1128/IAI.02601-14
Wang R, Khan BA, Cheung GYC, Bach T-HL, Jameson-Lee M, Kong K-F, et al. Staphylococcus epidermidis surfactant peptides promote biofilm maturation and dissemination of biofilm-associated infection in mice. J Clin Invest. 2011;121(1):238–48.
Periasamy S, Joo H-S, Duong AC, Bach T-HL, Tan VY, Chatterjee SS, et al. How Staphylococcus aureus biofilms develop their characteristic structure. Proc Natl Acad Sci U S A. 2012;109(4):1281–6.
Rice SA, Tan CH, Mikkelsen PJ, Kung V, Woo J, Tay M, et al. The biofilm life cycle and virulence of Pseudomonas aeruginosa are dependent on a filamentous prophage. ISME J. 2009;3(3):271–82.
Rossmann FS, Racek T, Wobser D, Puchalka J, Rabener EM, Reiger M, et al. Phage-mediated Dispersal of Biofilm and Distribution of Bacterial Virulence Genes Is Induced by Quorum Sensing. PLoS Pathog. 2015;11(2):e1004653.
Beloin C, Valle J, Latour-Lambert P, Faure P, Kzreminski M, Balestrino D, et al. Global impact of mature biofilm lifestyle on Escherichia coli K-12 gene expression. Mol Microbiol. 2004;51(3):659–74.
Dötsch A, Eckweiler D, Schniederjans M, Zimmermann A, Jensen V, Scharfe M, et al. The Pseudomonas aeruginosa Transcriptome in Planktonic Cultures and Static Biofilms Using RNA Sequencing. PLoS ONE. 2012;7(2):e31092.
Rumbo-Feal S, Gómez MJ, Gayoso C, Alvarez-Fraga L, Cabral MP, Aransay AM, et al. Whole transcriptome analysis of Acinetobacter baumannii assessed by RNA-sequencing reveals different mRNA expression profiles in biofilm compared to planktonic cells. PLoS One. 2013;8(8):e72968.
Balestrino D, Ghigo J-M, Charbonnel N, Haagensen JAJ, Forestier C. The characterization of functions involved in the establishment and maturation of Klebsiella pneumoniae in vitro biofilm reveals dual roles for surface exopolysaccharides. Environ Microbiol. 2008;10(3):685–701.
Schroll C, Barken KB, Krogfelt KA, Struve C. Role of type 1 and type 3 fimbriae in Klebsiella pneumoniae biofilm formation. BMC Microbiol. 2010;10:179.
Chua SL, Liu Y, Yam JKH, Chen Y, Vejborg RM, Tan BGC, et al. Dispersed cells represent a distinct stage in the transition from bacterial biofilm to planktonic lifestyles. Nat Commun. 2014;5:4462.
Chua SL, Hultqvist LD, Yuan M, Rybtke M, Nielsen TE, Givskov M, et al. In vitro and in vivo generation and characterization of Pseudomonas aeruginosa biofilm-dispersed cells via c-di-GMP manipulation. Nat Protoc. 2015;10(8):1165–80.
Lange R, Hengge-Aronis R. Identification of a central regulator of stationary-phase gene expression in Escherichia coli. Mol Microbiol. 1991;5(1):49–59.
Dong T, Schellhorn HE. Global effect of RpoS on gene expression in pathogenic Escherichia coli O157:H7 strain EDL933. BMC Genomics. 2009;10:349.
Borisov VB, Gennis RB, Hemp J, Verkhovsky MI. The cytochrome bd respiratory oxygen reductases. Biochim Biophys Acta. 2011;1807(11):1398–413.
VanOrsdel CE, Bhatt S, Allen RJ, Brenner EP, Hobson JJ, Jamil A, et al. The Escherichia coli CydX Protein Is a Member of the CydAB Cytochrome bd Oxidase Complex and Is Required for Cytochrome bd Oxidase Activity. J Bacteriol. 2013;195(16):3640–50.
Mika F, Hengge R. Small RNAs in the control of RpoS, CsgD, and biofilm architecture of Escherichia coli. RNA Biol. 2014;11(5):494–507.
MacKenzie KD, Wang Y, Shivak DJ, Wong CS, Hoffman LJL, Lam S, et al. Bistable Expression of CsgD in Salmonella enterica Serovar Typhimurium Connects Virulence to Persistence. Infect Immun. 2015;83(6):2312–26.
Waite RD, Paccanaro A, Papakonstantinopoulou A, Hurst JM, Saqi M, Littler E, et al. Clustering of Pseudomonas aeruginosa transcriptomes from planktonic cultures, developing and mature biofilms reveals distinct expression profiles. BMC Genomics. 2006;7:162.
Valle J, Da Re S, Schmid S, Skurnik D, D’Ari R, Ghigo J-M. The Amino Acid Valine Is Secreted in Continuous-Flow Bacterial Biofilms. J Bacteriol. 2008;190(1):264–74.
Hamilton S, Bongaerts RJ, Mulholland F, Cochrane B, Porter J, Lucchini S, et al. The transcriptional programme of Salmonella enterica serovar Typhimurium reveals a key role for tryptophan metabolism in biofilms. BMC Genomics. 2009;10:599.
Domka J, Lee J, Wood TK. YliH (BssR) and YceP (BssS) Regulate Escherichia coli K-12 Biofilm Formation by Influencing Cell Signaling. Appl Environ Microbiol. 2006;72(4):2449–59.
Kuczyńska-Wiśnik D, Matuszewska E, Laskowska E. Escherichia coli heat-shock proteins IbpA and IbpB affect biofilm formation by influencing the level of extracellular indole. Microbiol Read Engl. 2010;156:148–57.
Solano C, Echeverz M, Lasa I. Biofilm dispersion and quorum sensing. Curr Opin Microbiol. 2014;18:96–104.
Pereira CS, Thompson JA, Xavier KB. AI-2-mediated signalling in bacteria. FEMS Microbiol Rev. 2013;37(2):156–81.
Webb JS, Thompson LS, James S, Charlton T, Tolker-Nielsen T, Koch B, et al. Cell Death in Pseudomonas aeruginosa Biofilm Development. J Bacteriol. 2003;185(15):4585–92.
Brissette JL, Russel M, Weiner L, Model P. Phage shock protein, a stress protein of Escherichia coli. Proc Natl Acad Sci U S A. 1990;87(3):862–6.
Darwin AJ. Stress Relief during Host Infection: The Phage Shock Protein Response Supports Bacterial Virulence in Various Ways. PLoS Pathog. 2013;9(7):e1003388.
Roy AB, Petrova OE, Sauer K. The Phosphodiesterase DipA (PA5017) Is Essential for Pseudomonas aeruginosa Biofilm Dispersion. J Bacteriol. 2012;194(11):2904–15.
Wilksch JJ, Yang J, Clements A, Gabbe JL, Short KR, Cao H, et al. MrkH, a Novel c-di-GMP-Dependent Transcriptional Activator, Controls Klebsiella pneumoniae Biofilm Formation by Regulating Type 3 Fimbriae Expression. PLoS Pathog. 2011;7(8):e1002204.
Barraud N, Hassett DJ, Hwang S-H, Rice SA, Kjelleberg S, Webb JS. Involvement of Nitric Oxide in Biofilm Dispersal of Pseudomonas aeruginosa. J Bacteriol. 2006;188(21):7344–53.
Takada H, Morita M, Shiwa Y, Sugimoto R, Suzuki S, Kawamura F, et al. Cell motility and biofilm formation in Bacillus subtilis are affected by the ribosomal proteins, S11 and S21. Biosci Biotechnol Biochem. 2014;78(5):898–907.
Chacón KN, Mealman TD, McEvoy MM, Blackburn NJ. Tracking metal ions through a Cu/Ag efflux pump assigns the functional roles of the periplasmic proteins. Proc Natl Acad Sci U S A. 2014;111(43):15373–8.
Guilhen C, Taha M-K, Veyrier FJ. Role of transition metal exporters in virulence: the example of Neisseria meningitidis. Front Cell Infect Microbiol. 2013;3:102.
Yan Y, Su S, Meng X, Ji X, Qu Y, Liu Z, et al. Determination of sRNA Expressions by RNA-seq in Yersinia pestis Grown In Vitro and during Infection. PLoS ONE. 2013;8(9):e74495.
Ortega ÁD, Gonzalo-Asensio J, Portillo FG. Dynamics of Salmonella small RNA expression in non-growing bacteria located inside eukaryotic cells. RNA Biol. 2012;9(4):469–88.
Guilhen C, Iltis A, Forestier C, Balestrino D. Genome Sequence of a Clinical Klebsiella pneumoniae Sequence Type 6 Strain. Genome Announc. 2015;3(6):e01311-5.
Weiss Nielsen M, Sternberg C, Molin S, Regenberg B. Pseudomonas aeruginosa and Saccharomyces cerevisiae Biofilm in Flow Cells. J Vis Exp JoVE. 2011;(47):e2383.
Toledo-Arana A, Dussurget O, Nikitas G, Sesto N, Guet-Revillet H, Balestrino D, et al. The Listeria transcriptional landscape from saprophytism to virulence. Nature. 2009;459(7249):950–6.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Sherlock G. Analysis of large-scale gene expression data. Curr Opin Immunol. 2000;12(2):201–5.
We thank Caroline Vachias, Pierre Pouchin and Jean-Louis Couderc for their technical help in confocal imaging acquisition and data analyses. We thank Marine Rohmer and Stéphanie Rialle for their assistance in RNAseq data analysis and GEO submission. We thank Sylvie Miquel for helpful discussion and critical reading of the manuscript. Cyril Guilhen is supported by a fellowship from Ministère de l’Education Nationale, de l’Enseignement Supérieur et de la Recherche. This work was supported by a ‘Contrat Quinquennal Recherche, UMR CNRS 6023’ and ‘Nouveau Chercheur 2012, Région Auvergne’.
The authors declare that they have no competing interests.
CG, CF and DB conceived and designed the experiments. CG, NC, NP, NG and AI performed the experiments. CG, CF and DB analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
Biofilm development of K. pneumoniae CH1034. K. pneumoniae CH1034-gfp was cultivated in flow-cell at 37 °C with a constant flux of medium. Biofilm development and maturation were monitored by confocal microscopy. The biofilm structure evolved from a flat to a three-dimensional structure. (MPG 3450 kb)
Data relative to the 2 052 selected CDS. (XLSX 974 kb)
Data relative to the 19 selected ncRNA genes. (XLSX 18 kb)
Determination of the correlation index between RNAseq and RT-qPCR data. Relative expression levels of 20 randomly selected genes were determined in bacteria collected in the effluent compared to the 13 h-old biofilm. The RNAseq and RT-qPCR ratios were then log2 transformed and values were plotted against each other to evaluate their correlation. The correlation coefficient was deduced from a linear regression of the plotted values using Pearson’s correlation test in GraphPad Prism. RT-qPCRs were performed with three biological replicates of total RNA extracts. Data were normalized to the endogenous reference gene cpxR, whose expression did not show significant variation between the tested conditions according to the RNAseq data. (PDF 123 kb)
Representation of the transcriptomic profiles of planktonic, sessile and biofilm-dispersed cells. The heatmap represents the hierarchical clustering of the Z-score of each of the 2 052 genes differentially expressed in at least one of the 10 possible pairs of conditions. Each condition was composed of three biological replicates, which were clustered together. Columns were clustered with the hierarchical clustering. (PDF 926 kb)
Clusters of Orthologous Group (COG) affiliation of the genes of each K-means cluster. The circle size is proportional to the percentage of genes (indicated by numbers) affiliated to a COG category for one given cluster group. Percentages in bold characters correspond to the major part of each cluster. (PDF 311 kb)
Strategy used for signature gene identification. Two thresholds were used: an “Identity Threshold Fold-Change” and a “Differential Threshold Fold-Change. Their respective values are indicated below. As an example, here is presented the strategy employed to identify one signature gene of the 13 h-old biofilm condition. The absolute expression (baseMeans) of the gene is represented by a filled circle in the 13 h-old biofilm condition, and by empty circles in the other conditions. Signature gene is defined according to two characteristics: i) differential expression levels between the 13 h-old biofilm condition (filled circle) and the other conditions (empty circles) higher than 4 (Differential Threshold Fold-Change), and ii) differential expression levels between all other conditions (empty circles) less than 2.5 (Identity Threshold Fold-Change). BaseMeans correspond to the absolute expression values averaged for triplicates of a condition as calculated by the DESeq package. (PDF 147 kb)
List of primers used in this study. (XLSX 12 kb)
About this article
Cite this article
Guilhen, C., Charbonnel, N., Parisot, N. et al. Transcriptional profiling of Klebsiella pneumoniae defines signatures for planktonic, sessile and biofilm-dispersed cells. BMC Genomics 17, 237 (2016). https://doi.org/10.1186/s12864-016-2557-x