Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, March 2005, p. 931-944, Vol. 49, No. 3
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.3.931-944.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Institut für Medizinische Mikrobiologie und Hygiene, Technische Universität Dresden, Dresden,1 MWG Biotech AG, Ebersberg, Germany2
Received 13 March 2004/ Returned for modification 6 July 2004/ Accepted 29 October 2004
|
|
|---|
|
|
|---|
In 1997, the complete genome sequence of Escherichia coli K-12 strain MG1655, which consists of 4,639 Mb, was published (4). Four years later, the complete genomes of two serious bacterial pathogens, enterohemorrhagic E. coli (EHEC) O157:H7 strains EDL933 and RIMD0509952, were sequenced; and genome analysis revealed surprising results (15, 35). In comparison to E. coli K-12, their genomes were approximately 880 and 860 kb larger, respectively. For example, E. coli O157:H7 strain EDL933 contains 1,387 genes which are not present in E. coli K-12 (35). E. coli K-12 and EDL933 share a common backbone of 4.1 Mb which is predominantly colinear. The similarity between both chromosomes is discontinued by a number of so-called K islands (present in K-12 strain M1655) and O islands (present in EDL933). A total of 177 O islands ranging in size from 50 to 88 kb and comprising a total of 1.34 Mb and 234 K islands of >50 bp comprising 0.53 Mb have been detected (35). Among these O islands, 18 prophages and prophage-like elements have been identified in the EDL933 genome, including the Shiga toxin 1 (Stx1)-encoding prophage CP-933V and Stx2-encoding prophage BP-933W. Only phage BP-933W has yet been shown to produce infectious phage particles. The other phage sequences are mainly incomplete and cryptic (35). Hayashi et al. (15) also identified 18 prophages and 6 prophage-like elements in the chromosome of E. coli O157:H7 strain RIMD0509952. Phages represent a very mobile group of genetic elements which integrate into the bacterial genome at specific sites, the most prominent of which were found to be close to tRNA genes (6). Since prophage sequences comprise approximately 5% of the E. coli O157 genome, there is evidence that phages may play a role on genomic variability in E. coli (35).
From its first description as a food-borne pathogen in 1983 (40), E. coli O157:H7 has emerged as an important cause of gastrointestinal disease that has been detected throughout the world. E. coli O157:H7 is the prototype of EHEC, a subgroup of Stx-producing Escherichia coli (STEC) causing hemorrhagic colitis and the hemolytic-uremic syndrome (HUS). EHEC and STEC strains can be found among a large number of E. coli O serotypes and represent a genetically heterogeneous group of pathogens. STEC strains are now among the strains at the top of the list of bacterial causes of food-borne infections in the United States and Germany (26, 41).
Stx is the major pathogenicity factor of EHEC and is responsible for multiple damaging effects on eukaryotic cells, and the most convincing evidence that has been provided indicates that it is responsible for the clinical symptoms of bloody diarrhea and HUS (27, 38, 49).
The production of Stx by Stx-producing E. coli strains can be modulated by antibiotics (11, 18, 20, 54, 57). Karch et al. (18) have shown that subinhibitory concentrations of trimethoprim-sulfamethoxazole increased the yield of Stx produced by Shigella dysenteriae and STEC strains. In another study, E. coli O157:H7 strains expressing either one or two Stxs were treated with 13 antibiotics. It could be shown that Stx production was modulated differently in a strain-specific manner (11). Zhang et al. (57) demonstrated that ciprofloxacin induced Stx and Stx-encoding phages (Stx phages) in a mouse model, and Matsushiro et al. (25) induced Stx phages with norfloxacin. In a study by Köhler et al. (21), the simultaneous induction of Stx and Stx phages was triggered by particular growth promoters used in animal husbandry. In contrast, azithromycin showed strong in vitro bactericidal activity without inducing Stx (32).
The general location of the stx gene in the late-phase prophage region and the recA dependency of Stx induction led several investigators to suggest that Stx production is linked to the phage growth cycle. This was proven on a molecular level by Neely and Friedman (28) and Wagner and colleagues (51, 52). These findings shed new light on the controversial discussion about antibiotic treatment for HUS. Whereas some studies have reported on the beneficial effects of early antibiotic therapy, another study (21) showed that antibiotics trigger the progression to HUS.
The effects of antibiotics have hitherto been investigated only on individual Stx phages. In order to study the impact of prophages and prophage genes on the physiology of EDL933 at the transcriptional level, we investigated the influence of a low concentration of an antibiotic substance on the whole transcriptome of E. coli O157 by DNA microarray analysis, with particular emphasis on the expression of prophage genes. We have chosen norfloxacin as a model substance since it is a representative of the group of gyrase inhibitors, it is effective against Stx production, and it is used for the treatment of bacterial infections.
|
|
|---|
Growth curves. E. coli O157:H7 strain EDL933 was grown in five Erlenmeyer flasks to an optical density at 600 nm (OD600) of 0.3, as described above. Norfloxacin was added to four of the cultures to obtain final concentrations of 100, 200, 500, and 1,000 ng/ml. The fifth culture was used as a control and did not receive norfloxacin. The ODs were measured at 30-min intervals. Samples of 5 ml were taken 1 and 2 h after addition of norfloxacin. Following centrifugation at 3,740 x g for 5 min at 4°C, the supernatants were diluted in bacterial dilution buffer (Alexon Trend) and Stx production was measured by a ProspectT Shiga Toxin E. coli microplate assay (Alexon Trend), according to the recommendations of the manufacturer.
To analyze the time course of induction in more detail, tests for ODs, viable cell counts, Stx production, and plaque formation were performed in parallel with 200 ng of norfloxacin per ml. Viable cell counts were calculated after serial dilution in phosphate-buffered saline and plating on Luria-Bertani (LB) agar. Stx production was determined as described above. Concentrations of mature phage particles were determined by a double-layer plaque test, as described previously (43).
Induction experiments. EHEC O157:H7 strain EDL933 was grown overnight at 180 rpm in a rotary shaker at 37°C in 10 ml of LB broth (pH 7.5). Six 500-ml Erlenmeyer flasks, each of which contained 100 ml of LB broth, were inoculated with an overnight culture to an initial OD600 of 0.05. The bacteria were grown at 37°C at 180 rpm to an OD600 of 0.3. Subsequently, 0.8 µl of a 25-mg/ml norfloxacin (ICN) stock solution, prepared in glacial acetic acid, was added to three of the cultures (experimental cultures), giving a final concentration of 200 ng/ml. This amount of glacial acetic acid did not alter the pH of the medium. The other three cultures served as controls and were not supplemented with antibiotics. All bacterial suspensions (both experimental and control suspensions) were further incubated for 120 min at 37°C. RNA isolation was then performed at this time.
Isolation of total RNA. Total bacterial RNA was isolated from 50 ml of the EDL933 cultures on three different days by using RNeasy kits (Qiagen), according to the recommendations of the manufacturer, with minor modifications. Briefly, the bacteria were centrifuged at 3,740 x g for 5 min at 4°C. The pellets were resuspended in 1 ml of fresh lysis buffer (RLT) included in the RNeasy kit. The suspensions were transferred to Biopulverizer FastRNA-Blue tubes (Qbiogen) and mechanically treated in a Fast Prep instrument (Qbiogen) for 45 s at a speed setting of 6.5 to obtain crude extracts. After an additional centrifugation at 13,000 rpm for 3 min in a benchtop centrifuge (Eppendorf), the supernatants were processed as described in the RNeasy kit manual. To remove the DNA, we digested the sample with DNase I (Roche) and purified the RNA using the same kit. Elution of RNA from the columns was performed with diethyl pyrocarbonate (DEPC)-treated distilled water. DEPC-treated distilled water were prepared as described in the Bacteria MWG-Application Guide (http://ecom2.mwgdna.com/download/arrays/manuals/manual_bacteria.pdf). The RNA preparations from the three experiments and from the three controls were pooled. Total RNA was quantified by measuring the absorbance at 260 nm, and the quality and purity were confirmed by electrophoresis on 1.2% agarose gels containing 1.9% (vol/vol) formaldehyde. Morpholinopropanesulfonic acid (Sigma-Aldrich) was used as the electrophoresis buffer.
Synthesis and preparation of cDNA. Hybridization probes were prepared by incorporation of cyanine 3 (Cy3)-dCTP and Cy5-dCTP (NEN) into the cDNA during reverse transcription of the experimental and the control RNA, as described in the Bacteria MWG-Application Guide. For trials 1 and 2, the experimental and the control RNAs were also reverse labeled in order to control for differences in the labeling efficiencies between the fluorescence dyes.
Briefly, 6 µl of RNA solution containing 50 µg of total RNA in a volume of 6 µl DEPC-treated distilled water and 3 µl of random primer solution (3 µg/µl; Invitrogen) were mixed, incubated at 65°C for 10 min, subsequently incubated for 10 min at room temperature, and finally cooled on ice for 2 min. Two microliters of a reaction mixture containing 5 mM each dATP, dGTP, and dTTP and 2 mM dCTP (all from Invitrogen), 4 µl of fivefold-concentrated reaction buffer (Invitrogen), 2 µl of 100 mM dithiothreitol (Invitrogen), 2 µl of 1 mM Cy3-dCTP or 2 µl of 1 mM Cy5-dCTP, and 1 µl of Superscript II H (10 U; Invitrogen) was added to the mixture; and the mixture was incubated for 2 h at 42°C. To terminate the reaction and degrade the RNA, 5 µl of 1 N NaOH was added to the sample and the mixture was incubated at 65°C for 10 min. For neutralization, we added 5 µl of 1 N HCl and 200 µl of TE (Tris-EDTA) buffer (pH 7.5) to the mixture.
Labeled cDNA was then purified with a PCR purification kit (Qiagen) in accordance with the guidelines in the kit manual to remove the deoxynucleoside triphosphates that had not been incorporated, random primers, and fluorescent dyes. cDNA was eluted twice with 40 µl of nuclease-free water (Braun).
A total of 3,000 ng of Cy3-labeled cDNA from the experimental sample (Cy3 cDNAE) and 3,000 ng of Cy5-labeled cDNA from the control sample (Cy5 cDNAC), and vice versa for reverse labeling, were combined in one tube and dried in a Speed Vac apparatus (Heto Lab). The labeled cDNA was resolved in 250 µl of salt-based hybridization buffer (MWG Biotech AG).
Microarray hybridization, scanning, and analysis. The labeled cDNA was hybridized overnight with the oligonucleotides present on an E. coli O157 array (MWG Biotech AG) in a hybridization chamber (Corning). The slides were then washed three times and dried under standard conditions in accordance with the protocol of the manufacturer. For optimal utilization of spot intensities, the arrays were scanned at three different amplification settings (photomultiplier tube values of 35, 45, and 55) by using an Affymetrix 428 scanner.
After the slides were scanned, crude data from each slide were analyzed with BioDiscovery ImaGene (version 5.0) image-processing software (MWG Biotech AG). For optimal spot analysis, a grid pattern was adjusted for each slide, and awry spots (messed or fizzled spots) were tagged and excluded from the analysis. The spot intensities and local backgrounds around each spot were measured for both the Cy3 and the Cy5 channels.
Further analysis was performed with MAVI software, a proprietary product of MWG Biotech AG. This software calculates the slope for each spot and extrapolates the fluorescence intensity for faint and saturated spots scanned at three different amplification settings. Following this, the spot intensity of a specific amplification was determined by linear regression of the log signals. These values were normalized because of the different incorporation efficiencies of the Cy3 and Cy5 dyes (data not shown).
For each slide, the mean of the median background intensity values of each spot (MB) and the standard deviation (SD) of MB were calculated for both channels. Intensity cuts (Icut) were determined from these data by use of the formula MB + (2 x SD) for each channel (the Cy3- and Cy5-channels) (Table 1).
|
View this table: [in a new window] |
TABLE 1. Analysis of microarray data
|
The quotients of the fluorescence intensities of hybridized cDNA for each spot obtained by dividing the values of experiments and controls were calculated and represent a ratio. Spots with ratios of >2.0 were classified as upregulated genes, and spots with ratios of <0.5 were classified as downregulated genes. These numerical values for regulated genes have been adapted from other studies (12, 55), in which they were determined with confidence intervals >99.9%. P values were determined by the t test with Expressionist software (Genedata, Basel, Switzerland). Only data with P values <0.05 were included in the study.
E. coli O157 array. For our studies, we used the commercially available O157 array (MWG Biotech AG). The E. coli O157 array is composed of 6,176 50-mer oligonucleotides specific for 4,288 genes of E. coli K-12 strain MG1655, 5,358 genes of E. coli O157:H7 strain RIMD0509952, and 5,336 genes of E. coli O157:H7 reference strain EDL933. The array contains oligonucleotides specific for open reading frames (ORFs) of all three strains, two of the strains, or only one of them. Gene sequences with identities more than 92% (less than four mismatches) are represented by only one oligonucleotide. For example, oligonucleotide ecov2#4432 reports 18 homologous endolysin genes of different prophages (see https://ecom.mwgdna.com/services/cgx/cgx-search). On the other hand, certain genes, e.g., mosaic genes or truncated ORFs, may be specified by two oligonucleotides.
The prefixes Z (E. coli O157:H7 strain EDL933), B (E. coli K-12), and Ecs (E. coli O157:H7 strain RIMD0509952) describe oligonucleotides specific for each of the strains. In addition, 68 Arabidopsis control genes and 92 replicated spots increased the total number of spots to 6,336.
Real-time PCR. RNA isolation and cDNA synthesis for real-time PCR were performed as described above, but without cyanine dye labeling. Real-time PCR was accomplished in a final volume of 20 µl in glass capillaries with a LightCycler instrument (Roche Diagnostics). PCR was done with a FastStart DNA master kit SYBR Green I (Roche Diagnostics), according to the protocol of the manufacturer, with a final MgCl2 concentration of 4 mM. The cycling conditions were as follows: preincubation at 95°C for 10 min and 40 cycles of 95°C for 10 s, 55°C for 5 s, and 72°C for 10 s. Subsequently, a melting curve was calculated at temperatures ranging from 65 to 95°C. The specificity of the PCR for all amplicons was confirmed by agarose gel electrophoresis.
To obtain a standard curve, real-time reverse transcription (RT)-PCR was performed for serC, stxA2, stxB2, stxA1, stxB1, and recA with serial dilutions of the cDNA (Table 2). We defined a copy number for each dilution. Real-time PCR was also performed for experimental and control cDNA from E. coli O157:H7 strain EDL933 with an initial concentration of 2,800 pg of cDNA for each sample.
|
View this table: [in a new window] |
TABLE 2. Primer sequences and sizes of PCR products for real-time PCRa
|
|
|
|---|
At first, growth curve experiments (measurement of turbidity) were performed with 100, 200, 500, and 1,000 ng of norfloxacin per ml added 90 min (OD600 = 0.3) after inoculation of the culture and with a control without norfloxacin. Using norfloxacin concentrations of 500 and 1,000 ng/ml, we observed strong reductions in the OD600 30 min after the addition of the antibiotic (Fig. 1). With 100 and 200 ng of norfloxacin per ml, only a slight growth reduction was detected within 120 min, and this turned into a strong OD600 reduction after that time (Fig. 1). Using an induction time of 120 min and a norfloxacin concentration of 200 ng/ml, we noticed an increase in the level of Stx production (Fig. 2). Since we needed a marker for conditions under which particular EDL933 genes were expressed, we performed further experiments with 200 ng of norfloxacin per ml and an induction time of 120 min.
![]() View larger version (12K): [in a new window] |
FIG. 1. Semilogarithmic plot of growth curves of E. coli O157:H7 strain EDL933 in response to 100, 200, 500, and 1,000 ng of norfloxacin per ml. Control, culture grown without norfloxacin (see text for details); y axis, log OD600.
|
![]() View larger version (16K): [in a new window] |
FIG. 2. Growth curve, Stx production, and phage titer of E. coli O157:H7 strain EDL933 upon induction with 200 ng of norfloxacin per ml. N, viable cell counts; AU, absorbance units. The relative amount of Stx was measured by enzyme-linked immunosorbent assay. AUstx, OD450 of Stx determined by enzyme-linked immunosorbent assay upon induction with norfloxacin. The values for the control without norfloxacin were low (OD450 range, 1 to 30) and are therefore not included here. ODC and ODE, ODs for the control and the experimental samples, respectively; NC and NE, viable cell counts for the control and the experimental samples, respectively; pfuC and pfuE, PFU for the control and the experimental samples, respectively.
|
![]() View larger version (53K): [in a new window] |
FIG. 3. Scanning electron micrograph of E. coli EDL933 cells after 120 min of induction with 200 ng of norfloxacin per ml (A); (B) scanning electron micrograph of the control without norfloxacin. Preparation of samples was performed after fixation with 2.5% glutaraldehyde in 0.1 M phosphate-buffered saline. Microscopy was kindly performed by U. Nimmtschke and D. Streichert, Institute of Anatomy, Technical University of Dresden, and conducted with a S430 scanning microscope (LEO Electron Microscopy, Cambridge, United Kingdom). Magnification for both panels, x13,881.
|
Microarray experiments. RNA was prepared from three different norfloxacin-induced cultures (200 ng/ml, 120 min) and three cultures incubated without norfloxacin and subsequently reverse transcribed with incorporation of fluorescent cyanine dyes, as described below. The whole procedure was performed two times on two different days, and the two procedures were termed trial 2 and trial 3, respectively. For trial 1, the approach was conducted on a third day under the same conditions but without norfloxacin (Table 3). This additional trial was carried out to evaluate the expression profile under control conditions. The labeled cDNAs were hybridized with the MWG Biotech AG O157 arrays and analyzed as described below. In addition, cDNA from trials 1 and 2 was reverse labeled to control for the different incorporation efficiencies of the cyanine dyes.
|
View this table: [in a new window] |
TABLE 3. Experimental design and labeling scheme for the three trials
|
The three strains represented on the E. coli O157 array together contained 14,982 genes, which were specified by 6,167 oligonucleotides. A large part of these genes consisted of highly similar sequences, meaning that a given oligonucleotide may be specific for more than one gene. Therefore, we had to analyze a huge amount of raw data. Since we were primarily interested in EDL933 genes and to simplify the tables, we mention only genes with the prefix Z (EDL933) and omit those with the prefixes Ecs (RIMD0509952) and B (MG1655).
Due to the presence of a number of homologous phage sequences, even in the EDL933 genome, some problems arose with the assignment of spots to genes of particular phages. A reason for that is that those genes (e.g., lysis genes S, R, and Rz) consist of highly similar sequences in the different prophages present in EDL933. The oligonucleotides used for a given spot, e.g., for gene S, may be homologous to eight different S genes present in EDL933. Therefore, we used the designation "unequivocal" for the specific assignment of spots to a particular phage and "ambiguous" when the spots could be ascribed to homologous genes of distinct prophages.
In order to investigate the reproducibilities of trials 2 and 3, the ratios for slide 43 were plotted versus the ratios for slide 45 in a double-logarithmic plot (Fig. 4). The same was performed for slide 42 and slide 45 and for slide 42 and slide 43. These data demonstrate that from 84.6 to 96.7% of the ratios were in the twofold range, which indicates a high degree of reproducibility of the experiments described here (Table 4). The correlation coefficients show the variability between two experiments (slides), which ranged from 0.7 to 0.8 (Table 4).
![]() View larger version (33K): [in a new window] |
FIG. 4. Double-logarithmic plot (log-log plot) of the ratios for slide 43 versus those for slide 45. The plot shows up- and downregulated genes, indicated by ratios <0.5 and >2.0, respectively. Reproducible results were obtained when most of the ratios (genes) from different experiments are located in the twofold range (inner lines). In this range, ratios derived from different experiments do not differ by a factor greater than 2. The more that the ratios are located in the twofold range, the higher the concordance between the different trials is. The 5- and 10-fold ranges are indicated.
|
|
View this table: [in a new window] |
TABLE 4. Number of spots in twofold range and coefficients of correlation for trials 2 and 3
|
|
View this table: [in a new window] |
TABLE 5. Upregulated genes of Stx-encoding prophages BP-933W and CP-933V
|
|
View this table: [in a new window] |
TABLE 7. Upregulated prophage genes with more than two homologues sequences in the EDL933 genome
|
|
View this table: [in a new window] |
TABLE 6. Upregulated genes of non-Stx-encoding prophages present in the genome of E. coli O157:H7 strain EDL933
|
|
View this table: [in a new window] |
TABLE 8. Upregulated genes specific for E. coli O157:H7 strain RIMD0509952
|
|
View this table: [in a new window] |
TABLE 9. Upregulated genes of EDL933 not phage related
|
![]() View larger version (17K): [in a new window] |
FIG. 5. Genetic map and gene transcription profile of Stx2-converting phage BP-933W in response to 200 ng of norfloxacin per ml. The phage genome contains 71 ORFs. Each box shows a single gene. Gene sizes are drawn to scale. Gene designations were compiled from the NCBI database (accession no. AE005174). Genes with known functions are named in the respective boxes. Dark grey boxes, upregulated genes that could be unequivocally allocated to BP-933W; light grey boxes, genes that were ambiguously assigned to BP-933W; white boxes, nonregulated genes; boxes with hatched bars, genes that have not been included on the array; boxes with wavy lines, genes excluded from analysis.
|
Twenty-three upregulated genes were unequivocally assigned to the late region of BP-933W, starting downstream of gene Q and ending with the structural (head-tail) phage genes. The 99.0-fold upregulated gene Z1466 is homologous to the yihS gene of E. coli (37), the function of which is not known. The genes belonging to the lysis region, i.e., orf1467, S, R, and Rz, were upregulated but could be only ambiguously ascribed to the BP-933W genome. The tRNA genes ileZ, argN, and argO were not present on the array. Late-phase genes lom and bor appeared not to be regulated (Fig. 5), whereas the data for the late antiterminator gene Q, orf1490, orf1494, and orf1503 had to be excluded. Since bor is transcribed in the opposite direction from the other late-phase genes, this may be plausible. For lom, the situation is more difficult. The results of real-time RT-PCR with primers directed to lom were negative. Comparison of the restriction fragment length polymorphism of BP-933W derived from the EDL933 derivative used in our study and that of BP-933W from strain C600/933W revealed a slightly different pattern in this region, indicating that parts of lom have been deleted (data not shown). Since we used an E. coli EDL933 derivative that has been stored in Germany for >10 years, it may be possible that genetic rearrangements in the stab culture were responsible for this phenomenon. Such a phenomenon has been described by Papadopoulos et al. (34).
By taking a closer look at the ratios for the upregulated BP-933W genes, it is conspicuous that a number of ratios are approximately 10-fold higher than the others (Table 5). The genes with higher ratios are located mainly in the late region downstream of gene Q (Table 5). Interestingly, the ratios for the upregulated genes, which were ambiguously allocated to this region, i.e., lysis genes, were also in this range. This observation allow us to suggest that those genes could probably also be allocated to BP-933W.
The Stx1 prophage CP-933V contains 64 genes (NCBI database accession no. NC_002655). Eight upregulated genes could be unequivocally allocated to CP-933V, 19 upregulated genes appeared to be ambiguously located in the CP-933V genome, 36 genes were not regulated as determined by our approach, and 1 gene was excluded from analysis (P > 0,05). Most of the nonregulated genes occurred in the early and the beginning of the late-phase region (Fig. 6; Table 5). The upregulated genes of the lysis gene cassette appeared to be ambiguously assigned to CP-933V. In this case, however, the transcription ratios do not fit into the general ratio pattern for CP-933V, so we suggest that upregulation of the lysis gene spots was specified by BP-933W cDNA. Interestingly, the stx1 genes appeared not to be induced by norfloxacin under the conditions applied.
![]() View larger version (17K): [in a new window] |
FIG. 6. Genetic map and gene transcription profile of bacteriophage CP-933V. The map shows the gene transcription profile of CP-933V in response to 200 ng of norfloxacin per ml. Each box shows a single gene. Gene designations were compiled from a sequence from the NCBI database (accession no. AE005174). Symbols and shadings are as described in the legend to Fig. 5.
|
Upregulated genes specific for RIMD0509952. Twenty-two upregulated genes were assigned to prophages present in the genome of the Sakai outbreak strain RIMD0509952 (Table 8). We suggest that upregulation of RIMD0509952 phage-specific genes is caused by cross-hybridization of RIMD0509952 phage-specific oligonucleotides on the slide with EDL933 cDNA. One hint for this suggestion is the homology of EDL933 and RIMD0509952 phages (15). Most of the genes detected are yet uncharacterized (Table 8).
Upregulated genes outside prophages. We also found a category of regulated genes that are apparently not phage regulated (Table 9). These genes display ratios from 2.7 to 5.8 and were considered to be slightly induced. Among these were a putative phage integrase not present in any of the prophages, some stress-induced proteins, and DNA damage-inducible proteins associated with the SOS response. It is noteworthy that recA and uvrD appeared not to be upregulated under these conditions on the array.
Downregulated genes. According to the conditions applied, 122 genes were downregulated. Only two of these genes belonged to prophages. Most of the downregulated genes are responsible for metabolic functions such as energy, amino acid, fatty acid, protein, and cell wall metabolism (Table 10). We have omitted a detailed description of these genes because we turned our attention to prophage genes. Beside this basic cellular machinery, the seven locus of enterocyte effacement genes, espA, espB, espD, cesT, sepZ, sepL, and escV, were downregulated (Table 10). In addition, the genes for an EHEC factor for adherence cytotoxin (Z4332) and a Shet2-like toxin (Z4326) were also downregulated.
|
View this table: [in a new window] |
TABLE 10. Functional groups of downregulated genes of E. coli O157:H7 EDL933a
|
For the generation of standard curves, we analyzed the LightCycler runs by the fit-point method, as described in the manual accompanying the LightCycler instrument (LightCycler operator's manual, version 3.0, May 1999). Crossing points (CPs) were plotted against the input copy number to obtain standard curves and calculate the slopes with the software accompanying the LightCycler instrument. The corresponding real-time efficiency (E) of each standard curve was calculated by the equation 10(1/slope). To calculate the relative expression ratio of the regulated gene of interest (target [T]) versus that of the unregulated housekeeping gene (reference [R]) we used the formula described by Pfaffl et al. (36): E
CPT/E
CPR, where
CP is the change in CP (control CP sample CP). We provide a comparison of particular microarray data and real-time RT-PCR data in Table 11. The microarray-generated ratios for the control with the housekeeping gene serC and the stx genes are about 1, indicating nonregulated genes as expected. The means of the ratios in trial 2 and trial 3 also indicate nonregulation for serC. Array analysis has shown norfloxacin induced the stxA2 gene >157-fold but induced stxB2 only 40.7-fold. In general, the data obtained by real-time RT-PCR with the LightCycler instrument confirm these results (Table 11). However, the numerical values were higher than those obtained with the arrays, which is probably a result of the higher efficiency of the RT step (no cyanine dye incorporation) and the higher sensitivity of the real-time PCR. The same is true for stx1. Whereas the ratios for the stx1 genes were close to being below the threshold values for regulation on the array (1.8 for stxA1) and 1.5 (stxB1), real-time PCR indicated upregulation, with higher ratios of 7.0 (stxA1) and 5.9 (stxB1).
|
View this table: [in a new window] |
TABLE 11. Comparison of microarray analysis and real-time RT-PCR data
|
|
|
|---|
Little is known about the global genetic regulation of E. coli O157:H7 genes. Sperandio et al. (45) have shown by DNA microarray analysis that quorum sensing is a global regulatory mechanism for basic physiological functions as well as for the production of virulence factors in E. coli O157:H7. They found that 404 genes were upregulated by luxS at least fivefold. Moreover, an involvement of the SOS response was indicated in their study. Whereas the transcription of stx2 was 3-fold higher in strain 86-24 than in a luxS mutant, SOS genes such as recA, uvrA, and sulA were upregulated more than 20-fold. Dahan et al. (8) investigated the changes in gene expression upon binding of E. coli O157:H7 to plasma membranes. They found increased levels of stress-associated mRNA and decreased levels of mRNA encoding proteins involved in translation and type III secretion.
Bacteria are able to rapidly respond to environmental signals (24, 47). Besides signal transduction systems that forward stimuli from the outside of the cell through the cell wall to the inside, chemical substances may enter the cells and directly or indirectly affect gene transcription. Bacterial pathogens interact with a specific environment, the human host. Besides interaction with epithelial surfaces and competitive growth with the commensal flora, they come into contact with chemicals that are used for therapy. One group of chemicals used most frequently for the treatment of infectious diseases are antibiotics. Antibiotics are thought to interact specifically with bacterial metabolism and target specific structures. Targeting of these structures with therapeutic concentrations of antibacterial chemicals inhibits bacterial growth by interfering with cell wall biosynthesis, nucleotide synthesis, or protein biosynthesis. Previous studies from a number of investigators have shown that the mechanisms of antibiotic action are more complex and that subinhibitory concentrations of antibiotics may alter the bacterium's properties (2, 23). Since the use of antibiotics is common in veterinary and human medicine, as well as in animal husbandry and the food industry worldwide, bacteria frequently come into contact with antibiotics. Certain antibiotics may cause induction of bacteriophages, i.e., Stx-encoding bacteriophages (21, 57), and even subinhibitory concentrations may modulate Stx production (11, 21, 54, 56).
Norfloxacin is a member of the group of 4-quinolone antibiotics that interact with both DNA gyrase and topoisomerase IV (10). For E. coli, DNA gyrase is the primary target. Quinolones inhibit DNA synthesis by triggering the formation of stable complexes of DNA and topoisomerase IV as well as DNA gyrase. The quinolone-mediated cell death is thought to be initiated by the generation of double-stranded breaks (10). Norfloxacin is a good inducer of the SOS response, and it has been shown that genes of the SOS regulon are induced on E. coli K-12 microarrays upon incubation with norfloxacin in a dose-dependent manner (44). Those investigators used concentrations ranging from 0.03 to 8.0 µg of norfloxacin per ml for 30 min. The highest extent of gene regulation was observed with the highest norfloxacin concentration, which caused the upregulation of 22 genes involved in DNA metabolism. With a norfloxacin concentration comparable to that which we have used, only four genes of the SOS response were upregulated after 30 min of induction.
We have used 200 ng of norfloxacin per ml over an induction period of 120 min and detected by array analysis only two upregulated genes that were associated with the SOS response. One of these encodes a protein homologous to DinI of Serratia marcescens, and the other one encodes a DinD homologue (Table 9) (53). The ratios for both genes were 2.7 and 5.8, respectively, indicating moderate induction (Table 9). The real-time RT-PCR performed for determination of the expression of recA and uvrD confirmed the array data for uvrD and indicated a moderately increased level of recA expression. In general, RT-PCR is more sensitive because RT can be performed without incorporation of cyanine dyes. It is not clear why only three SOS-associated genes are switched on, but this may depend on the fact that particular DNA-damaging agents can preferentially induce particular genes of the SOS regulon (3), on the prolonged induction time, or on the low norfloxacin concentration.
The physiological response to norfloxacin probably cannot be compared between E. coli K-12 and EDL933, since EDL933 probably does not respond with a similar strength as E. coli K-12. However, 200 ng of norfloxacin per ml induced phage genes; and the genes most strongly affected were those of BP-933W in the late region, including stx2, which was induced more than 150-fold. It is intriguing that three SOS-related genes were induced. However, a low concentration of norfloxacin could probably switch on alternate DNA repair pathways. E. coli possesses multiple inducible DNA repair pathways; among these are recA-independent pathways, such as the UVM response (1, 17). However, the most plausible explanation would be that the low norfloxacin concentration is not able to fully induce the SOS regulon genes.
Regulation of gene transcription is an important response to environmental changes and is stringently regulated in bacteria. Our present understanding indicates that the coordination of gene transcription in E. coli involves a number of hierarchical levels, such as the local control of individual operons, the regional control of multiple operons within a regulon, or the regional control of multiple regulons in regulatory networks that are termed stimulons or modulons (29).
A high hierarchical level of gene regulation is represented by changes in the DNA conformation, such as DNA supercoiling or the presence of intrinsically curved DNA (14, 33). In studies with mutants it has been shown that the maximal abundance of proteins in a set of 88 proteins investigated occurred at supercoiling levels below that of the wild type, and others were most abundant with elevated levels of negative superhelices (46). Since the superhelical status is caused and influenced by the activities of DNA gyrase and topoisomerase, DNA gyrase inhibitors also influence supercoiling and gene transcription (10). It has been demonstrated that DNA gyrase inhibitors also lower the levels of DNA supercoiling. It may be hypothesized that norfloxacin at higher concentrations inhibits bacterial growth and probably modulates phage gene expression by causing changes in the DNA tertiary structure at lower levels. However, genes that have already been shown by other investigators to be regulated by supercoiling were not upregulated in our study. Future experiments will be needed to investigate this phenomenon.
It is conspicuous that the genes of prophage BP-933W are transcribed most strongly. In concordance with a temporal transcription scheme (48), genes of the late-phase region are upregulated manyfold (Table 5). However, some genes, such as int, cI, lom, bor, and some ORFs, do not fit in that scheme (Fig. 4). In order to obtain reproducible results, we used a strict system for the inclusion of spots in our analyses, as described above. For example, when the ratio for a spot on one of the three experimental slides was less than 2.0, the spot was precluded from further analyses. For the int gene, we obtained values of 4.0, 5.0, and 1.16 (data not shown). Although int appeared to be upregulated on two of the slides, it had to excluded since the ratio for one spot was less than 2.0. Therefore, from the slide data we cannot suggest exactly whether this gene is regulated or not. A similar situation was found for some other genes, such as cI. The bor gene is transcribed in the opposite direction and is not thought to be cotranscribed with the late phage genes, and lom was not present in our isolate in a fully intact form and therefore had to be excluded from the analysis.
Another interesting question to be addressed is the high level of expression of stx2 and the difference between the transcription efficiencies of stxA2 and stxB2. Up to now it has not been clear whether the relation of the A and B subunits of 1:5 is regulated on the transcriptional level or the translational level (30). By taking into account the fact that slide ratios cannot be considered absolute values, we have validated these data with real-time RT-PCR experiments. The A-subunit genes are transcribed more efficiently than the B-subunit genes. This was not expected and may be interpreted as evidence for regulation at the posttranscriptional level. This suggestion needs further research.
A number of phage genes are demonstratively regulated very strongly. These are mainly BP-933W genes. Whereas the stx2 gene is upregulated more than 150-fold, most genes of the late region displayed ratios from 10- to 99-fold. This is in contrast to CP-933V and non-Stx phages (Tables 5 to 7), which show lower regulation ratios, with a maximum of 20-fold. This points to the question of whether the different prophages possess different mechanisms of regulation of their growth cycle.
In older studies it was observed that toxin production is affected by the presence of a second toxin-encoding phage. The same could be true in the context of our study. Repressors of the BP-933W phage may alter gene transcription in CP-933V or the other prophages. Alternatively, these phages could contain other mechanisms of regulation. An interesting experiment would be investigation of whether the DNA supercoiling status of BP-933W may be responsible for that phenomenon.
Microarray analysis opens new dimensions in E. coli O157 research. The fact, however, that E. coli O157:H7 contains so many homologous DNA sequences aggravates the handling of E. coli O157:H7 arrays in this specific field of phage research. Prophage genes, present as functional or defective phage fragments, show gene transcription activity upon induction by environmental stimuli. Future research will be needed to clarify the impact of these observations for population biology and to determine regulation networks that enable the prophages to respond to and communicate with each other and alter transcription activity.
We thank Stefanie Müksch for skillful technical assistance. We are grateful to Matthias Frosch, Guido Dietrich, and Stefanie Theiss, University of Würzburg, who supported us with helpful discussions and technology transfer in the initial phase of the experiments. Andrea Hörster, MWG Biotech AG, supported us with valuable information on slide technology and the design of oligonucleotides for the real-time PCR.
|
|
|---|
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»