Previous Article | Next Article ![]()
Journal of Bacteriology, February 2009, p. 1044-1055, Vol. 191, No. 3
0021-9193/09/$08.00+0 doi:10.1128/JB.01270-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Ching-Seng Ang,1,
Paul D. Veith,1
Helen L. Mitchell,1
Alvin W. H. Lo,1
Christine A. Seers,1
Katrina A. Walsh,1
Nada Slakeski,1
Dina Chen,1
J. Patricia Lissel,1
Catherine A. Butler,1
Neil M. O'Brien-Simpson,1
Ian G. Barr,2 and
Eric C. Reynolds1*
Cooperative Research Centre for Oral Health Science, Melbourne Dental School, Bio21 Institute, The University of Melbourne, Victoria, Australia,1 WHO Collaborating Centre for Reference and Research on Influenza, Parkville, Victoria, Australia2
Received 10 September 2008/ Accepted 13 November 2008
|
|
|---|
|
|
|---|
P. gingivalis is a black-pigmented, asaccharolytic, gram-negative, anaerobic coccobacillus that relies on the fermentation of amino acids for energy production (51). Like most bacteria, P. gingivalis has an essential growth requirement for iron. It preferentially acquires iron in the form of heme, a molecule comprised of a protoporphyrin IX (PPIX) ring with a coordinated central ferrous atom (55). This utilization of heme as an iron source may reflect the inability of P. gingivalis to synthesize PPIX de novo (47). Heme is preferentially obtained from hemoglobin and is acquired through the activity of the cell surface Arg- and Lys-specific proteinase-adhesin complexes (16, 54, 55), possibly in conjunction with TonB-linked outer membrane receptors, such as HmuR, and accessory proteins, such as HmuY (24, 32, 39). Unlike aerobic or facultative bacteria, P. gingivalis does not produce siderophores to chelate environmental iron and lacks the ferric reductase activity usually associated with siderophore-mediated iron acquisition (7, 9). P. gingivalis stores heme on its surface in the form of µ-oxo bis-heme, which has inherent catalase activity that helps to protect the cell from oxidative stress (58). The initiation and progression of periodontal disease are associated with bleeding at the site of disease, thereby providing an elevated level of hemoglobin. Environmental heme availability has been reported to affect the virulence of P. gingivalis in animal models of disease, although the exact effects on virulence are still debated. Many studies have indicated an increase in P. gingivalis virulence when the bacterium is grown under conditions of heme limitation (8, 24), and many putative virulence traits of P. gingivalis, including extracellular Arg- and Lys-specific proteolytic activity, are reported to be upregulated by growth under these conditions (12, 18).
To better understand the role of P. gingivalis in initiating chronic periodontitis, it is necessary to investigate the mechanism by which P. gingivalis establishes and proliferates in subgingival plaque. In periodontally healthy individuals, most subgingival plaque bacteria are likely to be growing under conditions of heme limitation, while during active disease this will change to heme-excess growth conditions. The response of P. gingivalis to environmental heme availability can be mapped on a global scale by transcriptomic analysis using DNA microarrays, or the more direct approach of proteomics can be employed. The traditional approach for comparative proteomics involves analyzing protein extracts from two different cell conditions on two different two-dimensional gels, staining and imaging the gels, and overlaying and aligning the two images using sophisticated analysis software. An alternative and rapidly advancing approach is based on stable isotope labeling of protein or peptides followed by mass spectrometry (MS) or tandem MS (MS/MS). This procedure involves labeling with chemically identical but isotopically different tags for the two sample conditions. The combined samples are then processed and analyzed by MS. One such isotope tag was pioneered by Gygi et al. (25) and termed the isotope-coded affinity tag (ICAT). This methodology involves labeling of proteins via alkylation of cysteinyl residues with either an isotopically light version (12C) or isotopically heavy version (13C) of the ICAT reagent. Of the 1,988 protein-encoding genes of P. gingivalis, 88% code for a protein containing at least one cysteine residue, and hence this technique should be suitable for the vast majority of P. gingivalis proteins.
Although many P. gingivalis proteins have been associated with growth under conditions of heme limitation (7, 17), no extensive work on the changes to the P. gingivalis W50 proteome or transcriptome during heme limitation has been reported. In order to gain insight into the response of P. gingivalis to heme limitation, we compared the transcriptome and proteome of P. gingivalis grown in continuous culture under conditions of heme excess or limitation. This combined approach revealed many changes to the P. gingivalis proteome and transcriptome, in particular the expression of proteins involved in glutamate/aspartate catabolism. This shift in catabolism under heme limitation was confirmed by organic acid end product analysis. Also, a protein upregulated under heme limitation was shown, by insertional inactivation, to be involved in epithelial cell invasion.
|
|
|---|
Growth and harvesting of P. gingivalis. P. gingivalis W50 was grown in continuous culture in a Bioflo 110 fermentor/bioreactor (New Brunswick Scientific) with a 400-ml working volume. The growth medium was 37 g liter–1 brain heart infusion medium (Oxoid) supplemented with 5 mg ml–1 filter-sterilized cysteine hydrochloride and either 5.0 µg ml–1 hemin (heme excess) or 0.1 µg ml–1 hemin (heme limitation). Growth was initiated by inoculating the culture vessel with a 24-h batch culture (100 ml) of P. gingivalis grown in the same medium (heme excess). After 24 h of batch culture growth, the medium reservoir pump was turned on and the medium 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 was maintained at 7.4 ± 0.1. The culture was continuously gassed with 5% CO2 in 95% N2. Cellular dry weights were determined in triplicate essentially as described previously (14).
Analysis of culture fluid organic acids. The concentrations of short-chain fatty acids in cell-free continuous culture supernatants (from uninoculated, heme-excess, and heme-limited cultures) were determined by capillary gas chromatography based on the derivatization method described by Richardson et al. (44).
Construction of ECR312 mutant. A P. gingivalis pg1374 mutant was constructed by allelic exchange with plasmid pAL36.1, containing an ermF cassette flanked by the pg1374 gene. Plasmid pAL36.1 was constructed in a stepwise manner as follows. A 672-bp 5' fragment of the pg1374 gene with flanking ApaI and AatII restriction sites (underlined) was generated from W50 genomic DNA by PCR with primers 5'-AGAGGGCCCTAGCAATCATTGCATTGCT and 5'-TGCGACGTCGTGTTACCAATAGAGGATT. This was then cloned into ApaI- and AatII-digested pAL30 to create pAL36. Plasmid pAL30 (previously constructed in our laboratory) contained an ermF cassette that had been PCR amplified, using pVA2198 (23) as the template and the primers 5'-CGCGGATCCCCGATAGCTTCCGCTATT and 5'-GCCGGTACCTCCATCGCCAATTTGCCA, and then ligated into pGEM-T Easy (Promega). A 3' 565-bp fragment of the pg1374 gene with flanking PstI and NdeI restriction sites was also PCR amplified, using primers 5'-TGACTGCAGGCTTTCGACCTTGGATCTT and 5'-TCGCATATGAATAAATAAGTGCCGTCGG, and cloned into PstI- and NdeI-digested pAL36 to produce pAL36.1. Plasmid pAL36.1 was linearized using ScaI and transformed into P. gingivalis W50 as previously described (14). Recombinant colonies were selected after 7 days of incubation at 37°C under anaerobic conditions on horse blood agar containing 10 µg ml–1 erythromycin. Confirmation of DNA integration was performed by PCR analysis (not shown), and the resulting mutant was designated ECR312.
Antibiotic protection invasion assay.
To compare the invasion efficiencies of W50 and ECR312, antibiotic protection assays were carried out in a 24-well cell culture plate as described previously, with three biological replicates (29). The human epithelial cell line KB (ATCC CCL-17) was maintained in Earl's minimum essential medium (EMEM) (JRH Biosciences) supplemented with 25 mM L-glutamine, 10% (vol/vol) heat-inactivated fetal calf serum (JRH Biosciences), 100 IU/ml penicillin-streptomycin (JRH Biosciences), and 30 µg/ml gentamicin (JRH Biosciences) in a 5% CO2 atmosphere at 37°C. Cells were transferred to 24 wells in EMEM without antibiotics and grown overnight as a monolayer to a cell density of
105 cells. P. gingivalis was cultured in 37 g/liter brain heart infusion medium (Oxoid) supplemented with 5 mg/ml filter-sterilized cysteine hydrochloride, 5.0 µg/ml hemin, and 5.0 µg/ml erythromycin (for ECR312) until a cell density of
2.9 x 109 cells/ml. The cells were then washed three times with phosphate-buffered saline (PBS), adjusted to the required cell concentrations, and resuspended in EMEM without antibiotics.
Prior to the invasion assay, the KB cell monolayer was washed twice with PBS before being infected with 200 µl of P. gingivalis at the designated cell densities (without centrifugation) for 90 min at 37°C. Following incubation, the supernatant containing the KB cells and bacteria was transferred to a sterile 2-ml tube. The remaining bound cells were detached using a cell scraper and three 500-µl aliquots of PBS. External nonadherent bacteria were then removed by washing the cells three times with 1.8 ml PBS at 500 x g for 5 min. External adherent P. gingivalis cells and those that were not removed by washing were killed by incubation for 1 h with constant rotation in 1.5 ml of EMEM supplemented with 300 µg of gentamicin and 200 µg of metronidazole per ml. These concentrations of antibiotics were sufficient to kill 109 bacteria per ml in 1 h (29). After exposure to the antibiotics, the cells were washed a further three times with PBS at 500 x g for 5 min, and the supernatant was then used as a control to make sure that no bacteria survived. The KB cells with internalized bacteria were then lysed with 1 ml of sterile water for 20 min, and 100 µl of the lysed cells was immediately transferred to 900 µl sterile brain heart infusion medium to prevent lysis of P. gingivalis from prolonged exposure to water. Internalized bacteria were then enumerated on horse blood agar plates after incubation under anaerobic conditions at 37°C for 7 days. Invasion efficiency was expressed as a percentage of the initial inoculum recovered after antibiotic treatment and lysis of KB cells. The negative controls used were KB cells preinoculated with paclitaxel (Taxol) and cytochalasin D to stabilize microtubules and inhibit actin polymerization, a process required for invasion.
Cell binding assays.
Assays of P. gingivalis binding to KB cells were carried out essentially as described previously (41). Briefly, P. gingivalis was grown to mid-log phase (
2.9 x 109 cells ml–1), washed three times with PBS, and labeled with fluorescein isothiocyanate (FITC; Molecular Probes) at a final concentration of 100 µg ml–1, with gentle shaking, for 45 min at 37°C. After being labeled, the cells were washed twice with PBS and resuspended in EMEM (JRH Biosciences), and the P. gingivalis cell density was adjusted to give final P. gingivalis/KB cell ratios of 250:1, 500:1, 750:1, 1,000:1, 1,250:1 and 1,500:1.
Binding of wild-type P. gingivalis W50 and the ECR312 mutant was carried out in parallel in 24-well plates by incubating 200 µl of each bacterial cell suspension with KB cells in 5% CO2 at 37°C for 40 min. Following incubation, the supernatant containing detached KB cells and bacteria was transferred to a 1.5-ml tube. The remaining attached KB cells were then removed from the well by trypsin-EDTA treatment and pooled with the corresponding collected supernatant. Complete EMEM (500 µl) was then added to inactivate the trypsin. The KB cells were centrifuged (400 x g, 5 min) at room temperature and washed twice in PBS to remove the unbound bacteria. The attachment of FITC-labeled bacteria to KB cells was measured by flow cytometry (FC 500; Beckman Coulter). KB cells were distinguished from the small bacterial cells as large granular cells. FITC-labeled bacteria attached to KB cells were measured using a 525-nm filter (FL1), and the relative amount of P. gingivalis attached was determined by the change in mean fluorescence intensity (MF1).
Preparation of samples for quantitative ICAT analysis.
P. gingivalis cells were harvested from continuous culture during steady-state growth, washed three times with wash buffer (50 mM Tris-HCl, pH 8.0, 150 mM NaCl, 5 mM MgCl2) at 5,000 x g for 30 min, and disrupted with three passes through a French pressure cell (SLM; AMINCO) at 138 MPa. The lysed cells were then centrifuged at 2,000 x g for 30 min to remove unbroken cells, followed by ultracentrifugation at 100,000 x g, producing soluble (supernatant) and membrane fractions. All procedures were carried out on ice. Protein labeling and separation were based on the geLC-MS/MS approach (33), using the cleavable ICAT reagent (Applied Biosystems). Protein was first precipitated using trichloroacetic acid (16%) and solubilized with 6 M urea, 5 mM EDTA, 0.05% sodium dodecyl sulfate (SDS), and 50 mM Tris-HCl, pH 8.3. Protein concentration was determined using the bicinchoninic acid protein reagent and adjusted to 1 mg ml–1. Protein (100 µg) from cultures under each growth condition was individually reduced using 2 µl of 50 mM Tris(2-carboxy-ethyl)phosphine hydrochloride for 1 h at 37°C. Reduced protein from heme-limited growth conditions was then alkylated with the ICATheavy reagent and protein from heme-excess growth conditions was alkylated with the ICATlight reagent. The two samples were then combined and subjected to SDS-polyacrylamide gel electrophoresis (PAGE) on a precast Novex 10% Nupage gel (Invitrogen). The gel was stained for 5 min, using SimplyBlue SafeStain (Invitrogen), followed by destaining with water. The gel was then divided into 20 sections from the top of the gel to the dye front by use of a scalpel. Each section was approximately 3.5 mm wide and was further sliced into small pieces (
1 mm3). Gel pieces were washed three times with a solution of 50% ethanol and 25 mM ammonium bicarbonate (ABC) buffer, followed by dehydration with 100% ethanol. Gel pieces were reswollen in 40 µl of modified sequencing-grade trypsin (Promega) at a concentration of 5 µg/ml in 25 mM ABC buffer, and 1 mM CaCl2 was added and incubated at 4°C for 30 min. Excess trypsin solution was then removed, and 15 µl of 25 mM ABC buffer was added. Digestion was carried out overnight at 37°C and terminated by trifluoroacetic acid (TFA) addition to a final concentration of 0.1% (vol/vol). The supernatant was then transferred to a microcentrifuge tube. To the gel pieces, 50 µl of 50% ethanol in 0.1% TFA was then added and sonicated for 15 min. The process was repeated, and all supernatants were pooled and reduced in volume by use of a vacuum centrifuge. The extracted peptides were then mixed with 500 µl of affinity loading buffer before being loaded onto an affinity column per the manufacturer's instructions (Applied Biosystems). Eluted peptides were dried and the biotin tag cleaved with undiluted TFA at 37°C for 2 h, followed by drying under a reduced vacuum. The dried samples were suspended in 35 µl of 5% acetonitrile in 0.1% TFA.
LC-MS. MS was carried out using an Esquire HCT ion-trap mass spectrometer (Bruker Daltonics) coupled to an UltiMate nano-LC system (LC Packings-Dionex). Separation was achieved using an LC Packings reversed-phase column (C18 PepMap100) (75-µm inner diameter by 15 cm by 3 µm by 100 Å) and eluted in 0.1% formic acid with the following acetonitrile gradient: 0 to 5 min, 0%; 5 to 10 min, 0 to 10%; 10 to 100 min, 10 to 50%; 100 to 120 min, 50 to 80; and 120 to 130 min, 80 to 100%.
The LC output was directly interfaced to the nanospray ion source. MS acquisitions were performed under an ion charge control of 100,000 in the m/z range of 300 to 1,500, with a maximum accumulation time of 100 ms. MS/MS acquisition obtained over a mass range of 100 to 3,000 m/z was performed on three precursor ions for the most intense multiply charged ions, with an active exclusion time of 2 min.
Protein identification for ICAT.
Peak lists were generated using DataAnalysis 3.2 (Bruker Daltonics), using the Apex peak finder algorithm, with a compound detection threshold of 10,000 and a signal-to-noise threshold of 5. Global charge limitations of +2 and +3 were set for exported data. Protein identification was achieved using the MASCOT search engine (MASCOT 2.1; Matrix Science) on MS/MS data queried against the P. gingivalis database (obtained from The Institute for Genomic Research [TIGR] website [www.tigr.org] under accession no. NC_002950.2). Parameters for database searching are listed in Table S1 in the supplemental material. The matched peptides were further evaluated using the following criteria: (i) peptides with a probability-based Mowse score corresponding to a P value of at most 0.05 were regarded as positively identified, where the score is –log X 10(log(P)) and P is the probability that the observed match is a random event; (ii) where only one peptide was used in the identification of a specific protein and the MASCOT score was below 30, manual verification of the spectra was performed; (iii) the heavy and light peptides of an ICAT pair must have exhibited closely eluting peaks, as determined from their extracted ion chromatograms; (iv) for proteins with a single unique peptide, this peptide must have been identified more than once (e.g., in different SDS-PAGE fractions or in both the light and heavy ICAT forms); and (v) if a single peptide did not meet criterion iv, then the MASCOT score must have been
25, the expectation value must have been
0.01, and the MS/MS spectrum must have exhibited a contiguous series of b- or y-type ions.
Estimation of false-positive results. To independently estimate the level of false-positive assignments, we created a reverse database of P. gingivalis peptides by reversing the order of the amino acid sequence for each protein such that the database is identical in size to the normal database in terms of the protein number, size, and distribution of amino acids (42). The false-positive rate was thus estimated as NR/NF, where NR is the number of peptides identified with the reverse database (MASCOT score of peptides is above the threshold for the reverse database) and NF is the number of peptides identified with the normal database (MASCOT score of peptides is above the threshold for the normal database).
Quantification of relative abundances of proteins.
The ratio of isotopically heavy 13C to light 12C ICAT-labeled peptides was determined using a script from DataAnalysis (Bruker Daltonics) and verified manually based on measurement of the monoisotopic peak intensity (signal intensity and peak area) in a single MS spectrum. The minimum ion count of parent ions used for quantification was 2,000, although >96% of both heavy and light precursor ions had counts of >10,000. In the case of poorly resolved spectra, the ratio was determined from the area of the reconstructed extracted ion chromatograms of the parent ions or the entire isotopic envelope. Averages were calculated for multiple peptides derived from a single parent protein, and outliers were removed using Grubb's test, with an
level of 0.05.
Extraction of nucleic acids for transcriptomic analysis. RNAs were extracted from 5-ml samples of P. gingivalis cells harvested directly from the chemostat. To each sample, 0.2 volume of RNA stabilization reagent (5% [vol/vol] phenol in absolute ethanol) was added (6). Cells were pelleted by centrifugation (9,000 x g, 5 min, 25°C), immediately frozen in liquid nitrogen, and stored at –70°C for later processing. Frozen cells were suspended in 1 ml of Trizol reagent (Invitrogen) per 1 x 1010 cells and then disrupted using lysing matrix B glass beads (MP Biomedicals) and a Precellys 24 homogenizer (Bertin Technologies, France). The glass beads were removed by centrifugation, and the RNA fraction was purified according to the Trizol manufacturer's (Invitrogen) protocol, except that ethanol (at a final concentration of 35%) rather than isopropanol was added at the RNA precipitation stage and samples were then transferred to the spin columns from an Illustra RNAspin Mini RNA isolation kit (GE Healthcare). RNA was purified according to the manufacturer's instructions from the binding step onwards, including on-column DNase treatment to remove any residual DNA. RNA integrity was determined using an Experion automated electrophoresis station (Bio-Rad).
Genomic DNA was extracted from P. gingivalis cells growing in continuous culture by use of a DNeasy blood and tissue kit (Qiagen) in accordance with the manufacturer's instructions.
Microarray design, hybridization, and analysis. Microarray slides were printed by the Australian Genome Research Facility and consisted of 1,977 custom-designed 60-mer oligonucleotide probes for the predicted protein coding regions of the P. gingivalis W83 genome (GenBank accession no. NC_002950.2) (37). Gene identifiers used in the text are based on the locus tags for this database. In addition, the microarray design included probes for several predicted open reading frames found exclusively in the database of the Los Alamos National Laboratory Oralgen Project. These are identified by the suffix "-L." Microarray sample pool control probes were included to aid intensity-dependent normalization (65). The full complement of probes was printed three times per microarray slide onto Corning UltraGAPS coated slides.
Slides were hybridized using either heme-excess or heme-limited samples labeled with Cy3, combined with a universal genomic DNA reference labeled with Cy5 (GE Lifesciences). cDNA was synthesized from 10 µg of total RNA, using the SuperScript Plus indirect cDNA labeling system (Invitrogen), with 5 µg of random hexamers (Invitrogen) for priming of the cDNA synthesis reaction. cDNA was labeled with Cy3 by use of an Amersham CyDye postlabeling reactive dye pack (GE Lifesciences) and purified using the purification module of the Invitrogen labeling system. Cy5-dUTP-labeled genomic cDNA was synthesized in a similar manner from 400 ng of DNA, using the BioPrime Plus Array CGH indirect genomic labeling system (Invitrogen).
Prior to hybridization, microarray slides were immersed in blocking solution (35% formamide, 1% bovine serum albumin, 0.1% SDS, 5x SSPE [1x SSPE is 150 mM NaCl, 10 mM NaH2PO4, 1 mM EDTA]) for 1 h at 42°C. After being blocked, slides were washed briefly in H2O followed by 99% ethanol and then dried by centrifugation. Labeled cDNAs were resuspended in 55 µl of hybridization buffer (35% formamide, 5x SSPE, 0.1% SDS, 0.1 mg ml–1 salmon sperm DNA), denatured at 95°C for 5 min, and then applied to slides and covered with LifterSlips (Erie Scientific). Hybridization was performed at 42°C for 16 h. Following hybridization, slides were washed successively in 0.1% SDS plus 2x SSC (1x SSC is 150 mM NaCl plus 15 mM sodium citrate) (5 min at 42°C; all further washes were performed at room temperature), 0.1% SDS plus 0.1x SSC (10 min), and 0.1x SSC (4 washes; 1 min each) and then quickly immersed in 0.01x SSC and then 99% ethanol with centrifugation to dry the slides.
Slides were scanned using a GenePix 4000B microarray scanner, and images were analyzed using GenePix Pro 6.0 software (Molecular Devices). Three slides were used for each treatment (heme limitation or heme excess), representing three biological replicates.
Image analysis was performed using GenePix Pro 6.0 software (Molecular Devices), and "morph" background values were used as the background estimates for further analysis. To identify differentially expressed genes, the LIMMA software package (59-61) was used, with a cutoff of P values of <0.005. Within-array normalization was performed by fitting a global loess curve through the microarray sample pool control spots and applying the curve to all other spots. The Benjamini-Hochberg method was used to control the false discovery rate to correct for multiple testing (5). Operon prediction was carried out using the Microbesonline website (2; http://microbesonline.org).
Microarray data accession number. The microarray data presented in this paper have been entered into the NCBI GEO data bank (www.ncbi.nlm.nih.gov/projects/geo) under accession number GSE13375.
|
|
|---|
Response of P. gingivalis to heme limitation, as determined using ICAT analysis.
To carry out a quantitative ICAT analysis of the P. gingivalis response to heme limitation, a geLC-MS approach based on that of Li et al. (33) was chosen as a means to reduce sample complexity. The ICAT-labeled soluble and insoluble protein fractions were separated independently by SDS-PAGE, and each gel lane was divided into 20 sections for in-gel tryptic digestion followed by affinity purification and LC-MS. In total, 142 unique proteins were identified. No positive matches to the reverse database were obtained, indicating a low level of false-positive identification. Of the proteins detected in both the soluble and insoluble fractions, 53 (34.0%) were identified based on the presence of two or more unique peptides, with a probability-based Mowse score corresponding to a P value of at most 0.05, which typically equates to a MASCOT score of
20. Sixty proteins (38.5%) were identified based on the presence of one unique peptide identified from two or more different fractions or both ICAT labeling states (of those, 58 had MASCOT scores of
25). Forty-three of the proteins (27.5%) were identified on the basis of a single unique peptide with a MASCOT score of
25, an expectation value of
0.01, a contiguous series of b- or y-type ions, and the intense ions being accounted for when interpreted manually. An example of protein identification based on the analysis of a single unique peptide is shown in Fig. 1.
![]() View larger version (32K): [in a new window] |
FIG. 1. Identification of PG0390 based on detection of a single peptide. (A) Total ion chromatogram. (B) Mass spectrum at 54.8 min. (Inset) Enlarged ICAT peptide ion pair at 703.9 and 708.4 m/z at a ratio of 1:2 (light to heavy). (C) Product ion spectrum for precursor 703.9 m/z. This peptide ion was identified as having the sequence LVDLNC*FDIK (MASCOT score = 40; C* denotes an ICAT-modified cysteine).
|
![]() View larger version (13K): [in a new window] |
FIG. 2. Distribution of protein abundances based on ratio of heme limitation over heme excess (every unit on the log2 scale indicates a twofold change). In total, the change in abundance of 142 proteins (including 14 proteins found in both the soluble and insoluble fractions) was determined, with 70 of these exhibiting a twofold or greater change in abundance.
|
|
View this table: [in a new window] |
TABLE 1. Proteomic and transcriptomic analyses of gene products involved in glutamate/aspartate catabolism in P. gingivalis during growth under conditions of heme limitation and heme excessa
|
|
View this table: [in a new window] |
TABLE 2. Proteomic and transcriptomic analyses of proteins of P. gingivalis grown with heme limitation or with heme excessa
|
![]() View larger version (36K): [in a new window] |
FIG. 3. Proposed catabolic pathways for glutamate and aspartate in P. gingivalis (27, 63). Enzymes identified by the proteomic analysis are represented by their TIGR accession numbers (Table 1). Shaded accession numbers indicate enzymes that were identified in the transcriptomic analysis. Underlined accession numbers indicate enzymes that were detected in a separate qualitative proteomic analysis (data not shown), not the current ICAT study.
|
In addition to the gene products discussed in more detail below, transcription of several genes of interest was significantly upregulated, including the genes of a putative operon of two genes, PG1874 and PG1875, one of which encodes hemolysin A; nine concatenated genes, PG1634 to PG1641, of which PG1638 encodes a putative thioredoxin; and PG1043, which encodes FeoB2, a manganese transporter (15). PG1858, which encodes a flavodoxin, was the most highly upregulated gene, at 15.29-fold. Of the 152 significantly upregulated genes,
55 currently have no predicted function.
Glutamate/aspartate catabolism. Twenty proteins involved in glutamate/aspartate catabolism were identified by the proteomic analysis, and 14 of these proteins were encoded by four predicted operons (Fig. 3; Table 1). Enzymes catalyzing six of the eight steps directly involved in the catabolism of glutamate to butyrate were identified and found to have increased 1.8- to 4-fold under conditions of heme limitation (Fig. 3; Table 1). Although the other two catalytic enzymes (PG0690 [4-hydroxybutyrate coenzyme A-transferase {4-hydroxybutyrate CoA-transferase}] and PG1066 [butyrate-acetoacetate CoA-transferase]) were not detected using ICAT analysis, PG0690 was identified as significantly upregulated in the microarray analysis, as were the other genes in this operon (Table 1). The effects of heme limitation on the abundances of the enzymes of the aspartate catabolic pathway were mixed: enzymes catalyzing the breakdown of aspartate to oxaloacetate in the oxidative degradation pathway were unchanged in abundance and the enzymes involved in the conversion of pyruvate to acetate showed an increase of 2- to 4.4-fold.
Analysis of organic acid end products. The amounts of acetate, butyrate, and propionate in the spent culture medium of P. gingivalis grown with a heme excess were 6.00 ± 0.36, 6.51 ± 0.04, and 0.66 ± 0.07 mmol g–1 bacterial cellular dry weight, respectively. Levels of acetate, butyrate, and propionate in the spent culture medium of P. gingivalis grown with heme limitation were 13.1 ± 1.8, 7.77 ± 0.40, and 0.71 ± 0.05 mmol g–1 bacterial cellular dry weight, respectively, indicating that the relative level of acetate had more than doubled during heme limitation, while the levels of butyrate and propionate remained at similar levels.
KB cell binding and invasion by P. gingivalis W50 and ECR312. Two P. gingivalis surface proteins possibly involved in invasion of host cells, PG0350 and the immunoreactive 47-kDa protein (PG1374), were identified by the proteomic analysis to increase in abundance during heme limitation, and the transcription of these two genes (PG0350 and PG1374) was also shown to be upregulated by the microarray analysis (Table 2).
To determine the role of PG1374 in epithelial cell invasion, the rate of invasion of KB cells by wild-type W50 was compared with that of the ECR312 mutant lacking a functional PG1374 gene. In triplicate antibiotic protection assays, an average of 33,000 ± 2,600 CFU of internalized P. gingivalis W50 was recovered from the KB cells, compared with 16,000 ± 1,100 CFU of P. gingivalis mutant ECR312 cells (P < 0.01; Student's t test), showing that the ability of the mutant to invade KB cells was impaired relative to that of the wild type. In a separate assay measuring the binding of bacteria to the surfaces of the KB cells, there was no significant difference in the binding of the ECR312 mutant compared with that of the wild-type W50 strain (Fig. 4), indicating that the reduction in invasiveness was not due to a difference in adherence.
![]() View larger version (16K): [in a new window] |
FIG. 4. P. gingivalis W50 ( ) and ECR312 ( ) binding to KB cells. The assay was carried out with two biological replicates (n = 6). (Inset) Gating of live KB cells based on forward and side scattering properties (FSC and SSC, respectively) (top left), with five peaks representing FITC fluorescence of P. gingivalis W50 cells bound to KB cells at P. gingivalis-to-KB cell ratios of 250 to 1,250 (top right) and to ECR312 cells (bottom right). MOI, multiplicity of infection.
|
|
|
|---|
The principal source of energy for P. gingivalis is the catabolism of amino acids which are obtained in peptide form by the hydrolysis of host proteins (14, 62, 63). P. gingivalis preferentially utilizes aspartate/asparagine, glutamate/glutamine, threonine, serine, leucine, and valine (14, 62, 63). The major catabolic pathways of P. gingivalis involve glutamate and aspartate, where glutamate is metabolized to butyrate, propionate, and ammonia and aspartate is metabolized to butyrate, propionate, acetate, and ammonia (Fig. 3). Heme limitation caused major changes to the abundances of enzymes in the glutamate and aspartate catabolic pathways. Enzymes of the glutamate pathway, encoded by PG0687 to PG0692, which catalyze the conversion of succinyl-CoA to crotonyl-CoA were found to be upregulated by both the proteomic and transcriptomic analyses. Two P. gingivalis proteins involved in glutamate catabolism, acyl-CoA dehydrogenase (PG1076) and the
subunit of an electron transfer flavoprotein (PG1078), significantly increased in abundance during heme limitation. These proteins are encoded by genes arranged in a predicted large operon of 15 genes. This operon also includes two genes encoding a hypothetical protein and a conserved hypothetical protein whose abundances increased during heme limitation (PG1067 and PG1068, respectively).
The conversion of fumarate to succinate via the pathway from aspartate is catalyzed by a heterotrimeric succinate-quinone oxidoreductase complex consisting of two cytoplasmic enzymes, FrdA (PG1615) and FrdB (PG1614), and a transmembrane protein, FrdC (PG1616) (3). The two cytoplasmic fumarate reductase (Frd) proteins showed similar decreases in abundance in response to heme limitation (Fig. 3; Table 1). Previous studies with Bacteroides fragilis have suggested that heme is required for synthesis of the cytochrome b-dependent Frd complex (34). It is therefore not surprising to see lower levels of Frd during heme-limited growth of P. gingivalis considering that it is unable to synthesize PPIX de novo. Growth studies of both B. fragilis and P. gingivalis have shown that they require heme for growth and that exogenous succinate can partially substitute for this requirement (4, 36). This observation was supported in B. fragilis by use of Frd-deficient mutants whose growth was not stimulated by heme but was stimulated by the addition of succinate (3). Molar growth yield studies further showed that the B. fragilis Frd-deficient mutants had a similar ATP yield to that of the heme-limited wild type. Fumarate respiration is the most widespread type of anaerobic respiration, and fumarate is the only metabolic intermediate known to serve as an electron acceptor, yielding
0.5 ATP/electron to form succinate as the end product (28). These studies demonstrate the importance of the conversion of aspartate to succinate for energy and efficient growth. Under heme-excess conditions, a portion of the aspartate catabolized by P. gingivalis is initially reduced to succinate via the succinate-quinone oxidoreductase enzyme complex (FrdBAC) and then catabolized via the glutamate pathway to produce butyrate (Fig. 3) (63). The three- to fourfold decrease in abundance of these fumarate reductase enzymes during heme limitation suggests that less succinate derived from aspartate enters the glutamate catabolic pathway, and as a consequence, most of the aspartate catabolized is converted to acetate (Fig. 3). Organic acid analysis of the spent culture medium from P. gingivalis grown in continuous culture showed that there was a twofold higher level of acetate produced under heme limitation than the level produced with heme excess, while the levels of the other major end products, butyrate and propionate, were similar under both growth conditions. This is consistent with our hypothesis of a shift in the pathway used for aspartate catabolism from energy-efficient anaerobic respiration to a less efficient process. The increase in abundance of the enzymes that catalyze the conversion of pyruvate to acetate (acetate kinase [PG1081], phosphotransacetylase [PG1082], and pyruvate ferredoxin oxidoreductase [PG0548]) is also consistent with the increased amounts of acetate found in the spent culture medium (Table 1; Fig. 3).
The role of host cell invasion by P. gingivalis in chronic periodontitis is not well defined, but it is believed that the bacterium is able to escape stress by hiding inside host cells and reemerging when the stress has passed (30). Since 99.9% of total body iron is intracellular, P. gingivalis may be able to obtain more iron/heme inside epithelial cells during periods of heme limitation. It is therefore tempting to speculate that heme-limited P. gingivalis may upregulate host cell invasion-related proteins to promote internalization by epithelial cells. Our proteomic and transcriptomic analyses indeed showed that the levels of the cell surface proteins PG1374 and PG0350 increased under conditions of heme limitation (Table 2). Epithelial cell invasion assays of the PG1374 isogenic mutant (ECR312) showed that PG1374 is required for maximal epithelial cell invasion, as the mutant was internalized at less than half the level of the wild type.
We have designated PG1374 an internalin-related protein (Irp) and believe it belongs to a new class of proteins with cysteine-containing leucine-rich repeats (LRR), similar to the Listeria monocytogenes internalin protein InlJ (46). In L. monocytogenes, there are at least 15 members of the internalin family, and all have been found to share certain structural features, including a signal peptide and an N-terminal LRR domain followed by a conserved interrepeat region. Many of these proteins are involved in host cell invasion processes (35), and although it has not been demonstrated why multiple internalins exist, they are proposed to confer tropism toward different cell types (20). From the sequence information and predicted structure, it has been shown that more than half of the residues in the internalin LRR domain face outwards and that these residues are variable, suggesting that they are involved in protein-protein interactions specific to the different internalin classes (49). In addition to containing a signal peptide and the LRR domain, PG1374 and PG0350 are part of a novel class of up to 34 P. gingivalis cell surface-located outer membrane proteins that have no significant sequence similarities apart from a conserved C-terminal domain of approximately 80 residues (50). The C-terminal domain has been shown to be essential for secretion and attachment of these proteins to the P. gingivalis cell surface. In addition, PG1374 is strongly immunogenic when probed with sera from human periodontitis patients (45), which is consistent with a cell surface location.
The process of internalization of P. gingivalis into gingival epithelial cells is thought to involve a coordinated process of attachment and invasion mediated by fimbriae and a variety of cell surface proteins (11, 31, 40). A P. gingivalis 33277 mutant lacking PG0350 was shown to exhibit similar invasive characteristics but reduced biofilm formation capability compared to wild-type bacteria (10, 66). The similar invasiveness of the PG0350– mutant to that of the wild type was attributed to the presence of fimbriae, which play a role in epithelial cell invasion by strain 33277, although there is also a possibility that the similar invasiveness of this mutant may have been due to the presence of the second LRR-containing internalin-related protein PG1374 that we have identified in this study.
The significant reduction in epithelial cell invasion and similarity of cell binding by the P. gingivalis W50 PG1374 mutant compared to the wild type (Fig. 4) demonstrate that the observed reduction in invasion of epithelial cells can be attributed to a defect in the invasion process, not to lesser adherence. Bacterial invasion has been shown to be a highly complex process involving numerous proteins and receptors (13). The involvement of multiple factors, including haloacid dehalogenase (PG0653), an endopeptidase (PepO; PG0159), a cation-transporting ATPase (PG1642), and an ATP-binding cassette transporter (PG2206), in P. gingivalis invasion has been demonstrated previously (40, 64). It is interesting that PG0159 and PG1642 were also identified as being upregulated under heme limitation conditions in our study (Table 2). The characterization of the role of PG1374 in epithelial cell invasion adds to the list of P. gingivalis proteins involved in host cell invasion.
The upregulation of characterized and putative heme/iron/PPIX transport systems under heme limitation conditions was notable in this study. The genes of the hmu locus (PG1551 to PG1553), encoding an outer membrane heme/hemoglobin binding and transport system (32) and including the gene encoding the TonB-linked outer membrane receptor protein HmuR, were all significantly upregulated. For P. gingivalis, the Hmu system has been characterized as the major heme transporter (56) and, in common with the P. gingivalis Iht heme transport system (17), has an accessory outer membrane lipoprotein (HmuY and IhtB, respectively) that has been proposed to interact with the TonB-linked outer membrane receptor to mediate iron complex uptake across the outer membrane.
All five genes of the htr locus (PG0644 to -8), encoding the components of a predicted heme ABC transport system and a TonB-linked outer membrane receptor that has been linked to heme transport (57), were also significantly upregulated under conditions of heme limitation. In addition, a hypothetical protein, PG1019, was observed to be 25 times more abundant under heme limitation conditions, and both PG1019 and PG1020 were significantly upregulated in the transcriptomic analysis. Bioinformatic analyses suggest that PG1019 is a lipoprotein that is encoded by a gene located immediately upstream of a gene encoding a putative outer membrane receptor protein (PG1020) in a predicted operon. Homology searches of PG1020 at the Pfam protein domain website (www.sanger.ac.uk/software/Pfam/) and multiple alignments of PG1020 with known P. gingivalis TonB-linked outer membrane receptors, such as PG0668 and PG2008, show the presence of a putative TonB box, a conserved stretch of residues near the N termini of these TonB-linked receptors (48), and a conserved N-terminal region which Simpson et al. (56) refer to as the TonB box IV region (not shown). The high abundance of PG1019 is consistent with this protein being an accessory lipoprotein to a TonB-linked system involved in the transport of iron/heme/PPIX into the cell.
The largest change in P. gingivalis protein abundance from heme-excess to heme-limited growth was observed with an alkyl hydroperoxide reductase protein, AhpC (PG0618) (Table 2), a peroxide-scavenging enzyme that has been shown to play an important role in peroxide resistance in P. gingivalis (38). The increase in abundance of alkyl hydroperoxide reductase was supported by the transcriptomic analysis, where both PG0618 (ahpC) and PG0619 (ahpF) were significantly upregulated. In P. gingivalis, formation of the µ-oxo bis heme form of iron PPIX with oxygen on the cell surface is thought to act as an oxidative buffer due to its inherent catalase-like activity (58). This surface layer may also serve as storage for iron and PPIX (52). During heme limitation, depletion of this source of iron PPIX was shown to result in an increased susceptibility to oxidative stress (18). The substantial increase in abundance of alkyl hydroperoxide reductase during heme limitation could therefore be in response to the increased oxidative stress caused by the absence/reduction of the µ-oxo bis heme layer. OxyR, an oxygen-sensitive transcriptional activator, also plays a role in the expression of alkyl hydroperoxide during anaerobic growth (19), as a P. gingivalis OxyR– mutant showed a 16-fold decrease in ahpC gene expression. More recently, Duran-Pinedo et al. (21) demonstrated the positive regulation of ahpC expression by the RprY response regulator. The substantial increase in abundance thus suggests that heme availability may have a role in RprY- and OxyR-controlled gene expression. Interestingly, the very high codon adaptation index (CAI) value of the gene encoding this protein (0.84) suggests that it can be highly expressed in the cell for rapid induction in response to such stress. Highly expressed genes in many bacteria have a strong composition bias in terms of codon usage. The CAI can be used to predict the expression level of a gene based on its codon sequence, with a higher CAI value indicating increased expression (53).
This study is the first combined quantitative proteomic and transcriptomic analysis of the responses of P. gingivalis to a change in environmental heme availability. P. gingivalis responds to limitation of the essential micronutrient heme by increasing the abundance of a number of proteins linked to metabolism, oxidative stress response, virulence, and invasion of host cells. We have demonstrated the presence of an L. monocytogenes internalin-like protein that is involved in the invasion of human epithelial cells by P. gingivalis. A shift in catabolic pathways leading to increased acetate production was detected upon heme limitation, which was consistent with a change to a less energy-efficient process.
We thank Anthony R. Bird of CSIRO Human Nutrition, Adelaide, for his kind assistance in the short-chain fatty acid analysis and Steven Cleal for technical assistance with DNA microarray design.
Published ahead of print on 21 November 2008. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
S.G.D. and C.-S.A. contributed equally to this work. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»