Previous Article | Next Article ![]()
Infection and Immunity, May 2007, p. 2540-2547, Vol. 75, No. 5
0019-9567/07/$08.00+0 doi:10.1128/IAI.01957-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Department of Oral Biology and Center for Molecular Microbiology, College of Dentistry,1 Department of Molecular Genetics and Microbiology, College of Medicine, University of Florida, Gainesville, Florida2
Received 13 December 2006/ Returned for modification 28 January 2007/ Accepted 11 February 2007
|
|
|---|
|
|
|---|
More than 700 species, or phylotypes, of bacteria can inhabit the oral cavity (1). Temporally distinct patterns of microbial colonization result in biofilm formation on all surfaces in the oral cavity. On the tooth surfaces, the initial colonizers of the dental plaque biofilm are principally oral streptococci and actinomyces. Establishment of these organisms facilitates the subsequent colonization of additional actinomyces and related gram-positive rods along with gram-negative bacteria, such as Fusobacterium nucleatum. Further maturation is characterized by colonization by gram-negative anaerobes, such as Porphyromonas gingivalis (34, 52). Once colonization of the subgingival area has occurred, organisms shed from the plaque biofilm can interact with host epithelial cells that both have a barrier function and act as sensors of microbial infection (28). While many common oral organisms can adhere to gingival epithelial cells, only a subset of these organisms, including F. nucleatum, Aggregatibacter (Actinobacillus) actinomycetemcomitans, and P. gingivalis, can invade host cells (22, 23, 39, 59). Although it is well established that the bacterial inhabitants of the subgingival crevice are direct precursors of periodontal disease, the oral microbiota includes a spectrum ranging from commensals, such as Streptococcus gordonii, to aggressive pathogens, such as P. gingivalis. Some species, such as F. nucleatum, are located near the center of this spectrum as opportunistic commensals and are frequently found in individuals with good oral health but are also potentially able to contribute to disease (4, 15, 17, 44, 58, 62).
Transcriptional profiling using microarrays provides a way to monitor host cell responses to colonizing microorganisms on a global scale (8, 41). Numerous innate immune factors, for example, have been consistently found to be differentially regulated in host cells infected with pathogenic organisms compared to the regulation in uninfected controls (23, 29, 42). In the oral cavity, the pathogens P. gingivalis and A. actinomycetemcomitans induce widespread changes in the gingival epithelial cell transcriptome that are largely organism specific (23). There have been only a limited number of studies of the transcriptional responses to commensal organisms, and in these studies the focus is almost exclusively on nonoral mucosal ecosystems. For example, a study of the gastrointestinal tract commensal Bacteroides thetaiotaomicron showed it can modulate expression of genes involved in several important intestinal functions, including nutrient absorption, mucosal barrier fortification, xenobiotic metabolism, angiogenesis, and postnatal intestinal maturation (25). Commensal bacterial reconstitution of germfree mice has been shown to up-regulate expression of colonic epithelial cell genes associated with growth, apoptosis, and immune responses (19). In contrast, genes that may participate in extracellular oxidant defense and cellular metabolism were down-regulated by a nonpathogenic bacterial challenge (19). Hence, there is evidence that the host cell transcriptional response can be specific for the infecting organism and that commensals and pathogens can regulate distinct physiological functions in host cells.
In this study we utilized expression microarrays to investigate the transcriptional responses of oral epithelial cells to challenge with the commensal S. gordonii and the opportunistic commensal F. nucleatum. The transcriptional responses induced by these organisms were very similar to each other yet significantly different than the responses reported previously for oral pathogens (23). Among the biological processes altered most significantly in the host cells was the signal transduction pathway associated with mitogen-activated protein kinase (MAPK) and downstream effector molecules, including interleukins. Understanding how the host has adapted to commensals and how barrier cells respond to limit the impact of commensals should provide a mechanistic biological basis for health in the mixed bacterium-human ecosystem of the oral cavity.
|
|
|---|
Epithelial cells. Human immortalized gingival keratinocytes (HIGK) were originally generated by transfection of primary gingival epithelial cells with E6/E7 from human papillomavirus (45). HIGK were cultured in the presence of 5% CO2 in keratinocyte serum-free medium (Gibco/Invitrogen, Carlsbad, CA) supplemented with 0.05 mM calcium chloride, 200 mM L-glutamine, and an antibiotic-antimycotic (Invitrogen).
Microbe-host cell coculture. Bacteria in the mid-log phase were harvested and washed by centrifugation and then resuspended in antibiotic-free keratinocyte serum-free medium. HIGK (107 cells) were cocultured with bacteria to obtain a total association (adherent plus invading bacteria) of approximately 100 bacteria per epithelial cell. The numbers of adherent and invading organisms were confirmed in parallel experiments by plate counting (23). After 2 h of incubation at 37°C in the presence of 5% CO2, the HIGK were lysed with Trizol (Invitrogen) prior to RNA extraction. Cocultures were carried out in quadruplicate.
RNA isolation, cDNA-cRNA synthesis, and chip hybridization. RNA isolation, cDNA synthesis, labeled cRNA synthesis, and chip hybridization were performed as previously described (23). Briefly, total RNA was extracted from Trizol-lysed cells, treated with DNase I, purified, and quantified by using standard methods (QIAGEN, Valencia, CA, and Affymetrix, Santa Clara, CA). cDNA was synthesized by using the Affymetrix protocol (SuperScript double-stranded cDNA synthesis kit; Invitrogen); 5 to 8 µg of total cellular RNA was used as a template to amplify mRNA species for detection. Double-stranded cDNA was purified and used as a template for labeled cRNA synthesis. In vitro transcription was performed using a BioArray high-yield RNA transcript labeling kit (T7; Enzo Life Science, Farmingdale, NY) to incorporate biotinylated nucleotides. cRNA was subsequently fragmented and hybridized on Genechip human genome U133-A oligonucleotide arrays (Affymetrix) with proper controls. RNA samples were not pooled. The microarrays were hybridized for 16 h at 45°C, stained with phycoerythrin-conjugated streptavidin, and washed using the Affymetrix protocol (EukGE-WS2v4) with an Affymetrix fluidics station, and then they were scanned with an Affymetrix GeneChip 3000 scanner.
Microarray data analysis. Microarray data analysis was performed as previously described (23, 41). Briefly, expression filters were applied to remove Affymetrix controls and probe sets whose signals were not detected in all samples. The signal intensity values in the resulting data set were variance normalized, mean centered, and ranked by their coefficients of variation. Normalization was performed to give equal weight to all probe sets in the analysis, regardless of the order of magnitude of the raw signal intensity. To reduce the confounding effect of background signal variation on the analysis, only the half of the data set exhibiting the most variation across samples was used to perform an unsupervised hierarchical cluster analysis using the Cluster software (16). The resulting heat map and cluster dendrograms were visualized with the Treeview software (16) to determine the extent of characteristic host cell responses to each infection state, defined as identical treatments clustering together. Additional quality control data for the arrays are provided in the supplemental material.
Following the initial assessment of the host cell response to each condition, a supervised analysis was performed to investigate differences in gene regulation among experimental conditions. For this analysis, the raw signal intensities were log transformed for all probe sets that passed the initial expression filters and were correlated using BRB Array Tools (R. Simon and A. Peng-Lam, National Cancer Institute, Rockville, MD). In each supervised analysis, biological replicates were grouped into classes based on their infection states during coculture experiments, and probe sets significant at the P <0.001 level for the class were identified. To test the abilities of the significant probe sets to truly distinguish between the classes, leave-one-out cross-validation (LOOCV) studies were performed. In these LOOCV studies each array was left out in turn and a classifier was derived for the three groups by selecting probe sets significant at a P value of <0.001. The significant probe sets were then used with several prediction models (compound covariate predictor, nearest-neighbor predictor, and support vector machine predictor) to predict the class identity of the array that was left out and not included when the classification model was constructed. The ability of the classifier to correctly predict the class identity of the array that was left out was estimated using Monte Carlo simulations with 2,000 permutations of the data set.
Functional categorization by gene ontology and bioinformatics analyses. Gene ontology trees were populated using Pathway Express (33), available at http://vortex.cs.wayne.edu/projects.htm.
Immunoblotting.
HIGK were infected with F. nucleatum or S. gordonii as described above, using a time course of 1, 2, or 6 h. Cells were washed three times with phosphate-buffered saline and lysed in radioimmunoprecipitation buffer with proteinase inhibitors (Sigma, St. Louis, MO). Twenty-five micrograms of protein was separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis, electrophoretically transferred to polyvinylidene difluoride membranes (Millipore, Billerica, MA), and blocked with 5% skim milk in Tris-buffered saline-0.1% Tween 20. The membranes were incubated for 1 h with primary antibodies to GADD45
(1:200; Santa Cruz Biotechnology, Santa Cruz, CA), GADD45ß (1:200; Santa Cruz Biotechnology), DUSP4 (1:2,000; Abcam, Cambridge, MA), and ß-actin (1:10,000; Abcam). After three washes in Tris-buffered saline-0.1% Tween 20, each membrane was reacted with the species-appropriate peroxidase-coupled secondary antibody (1:1,000). Visualization was performed with the enhanced chemiluminescence system (Amersham, Buckinghamshire, United Kingdom). Band intensities were scanned and quantified using the Kodak 1D image analysis software (v.3.6.1).
Detection of cytokines. Supernatants of HIGK infected with F. nucleatum or S. gordonii were collected and filter sterilized. Interleukin-6 (IL-6) and IL-8 concentrations were determined by an enzyme-linked immunosorbent assay (Quantikine, Minneapolis MN) performed according to the manufacturer's protocol. Experiments were conducted in triplicate.
Microarray data accession numbers. The array results have been deposited in the GEO repository (http://www.ncbi.nlm.nih.gov/geo/index.cgi) under accession numbers GSM159477 to GSM159488, series GSE6927.
|
|
|---|
Comparison of the transcriptional profiles of HIGK infected with F. nucleatum and HIGK infected with S. gordonii. All samples of uninfected HIGK, F. nucleatum-infected cells, and S. gordonii-infected cells were used to determine the overall similarity of the transcriptional responses. After elimination of probe sets whose signals were not greater than background levels on all arrays, signal intensity data for the 12,125 probe sets that passed the initial expression filters were included in an unsupervised cluster analysis and supervised class prediction analysis. The unsupervised hierarchical cluster analysis revealed an infection state-dependent host cell transcriptional profile, as biological replicates clustered together (data not shown). In a supervised analysis based on the infection state, at a significance level of P < 0.001, 240 probe sets were differentially expressed. Assuming normality of the data set, the 240 significant genes is 20-fold higher than the 12 probe sets that would be expected by chance at a significance threshold of P < 0.001, given that 12,125 probe sets passed the expression filter. Treeview visualization of the 240 probe sets differentially expressed in the three classes is shown in Fig. 1. The major separation node occurred between uninfected cells and infected organisms, independent of bacterial species, and the differences between the uninfected and infected states exceeded the observable differences between F. nucleatum and S. gordonii, although organism-specific gene signatures could be discerned. Overall, 84% of the probe sets were regulated similarly by F. nucleatum and S. gordonii, which may reflect the degree of adaptation to the host of both these species.
![]() View larger version (24K): [in a new window] |
FIG. 1. Hierarchical clustering of variance-normalized gene expression data. The expression pattern of the cRNAs analyzed by using microarrays was determined by a supervised analysis of the variance-normalized data set of differentially expressed genes (P < 0.001, BRB ArrayTools) with the algorithm Cluster and was displayed with Treeview. Each row represents an individual gene element on the array, and each column represents the expression states of cRNAs for the challenge condition indicated. Each expression data point represents the relative fluorescence intensity of the cRNA from F. nucleatum-infected cells (columns Fn R1 to Fn R4) or S. gordonii-infected cells (columns Sg R1 to Sg R4) to the fluorescence intensity of the cRNA from uninfected cells (columns CTRL R1 to CTRL R4). The distance matrix used to determine the relatedness of samples through gene expression space was 1 Pearson's correlation coefficient. The cluster is subdivided into three groups consisting of genes that were repressed (green), genes that were induced (red), and genes whose expression did not change (black). The variation in gene expression for a given gene is expressed as the distance from the mean observation for that gene according to the color scale below the heat map.
|
|
View this table: [in a new window] |
TABLE 1. Ontology analysis of epithelial cell pathways impacted by infection with F. nucleatum or S. gordoniia
|
and GADD45ß were transcriptionally up-regulated following F. nucleatum infection (P = 0.0005 and P < 106, respectively), whereas expression of the gene encoding DUSP4 (one of several MAPK phosphatase family members) was down-regulated (P < 106). S. gordonii up-regulated GADD45ß gene expression (P < 106), down-regulated DUSP4 gene expression (P < 106), and had no detectable effect on GADD45
gene expression. The protein expression data obtained by Western blotting (Fig. 2) were consistent with the transcriptional data. Both F. nucleatum and S. gordonii decreased the levels of DUSP4; F. nucleatum up-regulated GADD45
and GADD45ß expression; and S. gordonii infection resulted in an increase in the GADD45ß level. In general, the amount of regulated protein continued to change over a 6-h period, corroborating the predictive power of the transcriptional "snapshot" for the longer-term phenotype of the epithelial cells, at least for these genes. The roles of these genes in signal transduction through the MAPK pathway and their impact on cell physiology are shown diagrammatically in Fig. 3.
![]() View larger version (32K): [in a new window] |
FIG. 2. Western immunoblots of HIGK infected with F. nucleatum (Fn) or S. gordonii (Sg) and uninfected controls (CTRL) for 1, 2, or 6 h. The blots were probed with antibodies to GADD45 (upper panel), GADD45ß (middle panel), and DUSP4 (lower panel) and then stripped and reprobed with antibodies to ß-actin. The graphs show the results of densitometric analyses of the ratio of test protein band intensity to ß-actin band intensity.
|
![]() View larger version (32K): [in a new window] |
FIG. 3. MAPK-related pathways containing F. nucleatum and S. gordonii differentially regulated genes at P < 0.05, adapted from Pathway Express and using the Kyoto Encyclopedia of Genes and Genomes nomenclature (see text for details). Red indicates up-regulation, green indicates down-regulation, and gray indicates no change in expression. +P indicates phosphorylation, and P indicates dephosphorylation. An arrow indicates a molecular interaction resulting in activation, and a line without an arrowhead indicates a molecular interaction resulting in inhibition. e, expression.
|
![]() View larger version (19K): [in a new window] |
FIG. 4. Enzyme-linked immunosorbent assay of IL-6 (A) and IL-8 (B) accumulation in HIGK supernatants following coculture with F. nucleatum or S. gordonii for 2, 4, 6, or 8 h. Ctrl, uninfected control. The error bars indicate standard deviations (n = 3). One asterisk, P < 0.05; two asterisks, P < 0.001 (as determined by a t test).
|
|
|
|---|
One pathway that was significantly impacted by F. nucleatum (P = 0.00825) and S. gordonii (P = 0.00437) involves MAPK signal transduction. This signaling pathway, which is evolutionarily conserved, connects cell surface receptors to regulatory targets within cells. MAPK signaling is involved in a vast array of physiological processes, including cell growth, migration, proliferation, differentiation, survival, development, and innate immunity (37, 49, 64). Consistent with these diverse roles, a variety of stimuli can activate MAPK pathways; these stimuli include growth factors, cytokines, ligands for G protein-coupled receptors, transforming agents, environmental stress, and viral and bacterial ligands. In the MAPK pathway, expression of GADD45
and GADD45ß was up-regulated by F. nucleatum infection, whereas expression of DUSP4 was down-regulated. Similarly, S. gordonii up-regulated GADD45ß expression and down-regulated DUSP4 expression. Regulation occurred at both the mRNA and protein levels. The GADD45 gene was originally identified as a gene that is rapidly induced by agents that cause DNA damage (18, 48). Transcriptional regulation of the GADD45 gene is mediated by both p53-dependent and p53-independent mechanisms (56), and GADD45 family members (
, ß, and
) are involved in the activation of the p38 and JNK pathways through MEKK4. Up-regulation of GADD45 expression may ultimately converge on growth arrest and on the activation of the nuclear transcription factor NF-
B (47, 56, 57, 63). The dual-specificity phosphatase DUSP4 is involved in the inactivation of the MAPKs by dephosphorylating both Thr and Tyr residues of ERK1 and ERK2 (30). One role of the DUSP family of phosphatases may be in the postinduction repression of MAPK activity (60). Hence, both up-regulation of GADD45 expression and down-regulation of DUSP4 expression by F. nucleatum or S. gordonii engender the concordant phenotype of enhanced information flow through the MAPK pathway. Together, these results provide a mechanistic framework for previous reports which showed that F. nucleatum and S. gordonii activated MAPK and NF-
B in gingival epithelial cells (7, 9, 26, 27, 35, 61).
One of the downstream targets of MAPK signaling is the production of cytokines. F. nucleatum induced IL-6 and IL-8 secretion, whereas S. gordonii repressed secretion of both of these proinflammatory cytokines. Differences in the cytokine secretion profiles may be related to differential regulation of individual components of the MAPK pathway, such as the ability of F. nucleatum to up-regulate both GADD45
expression and GADD45ß expression. In addition, F. nucleatum, but not S. gordonii, down-regulated expression of RSK2, whereas S. gordonii, but not F. nucleatum, down-regulated expression of Sap1a. RSK2 is a serine/threonine kinase that may play a role in mediating the growth factor- and stress-induced activation of the transcription factor CREB. SAP1a is a nuclear protein that stimulates transcription via the c-fos serum response element and also via an Ets binding site independent of the serum response factor. However, although RSK2 and Sap1a are involved in distinct aspects of MAPK signaling, their precise physiological roles are not known. The induction of expression of proinflammatory cytokines in gingival tissues and the subsequent inflammatory tissue damage are considered contributory factors in the pathogenesis of periodontal disease (10, 20, 53). Hence, stimulation of IL-6 and IL-8 by F. nucleatum may be one property that increases the pathogenicity of this organism compared to the pathogenicity of S. gordonii. Indeed, a general hyporesponsiveness to commensals may be advantageous in order to limit tissue destruction that might occur if a strong proinflammatory response were induced. Furthermore, it has been proposed that commensal species can "program" host cells to limit subsequent responses to more pathogenic organisms (2, 54) Conversely, in gingival epithelial cells, prior infection with the pathogen P. gingivalis can paralyze the local chemokine response to F. nucleatum (9). The extent to which challenge with S. gordonii can modulate subsequent host cell responses to other organisms is currently being investigated. Of possibly greater in vivo relevance is the potential of a complex microbial community to impact the transcriptional profile and phenotypic responses of host cells in a manner distinct from the summed activities of its constituents. Such community-based responses could partially explain apparent discrepancies that are observed between certain oral clinical manifestations and the putative pathogenic potentials of causative microbial species. For example, it is generally recognized that clinical cases of localized aggressive periodontitis that are associated with proapoptotic and proinflammatory A. actinomycetemcomitans do not result in significant inflammation or gingival destruction (3). It is tempting to speculate that the anti-inflammatory nature of certain commensal species has the potential to restrain the proinflammatory capability of pathogenic species.
Published ahead of print on 16 February 2007. ![]()
Supplemental material for this article may be found at http://iai.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»