- Research article
- Open Access
Comparative transcriptomic analysis of Porphyromonas gingivalisbiofilm and planktonic cells
BMC Microbiology volume 9, Article number: 18 (2009)
Porphyromonas gingivalis in subgingival dental plaque, as part of a mature biofilm, has been strongly implicated in the onset and progression of chronic periodontitis. In this study using DNA microarray we compared the global gene expression of a P. gingivalis biofilm with that of its planktonic counterpart grown in the same continuous culture.
Approximately 18% (377 genes, at 1.5 fold or more, P-value < 0.01) of the P. gingivalis genome was differentially expressed when the bacterium was grown as a biofilm. Genes that were down-regulated in biofilm cells, relative to planktonic cells, included those involved in cell envelope biogenesis, DNA replication, energy production and biosynthesis of cofactors, prosthetic groups and carriers. A number of genes encoding transport and binding proteins were up-regulated in P. gingivalis biofilm cells. Several genes predicted to encode proteins involved in signal transduction and transcriptional regulation were differentially regulated and may be important in the regulation of biofilm growth.
This study analyzing global gene expression provides insight into the adaptive response of P. gingivalis to biofilm growth, in particular showing a down regulation of genes involved in growth and metabolic activity.
The gram-negative obligate anaerobe Porphyromonas gingivalis, in subgingival dental plaque, has been strongly implicated in the onset and progression of chronic periodontitis, a disease characterized by the destruction of the tooth supporting (periodontal) tissues [1, 2]. There is increasing evidence that P. gingivalis is also associated with systemic diseases such as atherosclerosis [3, 4] and preterm birth . P. gingivalis is an asaccharolytic organism that relies on the catabolism of amino acids for energy production and growth . An array of virulence factors has been associated with P. gingivalis pathogenicity, including proteases, adhesins, fimbriae and capsular polysaccharide [6, 7]. The persistence of P. gingivalis in subgingival plaque for periods sufficiently long enough to elicit disease is inherently dependent on it surviving as part of a mature biofilm. Although mutational analyses have been employed to study genes associated with biofilm development by P. gingivalis [8–14], very little is known about the nature of P. gingivalis physiology and the crucial regulatory processes occurring in the mature P. gingivalis biofilm and how this relates to pathogenicity. In our laboratory we have devised a reproducible continuous culture method to grow biofilm and planktonic cells simultaneously in the same fermentor vessel. Using this approach we have compared the cell envelope proteome of P. gingivalis W50 biofilm and planktonic cells . In this current study, we have expanded our investigation of these cells, comparing the global gene expression within P. gingivalis biofilm and planktonic cells using microarray analysis.
Continuous culture conditions and biofilm formation
The growth and physical characterization of the biofilm and planktonic cells analysed in this study have been described in Ang et al. . The continuous culture system allows the simultaneous co-culture of planktonic cells and biofilm cells under identical growth conditions . Briefly, the methods used were as follows. To produce biofilm and planktonic cells for RNA harvest P. gingivalis was grown in continuous culture, in duplicate, using a Bioflo 110 fermentor with a total volume of 400 mL (New Brunswick Scientific, Edison, NJ, USA) in BHI medium supplemented with 5 mg mL-1 cysteine hydrochloride and 5.0 μg mL-1 haemin. Growth was initiated by inoculating the fermentor vessel with a 24 hour batch culture (100 mL) of P. gingivalis grown in the same medium. After a 24 h incubation the media reservoir pump was turned on and the media flow adjusted to give a dilution rate of 0.1 h-1(mean generation time of 6.9 h). The temperature of the vessel was maintained at 37°C and the pH at 7.4 ± 0.1. The culture was continuously gassed with 5% CO2 in 95% N2. Optical density readings (OD650 nm) and purity of the culture were analyzed daily. Biofilm could be seen to be forming on the fermentor vessel walls and on glass microscope slides that were fixed to the vessel walls. Each P. gingivalis W50 culture was maintained for 40 days until harvesting. Planktonic cells were harvested by rapidly pumping them out of the fermentor vessel. The microscope slides were then removed from the fermentor vessel for examination of biofilm thickness and cell viability. The biofilm was rinsed twice with cold PGA buffer  to remove contaminating planktonic cells and then removed by scraping with a spatula and suspended in cold PGA buffer in a 50 mL centrifuge tube. PGA buffer contained 10.0 mM NaH2PO4, 10.0 mM KCl, 2.0 mM citric acid, 1.25 mM MgCl2, 20.0 μM CaCl2, 25.0 μM ZnCl2, 50.0 μM MnCl2, 5.0 μM CuCl2, 10.0 μM CoCl2, 5.0 μM H3BO3, 0.1 μM Na2MoO4, 10 mM cysteine-HCl with the pH adjusted to 7.5 with 5 M NaOH.
The viability of cells comprising the biofilms that were on the glass microscope slides were determined using LIVE/DEAD® BacLight™ stain as per manufacturer instructions (Invitrogen) with visualized using confocal laser scanning microscopy (CLSM) essentially as described by Loughlin et al. . CLSM was done using an Axiovert 200 M inverted microscope (Carl Zeiss Pty Ltd. Germany) fitted with a Zeiss LSM 510 META Confocal scan head. Imaging was carried out using the 458/477/488 nm Argon and 543 nm HeNe laser lines and a 63× C-Apochromat® water immersion lens. Live and dead cells in the stained biofilms were quantified using COMSTAT software  with the viability of the biofilm obtained by averaging the number of live cells over the entire z-stack . Biofilm thickness was also measured using light microscopy .
Total RNA extraction
P. gingivalis W50 biofilm and planktonic samples (40 mL) were immediately added to 0.125 volume of ice-cold Phenol solution (phenol saturated with 0.1 M citrate buffer, pH 4.3, Sigma-Aldrich, Inc. Saint Louis, MO). The mixture was centrifuged and the pellet suspended in 800 μL of ASE lysis buffer (20 mM Na acetate, 0.5% SDS, 1 mM EDTA pH 4.2) and transferred into a 2 mL microcentrifuge tube. An equal volume of ice cold Phenol solution was added and the mixture was vortexed for 30 s before incubation at 65°C for 5 min. The mixture was then chilled on ice for 3 min after which of 200 μL of chloroform was added and mixed by brief vortexing. The mixture was centrifuged at 16,100 × g and the aqueous phase collected and extracted using a Phenol solution/chloroform (1:1 vol:vol) mix. The RNA in the aqueous phase was precipitated by addition of 700 μL of 4 M LiCl and incubated overnight at -20°C. Samples were then thawed and the total RNAs were pelleted by centrifugation. The pellet was washed with cold 70% ethanol, air dried and suspended in 50 μL of 0.1% diethylpyrocarbonate treated water. The samples were then treated with DNase I (Promega, Madison, WI) and purified using RNeasy Mini columns (Qiagen, Valencia, CA) according to protocols supplied by the manufacturer. The quality of the total RNA was verified by analytical agarose gel electrophoresis and the concentration was determined spectrophotometrically.
Reverse transcription reactions contained 10 μg of total RNA, 5 μg of random hexamers, the first strand buffer [75 mM KCl, 50 mM Tris-HCl (pH 8.3), 3 mM MgCl2], 0.63 mM each of dATP, dCTP, and dGTP, 0.31 mM dTTP (Invitrogen Life Technologies, Carlsbad, CA) and 0.31 mM aminoallyl dUTP (Ambion, Austin TX), 5 mM DTT, and 800 u of SuperScript III reverse transcriptase (Invitrogen). The reaction mixture was incubated at 42°C for 2 h. The RNA was hydrolysed by incubation with 0.5 M EDTA and 1 M NaOH at 65°C for 15 min and the sample neutralized with 1 M HCl before purification of the cDNA with QIAquick columns (Qiagen). The cDNAs were coupled with monoreactive Cy3 or Cy5 (40 nmol) (Amersham Biosciences, Piscataway, NJ) in the presence of 0.1 M NaHCO3 for 60 min at room temperature. The labeled cDNAs were purified using QIAquick columns (Qiagen), combined and vacuum dried. Samples were then suspended in hybridization buffer containing 50% formamide, 10× SSC (150 mM sodium citrate, pH 7.0 and 1.5 M NaCl), 0.2% SDS and 1 μg μL-1 salmon sperm DNA.
P. gingivalis microarrays were kindly provided by The Institute for Genomic Research (TIGR) (now The J. Craig Venter Institute). Each microarray consisted of 1907 70-mer oligonucleotides spotted in quadruplicate on a glass slide (CMT-GAPS; Corning, Corning, N.Y.). Detailed array information can be viewed at http://www.tigr.org and http://www.brop.org. A total of four slides were used for each planktonic-biofilm pair, where the cDNAs were labeled with the alternative dye and hybridized to the microarray slides using a dye-swapping design.
Slides were prehybridized at 42°C in 5× SSC, 0.1% SDS and 2% bovine serum albumin for 2 h and then briefly rinsed with distilled water and isopropanol. Slides were dried by centrifugation for 3 min at 1,500 × g. The labeled cDNAs hybridization mix was heated to 100°C for 2 min before adding to the DNA microarray. Each array was covered with a coverslip and placed inside a hybridization chamber (Corning Incorporated Life Sciences, Acton, MA). Hybridization was carried out in a 42°C water bath for approximately 16 h after which the coverslips were removed and the slides washed in 2× SSC, 0.1% SDS at 42°C. The arrays were washed at room temperature once with 0.1× SSC, 0.1% SDS for 10 min, four times for 1 min in 0.1× SSC, and then rinsed with distilled water followed by 100% ethanol. The arrays were dried immediately by centrifugation (3 min, 1,000 × g).
Image and data analysis
The hybridized arrays were scanned using an Agilent G2565AA microarray scanner system (Agilent Technologies, Santa Clara, CA). Imagene 6.0 software (Biodiscovery, Los Angeles, CA) was used for spot finding, signal-background segmentation, and intensity quantification. The intensity of each spot was local background corrected using GeneSight 4.1 (Biodiscovery) and the resultant data were log transformed such that the mean value for each channel (Cy3 and Cy5) had a log ratio of zero. The signal intensities for each dye swap hybridization were combined and the average log ratios were used for all further analysis. The data were normalized using intensity dependent Lowess normalization  per spot and per slide to remove the intensity-dependent deviation in the log2 (ratio) values. Identification of differentially regulated genes was performed using the GeneSight 4.1 confidence analyzer [based on an ANOVA approach of Kerr et al ]. This statistical analysis uses replicate spots to estimate an empirical distribution of noise. The constructed noise model is then used to determine the statistical measures for the likelihood of false positives above or below a certain expression ratio. The differentially regulated genes were identified at 99% confidence intervals with a cut-off value of log2 > 0.6 or log2 < -0.6. These values correspond to approximately 1.5 fold up- and down-regulated genes, respectively, a ratio considered biologically relevant [21, 22]. The stringent 99% confidence interval was selected to reduce the chance of false positive genes to 1% (P-value < 0.01). All DNA microarray work in this study was in compliance with MIAME guidelines and all data have been deposited under accession number E-TABM-467, in the ArrayExpress databases http://www.ebi.ac.uk/arrayexpress.
Validation of microarray data by real time, reverse transcription-PCR
Total RNA (1 μg) was reverse transcribed to cDNA using SuperScript III First Strand Synthesis Supermix (Invitrogen) in the presence of random primers (50 ng) according to the manufacturer's recommendations. Real time-PCR was carried out using a Rotor-Gene 3000 (Corbett Research, Sydney, Australia). The primers for the real-time analysis (Table 1) were designed using Primer3 software http://primer3.sourceforge.net/. The lengths of the primers were 18 to 20 nucleotides and the amplified products between 109 and 130-bp. The amplification efficiency of each primer set was determined empirically by using cDNA template dilutions over four orders of magnitude. The amplification efficiency for each primer set varied between 95.4% and 106.6%, showing that the amplicons were generated with comparable efficiency.
The real time-PCR reaction contained 12.5 μL of Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen), 0.2 μM of each gene-specific primer and 5 μL of cDNA template. The cycling conditions were 50°C for 2 min, 95°C for 2 min, then 40 cycles of 95°C for 15 s, 58°C for 30 s, and 72°C for 30 s. Negative controls of distilled water and total RNA samples were included in each run. All reactions were carried out in triplicate and melting curve analysis indicated that in each reaction a single product was amplified.
PG0347 encoding a putative UDP-glucose 4-epimerase, galE, was selected as normalizer for all reactions. The critical threshold cycle, CT for each gene was generated by the Rotor-Gene 6 software (Corbett Research) and the relative expression ratio of the selected genes calculated and analyzed using the relative expression software tool (REST) http://www.gene-quantification.info. Each real time-PCR reaction was performed using the biological replicate total RNA samples that were used for microarray analysis.
Results and Discussion
P. gingivalisW50 growth in continuous culture and biofilm formation
P. gingivalis is a slow growing anaerobe that even in rich media has a generation time of 4.65 h . In the continuous culture system we employed here P. gingivalis W50 replicated with a mean generation time of 6.9 h and reached steady state approximately 10 days after inoculation. The cell density of the culture remained constant, after it had reached steady state, at an OD650 nm of 2.69 ± 0.21 and 2.80 ± 0.52 for the first and second biological replicates respectively. Robust biofilm was obtained on the vertical surfaces of the fermentor vessel walls and at 40 days of culture the planktonic and biofilm cells from the fermentor vessel were harvested for analysis. The glass microscope slides that were fixed to the fermentor vessel walls were used for physical characterization of the biofilm. CLSM revealed that the surface of the biofilm featured variable structures and the average percentage of viable cells within the biofilm was 91.2 ± 7.3% . The biofilms were on average 240 ± 88 μm thick. Our continuous culture system allowed us to obtain a direct paired comparison of transcriptomic profiles of both the planktonic and biofilm grown cells that were cultivated in the same fermentor vessel and therefore were subjected to identical gross environmental influences (such as media composition and temperature).
Identification of genes differentially regulated during biofilm growth
Microarray hybridizations were conducted using the paired planktonic cell and biofilm total RNA samples obtained from the two independent continuous cultures. For each culture planktonic cell and biofilm pair, four technical replicates of array hybridizations were performed (2 array slides for each dye swap) yielding 16 measurements per gene as each gene was represented in quadruplicate on each slide. We designated all genes with an average expression ratio of 1.5-fold (up or down) differentially regulated, a threshold reported to be biologically significant [21, 22]. Moreover, we used the GeneSight 4.1 (Biodiscovery) confidence analyzer to discriminate genes that had a 99% likelihood of being differentially regulated at above or below the 1.5 threshold.
A total of 561 and 568 genes were identified to be differentially regulated (1.5 fold or more, P-value < 0.01) between the biofilm and planktonic cells of the first and second replicates respectively (data not shown). Of the identified genes, 377 belonged to a common data set (67% and 66% of the total genes identified for the first and second replicates respectively). Of the 377 genes in the common dataset 191 were up-regulated and 186 were down-regulated (see Additional files 1 and 2). This represents approximately 18% of the P. gingivalis genome.
To validate the microarray data real time-PCR of selected genes PG0158, PG0270, PG0593, PG0914, PG1055, PG1431 and PG1432 was performed. Six of the genes were selected from the up-regulated group and one from the down-regulated group in biofilm cells. The expression of galE was detected to remain unchanged during biofilm and planktonic growth (data not shown) and was used for normalization. There was a high correlation between the expression ratios determined by both methods (R2 = 0.9002) (Fig. 1).
Although in some studies the differential expression of only a small percentage of the genome has been suggested following comparison of gene expression in biofilm and planktonic cells [25–28] differential expression of a large number of genes has been demonstrated in other studies. For example, in Escherichia coli, using gene-fusion studies, 38% (out of 446 clones) underwent altered expression during biofilm development . Sauer and co-workers demonstrated that more than 50% (over 800 proteins) of the Pseudomonas aeruginosa proteome was differentially regulated between planktonic growth and the fully mature biofilm . Moreover, DNA microarray analysis indicated that up to 22% (a total of 580 genes) of the Staphylococcus aureus genome underwent expression changes during mature biofilm growth .
Factors shown to be relevant to P. gingivalis homotypic biofilm formation and heterotypic biofilm formation with Streptococccus gordonii include InlJ, an internalin family-related protein , the minor fimbrial protein MfaI , ClpXP  and the low molecular weight tyrosine phosphatase Ltp1 . In the sequenced P. gingivalis strain W83  and in our laboratory strain W50 (data not shown) the mfa1 gene encoding Mfa1 has been interrupted by an insertion of the mobile element ISPg4. The microarray data indicated that in strain W50 biofilm cells there was increased expression of PG0176 which is the 5-prime region of mfa1. Thus there is an indication that in P. gingivalis strains where mfa1 is intact and Mfa1 produced that the minor fimbrillin may be upregulated in a biofilm. P. gingivalis coaggregation with S. gordonii mediated by MfaI is suggested to be relevant to P. gingivalis host colonizaton . Increased Mfa1 production may in some strains improve host colonization, but for strains such as W50 it would not play a role in their pathogenesis. Differential expression of the genes encoding InlJ (PG0320) and ClpXP (PG0417, PG0418) was not observed in the current study.
The predicted cellular roles of the differentially regulated P. gingivalis gene products in this study encompass widespread functional categories (Fig. 2). However, 40% (77) of the up-regulated genes and 31% (58) of the down-regulated genes were annotated as encoding hypothetical or conserved hypothetical proteins. Genes encoding proteins with similarity to experimentally identified proteins with unknown functions accounted for about 10% of the differentially expressed genes.
The physiology of the biofilm
The down-regulation of many genes involved in cell envelope biogenesis, biosynthesis of cofactors, prosthetic groups and carriers and other cellular processes was observed in this study (Fig. 2). Similarly, many genes involved in energy production, DNA replication, fatty acid and phospholipid metabolism and central intermediary metabolism were also down-regulated. Taken together, these observations suggest a down-turn in cell replication and a slowed growth rate in biofilm cells.
The primary indication of the slowing of cell replication in the biofilm was the down-regulation of genes encoding proteins involved in chromosome replication such as DnaA (PG0001), the primosomal protein n' PriA (PG2032), single-stranded binding protein Ssb (PG0271), the DNA polymerase III alpha subunit DnaE (PG0035) and the DNA polymerase III beta subunit DnaN(PG1853). Also down-regulated in biofilm cells were genes encoding homologues of proteins involved in DNA repair and recombination, MutS  (PG0412), radA  (PG0227) and recN [39, 40] (PG1849). The biofilm cells also displayed up-regulation of a putative translational regulator, RecX (PG0157) that in E. coli has been shown to inhibit RecA activity which is important in homologous recombination and in the SOS pathway of DNA repair and mutagenesis .
The down-regulation of a significant number of genes associated with cell envelope biogenesis (see Additional files 1 and 2) also suggests that the growth rate was reduced in biofilm cells. The slower growth rate of cells in a biofilm has been previously attributed to restricted penetration of nutrients and helps explain the relative resistance of biofilms to antibiotics targeting growth [42, 43]. As biofilm cells exhibit a slower growth rate then the need for energy would decrease correspondingly. Indeed, the transcriptomic data showed that expression of seven genes involved in the glutamate catabolism pathway, one of the key sources of energy for P. gingivalis , were simultaneously down-regulated in biofilm cells. One down-regulated gene in this pathway was PG1812 which is predicted to encode the alpha subunit of 2-oxoglutarate oxidoreductase, an enzyme located at the branching point in this pathway between butyrate and propionate end-products. Three genes PG0690, PG1075 and PG1076 encoding 4-hydroxybutyrate CoA-transferase, the coenzyme A transferase beta subunit and acyl-CoA dehydrogenase (short-chain specific) respectively, that are in the pathway branch that produces butyrate, were down-regulated, as were a cluster of genes encoding a methylmalonyl-CoA decarboxylase (PG1608-1612) that is part of the pathway branch that produces propionate.
Signal transduction, regulatory and transcription genes
It has been well established that two-component signal transduction systems (TCSTSs) play an important role in biofilm formation in many bacteria, including E. coli , Enterococcus faecalis  and Streptococcus mutans . Interrogation of the P. gingivalis W83 ORFs revealed only 6 putative TCSTSs. The transcriptomic analysis indicated that one of these TCSTSs, comprising PG1431 and PG1432, that encode a DNA-binding response regulator of the LuxR family and a putative sensor histidine kinase respectively, was up-regulated in biofilm cells. To date, the involvement of signal transduction, transcriptional regulators and other transcription factors in P. gingivalis biofilm development has yet to be established.
Homologues of the TCSTSs PG1431 and PG1432 have been found in P. gingivalis strain ATCC 33277 and were designated fimR and fimS, respectively . FimR and FimS are known to regulate FimA-associated fimbriation . Comparative transcriptomic profiling of P. gingivalis ATCC 33277 and its fimR deficient mutant indicated only a limited number of genes were part of the fimR regulon including PG1974, PG0644 (tlr) and the first gene of the fim locus, PG2130 . Binding of FimR upstream of PG2130 initiates an expression cascade involving PG2131-34. The transcriptomic data presented here do concur with the possible positive regulation of PG1974 by PG1431, however, they are in conflict with the role of PG1431 in the positive regulation of the fim locus because in strain W50 biofilms we observed decreased expression of PG2133 and PG2134 with no differential expression of fimA. Thus, the role of PG1431 and PG1432 in P. gingivalis W50 biofilm growth may not be reflected in the earlier study of P. gingivalis ATCC 33277 FimS and FimR mutants.
It is predicted that there are 29 orphan transcriptional regulatory proteins in P. gingivalis but only 4 of these were differentially regulated in biofilm cells, one of which was the down-regulated PG0270, oxyR. The remaining three possible transcriptional regulators PG0173, PG0826 (of the AraC family of transcriptional regulators) and PG2186 were found to be up-regulated. Members of the AraC family of transcriptional regulators have been shown to be important in carbon metabolism, stress response and virulence in other species (for review see Gallegos),  and in the regulation of quorum sensing signaling in P. aeruginosa .
Sigma factors are the subunit of RNA polymerase responsible for the recognition of the specific sequence of the target gene promoter  and are involved in the regulation of diverse physiological processes, particularly virulence [53, 54] and biofilm formation [55, 56]. The array data indicated that three putative sigma factors of the σ70 family PG0594 (rpoD), PG1660 and PG1827 were differentially regulated in biofilm cells. Both PG0594 and PG1660 were up-regulated whilst PG1827 was down-regulated in biofilm cells. The observed differential expression of these sigma factors in biofilm cells may indicate that these proteins are important regulators of P. gingivalis during biofilm growth.
Genes encoding transport and binding proteins
Many genes predicted to encode transport and binding proteins were up-regulated in biofilm cells (Fig. 2). Six of these genes encode components of putative ABC transporter systems (PG0280, PG0281, PG1175, PG1663, PG2199 and PG2206). PG1175 and PG1663 are each predicted to be the inner membrane components of an ABC transporter complex, each having an N-terminal transmembrane domain and a C-terminal ABC ATPase domain. Interestingly, a RPSBLAST search based on the conserved domain database CDD  revealed that PG0280 and PG0281 encode putative permeases belonging to the family which includes LolC that has been shown to transport lipids across the inner membrane .
Potential virulence determinants and hypothetical genes
The complete P. gingivalis genome sequence has revealed a number of putative virulence determinants, several of which were highly up-regulated in biofilm cells. These include a putative sialidase (PG0352) and ADP-heptose-LPS heptosyltransferase (PG1155) with an average fold change of 3.22 and 2.58 respectively, a putative extracellular protease (PG0553) and thiol protease, tpr (PG1055)  with average fold changes of 6.22 and 12.28 respectively. We also observed an increased expression of the gene encoding HtrA, a putative periplasmic serine protease (htrA; PG0593) with an average fold change of 2.96. HtrA is known to play a role in biofilm formation of Streptococcus mutans  and virulence in a variety of bacterial species [61–63]. In P. gingivalis, HtrA has been shown to confer protection against oxidative stress and be involved in long term adaptation to elevated temperature [64, 65]. HtrA has also been implicated in the modulation of the activity of the gingipain cysteine proteinases at elevated temperature but it is not essential for the maturation or activation of the gingipains under normal conditions . Interestingly htrA occurs in a predicted operon upstream of rpoD. In Salmonella enterica serovar Typhimurium [66, 67] and Yersinia enterocolitica [68, 69] an alternative sigma factor RpoE has been implicated in the regulation of htrA and resistance to oxidative stress. Taken together, these results suggest that perhaps HtrA in concert with RpoD may be part of a stress response that is activated during P. gingvalis biofilm growth.
The majority of the differentially regulated P. gingivalis genes were of unknown or poorly characterized function. Three of the genes encoding the hypothetical proteins, PG0914, PG0844, and PG1630 were also amongst the most highly up-regulated genes in biofilm cells with an average fold change of 11.69, 9.35 and 8.21 respectively. RPSBLAST search indicated that some of the hypothetical P. gingivalis proteins do have similarities to proteins of known function such as HslJ, a heat shock protein (PG0706) and DegQ, a trypsin-like serine proteases (PG0840) (Table 2).
Comparison of our microarray results with the cell envelope proteome analysis of P. gingivalis W50 biofilm and planktonic cells performed by Ang et al. , using the same cells as in this study, indicates that 5 out of the 47 proteins that were of differential abundance in that study correlate with the protein abundances (up or down-regulated) that could be expected based on our microarray data. While this correlation is modest, it is important to bear in mind that protein cellular distribution, stability, post-translation modifications and/or turnover may result in measured protein abundances that differ from those expected from the transcriptomic data [70–72]. Some P. gingivalis proteins known to be associated with the outer membrane and virulence of the bacterium, such as the gingipains (RgpA and Kgp), HagA and CPG70, that were of differential abundance in the proteome study of Ang et al.  were not shown to be differentially expressed at the transcript level in this study. One of these proteins, the Lys-specific gingipain proteinase Kgp (PG1844) has been shown to be a major virulence factor for P. gingivalis in assimilating the essential nutrient haem . In this current study the Kgp transcript level was unchanged between planktonic and biofilm growth. However, in the Ang et al.  study significantly less of the Kgp protein was found on the cell surface in the biofilm relative to planktonic cells. Kgp, along with other surface proteins, is known to be released from the cell surface by as yet undefined mechanisms to be present in the extracellular environment. Hence the transcriptomic and proteomic data from the same cells suggests that a major virulence factor, Kgp, may be released from the surface of the biofilm cells with no reduction in expression. This mobilization of a major virulence factor involved in assimilation of an essential nutrient may be an important survival mechanism for P. gingivalis in a biofilm.
It must be noted that the study presented here is of P. gingivalis grown as a monospecies biofilm and not as part of a multispecies biofilm as in subgingival dental plaque. Nonetheless the study does provide useful insights into the global events occurring when the bacterium is grown as a biofilm for an extended period, reflective of the chronic infection of the host. Analyses of P. gingivalis gene expression when it is grown as part of a multispecies biofilm are currently underway in our laboratory.
In this study, we have shown 18% of the P. gingivalis W50 genome exhibited altered expression upon mature biofilm growth. Despite the intrinsic spatial physiological heterogeneity of biofilm cells we were able to identify a large subset of genes that were consistently differentially regulated within our biofilm replicates. From the downturn in transcription of genes involved in cell envelope biogenesis, DNA replication, energy production and biosynthesis of cofactors, prosthetic groups and carriers, the transcriptomic profiling indicated a biofilm phenotype of slow growth rate and reduced metabolic activity. The altered gene expression profiles observed in this study reflect the adaptive response of P. gingivalis to survive in a mature biofilm.
Lamont RJ, Jenkinson HF: Life below the gum line: pathogenic mechanisms of Porphyromonas gingivalis. Microbiol Mol Biol Rev. 1998, 62: 1244-1263.
Griffen AL, Becker MR, Lyons SR, Moeschberger ML, Leys EJ: Prevalence of Porphyromonas gingivalis and periodontal health status. J Clin Microbiol. 1998, 36: 3239-3242.
Chun YH, Chun KR, Olguin D, Wang HL: Biological foundation for periodontitis as a potential risk factor for atherosclerosis. J Periodontal Res. 2005, 40: 87-95. 10.1111/j.1600-0765.2004.00771.x.
Offenbacher S, Jared HL, O'Reilly PG, Wells SR, Salvi GE, Lawrence HP, Socransky SS, Beck JD: Potential pathogenic mechanisms of periodontitis associated pregnancy complications. Ann Periodontol. 1998, 3: 233-250.
Shah H, Gharbia S: Batch culture and physiological properties. Biology of the species Porphyromonas gingivalis. Edited by: Shah HN, Mayrand D, Genco RJ. 1993, Florida: Boca Raton CRC Press Inc, 85-103.
Holt SC, Kesavalu L, Walker S, Genco CA: Virulence factors of Porphyromonas gingivalis. Periodontol 2000. 1999, 20: 168-238. 10.1111/j.1600-0757.1999.tb00162.x.
O'Brien-Simpson N, Veith PD, Dashper SG, Reynolds EC: Porphyromonas gingivalis gingipains: the molecular teeth of a microbial vampire. Curr Protein Pept Sci. 2003, 4: 409-426. 10.2174/1389203033487009.
Chen W, Palmer RJ, Kuramitsu HK: Role of polyphosphate kinase in biofilm formation by Porphyromonas gingivalis. Infect Immun. 2002, 70: 4708-4715. 10.1128/IAI.70.8.4708-4715.2002.
Davey ME, Duncan MJ: Enhanced biofilm formation and loss of capsule synthesis: deletion of a putative glycosyltransferase in Porphyromonas gingivalis. J Bacteriol. 2006, 188: 5510-5523. 10.1128/JB.01685-05.
Kuramitsu HK, Chen W, Ikegami A: Biofilm formation by the periodontopathic bacteria Treponema denticola and Porphyromonas gingivalis. J Periodontol. 2005, 76: 2047-2051. 10.1902/jop.2005.76.11-S.2047.
Lin X, Wu J, Xie H: Porphyromonas gingivalis minor fimbriae are required for cell-cell interactions. Infect Immun. 2006, 74: 6011-6015. 10.1128/IAI.00797-06.
Nakao R, Senpuku H, Watanabe H: Porphyromonas gingivalis galE is involved in lipopolysaccharide O-antigen synthesis and biofilm formation. Infect Immun. 2006, 74: 6145-6153. 10.1128/IAI.00261-06.
Capestany CA, Kuboniwa M, Jung IY, Park Y, Tribble GD, Lamont RJ: Role of the Porphyromonas gingivalis InlJ protein in homotypic and heterotypic biofilm development. Infect Immun. 2006, 74: 3002-3005. 10.1128/IAI.74.5.3002-3005.2006.
Chen W, Honma K, Sharma A, Kuramitsu HK: A universal stress protein of Porphyromonas gingivalis is involved in stress responses and biofilm formation. FEMS Microbiol Lett. 2006, 264: 15-21. 10.1111/j.1574-6968.2006.00426.x.
Ang CS, Veith PD, Dashper SG, Reynolds EC: Application of 16O/18O reverse proteolytic labeling to determine the effect of biofilm culture on the cell envelope proteome of Porphyromonas gingivalis W50. Proteomics. 2008, 8: 1645-1660. 10.1002/pmic.200700557.
Dashper SG, Brownfield L, Slakeski N, Zilm PS, Rogers AH, Reynolds EC: Sodium ion-driven serine/threonine transport in Porphyromonas gingivalis. J Bacteriol. 2001, 183: 4142-4148. 10.1128/JB.183.14.4142-4148.2001.
Loughlin PM, Cooke TG, George WD, Gray AJ, Stott DI, Going JJ: Quantifying tumour-infiltrating lymphocyte subsets: a practical immuno-histochemical method. J Immunol Methods. 2007, 321: 32-40. 10.1016/j.jim.2007.01.012.
Heydorn A, Nielsen AT, Hentzer M, Sternberg C, Givskov M, Ersboll BK, Molin S: Quantification of biofilm structures by the novel computer program COMSTAT. Microbiology. 2000, 146 (Pt 10): 2395-2407.
Cleveland W: Robust locally weighted regression and smoothing scatterplots. J Am Stat Assoc. 1974, 74: 829-836. 10.2307/2286407.
Kerr MK, Martin M, Churchill GA: Analysis of variance for gene expression microarray data. J Comput Biol. 2000, 7: 819-837. 10.1089/10665270050514954.
Hughes TR, Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD, et al: Functional discovery via a compendium of expression profiles. Cell. 2000, 102: 109-126. 10.1016/S0092-8674(00)00015-5.
Smoot LM, Smoot JC, Graham MR, Somerville GA, Sturdevant DE, Migliaccio CA, Sylva GL, Musser JM: Global differential gene expression in response to growth temperature alteration in group A Streptococcus. Proc Natl Acad Sci USA. 2001, 98: 10416-10421. 10.1073/pnas.191267598.
Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.
Milner P, Batten JE, Curtis MA: Development of a simple chemically defined medium for Porphyromonas gingivalis: requirement for alpha-ketoglutarate. FEMS Microbiol Lett. 1996, 140: 125-130.
Beloin C, Valle J, Latour-Lambert P, Faure P, Kzreminski M, Balestrino D, Haagensen JA, Molin S, Prensier G, Arbeille B, et al: Global impact of mature biofilm lifestyle on Escherichia coli K-12 gene expression. Mol Microbiol. 2004, 51: 659-674. 10.1046/j.1365-2958.2003.03865.x.
Schembri MA, Kjaergaard K, Klemm P: Global gene expression in Escherichia coli biofilms. Mol Microbiol. 2003, 48: 253-267. 10.1046/j.1365-2958.2003.03432.x.
Shemesh M, Tam A, Steinberg D: Differential gene expression profiling of Streptococcus mutans cultured under biofilm and planktonic conditions. Microbiology. 2007, 153: 1307-1317. 10.1099/mic.0.2006/002030-0.
Whiteley M, Bangera MG, Bumgarner RE, Parsek MR, Teitzel GM, Lory S, Greenberg EP: Gene expression in Pseudomonas aeruginosa biofilms. Nature. 2001, 413: 860-864. 10.1038/35101627.
Prigent-Combaret C, Vidal O, Dorel C, Lejeune P: Abiotic surface sensing and biofilm-dependent regulation of gene expression in Escherichia coli. J Bacteriol. 1999, 181: 5993-6002.
Sauer K, Camper AK, Ehrlich GD, Costerton JW, Davies DG: Pseudomonas aeruginosa displays multiple phenotypes during development as a biofilm. J Bacteriol. 2002, 184: 1140-1154. 10.1128/jb.184.4.1140-1154.2002.
Beenken KE, Dunman PM, McAleese F, Macapagal D, Murphy E, Projan SJ, Blevins JS, Smeltzer MS: Global gene expression in Staphylococcus aureus biofilms. J Bacteriol. 2004, 186: 4665-4684. 10.1128/JB.186.14.4665-4684.2004.
Park Y, Simionato MR, Sekiya K, Murakami Y, James D, Chen W, Hackett M, Yoshimura F, Demuth DR, Lamont RJ: Short fimbriae of Porphyromonas gingivalis and their role in coadhesion with Streptococcus gordonii. Infect Immun. 2005, 73: 3983-3989. 10.1128/IAI.73.7.3983-3989.2005.
Capestany CA, Tribble GD, Maeda K, Demuth DR, Lamont RJ: Role of the Clp system in stress tolerance, biofilm formation, and intracellular invasion in Porphyromonas gingivalis. J Bacteriol. 2008, 190: 1436-1446. 10.1128/JB.01632-07.
Maeda K, Tribble GD, Tucker CM, Anaya C, Shizukuishi S, Lewis JP, Demuth DR, Lamont RJ: A Porphyromonas gingivalis tyrosine phosphatase is a multifunctional regulator of virulence attributes. Mol Microbiol. 2008, 69: 1153-1164. 10.1111/j.1365-2958.2008.06338.x.
Nelson KE, Fleischmann RD, DeBoy RT, Paulsen IT, Fouts DE, Eisen JA, Daugherty SC, Dodson RJ, Durkin AS, Gwinn M, et al: Complete genome sequence of the oral pathogenic bacterium Porphyromonas gingivalis strain W83. J Bacteriol. 2003, 185: 5591-5601. 10.1128/JB.185.18.5591-5601.2003.
Lamont RJ, El-Sabaeny A, Park Y, Cook GS, Costerton JW, Demuth DR: Role of the Streptococcus gordonii SspB protein in the development of Porphyromonas gingivalis biofilms on streptococcal substrates. Microbiology. 2002, 148: 1627-1636.
Kunkel TA, Erie DA: DNA mismatch repair. Annu Rev Biochem. 2005, 74: 681-710. 10.1146/annurev.biochem.74.082803.133243.
Beam CE, Saveson CJ, Lovett ST: Role for radA/sms in recombination intermediate processing in Escherichia coli. J Bacteriol. 2002, 184: 6836-6844. 10.1128/JB.184.24.6836-6844.2002.
Picksley SM, Attfield PV, Lloyd RG: Repair of DNA double-strand breaks in Escherichia coli K12 requires a functional recN product. Mol Gen Genet. 1984, 195: 267-274. 10.1007/BF00332758.
Sanchez H, Alonso JC: Bacillus subtilis RecN binds and protects 3'-single-stranded DNA extensions in the presence of ATP. Nucleic Acids Res. 2005, 33: 2343-2350. 10.1093/nar/gki533.
Stohl EA, Brockman JP, Burkle KL, Morimatsu K, Kowalczykowski SC, Seifert HS: Escherichia coli RecX inhibits RecA recombinase and coprotease activities in vitro and in vivo. J Biol Chem. 2003, 278: 2278-2285. 10.1074/jbc.M210496200.
Gilbert P, Collier PJ, Brown MR: Influence of growth rate on susceptibility to antimicrobial agents: biofilms, cell cycle, dormancy, and stringent response. Antimicrob Agents Chemother. 1990, 34: 1865-1868.
Walters MC, Roe F, Bugnicourt A, Franklin MJ, Stewart PS: Contributions of antibiotic penetration, oxygen limitation, and low metabolic activity to tolerance of Pseudomonas aeruginosa biofilms to ciprofloxacin and tobramycin. Antimicrob Agents Chemother. 2003, 47: 317-323. 10.1128/AAC.47.1.317-323.2003.
Takahashi N, Sato T, Yamada T: Metabolic pathways for cytotoxic end product formation from glutamate- and aspartate-containing peptides by Porphyromonas gingivalis. J Bacteriol. 2000, 182: 4704-4710. 10.1128/JB.182.17.4704-4710.2000.
Ferrieres L, Clarke DJ: The RcsC sensor kinase is required for normal biofilm formation in Escherichia coli K-12 and controls the expression of a regulon in response to growth on a solid surface. Mol Microbiol. 2003, 50: 1665-1682. 10.1046/j.1365-2958.2003.03815.x.
Hancock LE, Perego M: The Enterococcus faecalis fsr two-component system controls biofilm development through production of gelatinase. J Bacteriol. 2004, 186: 5629-5639. 10.1128/JB.186.17.5629-5639.2004.
Li YH, Lau PC, Tang N, Svensater G, Ellen RP, Cvitkovitch DG: Novel two-component regulatory system involved in biofilm formation and acid resistance in Streptococcus mutans. J Bacteriol. 2002, 184: 6333-6342. 10.1128/JB.184.22.6333-6342.2002.
Hayashi J, Nishikawa K, Hirano R, Noguchi T, Yoshimura F: Identification of a two-component signal transduction system involved in fimbriation of Porphyromonas gingivalis. Microbiol Immunol. 2000, 44: 279-282.
Nishikawa K, Yoshimura F, Duncan MJ: A regulation cascade controls expression of Porphyromonas gingivalis fimbriae via the FimR response regulator. Mol Microbiol. 2004, 54: 546-560. 10.1111/j.1365-2958.2004.04291.x.
Gallegos MT, Schleif R, Bairoch A, Hofmann K, Ramos JL: Arac/XylS family of transcriptional regulators. Microbiol Mol Biol Rev. 1997, 61: 393-410.
Dong YH, Zhang XF, Xu JL, Tan AT, Zhang LH: VqsM, a novel AraC-type global regulator of quorum-sensing signalling and virulence in Pseudomonas aeruginosa. Mol Microbiol. 2005, 58: 552-564. 10.1111/j.1365-2958.2005.04851.x.
Raivio TL, Silhavy TJ: Periplasmic stress and ECF sigma factors. Annu Rev Microbiol. 2001, 55: 591-624. 10.1146/annurev.micro.55.1.591.
Libby SJ, Lesnick M, Hasegawa P, Weidenhammer E, Guiney DG: The Salmonella virulence plasmid spv genes are required for cytopathology in human monocyte-derived macrophages. Cell Microbiol. 2000, 2: 49-58. 10.1046/j.1462-5822.2000.00030.x.
Lingnau A, Domann E, Hudel M, Bock M, Nichterlein T, Wehland J, Chakraborty T: Expression of the Listeria monocytogenes EGD inlA and inlB genes, whose products mediate bacterial entry into tissue culture cell lines, by PrfA-dependent and -independent mechanisms. Infect Immun. 1995, 63: 3896-3903.
Adams JL, McLean RJ: Impact of rpoS deletion on Escherichia coli biofilms. Appl Environ Microbiol. 1999, 65: 4285-4287.
Kojic M, Venturi V: Regulation of rpoS gene expression in Pseudomonas: involvement of a TetR family regulator. J Bacteriol. 2001, 183: 3712-3720. 10.1128/JB.183.12.3712-3720.2001.
Marchler-Bauer A, Anderson JB, Cherukuri PF, DeWeese-Scott C, Geer LY, Gwadz M, He S, Hurwitz DI, Jackson JD, Ke Z, et al: CDD: a Conserved Domain Database for protein classification. Nucleic Acids Res. 2005, 33: D192-196. 10.1093/nar/gki069.
Narita S, Matsuyama S, Tokuda H: Lipoprotein trafficking in Escherichia coli. Arch Microbiol. 2004, 182: 1-6. 10.1007/s00203-004-0682-4.
Bourgeau G, Lapointe H, Peloquin P, Mayrand D: Cloning, expression, and sequencing of a protease gene (tpr) from Porphyromonasgingivalis W83 in Escherichia coli. Infect Immun. 1992, 60: 3186-3192.
Biswas S, Biswas I: Role of HtrA in surface protein expression and biofilm formation by Streptococcus mutans. Infect Immun. 2005, 73: 6923-6934. 10.1128/IAI.73.10.6923-6934.2005.
Cortes G, de Astorza B, Benedi VJ, Alberti S: Role of the htrA gene in Klebsiella pneumoniae virulence. Infect Immun. 2002, 70: 4772-4776. 10.1128/IAI.70.9.4772-4776.2002.
Stack HM, Sleator RD, Bowers M, Hill C, Gahan CG: Role for HtrA in stress induction and virulence potential in Listeria monocytogenes. Appl Environ Microbiol. 2005, 71: 4241-4247. 10.1128/AEM.71.8.4241-4247.2005.
Ibrahim YM, Kerr AR, McCluskey J, Mitchell TJ: Control of virulence by the two-component system CiaR/H is mediated via HtrA, a major virulence factor of Streptococcus pneumoniae. J Bacteriol. 2004, 186: 5258-5266. 10.1128/JB.186.16.5258-5266.2004.
Roy F, Vanterpool E, Fletcher HM: HtrA in Porphyromonas gingivalis can regulate growth and gingipain activity under stressful environmental conditions. Microbiology. 2006, 152: 3391-3398. 10.1099/mic.0.29147-0.
Yuan L, Rodrigues PH, Belanger M, Dunn WA, Progulske-Fox A: Porphyromonas gingivalis htrA is involved in cellular invasion and in vivo survival. Microbiology. 2008, 154: 1161-1169. 10.1099/mic.0.2007/015131-0.
Johnson K, Charles I, Dougan G, Pickard D, O'Gaora P, Costa G, Ali T, Miller I, Hormaeche C: The role of a stress-response protein in Salmonella typhimurium virulence. Mol Microbiol. 1991, 5: 401-407. 10.1111/j.1365-2958.1991.tb02122.x.
Humphreys S, Stevenson A, Bacon A, Weinhardt AB, Roberts M: The alternative sigma factor, sigmaE, is critically important for the virulence of Salmonella typhimurium. Infect Immun. 1999, 67: 1560-1568.
Heusipp G, Nelson KM, Schmidt MA, Miller VL: Regulation of htrA expression in Yersinia enterocolitica. FEMS Microbiol Lett. 2004, 231: 227-235. 10.1016/S0378-1097(03)00962-5.
Li SR, Dorrell N, Everest PH, Dougan G, Wren BW: Construction and characterization of a Yersinia enterocolitica O:8 high-temperature requirement (htrA) isogenic mutant. Infect Immun. 1996, 64: 2088-2094.
Corbin RW, Paliy O, Yang F, Shabanowitz J, Platt M, Lyons CE, Root K, McAuliffe J, Jordan MI, Kustu S, et al: Toward a protein profile of Escherichia coli: comparison to its transcription profile. Proc Natl Acad Sci USA. 2003, 100: 9232-9237. 10.1073/pnas.1533294100.
Eymann C, Homuth G, Scharf C, Hecker M: Bacillus subtilis functional genomics: global characterization of the stringent response by proteome and transcriptome analysis. J Bacteriol. 2002, 184: 2500-2520. 10.1128/JB.184.9.2500-2520.2002.
Pratt JM, Petty J, Riba-Garcia I, Robertson DH, Gaskell SJ, Oliver SG, Beynon RJ: Dynamics of protein turnover, a missing dimension in proteomics. Mol Cell Proteomics. 2002, 1: 579-591. 10.1074/mcp.M200046-MCP200.
This work was supported by the Australian National Health and Medical Research Council (Project Grant No. 300006) and Australian Government's Cooperative Research Centres program, through the Cooperative Research Centre for Oral Health Science. Microarray slides were kindly provided by TIGR and NIDCR. We also thank Rebecca Fitzgerald for helpful discussions on real time reverse transcription-PCR analysis.
The following material was obtained through NIAID's Pathogen Functional genomics Resource Center, managed and funded by Division of Microbiology and Infectious Diseases, NIAID, NIH, DHHS and operated by the J. Craig Center Institute.
CAS, SGD, NS and ECR designed the study. AWL performed the array and real time PCR analyses and wrote the initial draft of the manuscript. JPL carried out the continuous culture of P. gingivalis planktonic and biofilm cells. CAS, JB, SGD, NS and ECR revised the draft critically for important intellectual content. All authors have read and approved the manuscript.
Electronic supplementary material
Additional file 1: Genes differentially expressed in both P. gingivalis biofilm biological replicates arranged by functional category. The data provided represent the genes differentially expressed in P. gingivalis strain W50 biofilm grown cells relative to planktonic cells, arranged in order of predicted functional role of the gene product. (DOC 790 KB)
Additional file 2: Genes differentially expressed in both P. gingivalis biofilm biological replicates arranged by ORF number. The data provided represent the genes differentially expressed in P. gingivalis strain W50 biofilm grown cells relative to planktonic cells, arranged in order of TIGR ORF annotation. (DOC 703 KB)
About this article
Cite this article
Lo, A.W., Seers, C.A., Boyce, J.D. et al. Comparative transcriptomic analysis of Porphyromonas gingivalisbiofilm and planktonic cells. BMC Microbiol 9, 18 (2009) doi:10.1186/1471-2180-9-18
- Planktonic Cell
- Major Virulence Factor
- Average Fold Change
- Continuous Culture System
- Gingivalis Strain