Previous Article | Next Article ![]()
Infection and Immunity, October 2005, p. 6791-6802, Vol. 73, No. 10
0019-9567/05/$08.00+0 doi:10.1128/IAI.73.10.6791-6802.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Microbiology & Immunology, New York Medical College, Valhalla, New York 10595
Received 6 January 2005/ Returned for modification 16 February 2005/ Accepted 16 June 2005
|
|
|---|
|
|
|---|
Molecular analysis of B. burgdorferi isolates has resulted in the differentiation of B. burgdorferi sensu lato into 10 distinct species (77). In North America, all known cases of Lyme disease are caused by infection with B. burgdorferi sensu stricto (77). Previous studies have shown that B. burgdorferi sensu stricto isolates obtained from patients with early Lyme disease can be distinguished on the basis of several genotypic parameters (77), and isolates have been classified into several rRNA gene spacer types (RST) (29, 38, 39). We have previously shown that the distributions of these genotypes differ significantly between isolates cultured from the skin of patients with early Lyme disease and those cultured from the blood of such patients (80) and that there is a highly significant association between the frequency of hematogenous dissemination in patients with early Lyme disease and the genotype of the infecting B. burgdorferi strain. These findings were confirmed in a mouse model demonstrating that B. burgdorferi RST1 clinical isolates were more invasive and produced more-severe pathology in experimentally infected animals than certain RST3 isolates did (75, 76). This suggests that the genotypic heterogeneity of B. burgdorferi clinical isolates results in differing potentials for bloodstream dissemination. This variability in virulence properties could be the result of differences in genetic content and/or alterations in gene expression between isolates. It is therefore of interest to identify genes which may be responsible for the observed differences in pathogenic potentials.
Genome array technology has proven to be a powerful tool for the assessment of bacterial gene expression under differing physiological conditions, and comparative transcriptional profiling has been employed to identify virulence genes (17, 18, 61). We recently reported the preparation and use of B. burgdorferi genome arrays to assess gene expression in B. burgdorferi B31 grown at 23 and 35°C (46, 47), and these arrays have been used for several other global transcriptional analyses of B. burgdorferi (3, 10, 44, 74, 83). In the present study, these genome arrays were employed to assess differences in gene expression between two B. burgdorferi clinical isolates with differing capacities for hematogenous dissemination.
|
|
|---|
To determine B. burgdorferi growth characteristics, wild-type and mutant strains were grown in modified BSK medium (79) without antibiotics at 34°C to densities of 1 x 108 spirochetes/ml. Dilutions (1:1,000) were inoculated into triplicate 10-ml cultures and cultivated at either 25o or 34°C. Spirochetes were enumerated daily by dark-field microscopy. These studies were repeated a minimum of two times, yielding at least five values for each time point.
RNA and DNA isolation. RNA was prepared from four separate 50-ml mid-exponential-phase cultures by using a Purescript RNA isolation kit according to the manufacturer's instructions (Gentra Systems). Genomic DNA was removed by DNase I treatment (DNA-free kit; Ambion). The final RNA was resuspended in RNase-free water, and the RNAs from the four separate cultures were pooled and quantitated. DNA was prepared from 10-ml cultures of each isolate by use of a commercial nucleic acid extraction kit (Isoquick; Orca Research).
cDNA probe synthesis.
Hybridization probes for array experiments were generated from each RNA sample by standard cDNA synthesis as described previously (46, 47). For each isolate, two cDNA synthesis reactions were performed, with the first using a mixture of all 3' B. burgdorferi open reading frame (ORF)-specific primers (Sigma-Genosys Biotechnologies, Inc.) and the second using random hexamer oligonucleotides (Promega). In each case, 5 µg of RNA was mixed with 0.4 mM each of dTTP, dGTP, and dCTP, 200 units of Superscript II reverse transcriptase (Invitrogen, Life Technologies), 10 U of RNase inhibitor, and 20 µCi of [
-33P]dATP (Amersham) in a final volume of 20 µl of 50 mM Tris-HCl (pH 8.3), 75 mM KCl, 3 mM MgCl2 (first-strand buffer). The cDNA synthesis reaction mixtures were incubated at 42°C for 2 h. Trichloroacetic acid precipitation was used to determine the labeling efficiency of each reaction, which was typically 2 x 106 to 5 x 106 cpm/µg. Unincorporated nucleotides were removed from the labeling reaction by passage through a Sephadex G-25 spin column. After purification, both labeling reactions for each isolate were pooled and used as a probe.
Genomic DNA labeling. Five µg of sheared genomic DNA isolated from strain BL206 or strain B356 was labeled radioactively by random priming using the Klenow fragment of DNA polymerase (21).
Hybridization and array data acquisition. The B. burgdorferi genome arrays have been described previously (46, 47). Prior to hybridization, a pair of the membrane arrays was placed in a prehybridization solution (ExpressHyb; Clontech) and incubated at 55°C for 1 hour. The labeled cDNA probes were denatured at 95°C for 5 minutes, added to 10 ml of ExpressHyb solution, and incubated with the membrane for 14 to 16 h at 55°C. After hybridization, the membrane was washed twice with 2x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), 0.05% sodium dodecyl sulfate (SDS) for 30 min at room temperature and then with 0.1x SSC, 0.1% SDS at 65°C for 30 min. The arrays were exposed to a phosphorimager screen (Molecular Dynamics, California) for 48 h. Array hybridizations for each isolate were performed four times, and the same membranes were probed reciprocally (i.e., by switching the membranes) after stripping in 0.4 M NaOH for 20 min at 42°C, a procedure followed by two 10-min incubations in 200 mM Tris-HCl (pH 7.0), 0.1x SSC, 0.2% SDS at room temperature.
The exposed phosphorimager screen was scanned with a pixel size of 100 µm on a Storm 840 PhosphorImager (Molecular Dynamics). The image files were analyzed using ArrayVision software (Imaging Research, Canada), and the intensities for each spot were exported into ArrayStat (Imaging Research, Canada) for statistical analysis as described in detail elsewhere (47). Briefly, the background-corrected pixel intensity values for each spot were log transformed, common error was determined, and outliers were removed. The values were then normalized by the mean across conditions by an iterative process, and a false-positive error correction was achieved by an application of the false-discovery rate method (6). The z test for four independent conditions (eight independent spots for each ORF) using the false-discovery rate method with a nominal alpha value for P being
0.01 was used as the basis for computing confidence intervals.
An alternate method for identifying differentially expressed ORFs was employed (3). ORFs for which the raw pixel values were not >2 standard deviations (SD) above the pixel values for blank spots were eliminated from further analysis. Individual pixel values for the remaining ORFs were divided by the total pixel values for all ORFs, yielding the fractional intensity for each ORF. The ratios of fractional intensity values for each ORF for the two isolates were calculated and converted to log10. ORFs with log ratio values that were greater or smaller than the mean log ratio for all ORFs by >2 SD were considered to be differentially expressed.
Standard and real-time quantitative RT-PCR. The presence or the absence of BBA74 transcript in wild-type and mutant B. burgdorferi was determined by standard reverse transcription (RT)-PCR. RNA was isolated from spirochetes grown to late log phase in 10 ml of modified BSK (wild-type) or modified BSK containing 200 µg/ml kanamycin (mutant) as described above, and total RNA was treated with DNase (AMBION DNase kit). cDNA was generated by reverse transcription. Each reaction contained 0.4 mM of each dNTP, 10 units of RNase inhibitor, 500 ng of random hexamer primers, 500 ng of total RNA, and 10 units of Superscript II reverse transcriptase (Invitrogen) in first-strand buffer. Following a 1-h incubation at 42°C, 1 µl of cDNA was used as a template to perform PCR using BBA74-specific primers (Table 1). Each RT-PCR run included controls containing no reverse transcriptase to rule out DNA contamination.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotide primers used in this study
|

Ct) was used to determine the relative expression of each gene between isolates BL206 and B356 (36).
![]() View larger version (8K): [in a new window] |
FIG. 1. BBA74 transcript levels are higher in attenuated B. burgdorferi isolates. BBA74 mRNA and flaB mRNA were quantitated for five B. burgdorferi isolates by real-time RT-PCR. Data for each isolate are the average for two separate cultures, and for each culture the RT-PCR was performed in triplicate.
|
Protein lysates (0.5 to 5 µg) were separated on 10- to 15%-gradient SDS-PAGE gels. Proteins were transferred to polyvinylidene difluoride membranes using a semidry blotting apparatus and probed with either preimmune rabbit serum (1:1,000 dilution), polyclonal anti-BBA74 rabbit serum (1:1,000 dilution), polyclonal anti-OspA/B rabbit serum (1:10,000 dilution), or anti-OspA mouse monoclonal antibody (1:2,000 dilution), each diluted in Tris-buffered saline containing 0.05% Tween 20 (TBS-T), 5% nonfat dry milk. Alkaline phosphatase-conjugated goat anti-rabbit or anti-mouse antibodies (Kirkegaard & Perry Laboratories, Inc., Gaithersburg, MD) were used as the secondary antibodies at dilutions of 1:1,000. The blots were washed three times with TBS-T to remove the unbound conjugated antibody. Visualization was achieved by treating the blots with the substrate BCIP/NBT 1 (5-bromo-4-chloro-3-indoyl-phosphate/nitroblue tetrazolium) solution (Kirkegaard & Perry Laboratories, Inc.) for 3 min at room temperature. The color reaction was stopped by washing the blots thoroughly with deionized water. Blots were dried overnight at room temperature and scanned on an Alpha Imager 2200 instrument (Alpha Innotech Corporation, San Leandro, CA), and band densitometry was performed with the provided imaging software (AlphaEase FC).
Mutagenesis of BBA74.
The mutation of BBA74 was accomplished by allelic exchange of the wild-type gene with an inactivated copy. A schematic diagram showing the relevant regions of plasmid lp54 is presented in Fig. 3A. A 2,327-bp DNA fragment spanning BBA73 to BBA76 (nucleotides 50601 to 52908 on plasmid lp54; GenBank accession number, AE000790) was amplified by PCR using primers pA73F and pA76R (Table 1) as forward and reverse primers, respectively, and cloned into the plasmid vector pGEM-T (Promega). A kanamycin resistance cassette from Streptococcus faecalis (65) was PCR amplified using primers pKNF and pKNR (Table 1) as forward and reverse primers, respectively. The recombinant pGEM vector was linearized by cleavage at a unique BstEII site (nucleotide 195 of BBA74) and treated with T4 DNA polymerase to produce blunt ends, and the kanamycin cassette was inserted. The pGEM-T ampicillin resistance locus was inactivated by the insertion of an irrelevant 248-bp DNA fragment into the unique ScaI site in the ß-lactamase gene so as to avoid the introduction of ampicillin resistance into B. burgdorferi, and the plasmid phenotype (Knr and Amps) was confirmed in Escherichia coli DH5
. Plasmid DNA was purified from host E. coli by standard methods.
![]() ![]() View larger version (44K): [in a new window] |
FIG. 3. Disruption of BBA74 by insertion of kanamycin resistance cassette. (A) Schematic diagram of mutant construct. The region of lp54 between BBA73 and BBA76 was amplified by PCR using primers p73F/p76R and cloned into pGEM-T. The insertion of the kanamycin cassette (Kn) in the unique BstEII site in BBA74 is shown. Locations of forward and reverse primers used for confirmation of the disruption are indicated. (B) Confirmation of BBA74 gene disruption. Wild-type BL206 (lanes 1, 4, and 7), mutant BBA74 (lanes 2, 5, and 8), and pGEM-T vector construct containing disrupted BBA74 (lanes 3, 6 and 9) were subjected to PCR with primers p74F/p74R (lanes 1 to 3), p74F/pKNR (lanes 4 to 6), and pKNF/p74R (lanes 7 to 9). Products were separated by electrophoresis on a 1.5% agarose gel and stained with ethidium bromide. Lane M, DNA molecular size markers.
|
Plasmid profiles. The plasmid content of B. burgdorferi mutants and transformants was determined by plasmid-specific PCR as previously described (28).
|
|
|---|
Transcript analysis and comparison between B. burgdorferi clinical isolates BL206 and B356. An underlying hypothesis of this study is that the differences in the pathogenic potentials of the two B. burgdorferi isolates are to some measure the result of differences between them in terms of gene expression. By characterizing isolate-specific patterns of gene expression, it should be possible in principle to identify genes and proteins responsible for the differential pathogenesis among isolates. As a first step in this process, the gene expression profiles of a virulent isolate (BL206) and of an attenuated strain (B356) cultivated at a temperature of 34°C were determined using B. burgdorferi genome arrays. Probe synthesis and array hybridization were performed in quadruplicate for each isolate. Since each ORF is spotted in duplicate on a given array, eight data points were analyzed for each ORF. The intensities for each spot were analyzed using ArrayVision and statistically analyzed with ArrayStat software. In a second approach, the expression of each ORF was determined as the percentage of the total expression for all ORFs, and the expression ratios for each ORF between the isolates were calculated. Expression ratios were converted to log10 values, the mean log ratio for all ORFs was determined, and those ORFs with logged expression ratios of >2 SD above or below the mean were identified. Only those ORFs with significant differences in expression ratios between isolates by both analytical methods were considered to be differentially expressed. These ORFs are listed in Table 2.
|
View this table: [in a new window] |
TABLE 2. Genes with significantly differing expression in isolates B356 and BL206
|
To independently validate the microarray data, the relative expressions of selected ORFs were determined by real-time RT-PCR. As shown in Table 2, the array results were well supported by the real-time RT-PCR findings.
Comparative expression of selected genes. (i) Linear plasmids lp54 and lp36. A striking observation is the distinctive natures of gene expression for the two isolates with respect to two sets of plasmid-encoded genes. Eleven genes, each on both plasmids lp54 and lp36, were differentially expressed. Surprisingly, all lp54-encoded genes had significantly greater expression in isolate B356. By contrast, the lp36-encoded ORFs all had higher expression levels in isolate BL206. Whereas the observed differences in transcript levels between the isolates for lp54-encoded ORFs appear to reflect true differential expression (i.e., substantial transcript is produced in both isolates), this may not be the case for lp36-encoded genes, since none of the latter genes appear to have substantial expression in isolate B356.
(ii) Possible operons. Examination of the differentially expressed ORFs revealed two regions that may be expressed as operons. The chromosomal region between BB0328 and BB0335 contains genes encoding components of an oligopeptide permease (Opp). These ABC transporter systems are composed of transmembrane components (OppB and OppC), ATP-binding proteins (OppD and OppF), and peptide-binding proteins (OppA). BB0328 to BB0330 encode binding proteins that have been designated OppA-1, OppA-2, and OppA-3, respectively; BB0332 to BB0335 encode the transmembrane and ATP-binding components of the transporter. BBB16 (on plasmid cp26) and BBA34 (on lp54) encode two additional OppA proteins. BB0329 and BB0330 showed significantly greater expression in B356 than in BL206 (expression ratios of 4.1 and 7.1, respectively). BB0328, BB0333 to BB0335, and BBB16 also showed elevated expression in B356 by both array and RT-PCR analysis, although the differences measured by array were significant by only one of the two statistical methods (data not shown). The data are consistent with earlier studies suggesting that BB0328 to BB0330 and BB0332 to BB0335 may be cotranscribed (7). Alternatively, the three OppA genes may be individually transcribed but coordinately upregulated in B356.
BB0767 to BB0769 were expressed at significantly higher levels in isolate B356 than in BL206 (expression ratios of 3.9, 5.7, and 4.5, respectively). These genes code for three seemingly unrelated proteins: MurG [UDP-N-acetylglucosamine-N-acetylmuramyl-(pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase], pyridoxal kinase, and sialoglycoprotease. They are oriented in the same direction (opposite to the gene numbering), and the coding regions are tightly packed; the end of BB0769 overlaps BB0768 by 25 nucleotides, and BB0768 and BB0767 are separated by 149 nucleotides. Interestingly, the start codon of BB0769 overlaps the termination codon of the upstream BB0770, as do the initiation and termination codons of BB0766 and BB0767, but neither BB0766 nor BB0770 is expressed at significantly different levels in the two isolates. The data suggest that BB0767 to BB0769 may be cotranscribed.
(iii) Outer surface protein genes. Several well-studied genes encoding outer surface proteins showed significantly higher transcript levels in isolate B356. These include ospC, ospD, BB0603 (encoding P66), and BBA74 (expression ratios of 6.4, 5.2, 4.3, and 27.9, respectively).
Differential expression of BBA74.
The transcript levels of BBA74 were >20-fold higher in B356 than in BL206 (Table 2), and this gene was chosen for more-careful study. To explore whether the differences in expression were due to sequence variation in either upstream or coding regions of BBA74, the region between BBA73 and BBA75 from B. burgdorferi BL206 and B356 was amplified by PCR, and the products were sequenced. The BL206 sequence was identical to that of B31 MI (GenBank accession number, AE000790). The B356 sequence showed only two nucleotide changes: A
G at position 83 and T
A at position 375 relative to the BBA74 initiation codon. Whether either of these changes contributes to the differential expression of BBA74 remains to be determined.
BBA74 transcript levels were measured by real-time RT-PCR in two additional clinical isolates with differing capacities for dissemination. Results are presented in Fig. 1. In strain B515, which disseminates in mice (76), the ratio of mRNA levels of BBA74 and flaB was 0.94; this is comparable to that observed for two other disseminating isolates, BL206 and B31 MI (ratios of 1.06 and 0.81, respectively). In contrast, the BBA74 transcript level was substantially higher (BBA74-to-flaB ratio of 7.22) in strain B331, which does not disseminate (76).
BBA74 protein levels are higher in B. burgdorferi B356. Given the >20-fold increase in the level of BBA74 transcript in isolate B356, it was of interest to determine whether the level of BBA74 protein was also elevated in this isolate relative to BL206. This determination was accomplished by semiquantitative immunoblotting using anti-BBA74 polyclonal antibodies. The densitometry of the blots shown in Fig. 2 revealed that the amount of BBA74 protein is three- to fivefold higher in B356 than in BL206. Similarly, the densitometry of BBA74 spots on two-dimensional gels also yielded a four- to fivefold-higher level for the protein in isolate B356 (data not shown). These results correlate with the increased mRNA levels in B356 as determined by array and RT-PCR analysis.
![]() View larger version (37K): [in a new window] |
FIG. 2. BBA74 protein levels are higher in B. burgdorferi B356. Increasing amounts of whole-cell lysates of strains B356 and BL206 were subjected to SDS-PAGE followed by immunoblotting with either BBA74-specific antibodies (top panel) or OspA-specific antibodies (bottom panel). Lanes 1 to 4, strain B356; lanes 5 to 8, strain BL206.
|
The inactivation of BBA74 in the mutants was confirmed by RT-PCR and Western blotting. As shown in Fig. 4, the mutants did not produce BBA74 transcript (Fig. 4A) or BBA74 protein (Fig. 4B) in any of the B. burgdorferi strain backgrounds.
![]() View larger version (32K): [in a new window] |
FIG. 4. Transcription and translation of BBA74 in mutant strains. (A) RT-PCR. Lanes 1 and 2, wild-type B31; lane 3, mutant B31; lanes 4 and 5, wild-type BL206; lane 6, mutant BL206; lanes 7 and 8, wild-type B356; lane 9, mutant B356. Reactions depicted in lanes 1, 4, and 7 lacked reverse transcriptase. Amplified products were separated by electrophoresis on a 1.5% agarose gel and stained with ethidium bromide. Lane M, DNA molecular size markers. (B) Immunoblotting. Total cell lysate (3.5 µg) was separated by 12.5% SDS-PAGE and transferred to polyvinylidene difluoride membranes, and blots were developed with either OspA/B-specific or BBA74-specific antibodies. P, wild type; M, mutant.
|
![]() View larger version (10K): [in a new window] |
FIG. 5. Growth characteristics of wild-type () and BBA74 mutant ( ) BL206 at 25 ° and 34°C. Aliquots of starter cultures were inoculated into fresh media at densities of 1 x 105 cells/ml, and growth was monitored daily by dark-field microscopy.
|
|
|
|---|
The majority of the genes with expressions that differ significantly between the two isolates studied here are plasmid encoded. Interestingly, genes on two linear plasmids, lp36 and lp54, accounted for 22 of these 57 genes. Eleven lp36-encoded genes had substantially greater expression in BL206 than in B356. Little is known regarding the 54 putative ORFs encoded on this plasmid. BBK32 was identified as a fibronectin-binding protein (53), and BBK17 is annotated as adenine deaminase (24). BBK32 and BBK50 mRNA and protein were detected in engorged ticks and infected mouse tissue, indicating that the genes are expressed in vivo (22). Both of these genes have measurable and equivalent expression levels in both isolates studied here. In contrast, BBK18, BBK34, BBK40, BBK45 to -47, BBK51, and BBK52 had at least fivefold-greater expression levels in BL206. The expression of BBK45 and BBK47 had been shown previously to be downregulated by between 2.8- and 5-fold in B. burgdorferi B31 growing in dialysis membrane chambers (DMC) that mimic the host-adapted state (10, 56). The observed differences in expression of these lp36-encoded genes in the two B. burgdorferi isolates may be the result of true differences in transcript levels. Alternatively, the differences may result from sequence variation in these genes or from their complete absence in the B356 genome. Since the array is composed of ORFs amplified by PCR from strain B31 DNA, substantial genome sequence differences between B31 and B356 would be expected to result in a variation in hybridization efficiency. Although there was only a minor difference in the global hybridization efficiencies of BL206 and B356 genomic DNAs to the array, genomic sequences located on lp36 hybridized almost sevenfold more intensely with BL206 than with B356 genomic DNA. It is, therefore, likely that the differences observed for lp36-encoded genes result from substantial genomic differences between the two isolates relative to this plasmid. The absence of one or more of these genes may contribute to the attenuated phenotype of B356.
In contrast, 11 lp54-encoded genes had significantly greater expression in isolate B356. Numerous studies have shown that the expression of lp54-encoded genes is regulated in response to physiological and environmental conditions (10, 25, 47, 56, 62). These include genes encoding the major outer surface proteins OspA and OspB and decorin-binding proteins DbpA and DbpB. Eight of the 11 genes that exhibited differential expression in this study have been shown previously to have regulated expression. Interestingly, six of these (BBA52, BBA59, BBA67 to -69, and BBA74) were shown to have reduced expression in DMC (10, 56). lp54 DNA from both BL206 and B356 hybridized equally well to the array, so the observed differences in RNA hybridization reflect a true alteration in transcript level.
Based on B31 genome annotation (24), nearly one-fourth (19/78) of the differentially expressed genes are classified as encoding cell envelope or transport proteins. This is consistent with other global analyses of B. burgdorferi transcription (10, 47, 56, 74) and with the notion that the spirochete remodels its surface in response to changes in environmental conditions (35, 52, 81). Of particular interest is the fact that the expression levels in B356 of BB0328 and BB0329, which encode OppA-2 and OppA-3, were increased by 4.1-fold and 7.1-fold, respectively. As noted above, these genes are part of an opp operon spanning from BB0328 to BB0335, and other genes in this region also displayed expression levels elevated by approximately two- to sevenfold, although this increase was significant by only one of the two statistical tests. Previous studies by Bono et al. (7) and Wang et al. (78) demonstrated that (i) BB0328 to BB0330 can be transcribed as individual mRNAs or as a polycistronic transcript; (i) only BBA34 (encoding OppA-5) is subject to temperature regulation; (iii) BB0328, BB0330, and BBA34 are relatively well expressed in vitro at 37°C, but the expression of BB0329 and BBB16 is weak; (iv) BB0329 is not appreciably expressed in ticks, and expression is moderately enhanced in mice; and (v) BB0330 transcription is increased in both ticks and mice relative to in vitro-cultivated spirochetes. Taken together, these studies suggested that different OppA alleles are subject to alternative environmental signals. In the present study, only spirochetes grown at 34°C were analyzed, and while it is clear that oligopeptide permease-related genes have significantly altered expression in isolate B356 relative to that in BL206, whether the regulation of expression under other conditions is also changed remains to be elucidated.
In bacteria, oligopeptide permeases are membrane-associated complexes belonging to the ABC transporter family. Although the major role for most bacterial oligopeptide permeases is in nutrient uptake, they may also be involved in sensing environmental signals in the form of peptides, culminating in diverse responses, such as the expression of virulence determinants, quorum sensing, chemotaxis, sporulation, conjugation, and competence (33, 51, 57). Peptide transport systems have been shown to be associated with the virulences of Salmonella enterica serotype Typhimurium (50), Streptococcus pneumoniae (19), and Enterococcus faecalis (20). Based on these observations, it has been proposed that the OppA proteins may also have a role in the adaptation of B. burgdorferi to different environmental niches, such as the tick or the mammalian host (7, 78). This family of peptide-binding proteins also may be important in regulating spirochetal chemotaxis and/or gene expression for the targeting and penetration of specific tissues in invertebrate and vertebrate hosts.
Outer membrane proteins can serve a role in solute and protein transport, and substances can be taken up by bacteria by diffusion through porin proteins embedded in the outer membrane (1). Three porin activities have been identified for B. burgdorferi. These are encoded by BB0034 (P13), BB0603 (P66 or Oms66), and BBA74 (Oms28) (49, 66, 67). The latter two genes were transcribed in isolate B356 at levels significantly greater than in BL206. P66 localizes to the outer membrane (11, 12) and has been shown to bind ß3-chain integrins (15).
BBA74 encodes a product of 257 amino acids with a putative peptidase I signal sequence (66). Skare et al. demonstrated that native BBA74 protein possesses porin activity and designated it Oms28. The protein is associated with B. burgdorferi outer membrane vesicles, although only a fraction of recombinant BBA74 expressed in Escherichia coli is associated with the outer membrane (66). BBA74 expression is regulated by temperature (47, 74), the host environment (56), and the redox state (J. Skare, personal communication). In the present study, BBA74 expressions were found to be 7- to 28-fold higher in two nondisseminating strains (B356 and B331) than in disseminating strains (BL206 and B515). Furthermore, BBA74 protein levels were also higher in B356. The differences in transcript levels are unlikely to be the result of variation in a cis-acting sequence element, since the intergenic region between BBA73 and BBA74 contained only two nucleotide differences in B356 (at positions 83 and 375). Although a role for either of these nucleotides in transcriptional regulation cannot be completely ruled out, it is more likely that differences in the functions or levels of a trans-acting factor are responsible for the observed variation in BBA74 transcription. Interestingly, the mutagenesis of the C-terminal protease, CtpA, results in a dramatic increase in BBA74 protein in strain B31 cultivated at 35°C (48). In the present study, a 3.5- to 5-fold difference in BBA74 protein levels between B356 and BL206 grown at 34°C was noted. This suggests that the levels of BBA74 protein are subject to both transcriptional and posttranslational regulation.
Mutagenesis of BBA74 did not impair the growth of B. burgdorferi in vitro at either 25° or 34°C, and there were no significant differences in growth rates or final cell densities between the mutant and the wild type. However, the BBA74 mutant displayed a 2- to 5-day growth lag at both temperatures. Given the demonstrated porin activity of BBA74 (66), it is tempting to speculate that inactivation of the protein interferes with the nutrient uptake that may be required for the entry of cells into log-phase growth. Similar transitions may be required for the adjustment to growth in different environments during the natural tick-mammal cycle. It is of interest that the mutagenesis of P66, which has also been shown to possess porin activity (67), did not result in any growth impairment in vitro (16). In the present study, all mutants lacked plasmids lp25 and lp28-1 and are, therefore, noninfectious (26, 32, 54, 55). Thus, a definitive role for this protein in vivo has yet to be established.
The present study has several potential limitations. One is the nature of the B. burgdorferi genome itself. Nearly two-thirds of the annotated B31 ORFs have no known function (24), and more than half of the differentially expressed genes identified in this investigation are in this category. To the extent that some of these ORFs may play roles in dissemination, the lack of information regarding their functions is a drawback. A further possible limitation is the use of gene arrays that are based on the genome sequence of strain B31 for analysis of other isolates whose genomes have not been sequenced. Thus, genes with possible relevance to virulence that may be present in the other isolates (e.g., BL206 and B356 in the present investigation) but not in B31 would not be subject to analysis. In principle, substantial sequence differences between B31 and the tested strains could result in reduced signals for certain ORFs. Indeed, signals were somewhat lower when arrays were hybridized with B356 genomic DNA. This problem was addressed in our experimental design by scoring the expression of each ORF as a percentage of the total expression for all ORFs.
The present study was carried out with in vitro-cultivated organisms. It may be argued that certain virulence genes may be expressed only in the infected host. As noted earlier, however, this ideal situation has rarely been achieved in microarray studies of bacterial virulence, due to the low numbers of bacteria typically present in infected tissues (60). As a result, most studies of this kind have relied on the use of in vitro cultivation conditions, e.g., low iron or pH, designed to mimic the in vivo state (9, 17, 60). Two B. burgdorferi global expression studies have been performed with organisms cultivated in dialysis membrane chambers implanted in rats (10, 56). Unfortunately, a number of methodological differences between these two investigations resulted in a low concordance between the two data sets. Interestingly, nearly 30% of the ORFs listed in Table 2 (including BB0603, BBA74, ospC, and ospD) were identified as having significantly altered expression in DMC in at least one of these studies (10, 56). Narasimhan et al. have used a preamplification approach to study B. burgdorferi gene expression in feeding ticks (45) and in heart or central nervous system of nonhuman primates (44). Methodological differences between those reports and the present study make comparisons difficult. One of these studies utilized membrane arrays the same as those employed here. It is worthy of note that expression was detected for only 470 of the ORFs spotted on the membrane (44), suggesting that this preamplification method may alter the mRNA distribution or underestimate the levels of mRNAs with low expression.
Only two clinical isolates, a virulent strain (BL206) and an attenuated strain (B356), were analyzed by transcriptional profiling in this study. Numerous studies have demonstrated substantial genotypic variability among B. burgdorferi isolates (23, 39, 40, 77). Thus, the extent to which these two isolates are representative of all virulent and attenuated strains is unclear, although increased expression of BBA74 was observed in the two attenuated strains studied thus far (B356 and B331). BL206 belongs to a monophyletic clade based on rRNA gene and ospC sequence analysis, and all isolates of this genotype that have been tested in vivo produce severe disseminated infections in mice (76). Thus, it is reasonable to expect that BL206 is representative of other RST1 isolates. The question of whether the observed differences in global gene expressions between BL206 and B356 are reflective of all attenuated isolates characterized to date, which are more genotypically diverse (76, 80), warrants further study. Regardless, this study was designed to explore differences in gene expression for isolates that differ in their potentials for hematogenous dissemination, as is the case for BL206 and B356.
Despite these limitations, several of the present findings are of interest and may provide insights into B. burgdorferi pathogenesis. Evolution of a pathogenic lineage from a mildly pathogenic or commensal ancestor most frequently results from the acquisition of a virulence gene or genes. However, examples of enhanced pathogenicity arising from the loss of a gene or a change in its function have been described for Pseudomonas aeruginosa, E. coli, Shigella spp., and other bacterial pathogens (8, 41, 43, 68). Such alterations have been called pathoadaptive mutations (68). This may be the case for the disseminating lineage of B. burgdorferi. Differences in expression levels of lp54- and lp36-encoded genes in disseminating and nondisseminating strains resulting from either differential transcription or genomic variation were noted. Either mechanism would produce spirochetes with different protein contents. Additionally, a number of genes with potential involvement in nutrient uptake (BB0603, BBA74, BB0329, BB0330, and BBB29) have significantly higher expression in isolate B356. Moreover, nearly 25% of the differentially expressed genes are predicted to be localized on the cell surface. This latter observation implies that these two isolates have considerably different cell surface properties. Mouse infection with isolate BL206 resulted in dissemination to all peripheral tissues tested and a relatively high pathogen burden. In contrast, infection with B356 yielded no disseminated infection, and spirochetes could not be recovered from the inoculation site within 24 h after infection. Similar results were obtained with immunodeficient SCID mice (75, 76). This suggests that isolate B356 is rapidly cleared by components of the innate immune system. The recognition of pathogen-associated molecular patterns plays a central role in the innate immune response to infection (30, 42). Thus, the distinctive natures of the surfaces of BL206 and B356 may result in differential recognition by the innate immune system, leading to the observed differences in spirochete dissemination.
The ORFs listed in Table 2 were scored as having significantly different expression levels in the two tested isolates. These differences can be the result of either alterations in the transcriptional regulation of these genes or differences between the two isolates in terms of genetic content. In either case, these 78 ORFs are a good starting point for the investigation of factors involved in the hematogenous dissemination of B. burgdorferi. Ultimately, the possible roles of these putative virulence genes must be confirmed by mutagenesis and in vivo studies in the rodent-tick infectious cycle. The inactivation of BBA74 is a first step in this process, and more extensive studies with this and other genes are currently under way.
This work was supported by NIH grant AI45801 and the G. Harold and Leila Y. Mathers Charitable Foundation.
|
|
|---|
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»