This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tucker, D. L.
Right arrow Articles by Fox, G. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tucker, D. L.
Right arrow Articles by Fox, G. E.
Agricola
Right arrow Articles by Tucker, D. L.
Right arrow Articles by Fox, G. E.

 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

Effect of an Artificial RNA Marker on Gene Expression in Escherichia coli

Don L. Tucker,1 Fathi Karouia,1 Jim Wang,2 Yi Luo,2 Tong-Bin Li,2 Richard C. Willson,1,3 Yuriy Fofanov,2 and George E. Fox1,3*

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


arrow
ABSTRACT
 
Transcriptional analysis was used to examine the effect of a genomically encoded artificial RNA on Escherichia coli in rich and minimal media. Only the expression of a single gene, deoC, was unequivocally affected under both conditions. E. coli marker strains of this type may be useful in monitoring the fate and transport of bacteria in various applications.


arrow
INTRODUCTION
 
A previously constructed Escherichia coli tracking strain known as PCPHR expresses a genomic gene encoding an artificial RNA (aRNA) which contains a unique 17-nucleotide identifier sequence (2). The aRNA gene was originally derived from Vibrio proteolyticus 5S rRNA by deleting a key functional region and replacing it with the identifier sequence (21). The aRNA product is expressed at high levels from an rRNA promoter and accumulates in the cell without an obvious effect on growth rate (2). Although it resembles 5S rRNA, the aRNA is inactive and not integrated into the E. coli ribosomes (2, 11, 14, 21). The high accumulation level of the aRNA allows its presence to be readily detected by many routine molecular methods, including RNA size, hybridization, fluorescence, and immunohybridization (15, 16, 22). Large numbers of alternative but equally unique identifier sequences could be introduced into the aRNA (12), thereby making it possible to simultaneously differentiate between multiple sources of E. coli. These aRNA strains are, however, genetically modified organisms whose use in field applications, such as source tracking or bioremediation, would require regulatory approval.

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 this window]
[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

When rich medium is considered alone (Table 1), the fim genes (fimI and fimA) also met all three criteria for increased gene expression in PCPHR compared to the wild-type strain EMG2 described in this report. When gene transcription in rich medium is analyzed using only the Student t test (P value, <0.05), all genes of the fimA operon (fimAICDFGH) were more highly transcribed in PCPHR (data not shown). FimI is involved in type 1 pili biosynthesis, while a repeating monomer of FimA, arrayed helically to form a hollow core fiber, serves as a major component of type 1 pili (20, 26). Type I pili are filamentous bacterial surface attachment organelles that have been associated with virulence (20). The reason for the increased fimA operon gene transcription is unclear.

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.


arrow
ACKNOWLEDGMENTS
 
This work was supported in part by grants from the Gulf Coast Hazardous Waste Research Center (081UHH2793) and Environmental Institute of Houston to G.E.F., the Institute of Space Systems Operations to D.L.T., G.E.F., and R.C.W., and the Texas Learning and Computational Center to Y.F., R.C.W., and B. Montgomery Pettitt.


arrow
FOOTNOTES
 
* Corresponding author. Mailing address: Dept. Biology and Biochemistry, University of Houston, Houston, TX 77204-5001. Phone: (713) 743-8363. Fax: (713) 743-8351. E-mail: fox{at}uh.edu. Back


arrow
REFERENCES
 
    1
  1. Ahmad, S. I., and R. H. Pritchard. 1969. A map of four genes specifying enzymes involved in catabolism of nucleosides and deoxynucleosides in Escherichia coli. Mol. Gen. Genet. 104:351-359.[CrossRef][Medline]
  2. 2
  3. Ammons, D., J. Rampersad, and G. E. Fox. 1998. A genomically modified marker strain of Escherichia coli. Curr. Microbiol. 37:341-346.[CrossRef][Medline]
  4. 3
  5. Ammons, D., J. Rampersad, and G. E. Fox. 1999. 5S rRNA gene deletions cause an unexpectedly high fitness loss in Escherichia coli. Nucleic Acids Res. 27:637-642.[Abstract/Free Full Text]
  6. 4
  7. Amouyal, M., L. Mortensen, H. Buc, and K. Hammer. 1989. Single and double loop formation when deoR repressor binds to its natural operator sites. Cell 58:545-551.[CrossRef][Medline]
  8. 5
  9. Arfin, S. M., A. D. Long, E. T. Ito, L. Tolleri, M. M. Riehle, E. S. Paegle, and G. W. Hatfield. 2000. Global gene expression profiling in Escherichia coli K12. The effects of integration host factor. J. Biol. Chem. 275:29672-29684.[Abstract/Free Full Text]
  10. 6
  11. Arnold, C. N., J. McElhanon, A. Lee, R. Leonhart, and D. A. Siegle. 2001. Global analysis of Escherichia coli gene expression during the acetate-induced acid tolerance response. J. Bacteriol. 183:2178-2186.[Abstract/Free Full Text]
  12. 7
  13. Barbosa, T. M., and S. B. Levy. 2000. Differential expression of over 60 chromosomal genes in Escherichia coli by constitutive expression of MarA. J. Bacteriol. 182:3467-3474.[Abstract/Free Full Text]
  14. 8
  15. Buxton, R. S., H. Albrechtsen, and K. Hammer-Jespersen. 1977. Overlapping transcriptional units in the deo operon of Escherichia coli K-12. J. Mol. Biol. 114:287-300.[CrossRef][Medline]
  16. 9
  17. Condon, C., C. Squires, and C. L. Squires. 1995. Control of rRNA transcription in Escherichia coli. Microbiol. Rev. 59:623-645.[Abstract/Free Full Text]
  18. 10
  19. Conway, T., and G. K. Schoolnik. 2003. Microarray expression profiling: capturing a genome-wide portrait of the transcriptome. Mol. Microbiol. 47:879-889.[CrossRef][Medline]
  20. 11
  21. D'Souza, L. M., R. C. Willson, and G. E. Fox. 2000. Expression of marker RNAs in Pseudomonas putida. Curr. Microbiol. 40:91-95.[CrossRef][Medline]
  22. 12
  23. D'Souza, L. M., M. Larios-Sanz, R. A. Setterquist, R. C. Willson, and G. E. Fox. 2003. Small RNA sequences are readily stabilized by inclusion in a carrier rRNA. Biotech. Prog. 19:734-738.[CrossRef][Medline]
  24. 13
  25. Hammer-Jespersen, K., A. Munch-Petersen, P. Nygaard, and M. Schwartz. 1971. Induction of enzymes involved in the catabolism of deoxyribonucleosides and ribonucleosides in Escherichia coli K 12. Eur. J. Biochem. 19:533-538.[Medline]
  26. 14
  27. Hedenstierna, K. O. F., H. Y. Lee, Y. Yang, and G. E. Fox. 1993. A prototype stable RNA identification cassette for monitoring plasmids of genetically engineered microorganisms. Syst. Appl. Microbiol. 16:280-286.[Medline]
  28. 15
  29. Kourentzi, K. D., G. E. Fox, and R. C. Willson. 2001. Microbial detection with low molecular weight RNA. Curr. Microbiol. 43:444-447.[CrossRef][Medline]
  30. 16
  31. Kourentzi, K. D., G. E. Fox, and R. C. Willson. 2002. Microbial identification by immunohybridization assay of artificial RNA labels. J. Microbiol. Methods 49:301-306.[CrossRef][Medline]
  32. 17
  33. Mori, H. 2004. From the sequence to cell modeling: comprehensive functional genomics in Escherichia coli. J. Biochem. Mol. Biol. 37:83-92.[Medline]
  34. 18
  35. Mortensen, L., G. Dandanell, and K. Hammer. 1989. Purification and characterization of the deoR repressor of Escherichia coli. EMBO J. 8:325-331.[Medline]
  36. 19
  37. Neidhardt, F. C., P. L. Bloch, and D. F. Smith. 1974. Culture medium for enterobacteria. J. Bacteriol. 119:736-747.[Abstract/Free Full Text]
  38. 20
  39. Orndorff, P. E., A. Devapali, S. Palestrant, A. Wyse, M. L. Everett, R. R. Bollinger, and W. Parker. 2004. Immunoglobulin-mediated agglutination of and biofilm formation by Escherichia coli K-12 require the type 1 pilus fiber. Infect. Immun. 72:1929-1938.[Abstract/Free Full Text]
  40. 21
  41. Pitulle, C., K. O. F. Hedenstierna, and G. E. Fox. 1995. A novel approach for monitoring genetically engineered microorganisms by using artificial, stable RNAs. Appl. Environ. Microbiol. 61:3661-3666.[Abstract]
  42. 22
  43. Pitulle, C., L. Dsouza, and G. E. Fox. 1997. A low molecular weight artificial RNA of unique size with multiple probe target regions. Syst. Appl. Microbiol. 20:133-136.[Medline]
  44. 23
  45. Richmond, C. S., J. D. Glasner, R. Mau, H. Jin, and F. R. Blattner. 1999. Genome-wide expression profiling in Escherichia coli K-12. Nucleic Acids Res. 27:3821-3835.[Abstract/Free Full Text]
  46. 24
  47. Tao, H., C. Bausch, C. Richmond, F. R. Blattner, and T. Conway. 1999. Functional genomics: expression analysis of Escherichia coli growing on minimal and rich media. J. Bacteriol. 181:6425-6440.[Abstract/Free Full Text]
  48. 25
  49. Tucker, D. L., N. Tucker, and T. Conway. 2002. Gene expression profiling of the pH response in Escherichia coli. J. Bacteriol. 184:6551-6558.[Abstract/Free Full Text]
  50. 26
  51. Valenski, M. L., S. L. Harris, P. A. Spears, J. R. Horton, and P. E. Orndorff. 2003. The product of the fimI gene is necessary for Escherichia coli type 1 pilus biosynthesis. J. Bacteriol. 185:5007-5011.[Abstract/Free Full Text]
  52. 27
  53. Valentin-Hansen, P., K. Hammer, J. E. Love Larsen, and I. Svendsen. 1984. The internal regulated promoter of the deo operon of Escherichia coli K-12. Nucleic Acids Res. 12:5211-5224.[Abstract/Free Full Text]


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.





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tucker, D. L.
Right arrow Articles by Fox, G. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tucker, D. L.
Right arrow Articles by Fox, G. E.
Agricola
Right arrow Articles by Tucker, D. L.
Right arrow Articles by Fox, G. E.