Previous Article | Next Article ![]()
Infection and Immunity, December 2006, p. 6642-6655, Vol. 74, No. 12
0019-9567/06/$08.00+0 doi:10.1128/IAI.01250-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Microbiology and Immunology, Stanford University School of Medicine, Stanford, California 94305,1 Department of Biochemistry and Biophysics, University of California, San Francisco, San Francisco, California 94143,2 Biosciences Directorate, Lawrence Livermore National Laboratory, Livermore, California 945503
Received 4 August 2006/ Returned for modification 25 August 2006/ Accepted 1 September 2006
|
|
|---|
|
|
|---|
Francisella is found primarily within macrophages in the host (2). Therefore, Francisella-macrophage interactions have been studied in vitro as a model for in vivo interactions. Francisella readily enters unactivated macrophages and is initially contained within a vacuole. The Francisella-containing vacuole has been reported to avoid phagosome-lysosome fusion before the bacteria escape the phagocytic vacuole and replicate to high numbers in the cytosol (10, 17). It has also been reported that Francisella induces apoptosis (26) and abrogates Toll-like receptor signaling (44). However, when Francisella enters activated macrophages, cytosolic Francisella induces a host defense response, and the macrophage undergoes rapid caspase-1-dependent cell death (33). Several mutants that do not induce macrophage death in activated macrophages have been identified (33). Unlike wild-type bacteria, these mutants reside within a vacuole at 6 h postinfection and are attenuated for growth in macrophages and in mice (6, 35). Although it is known that replication in macrophages is important in the pathogenesis of tularemia, very little is known about the genetic and molecular mechanisms that govern Francisella interactions with macrophages.
Baron and Nano previously isolated a spontaneous Francisella mutant with altered colony morphology on plates containing the chromogenic phosphatase substrate XP, which indicated that this strain may have been a spontaneous phosphatase mutant or altered in the expression of an exported phosphatase (6). This mutant was unable to replicate in macrophages, and the mutated locus was identified by a complementation strategy and named mglAB, for macrophage growth locus (6). MglA and MglB are homologs of SspA and SspB from Escherichia coli, which are global regulators of the stringent starvation response. SspA is a transcription factor involved in the regulation of genes contributing to the stringent starvation response, acid resistance, and phage P1 late genes (19, 20, 49, 50). SspB is an accessory factor that determines substrate specificity for the ClpXP protease (31).
Additional studies have shown that the mglA transcript is induced inside macrophages (5), and MglA has subsequently been found to regulate the expression of several genes necessary for replication in macrophages and for virulence, including the intracellular growth locus genes iglA, iglC, and iglD and the pathogenicity determinant protein genes pdpD and pdpA (29). These genes lie within a
30-kb stretch of the genome that has characteristics of a pathogenicity island (13), including a lower G+C content in the DNA of this region as well as the presence of virulence genes (35). This region has been designated the Francisella pathogenicity island (FPI).
The sequencing of the Francisella genome (27, 35) allowed us to construct a Francisella microarray to examine global transcriptional responses. Since MglA was identified as a regulator of virulence genes, we began by investigating the transcriptional differences between wild-type Francisella and an mglA mutant during growth in rich media, conditions under which MglA-regulated genes are expressed (5). We identified 102 MglA-regulated genes, the majority of which are positively regulated, including the entire FPI. We mutated eight MglA-regulated genes both within and outside the FPI and show that five of these genes encode novel virulence factors utilized by Francisella during in vitro infections of macrophages and during in vivo infection in mice.
|
|
|---|
Microarray design. Open reading frames (ORFs) for the Francisella tularensis LVS (live vaccine strain) and SchuS4 genome sequences were predicted by using Glimmer2 (11, 40). Corresponding U112 genes are referred to by the SchuS4 locus tag. Oligonucleotides of 70 bp were selected for each ORF using ArrayOligoSelector (9). Oligonucleotides were designed to represent genes present in (i) SchuS4 and LVS, (ii) LVS and not SchuS4, or (iii) SchuS4 and not LVS based on a cutoff of 70% identity over 80% of the sequence length. Whenever possible, two 70-mer oligonucleotides were selected for each ORF (88%); however, this was not possible in some cases due to predicted cross-hybridization. The specificity of the designed oligonucleotides was verified in silico by an examination of BLAST similarity and the free energies of the oligonucleotides against their second-best matches in the genome. Oligonucleotides were purchased from QIAGEN (Valencia, CA) and resuspended in 3x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate). Glass slide preparation, printing, and processing after printing were performed as described previously (14). The final array contains 4,219 independent spots representing 2,242 ORFs.
Time courses and RNA isolation. Liquid cultures of U112 and mglA mutant strain GB2 were grown overnight at 37°C with shaking and subcultured to an optical density at 600 nm (OD600) of 0.006 in 300 ml. Samples were taken for serial dilutions to verify the OD600 and to determine the CFU for time zero, and the cultures were then placed into a shaking 37°C water bath to minimize changes in temperature. Samples were taken throughout the growth curve for RNA isolation and to determine the OD600 and CFU counts. Ten samples were taken during the first growth curve at 0.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, and 7.5 h. Seven samples were taken during the second growth curve at 1, 3, 5, 6, 7, 8, and 10 h. RNA was harvested as previously described (45). Briefly, bacteria were collected onto 0.45-µm nitrocellulose filters (Millipore, Bedford, Mass.) by vacuum filtration, placed into a 50-ml conical tube, and then immediately frozen in liquid N2 and stored at 80°C. RNA was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA) as previously described (45). The aqueous phase was used in the RNeasy Clean-Up protocol (QIAGEN, Chatsford, CA) with an on-column DNase digestion step for 40 min (QIAGEN, Chatsford, CA). The concentration and purity of RNA were assessed by measuring the OD260 and the OD280, and the integrity of the RNA was verified on a 1% agarose gel.
Preparation of cDNA probes, hybridization, and analysis.
For each of the time courses, equal quantities of RNA from each of the U112 and GB2 samples were combined to create a common reference for the experiment. For each time point or reference sample, cDNA was synthesized from 1 µg of total RNA in a standard reverse transcription (RT) reaction using 5 µg of random hexamers (Amersham Biosciences, Piscataway, NJ) and Superscript III (Invitrogen). Amino allyl dUTPs were incorporated at this step (2.5 mM dATP, 2.5 mM dCTP, 2.5 mM dGTP, 1 mM amino allyl dUTP, and 1.5 mM dTTP) (Invitrogen). cDNA was purified using Zymo DNA purification columns as specified by the manufacturer (Zymo Research Corp., Orange, CA), and samples from each time point were labeled with Cy5 (red) fluorophores and the reference was labeled with Cy3 (green) fluorophores as previously described (Amersham Biosciences) (45). Unincorporated fluorophores were quenched using 5 µl of 4 M hydroxylamine and incubated in the dark for 15 min. The samples from each time point and reference samples were then combined, and unincorporated dyes were removed using Zymo DNA purification columns. For hybridization, cDNA was eluted in 19 µl of Tris-EDTA buffer, and 2 µl of 20 mg/ml yeast tRNA (Invitrogen), 4.25 µl of 20x SSC, and 0.75 µl of 10% sodium dodecyl sulfate were added to the probe. Probes were denatured for 2 min at 99°C, spun at 17,900 x g, and cooled at room temperature before they were added to the arrays. The samples and arrays were incubated at 60°C for
14 h, and stringency washes were performed as described previously (14). Briefly, the arrays were washed for 2 min each in a series of four washes with increasing stringency: (i) 2x SSC-0.03% sodium dodecyl sulfate, (ii) 2x SSC, (iii) 1x SSC, and (iv) 0.2x SSC. The microarrays were scanned and analyzed using a Gene Pix 4000A scanner and GENEPIX5.1 software (both from Axon Instruments, Redwood City, CA).
Normalized data were collected using the Stanford Microarray Database (43). Spots were excluded from analysis due to obvious abnormalities, a regression correlation of <0.6, or a Cy3 net mean intensity of <100. Only spots with at least 70% good data across the experiment were included for analysis. The ratios of the red (time point sample) channels to green (reference) channels for each spot were expressed as log2(red/green) and used for hierarchical clustering using the CLUSTER program (15). Results were visualized using the TREEVIEW program (15). Using data from all of the microarrays and both growth curves, significant differences between the wild-type and mutant growth curves were determined using the Significance Analysis for Microarrays (SAM) program, v. 1.21 (46). Using the two-class statistical analysis tool in the SAM program, a list of genes whose transcript levels were significantly increased or decreased between the two groups was calculated. A calculated false discovery rate of <1% was used to assign significance, and a twofold cutoff in the change in expression level was imposed. The genes identified by SAM analysis are included in Table 1.
|
View this table: [in a new window] |
TABLE 1. List of MglA-regulated genesaa
|
|
View this table: [in a new window] |
TABLE 2. Primers used in this study
|
To select for a spontaneous streptomycin-resistant strain of GB2 for competition experiments, 100 µl of a culture grown overnight was plated onto modified MH agar supplemented with streptomycin (400 µg/ml). Individual streptomycin-resistant clones were passaged on modified MH agar with streptomycin to ensure that the clones were streptomycin resistant and tested for in vitro growth defects in rich media.
Complementation of mutants. The chloramphenicol cassette from pACYC184 (NEB) was amplified with primers CmF and CmR (Table 2), introducing ClaI and NotI restriction sites, respectively. The omp26 promoter was amplified from pLG62a (L. Gallagher and C. Manoil, unpublished data) with primers omp26f and omp26r (Table 2), introducing NotI and ClaI sites, respectively. The resulting PCR products were digested with ClaI and ligated to create a chloramphenicol resistance cassette driven by the omp26 promoter. This cassette was cloned into the NotI site of pMTM-T3 (34). Additionally, a chloramphenicol version of pDSK519 (23) was made by excising the kanamycin cassette using PvuII. The omp26 promoter-driven chloramphenicol cassette was amplified using omp26F_PvuII and cmR_PvuII, cut with PvuII, and ligated to the pDSK519 backbone. The resulting plasmid was called pDSC.
oppB, FTT1209c, and cds2 mutants were complemented in cis using a strategy similar to the mutagenesis strategy described above. In this strategy, the PCR of the 5' region included a full-length version of the gene of interest and sufficient upstream and downstream regions to allow for homologous recombination. The omp26 promoter-driven chloramphenicol cassette described above was cut with NotI and, along with the 5' and 3' digested regions, was ligated to SalI- and XbaI-digested pMPM-T6 omega (34). The plasmid was then transformed into the corresponding mutant, and allelic replacement was confirmed by screening for chloramphenicol resistance and kanamycin sensitivity. The FTT0989 mutant was complemented in trans by amplifying the FTT0989 gene and promoter using FTT0989compF1 and FTT0989compR1 (Table 2), introducing SpeI sites. pDSC was cut using XbaI and ligated to SpeI-cut FTT0989. Colonies were screened for the insertion of FTT0989 by PCR. Plasmid DNA was isolated from colonies with insertions and transformed into a chemically competent FTT0989 mutant as described previously (3). Complemented mutants were selected on Mueller-Hinton agar (Difco Laboratories) supplemented with chloramphenicol (7.5 µg/ml). All complemented clones were then sequenced and tested for their ability to grow in tryptic soy broth plus 0.2% cysteine by measuring the OD600 and counting CFU.
Bone-marrow-derived macrophages. Bone-marrow-derived macrophages were prepared as described elsewhere previously (41). Briefly, bone marrow was collected from the femurs of mice. Bone marrow cells were plated onto sterile petri dishes and incubated in Dulbecco's modified Eagle medium with high glucose, L-glutamine, 110 mg/liter sodium pyruvate, and pyroxidine hydrochloride (Invitrogen Life Technologies) and supplemented with 10% fetal calf serum (Invitrogen Life Technologies), 10 mM HEPES (Invitrogen Life Technologies), and 30% macrophage colony-stimulating factor (M-CSF)-conditioned medium. M-CSF-conditioned medium was collected from an L929 M-CSF cell line. Bone marrow cells were incubated at 37°C under 5% CO2, and macrophages were harvested after 6 days. All assays were performed using standard tissue culture plates at 37°C under 5% CO2 in similar media, excluding M-CSF-conditioned medium.
Macrophage replication assay. Bone-marrow-derived macrophages were seeded in 24-well plates at a density of 2.5 x 105 macrophages per well and incubated at 37°C and 5% CO2 overnight. The macrophages were then infected at a multiplicity of infection (MOI) of 25:1 (bacteria-to-macrophage ratio) and centrifuged at 730 x g for 15 min to mediate attachment. The macrophages were incubated (time zero) for 1 h and washed three times with warm media. Macrophages were lysed at 2 and 10 h postinfection with 1% saponin in water for 5 min and then diluted in phosphate-buffered saline and plated onto tryptic soy agar with 0.1% cysteine. Plates were incubated at 37°C and 5% CO2 for 2 days, and colonies were enumerated. Statistical analysis was performed using Student's two-tailed t test for independent means.
Cytotoxicity. Bone-marrow-derived macrophages were seeded in 96-well plates at 5 x 104 cells per well and incubated overnight with heat-killed F. novicida, with approximately 100 heat-killed bacteria per macrophage. Cells were infected at a multiplicity of infection of 100:1 or 10:1 (bacteria-to-macrophage ratio) and centrifuged at 730 x g for 15 min to mediate attachment. Medium was removed after 1 h and replaced with phenol red-free Dulbecco's modified Eagle medium with high glucose, L-glutamine, 110 mg/liter sodium pyruvate, and pyroxidine hydrochloride (Invitrogen Life Technologies) supplemented with 10% fetal calf serum. Cell death was quantified colorimetrically using the CytoTox96 lactate dehydrogenase (LDH) release kit according to the manufacturer's instructions (Promega, Madison, WI).
Mouse infections. Six- to eight-week-old female wild-type C57BL/6 mice (Jackson Laboratories, Bar Harbor, ME) were kept under specific-pathogen-free conditions in filter-top cages. Mice were provided with sterile water and food ad libitum. In competitive index (CI) experiments, mice were infected subcutaneously with equal amounts of the mutant strain (5 x 104 cells) and wild-type strain U112 (5 x 104 cells) for a total of 105 bacteria in a 0.05-ml volume. The mice were monitored for signs of morbidity and mortality twice daily for the entirety of the study. Spleens were harvested at 2 days postinfection and homogenized, and dilutions were plated to determine CFU/g tissue on both MH agar and MH agar supplemented with the appropriate antibiotic to select for mutant bacteria. Competitive indices were calculated as the ratio of mutant to wild type of the output, normalized for the input, and significance was calculated by comparing the log of the CI to 0. Standard errors were calculated and significance of results was determined by the application of the Mann-Whitney statistical test to the log10 value of the CI. All animal infection experiments were approved by the Institutional Animal Care and Use Committee and the Institutional Biosafety Committee of Stanford University.
Microarray accession numbers. Results of the SAM analyses are available at the Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/) under accession numbers GSE 5468.
|
|
|---|
100 genes not previously known to be regulated by MglA, suggesting that MglA regulates the expression of roughly 100 genes.
![]() View larger version (7K): [in a new window] |
FIG. 1. Wild-type and mglA mutant Francisella bacteria replicate with similar kinetics in rich media at 37°C. Shown is a representative growth curve for wild-type and mglA mutant Francisella strains. x, time points when RNA was isolated from each strain.
|
![]() View larger version (45K): [in a new window] |
FIG. 2. Transcriptional profiles of wild-type and mglA mutant Francisella strains in rich media. (A) Hierarchical gene cluster of microarrays from two independent wild-type and mglA mutant growth curves. Columns represent individual time points, increasing from left to right during the growth curve. Rows represent individual genes. Genes were clustered into groups that were expressed in lag/exponential phase or stationary phase or that were constitutively expressed. One distinct cluster was expressed in wild-type bacteria but not in the mglA mutant. (B) Expression patterns of the mglAB operon and the known MglA-regulated genes pdpD, iglA, iglC, iglD, and pdpA. Duplicate gene spots that passed the filtering criteria are included.
|
2-fold between the wild type and mutant with a false discovery rate of <1 are included in Table 1. Ten genes were expressed at higher levels in the mglA mutant than in the wild type, suggesting negative regulation of these genes by an MglA-dependent mechanism. However, the majority of the genes (92 genes) were downregulated in the mglA mutant compared to the wild type, suggesting positive regulation of these genes by an MglA-dependent mechanism (Fig. 3). These genes include the previously identified MglA-regulated genes in the FPI: pdpD, pdpA, iglA, iglC, and iglD (29). Gene FTT1091, which encodes a protein with homology to an amidase from Nostoc punctiforme, has previously been shown by proteomic analysis to be regulated by MglA (29). This gene was also transcriptionally regulated in an MglA-dependent manner in our experiments, further validating our results. We also identified many genes that lie within the FPI that were not previously known to be MglA regulated. Indeed, pdpB, pdpC, and all of the genes located between these open reading frames, CDS1 through CDS9, are also positively regulated by MglA (Table 1), showing that the entire FPI is under the control of MglA.
![]() View larger version (52K): [in a new window] |
FIG. 3. Clustering of the 102 MglA-regulated genes as determined by SAM. Seventeen spots representing 10 genes were expressed at higher levels in the mutant than in the wild type, and 183 spots representing 92 genes were downregulated in the mutant compared to the wild-type. Columns represent individual time points, increasing from left to right during the growth curve, and rows represent individual genes.
|
![]() View larger version (25K): [in a new window] |
FIG. 4. Quantification of selected MglA-regulated genes by real-time reverse transcription-PCR. Differences in transcript abundance of MglA-regulated genes were confirmed by real-time RT-PCR using RNA from the second growth curve. Gene-specific primers for (A) pdpB, (B) FTT0989, (C) FTT0612, (D) oppB, (E) FTT1209c, and (F) mglA were used to amplify either wild-type (wt) RNA, represented by a solid line, or mglA mutant RNA, represented by a dashed line. Testing of samples was performed in triplicate. Experiments were performed at least three times. Means and standard deviations from a representative experiment are shown. *, P < 0.05; **, P < 0.01; ***, P < 0.001 (compared with the wild type).
|
The F. tularensis SchuS4 genome sequence has recently been added to the Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.ad.jp/dbget-bin/www_bfind?F.tularensis). KEGG is a bioinformatics database that uses information about protein interaction networks in well-studied organisms to predict the cellular pathway(s) and function of proteins in other organisms (7, 22). The majority of the MglA-regulated genes were not assigned to a pathway (76 genes), which is consistent with the large number of hypothetical genes and pseudogenes. However, 26 genes were grouped into the categories of metabolism (20 genes), environmental information processing (3 genes), or genetic information and processing (3 genes). Within the category of metabolism, amino acid (10 genes) and carbohydrate (8 genes) metabolism genes appeared to be widely represented. For example, ivlE is involved in the biosynthesis and degradation of leucine, isoleucine, and valine, and galM is involved in glycolysis. The variety of categories identified by our analysis suggests that MglA regulates several different cellular processes during the growth conditions described here.
Construction of mutants in MglA-regulated genes. We selected eight MglA-regulated genes for mutagenesis and further study. Two genes within the FPI, pdpB and cds2, were chosen because this region has already been shown to contain virulence factors (35). We chose four genes outside the FPI based upon homology to known virulence factors as well as whether they encode either possible secretion signals or conserved catalytic domains (Table 3). These mutated MglA-regulated genes include FTT0612, which encodes a protein similar to OspD3 (SenA), a putative enterotoxin from Shigella flexneri (36); FTT0989, which encodes a protein with similarity to transglutaminase-like enzymes; FTT0123 (oppB), which encodes a protein similar to oligopeptide permease subunit B (OppB) from Pasteurella multocida; and FTT1209c, which encodes a protein with similarity to endothelin-converting enzyme (ece-1) from Shewanella putrefaciens. Lastly, two mutants were made using the hypothetical genes FTT1032 and FTT1771. Isogenic mutants were created for each gene by allelic exchange and replacement of the majority of the open reading frame with a kanamycin resistance cassette. All the mutants and complemented mutants grew similarly to the wild-type strain at 37°C in rich media (data not shown).
|
View this table: [in a new window] |
TABLE 3. Genes chosen for further studies
|
![]() View larger version (15K): [in a new window] |
FIG. 5. MglA-regulated genes contribute to Francisella replication in macrophages. Bone-marrow-derived macrophages were infected at an MOI of 25:1. Samples were taken at 2 h, 6 h, and 10 h postinfection. (A) Percent bacterial entry into macrophages was calculated by dividing the 2-h counts by the input. Macrophages were infected with (B) pdpB and cds2 mutants, (C) FTT0612 and FTT0989 mutants, and (D) oppB and FTT1209c mutants. Bacterial replication (n-fold) was calculated by dividing the 10-h bacterial counts by the averages of the 2-h bacterial counts for (E) pdpB and cds2 mutants and (F) FTT0612, FTT0989, oppB, and FTT1209c mutants. The wild type, the mglA mutant, and pdpA mutants were included as controls. Samples were performed in triplicate. Experiments were performed at least three times. Means and standard deviations from a representative experiment are shown. *, P < 0.05; **, P < 0.01; ***, P < 0.001 (compared with the wild type).
|
![]() View larger version (16K): [in a new window] |
FIG. 6. MglA-regulated genes are involved in macrophage cytotoxicity. Macrophages were prestimulated with heat-killed, wild-type Francisella cells 18 h prior to infection. Cytotoxicity was measured by LDH release at 8 h postinfection. The mglA and pdpA mutants were included as controls. Macrophages were infected at an MOI of 100:1 with (A) pdpB and cds2 mutants and (B) FTT0612, FTT0989, oppB, and FTT1209c mutants. Mutants were complemented, and cytotoxicity was measured at (C) 8 h postinfection with an MOI of 100:1; (D) 10 h postinfection with an MOI of 20:1; (E) 8 h postinfection with an MOI of 100:1; and (F) 10 h postinfection with an MOI of 10:1. Testing of samples was performed in triplicate in each experiment. Experiments were performed at least three times. Means and standard deviations from a representative experiment are shown. *, P < 0.05; **, P < 0.01; ***, P < 0.001 (compared with wild-type for A and B or mutant compared to the complemented mutant by the bar in C to F). wt, wild type.
|
The two mutants that had a higher rate of replication in macrophages than wild-type Francisella, oppB and FTT1209c (ece-1 homolog), induced significantly higher levels of cell death than the wild-type at 8 h (Fig. 6B). The oppB and FTT1209c mutants were complemented in cis, and both constructs were under the control of their native promoters and in single copy to avoid possible problems of overexpression. These complemented mutants displayed lower levels of macrophage cytotoxicity that were similar to the levels induced by wild-type bacteria (Fig. 6E and F). These data show that we have identified both attenuated and hypervirulent mutants and, together with the replication data, highlight the striking correlation between the ability of bacteria to replicate in unactivated macrophages and the induction of cytotoxicity in activated macrophages.
MglA-regulated genes contribute to Francisella virulence in mice. We identified mutations in MglA-regulated genes that were either avirulent or hypervirulent in in vitro assays and sought to further investigate the fitness of these mutants in Francisella infection in vivo. We performed competitive index experiments where the relative fitness of an individual Francisella mutant could be compared to that of the wild-type strain within a single mouse. The CI is expressed as the ratio of the number of mutant bacteria recovered from the tissue divided by that of the wild type and is normalized to the ratio of the mutant to the wild type from the input. Therefore, a competitive index value that is less than 1 indicates that the mutant is less fit than the wild-type strain in vivo. Mice were infected subcutaneously with an equal mixture of wild-type and mutant bacteria, and spleens were harvested at 2 days postinfection to determine the levels of bacterial colonization of each strain. Consistent with our in vitro results, we found that the pdpB and cds2 mutants showed dramatic decreases in fitness (CI of <1 x 106 [lower limit of detection]) (Fig. 7). In fact, the mglA mutant and all the FPI mutants tested here had similar CI values (Fig. 7). The FTT0989 (transglutaminase homolog) mutant, which showed a defect in intramacrophage replication and cytotoxicity, also showed a decrease in fitness in vivo (CI of 0.051). The FTT0612 (ospD3 homolog), FTT1032, and FTT1771 mutants, which behaved similarly to the wild-type strain in the macrophage assays, were not attenuated in vivo (CI of 0.997) (data not shown). The oppB and FTT1209c (ece-1 homolog) mutants, which replicated better than the wild-type strain in macrophages and induced higher levels of host cell death, were less fit in vivo than the wild type, with CI values of 0.19 and 0.0003, respectively. These results suggest that higher rates of replication seen in vitro are not necessarily beneficial in vivo. Taken together, these data highlight the role of MglA-regulated genes in Francisella pathogenesis in vivo.
![]() View larger version (13K): [in a new window] |
FIG. 7. MglA-regulated genes contribute to Francisella virulence in mice. Mice were infected subcutaneously with equal amounts of mutant and wild-type Francisella cells. At 2 days postinfection, spleens were taken for counts, and the competitive index was calculated. (A) Competitive indices of pdpB and cds2 FPI mutants. A spontaneous streptomycin-resistant mglA mutant (see Materials and Methods for a description) and the pdpA mutants were included as controls. (B) Competitive indices of FTT0612, FTT0989, oppB, and FTT1209c mutants. At least five mice were used to calculate the CI value for each mutant strain. Each spot represents the CI value from an individual mouse, and bars represent the geometric means. Significance was calculated by comparing the log of the CI value to 0. *, P < 0.05; **, P < 0.01; ***, P < 0.001 (compared with 0).
|
|
|
|---|
We mutated pdpB and cds2, two FPI genes that have not yet been characterized. Similar to the pdpA, pdpD, and iglC FPI mutants, the pdpB and cds2 mutants were avirulent, as determined by bacterial replication, macrophage cytotoxicity, and virulence in mice (Fig. 5, 6, and 7) (17, 28, 29, 33, 35). It is unclear how the FPI contributes to Francisella replication and virulence; however, all the FPI mutants tested here exhibited similar phenotypes in every assay. This is strongly suggestive of a system comprised of multiple components whereby mutating one component results in the inactivation of the entire system. Interestingly, PdpB has limited similarity to a domain in IcmF, a protein that may be involved in type IV secretion in Legionella pneumophila (42, 47). The FPI proteins IglA and IglB share homology with Edwardsiella tarda EvpA (22% identity) and EvpB (29% identity), respectively, which are involved in the secretion of virulence proteins (39). The fact that mutations in several FPI genes have similarly dramatic phenotypes, together with the homology of some of the FPI genes to secretion system genes in other bacteria, suggests that the FPI may encode a secretion system that is essential for Francisella intracellular replication, macrophage cytotoxicity, and virulence in vivo.
In addition to the FPI genes, we identified several novel genes that appear to be important for intracellular replication, macrophage cytotoxicity, and virulence in mice. For example, similar to FPI mutants, an FTT0989 mutant is attenuated for replication inside macrophages, does not induce macrophage death, and is outcompeted by wild-type bacteria in mice. The protein encoded by FTT0989 has a transglutaminase domain containing the catalytic triad characteristic of this domain, suggesting that the protein may have transglutaminase activity. The protein also contains a putative secretion signal and may therefore be secreted by Francisella, making it available to interact with host cells and thereby contribute to virulence. Indeed, FTT0989 plays an important role in Francisella pathogenesis, since it was attenuated for growth in macrophages and in mice (Fig. 5 and 7). Cytotoxic necrotizing factor type 1 from Escherichia coli and dermonecrotic toxin from Bordetella pertussis are two characterized transglutaminases that activate Rho GTPases of the host cells and are involved in virulence (25).
We identified two additional mutants, FTT1209c and oppB (FTT0123), that had similar phenotypes to one another in macrophages. Both of these mutants replicated to higher levels in macrophages than wild-type Francisella and induced higher levels of cytotoxicity in activated macrophages (Fig. 5 and 6). However, these mutants were attenuated in vivo (Fig. 7). FTT1209c encodes a predicted M13 metallopeptidase with striking similarity to mammalian endothelin-converting enzyme (ECE-1). FTT1209c has several homologs, including PepO from Porphyromonas gingivalis, which plays a role in host cell invasion (1, 4). FTT1209c has been annotated as a pseudogene in the sequenced F. tularensis strain SchuS4 (35). This gene does not appear to be a pseudogene in F. novicida strain U112, the strain chosen for our experiments (unpublished data). Both of the genes from the two species contain full-length N-terminal and C-terminal M13 peptidase domains as well as the HEXXH metal binding motif, suggesting that both gene products may act as metallopeptidases. However, FTT1209c in U112 has a predicted signal peptide that is absent in SchuS4. It is unknown if the signal sequence is necessary for the phenotypes described here.
FTT0123 (oppB) encodes a protein similar to OppB from Pasteurella multocida. The Opp complex has been identified as a virulence factor for several pathogens, including L. monocytogenes, group A Streptococcus, Bacillus cereus, and Bacillus thuringiensis (8, 18, 38, 48). The opp operon (oppABCDF) encodes an ATP-dependent transport system involved in the uptake of small peptides, which can serve as carbon and nitrogen sources and may have signaling properties (30). For group A Streptococcus and B. cereus, mutations in the opp operon result in an inability to express virulence factors, suggesting a regulatory role for Opp in these two species. oppB has also been annotated as a pseudogene in SchuS4 due to the presence of a nonsense mutation that results in a truncated protein. This mutation is not present in oppB from U112, where it appears to encode a full-length protein (unpublished data). Perhaps SchuS4 has evolved other ways to make it competitive in nature despite the alteration of this locus. Future experiments will be designed to investigate this possibility. Taken together, the phenotypes of the FTT1209c and OppB mutants suggest that limiting intracellular growth is a critical aspect of Francisella pathogenesis.
To our knowledge, this is the first description of Francisella mutants that replicate to higher levels than wild-type bacteria within macrophages, and these data suggest that oppB and FTT1209c may be negative regulators of Francisella growth. This suggests that Francisella encodes factors that limit its growth in macrophages and that controlling its growth rate is important for pathogenesis. Controlling replication within a macrophage may allow Francisella to grow intracellularly without triggering macrophage cell death and exposing the bacteria to complement, antibodies, etc. In addition, Francisella may use the macrophage for dissemination from the skin to the draining lymph nodes, and uncontrolled replication of Francisella could result in rapid macrophage death and elimination of the ability to disseminate efficiently.
Not all MglA-regulated genes are expected to play a role in virulence. In keeping with this assumption, we have presented data for FTT0612 (OspD3 homolog), FTT1032, and FTT1771 mutants that behaved similar to the wild type in the macrophage assays and did not contribute to virulence in a mouse model of infection. This suggests that these genes are dispensable for virulence in Francisella.
While the present work describes many new MglA-regulated genes, the mechanism of MglA regulation remains unknown. However, it is interesting that while MglA-regulated genes are maximally expressed during the stationary phase, the mglA transcript is maximally expressed in the lag and early exponential phases (our work and see reference 5). MglA is closely related to E. coli SspA, a stringent starvation protein. The SspA regulon overlaps that of another transcription factor (20), and SspA works with an additional cofactor to transcribe phase P1 late genes (19). It is possible that the MglA regulon presented here represents the combined action of MglA and an additional cofactor and/or multiple layers of regulation. Additional studies will be required to investigate this possibility.
We have identified five novel MglA-regulated virulence genes that either facilitate or limit replication within macrophages and play a role in virulence in vivo. This suggests that in addition to the bacterium's ability to grow intracellularly, the bacterium's ability to keep its own intracellular replication in check is also essential to Francisella pathogenesis. The growing number of MglA-regulated virulence factors underscores the idea that MglA is an important regulator of virulence factors. Nonetheless, the majority of MglA-regulated genes remain uncharacterized. Twenty out of 102 MglA-regulated genes are predicted to play a role in metabolism, particularly in amino acid metabolism (Table 1). Work with Listeria monocytogenes has shown that amino acid metabolism genes are essential for intracellular growth (21). This is particularly relevant, since L. monocytogenes escapes the vacuole and replicates in the cytosol of infected cells, similar to Francisella. However, the most interesting class of MglA-regulated genes may be the 43 hypothetical genes with no known or characterized homologs. Characterization of these genes may lead to the elucidation of new processes in bacterial pathogenesis that are shared with other pathogens as well as some that may be unique to Francisella.
This work was supported by grants RO1 AI063302 and AI-65359 from the NIAID, NIH, to D.M.M. A.B. was supported by a Cellular and Molecular Biology Training Program grant (NIH 5 T32 GM007276). D.S.W. was supported by a postdoctoral fellowship from the Giannini Family Foundation.
Published ahead of print on 25 September 2006. ![]()
|
|
|---|
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»