Previous Article | Next Article ![]()
Applied and Environmental Microbiology, July 2005, p. 4156-4159, Vol. 71, No. 7
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.7.4156-4159.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
| SHORT REPORT |
Department of Biology and Biochemistry,1 Department of Computer Science,2 Department of Chemical Engineering, University of Houston, Houston, Texas 77204-50013
Received 25 August 2004/ Accepted 4 January 2005
|
|
|---|
|
|
|---|
The issue, which was addressed here, is whether or not the presence of the aRNA tracking gene will significantly modify the behavior of the cells carrying it. E. coli has eight genes for 5S rRNA that are found in seven operons. The synthesis of the rRNAs is strongly linked with cell growth and subject to sophisticated regulation (9). For example, it has been shown that if the levels of 16S and 23S rRNA are artificially reduced by deletions in one or more rRNA operons, the cell will up-regulate its remaining rRNA operons to compensate for the loss (9). In contrast, the successive deletion of 5S rRNA genes revealed that there is apparently no ability to directly compensate for underexpression of 5S rRNA (3). These results suggest that the addition of an artificial RNA to the genome would not have a significant impact on other cellular processes. Nevertheless, the safety and utility of tracking strains expressing aRNAs that resemble 5S rRNA will depend on whether their presence affects other components of the cellular machinery in a significant or unexpected way.
In order to address this issue, arrays of PCR products targeting each known open reading frame in the E. coli genome were used to determine what transcriptional changes are induced in the PCPHR strain compared to the wild-type (WT) E. coli EMG2 strain. Such arrays have been widely employed in the study of E. coli gene expression with considerable success (5-7, 17, 23-25). These experiments have generated extensive lists of genes that respond to various experimental variables (10).
Cultures of the WT EMG2 and PCPHR strains were grown aerobically in either rich YT medium (8 g Bacto tryptone, 5 g yeast extract, 5 g NaCl per liter) or defined minimal MOPS (morpholinepropanesulfonic acid) medium supplemented with glucose (19). Total RNA was harvested from the mid-log phase of growth. RNA samples from both strains were isolated twice from YT medium and three times from M9 minimal medium. Hybridizations were conducted using matched pairs of Panorama E. coli gene arrays (Sigma-Genosys Biotechnologies Inc., Houston, TX). These arrays were spotted with PCR products in duplicate by the manufacturer. They contain 4,290 PCR products coinciding with all known E. coli MG1655 open reading frames. The arrays were rinsed, prehybridized (65°C for 1 to 4 h), hybridized (65°C for 16 to 18 h), washed, and wrapped in Saran Wrap as described previously (24) with minor modifications (25). Exposed phosphor screen images (Fuji BAS-IP MS phosphor screen from FUJIFILM Medical Systems USA Inc., Stafford, TX) were visualized at a pixel density of 50 µm (40,000 dots/cm2) with a Storm 860 PhosphorImager (Molecular Dynamics, Sunnyvale, CA). The array membranes used in these experiments were subsequently stripped of hybridized probe by boiling inverted arrays in a microwave (100°C) in 200 ml of stripping solution (10 mM Tris, pH 7.5, 1 mM EDTA, 1% sodium dodecyl sulfate) and reused as many as eight times. Two hybridization probes were created from each RNA sample for hybridization to each array in an array pair resulting in four hybridizations for each strain grown in rich YT medium and six hybridizations each in M9 minimal medium.
A software package, ImageAnalyzer, was developed to quantify the hybridization signals. ImageAnalyzer utilizes specific metrics that take into consideration the nonuniformity of pixels inside a spot. These metrics exclude the highest-intensity pixels, which may reflect spike noises, and lowest-intensity pixels resulting from the surface tension effect of the probe printing process. The results obtained with ImageAnalyzer were compared to results obtained with the demonstration version of the commercially available ArrayVision Software (Imaging Research Inc., St. Catharines, Ontario, Canada), and ImageAnalyzer was found to produce spot pixel density data of similar or better quality (data not shown). Analysis of variance and the Student t test were performed between the groups of arrays obtained for each probe/gene, and the genes were sorted according to their P values, which represent the level of significance in differential expression, and according to the changes (n-fold) in gene expression between samples. The spot finding and statistical analysis software can be obtained upon request from Yuriy Fofanov (yfofanov@bioinfo.uh.edu). A descriptive user's manual is available at http://www.bioinfo.uh.edu/ImageAnalyzer/.
Genes were determined to have significant changes in gene expression based on three previously identified criteria: (i) an overall P value of <0.05, which implies a 95% probability that a change in expression between strains or media was significant (5); (ii) an overall ratio of gene expression of >2.0 between strains or media (7); and (iii) a gene expression ratio of >2.0 in each of the individual comparisons. Combining statistical methods in this manner increases the probability that genes remaining after statistical analysis are in fact changing significantly under the experimental conditions (23, 25). The raw image data and results of the analysis of the expression studies are available at http://www.bioinfo.uh.edu/publications/aRNA_marker/.
Repeated growth studies revealed at most minor differences in lag phase, growth rate, and final cell density between WT EMG2 and PCPHR. The average doubling times from four high-quality growth curves in minimal media were 75 min in EMG2 and 73 min in PCPHR with standard deviations of 13 and 12 min, respectively. In accordance with previously published data (2), the average doubling times based on six high-quality growth curves in rich media for EMG2 and PCPHR were 25.1 and 27.7 min with respective standard deviations of 1.7 and 2.6 min. As expected, rich-medium cultures grew faster and reached a higher final cell density than cultures grown in M9 minimal medium. These data suggest that the stable aRNA present in PCPHR does not appreciably affect growth rate or energy consumption in either rich or minimal medium.
The array studies revealed that the expression level of only one gene, deoC, which encodes a 2-deoxyribose-5-phosphate aldolase, in the PCPHR strain was incontrovertibly elevated in both minimal and rich media. To ensure that all putative aRNA-impacted genes in minimal media had been identified, the 10 genes with the greatest overall changes in gene expression ratios between EMG2 and PCPHR were also examined (Table 1). The role that many of these genes could play in response to the stable aRNA is unclear, but it is noteworthy that the genes deoB (which encodes phosphopentomutase) and deoD (which encodes purine-nucleoside phosphorylase) were present in this less-stringent analysis. Along with deoA, whose expression was not elevated, and deoC, these genes belong to the deo operon, which is involved in deoxyribonucleoside and ribonucleoside catabolism (1, 4, 7, 8, 13, 18, 27). It therefore seems likely that the primary response to the presence of the aRNA gene in the PCPHR strain is an increase in ribonucleotide catabolism as the accumulating aRNA is gradually degraded. However, whereas deoB and deoD operate on ribonucleotides and deoxyribonucleotides, prior work suggests that the deoC gene product is involved only in deoxynucleoside catabolism (7, 13), and hence, its up-regulation was not actually expected. The PCPHR strain may have significant unexpected value in future work aimed at fully understanding the regulation and function of the deo operon.
|
View this table: [in a new window] |
TABLE 1. Genes expressed in E. coli PCPHR and WT EMG2 in M9 (minimal) and YT (rich) media, sorted by fold change in gene expression
|
None of the genes whose expression was elevated in EMG2 (i.e., expression decreased in PCPHR) passed the original >2.0-fold-change test for significance. Decreasing the stringency of the change (n-fold) in gene expression of YT medium-grown cultures from >2.0 to >1.5 revealed six additional genes with increased expression in PCPHR and eight genes with increased expression in EMG2. This slight increase in expression of some YT medium-grown PCPHR genes could be caused by the medium composition or the presence of additional false positives in the smaller YT medium data set. Thus, in conclusion, gene expression in the PCPHR strain is extremely similar to that in the WT EMG2 strain lacking the aRNA identifier in both minimal and rich media.
Quantitative methods of tracking specific microbial populations are essential to basic research on complex microbial ecosystems. Such capability is also needed in practical applications, such as source tracking and microbial monitoring in bioremediation. For example, the release of naturally occurring or genetically engineered organisms as part of a bioremediation effort raises considerable public concern regarding the possibility of organisms escaping to nearby ecosystems. Groups participating in such activities will require incontrovertible evidence as to where their organisms did or did not end up. Appropriately designed tracking strains possibly based on aRNA identifiers can meet these requirements.
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»