5-azacytidine induces transcriptome changes in Escherichia coli via DNA methylation-dependent and DNA methylation-independent mechanisms

Background Escherichia coli K-12 strains contain DNA cytosine methyltransferase (Dcm), which generates 5-methylcytosine at 5′CCWGG3′ sites. Although the role of 5-methylcytosine in eukaryotic gene expression is relatively well described, the role of 5-methylcytosine in bacterial gene expression is largely unknown. Results To identify genes that are controlled by 5-methylcytosine in E. coli, we compared the transcriptomes of cells grown in the absence and presence of the DNA methylation inhibitor 5-azacytidine. We observed expression changes for 63 genes. The majority of the gene expression changes occurred at early stationary phase and were up-regulations. To identify gene expression changes due to a loss of DNA methylation, we compared the expression of selected genes in a wild-type and dcm knockout strain via reverse transcription quantitative PCR. Conclusions Our data indicate that 5-azacytidine can influence gene expression by at least two distinct mechanisms: DNA methylation loss and a mechanism that is independent of DNA methylation loss. In addition, we have identified new targets of 5-methylcytosine-mediated regulation of gene expression. In summary, our data indicate that 5-azacytidine impacts the composition of the bacterial transcriptome, and the primary effect is increased gene expression at early stationary phase. Electronic supplementary material The online version of this article (doi:10.1186/s12866-016-0741-4) contains supplementary material, which is available to authorized users.


Background
The modified DNA base 5-methylcytosine (5-MeC) plays an important role in transcriptional regulation in higher eukaryotes. The presence of 5-MeC in eukaryotic promoters and CpG islands is generally repressive for transcription, whereas 5-MeC in gene bodies is positively correlated with transcription [1]. Bacteria, such as E. coli, contain 5-MeC [2]. In E. coli K-12 strains, the only known cytosine-5 DNA methyltransferase is DNA cytosine methyltransferase (Dcm) [3,4]. Dcm methylates the second cytosine in 5′CCWGG3′ sequences [3]. The dcm gene is in an operon with the vsr gene which codes for a protein that repairs T:G mismatches caused by deamination of 5-MeC [5][6][7]. The original function elucidated for Dcm was in restriction enzyme biology where Dcm promotes the loss of plasmids containing the EcoRII restriction enzyme gene (which cleaves 5′CCWGG3′ sites) and protects cells from post-segregational killing by the EcoRII restriction enzyme [8,9]. In addition, Dcm protects phage lambda against DNA cleavage when EcoRII is introduced into the cell [10]. However, Dcm is a solitary methyltransferase without a cognate restriction enzyme in K-12 cells. Other roles for Dcm are certainly possible.
Based on the important role of 5-MeC in eukaryotic transcription and the fact that there is little known about the relationship between 5-MeC and gene expression in bacteria, Dcm has been recently evaluated for an impact on the composition of the E. coli transcriptome. Our group has demonstrated that two ribosomal protein genes and the drug resistance transporter gene sugE are upregulated in the absence of the dcm gene at early stationary phase via reverse-transcription quantitative PCR (RT-qPCR) [11,12]. Kahramanoglou et al. demonstrated that there are gene expression changes in dcm knockout cells using DNA microarrays, and most changes are at stationary phase [13]. Taken together, these data suggest that Dcm influences the transcriptome. As the only known function of Dcm is cytosine DNA methylation, the simplest model is that Dcm mediates gene expression changes via the generation of 5-MeC. It is noteworthy that some DNA methyltransferases can methylate tRNA and influence gene expression via a DNA-methylation independent mechanism [14][15][16].
In order to test the model that Dcm-mediated cytosine DNA methylation directly influences gene expression in E. coli and identify new genes impacted by DNA methylation, we analyzed the E. coli transcriptome in the absence and presence of 5-azacytidine (5-azaC) treatment. 5-azaC is a nucleoside analog that is used clinically to treat myelodysplastic syndromes [17]. 5-azaC is phosphorylated upon cell entry and incorporated into both RNA and DNA [18,19]. When 5-azaC is incorporated into DNA, cytosine-5 DNA methyltransferases become covalently trapped on the DNA and are degraded, and this limits the amount of enzyme available for the generation of 5-MeC [18,19]. Thus, 5-azaC is a cytosine DNA methylation inhibitor. It is important to note that 5-azaC has effects on the cell beyond blocking DNA methylation. For example, 5-azaC can induce the SOS response [20,21], induce DNA mutations [22], block translation [23], and block RNA methylation [24]. Thus, the physiology of 5-azaC treated cells is not identical to cells lacking cytosine DNA methyltransferases. Although 5-azaC has been routinely used to demethylate DNA in a variety of eukaryotes to assess the consequences of cytosine DNA methylation loss [25,26], this is the first report of the response of the entire transcriptome to 5-azaC in a bacterial organism.

Effects of 5-azaC on global DNA methylation levels
First, we determined the concentration dependence of DNA methylation inhibition by 5-azaC using digestion of DNA with the restriction enzyme isoschizomers BstNI and PspGI (Fig. 1). Both enzymes cut DNA at Dcm recognition sites (5′CCWGG3′), but PspGI is blocked by Dcm-mediated methylation of the second cytosine. In the absence of 5-azaC, DNA from early stationary phase cells was largely resistant to PspGI indicating that the DNA is heavily methylated at this stage. At early logarithmic stage, DNA was slightly sensitive to PspGI, indicating that most but not all 5′CCWGG3′ sites are methylated. These data are consistent with our current Fig. 1 The effect of 5-azaC on global DNA methylation levels in Escherichia coli. E. coli wild-type strain BW25113 was grown in LB at 37°C in the absence of 5-azaC, the presence of 0.5 μg/mL 5-azaC, 5 μg/mL 5-azaC, and 50 μg/mL 5-azaC. At two hours (logarithmic phase, L) and eight hours (early stationary phase, S), genomic DNA was isolated. DNA was incubated in the absence of restriction enzyme, presence of BstNI (B) (not blocked by Dcm methylation), or presence of PspGI (P) (blocked by Dcm methylation). DNA samples were analyzed by agarose gel electrophoresis and ethidium bromide staining results ( Fig. 3) indicating that there is 100 % methylation of a single osmE 5′CCWGG3′ site analyzed by sodium bisulfite sequencing from early stationary phase DNA and 93.3 % methylation of a single osmE 5′CCWGG3′ site analyzed by sodium bisulfite sequencing from logarithmic phase DNA. The data are also consistent with previously published results indicating stationary phase DNA is heavily methylated, whereas logarithmic phase DNA is methylated but has some unmethylated sites [13,27]. It is currently unclear if the lack of 100 % methylation at logarithmic phase is due to active DNA replication where newly synthesized DNA is unmethylated or regulation of DNA methylation. When DNA was isolated from cells incubated in the presence of 5-azaC, the DNA was more sensitive to digestion with PspGI, indicating 5-azaC inhibits DNA methylation at all concentrations used at both phases. There was an increase in PspGI sensitivity as the concentration moved from 0.5 μg/mL 5-azaC to 5 μg/mL 5-azaC, especially in early stationary phase DNA.

Effects of 5-azaC on bacterial growth
We conducted experiments to determine the effects of 5-azaC on the kinetics of E. coli growth (Fig. 2). Bacterial cells were grown in the absence of 5-azaC and three concentrations of the drug as described previously. There was little effect of 5-azaC on bacterial growth for the first few hours. At two hours, which represents logarithmic phase and one time point for the DNA microarray analysis, ANOVA indicated that there is not a statistically significant difference in A 600 readings between the four cultures (p = 8.3E-2). This is likely due to the fact that 5-azaC has to be transported into the cells, phosphorylated, and incorporated into nucleic acids before it is active. After 8 h of growth in 5-azaC (early stationary phase), ANOVA indicated there was a statistically significant difference in A 600 readings (p = 9.9E-5). Post hoc analysis indicated that there was no difference in A 600 readings between untreated cultures and cultures grown in 0.5 μg/mL 5-azaC (p = 1.4E-1). However, there was a statistically significant difference between untreated cultures and cultures grown in 5 μg/ mL 5-azaC (p = 4.7E-3) and 50 μg/mL 5-azaC (p = 7.6E-5). In summary, the data are consistent with previous experiments indicating that high concentrations of 5-azaC can inhibit bacterial growth [28][29][30]. Based on the PspGI digestion data and the growth curve data above, we chose 5 μg/mL as our working concentration of 5-azaC. At 5 μg/mL, there is an obvious reduction of DNA methylation, and the reduction is greater than that Fig. 2 The effect of 5-azaC on Escherichia coli growth in liquid media. E. coli wild-type strain BW25113 was grown in LB at 37°C in the absence of 5-azaC (black squares), the presence of 0.5 μg/mL 5-azaC (red circles), 5 μg/mL 5-azaC (green triangles) and 50 μg/mL 5-azaC (blue diamonds) for eight hours. At every hour, each culture was analyzed by spectrophotometry at 600 nm. The data are from three biological replicates (n = 3). Error bars representing standard deviation are not shown to maintain clarity of the four lines. Brackets indicate time points of RNA isolation for microarray analysis, and data analyzed by ANOVA. The asterisk indicates a p-value of <0.05 observed using 0.5 μg/mL 5-azaC. With respect to growth defects, the growth defect with 5 μg/mL 5-azaC is less than that observed using 50 μg/mL 5-azaC. A severe growth defect could compromise the interpretation of downstream DNA microarray experiments.
Effects of 5-azaC on site specific methylation The PspGI assay described above to monitor DNA methylation is primarily qualitative. To better quantify the ability of 5-azaC to block DNA methylation in Escherichia coli, we analyzed the osmE promoter region for 5-MeC from cells grown in the absence and presence of 5 μg/mL 5-azaC using sodium bisulfite sequencing. Treatment of DNA with sodium bisulfite causes the deamination of cytosines to uracils, but 5-MeCs are not deaminated and remain as cytosines in DNA sequencing reactions [31]. OsmE was chosen as it contains one 5′CCWGG3′ site in its promoter. The sodium bisulfite sequencing data are shown in Fig. 3. At the second cytosine in the single 5′CCWGG3′ site in the DNA analyzed, 93.3 % of the logarithmic phase cytosines and 100 % of the early stationary phase cytosines were resistant to bisulfite conversion. At the remaining 34 non-5′CCWGG3′ sites, only 3.3 % of the logarithmic phase cytosines and 2.6 % of the early stationary phase cytosines were resistant to bisulfite conversion. The data are consistent with the model that Dcm methylates 5′CCWGG3′ sites. Low levels of non-5′CCWGG3′ methylation are most easily explained by a very low level of non-conversion of cytosines to uracil, which is always observed in bisulfite sequencing reactions. However, it is not possible to rule out low levels of non-5′CCWGG3′ methylation. In the presence of 5-azaC, the level of methylation at the 5′CCWGG3′ site was reduced from 93.3 % to 28.6 % at logarithmic phase and from 100 % to 50 % at early stationary phase. The potency of 5-azaC at early logarithmic phase is likely due to better incorporation of the drug into DNA during DNA replication at this stage, yet there is inhibition at early stationary phase as well. In summary, the data demonstrate that 5 μg/mL 5-azaC blocks, but does not eliminate, DNA methylation at 5′CCWGG3′ sites in E. coli.

Transcriptome changes in 5-azaC treated cells
We compared E. coli transcriptomes in the absence and presence of 5-azaC using DNA microarray analysis of logarithmic phase and early stationary phase RNAs (Additional files 1 and 2). Overall, we detected six upregulated genes and one down-regulated gene at logarithmic phase (Table 1), and 44 up-regulated genes and 12 down-regulated genes at early stationary phase ( Table 2). To confirm the validity of the microarray experiments, we analyzed expression of the recA and rsmI genes in 5-azaC treated cells via RT-qPCR (Table 3). The RT-qPCR experiments confirmed the direction of the changes observed in the microarray experiments; the RT-qPCR fold-changes were slightly higher.
At logarithmic phase, the most obvious changes are genes that respond to DNA damage such as recN, dinD, dinB, and dinG. DAVID analysis indicates that the DNA repair and SOS response functional categories are overrepresented in this list (Table 4). We used a moderately stringent cut off for gene expression changes and there are other up-regulated SOS genes that have a p value <0.05, but did not make the FDR cutoff. One example is recA, which was shown to be up-regulated by RT-qPCR (Table 3). To determine if the gene expression changes are due to a direct effect of 5-azaC or loss of DNA methylation, we analyzed expression of recA and rsmI in a dcm knockout strain. RecA expression did not change in the dcm knockout strain indicating 5-azaC-dependent recA gene expression changes are not due to a loss of DNA methylation. However, expression of rsmI, a 16S rRNA methyltransferase, was increased in the dcm knockout strain, indicating rsmI expression is influenced by Dcm.
We further inspected the data for expression of DNA repair pathway genes in the presence of 5-azaC. The need for DNA repair after 5-azaC treatment is predicted to be due to the formation of DNA-Dcm crosslinks, as Dcm overexpression strains display increased sensitivity to 5-azaC induced killing [21,28,29,32]. In addition, DNA repair may result in synthesis of new unmethylated DNA and explain methylation-dependent gene expression changes. At logarithmic phase, there is little evidence for 5-azaC induced expression of the mismatch repair genes mutS, mutH, and mutL, nucleotide excision repair pathway genes uvrA, uvrB, and uvrC, and pyrimidine base excision repair genes ung, nth, mug, and mutM (Table 1). We cannot rule out induction under conditions that were not evaluated. We did observe upregulation of one gene that functions in translesion synthesis (dinB). The significance of increased dinB expression is unclear. Translesion synthesis is not thought to be required for 5-azaC induced DNA-protein crosslink repair as translesion mutants do not display increased sensitivity to the drug [21,30,32,33]. Interestingly, three of the seven genes in Table 1 have a defined or predicted role in homologous recombination (recN, rmuC, dinG). Homologous recombination is thought to be the main pathway required for repair of 5-azaC induced damage as homologous recombination mutants are more sensitive to 5-azaC than wild-type strains [21,28,30,[32][33][34]. Thus, our data are consistent with a model where 5-azaC induced damage is repaired by homologous recombination, and components of this pathway are upregulated in 5-azaC treated cells. At early stationary phase, there are 56 differentially expressed genes. Interestingly, the early stationary phase gene list is completely different from the logarithmic phase list (Tables 1 and 2). DAVID analysis identified functional categories ribosome and organelle inner membrane as overrepresented in the early stationary phase list (Table 4). Ribosomal protein genes have previously been reported as up-regulated in the dcm knockout strain [12,13] and the organelle inner membrane list contains sugE, which has previously been reported to be up-regulated in the absence of dcm [11]. To identify new target genes impacted by cytosine DNA methylation, we analyzed a few targets in the dcm knockout strain via RT-qPCR (Table 3). We investigated osmE via RT-qPCR as it has a 5′CCWGG3′ site in its promoter (Fig. 3). OsmE expression decreases in the dcm knockout strain at early stationary phase, indicating Dcm influences osmE expression. OsmE was also down-regulated in the dcm knockout strain at stationary phase via RT-qPCR. Interestingly, osmE is the first example of a gene that is dependent upon Dcm for maximal expression. The expression of atpH, a subunit of ATP synthase, did not change in the dcm knockout strain, indicating the 5-azaC-dependent gene expression change is not due to a loss of DNA methylation. FimA, a fimbria subunit, was upregulated in the dcm knockout strain indicating control by Dcm. The expression of the cspB gene, a cold shock protein, increased in the dcm knockout strain, but not to the extent observed with 5-azaC treatment. Thus, cspB may respond to both the loss of DNA methylation and a 5-azaC-dependent DNA methylationindependent mechanism.

Discussion
At logarithmic phase, the dcm gene is expressed as detected by microarray analysis and the DNA is methylated at 5′CCWGG3′ sites ( Figs. 1 and 3). Yet there are very few 5-azaC-induced gene expression changes in the wild-type strain. The majority of 5-azaC changes are in DNA repair protein genes and are likely due to DNA damage and activation of the SOS response rather than a loss of DNA methylation (Table 4). 5-azaC-dependent activation of the SOS response has been reported [20,21], and the pathway promotes survival in the presence of 5-azaC [20,21,28,30,[32][33][34]. Although it is possible that we are missing dcm-dependent changes due to a moderately stringent statistical cutoff and/or the lack of 100 % methylation inhibition via 5-azaC, a simple model based on our data and three recent articles is that Dcm has only a minor effect on the logarithmic stage transcriptome [11][12][13]. It is noteworthy that the effect of Dcm on the transcriptome has only been tested under standard conditions (rich medium, 37°C, aerobic, laboratory strain), and future work will evaluate other conditions.
At early stationary phase, the dcm gene is expressed as detected by microarray analysis, the DNA is methylated at 5′CCWGG3′ sites, and there are numerous 5-azaCdependent gene expression changes. Not all of the genes induced by 5-azaC at stationary phase are predicted to  Fig. 3 5-azaC reduces cytosine DNA methylation of the osmE promoter in Escherichia coli. E. coli wild-type strain BW25113 was grown in LB at 37°C in the absence of 5-azaC and presence of 5 μg/mL 5-azaC. DNA was isolated after two hours (logarithmic phase) and eight hours (early stationary phase). DNA samples were treated with sodium bisulfite. The osmE locus was amplified by PCR using sodium bisulfite-treated DNA as a template. PCR products were inserted into the plasmid pGEM-T Easy, and analyzed by Sanger sequencing.  be due to a loss of DNA methylation and future work will elucidate the mechanism of these gene expression changes not described in this article. A simple model based on our work and three recent studies is that Dcm has a significant impact on the early stationary phase transcriptome [11][12][13]. The mechanism of Dcmdependent control of gene expression is not known and remains an enigma. Kahramanoglou et al. provided evidence that Dcm controls expression of the stationary phase specific RNA polymerase subunit rpoS gene and a loss of Dcm generates an increase in RpoS and RpoSdependent gene products [13]. However, we did not observe up-regulation of rpoS expression in 5-azaC treated cells or dcm knockout cells via DNA microarrays and RT-qPCR (Tables 1, 2 and 3). What could explain the differences in rpoS expression in the absence of the dcm gene in the two studies? It is possible that genetic differences between the MG1655 strain used in the Kahramanoglou et al. study and the BW25113 strain used in this study could explain the differences in rpoS expression in the absence of dcm. The genome sequence of BW25113 was released in 2014 [35]. BLASTN searches using the MG1655 rpoS gene and 1000 basepair upstream region as a query indicated the same region is  100 % identical in BW25113 (data not shown). We also sequenced the main rpoS promoter responsible for stationary phase specific expression [36,37], and there are no sequence differences between MG1655 and BW25113 (Additional file 3). We also sequenced four rpoS transcriptional start sites identified by Mendoza-Vargas et al. [38], and found no genetic differences between the two strains (Additional file 3). In summary, we have no evidence for genetic differences that potentially influence rpoS expression in MG1655 and BW25113. We certainly cannot rule out genetic changes outside of the rpoS loci in the two strains, epigenetic differences in the two strains, or posttranscriptional regulation of rpoS expression via Dcm. Therefore, our data are not consistent with a model of Dcm-mediated regulation via RpoS and we must consider alternative models. What other molecules are stationary phase specific that could be effector molecules? One possibility is small RNAs, as some small RNAs are up-regulated in E. coli during stationary phase [39], and small RNAs are known to guide DNA methylation in plants [40]. Also, two of the best described Dcm-influenced genes, sugE and osmE, contain 5′CCWGG′3 sites near the -10 region of the promoter. Therefore, a model where Dcm directly influences gene expression by direct methylation of target genes still must be evaluated. DNA methylation could influence gene expression via numerous mechanisms including the alteration of DNA structure and influencing transcription factor binding. With respect to osmE, Regulon DB indicates the transcription factors FIS and IHF bind to the osmE promoter, but the binding sites do not overlap the 5′CCWGG3′ site [41]. It is important to note that~11 % of E. coli promoters contain 5′CCWGG3′ sites [12]. Yet, most of these genes are not changing in response to DNA methylation loss, and therefore promoter methylation is not likely a common mechanism to control gene expression.
Since we observed gene expression changes in the absence of Dcm, another question concerns the biological consequences of the gene expression changes. As the vast majority of gene expression changes occur at early stationary phase, it is possible that Dcm has a role in stationary phase biology. During stationary phase of cells grown in liquid culture,~99 % of the cells die within 7-10 days [42]. Subsequently, cells persist for long periods of time where waves of mutations can generate cells that can outcompete their predecessors (growth advantage at stationary phase, GASP). We are currently working to determine if Dcm promotes or reduces GASP.

Conclusions
In this report, the E. coli transcriptomes of cells grown in the absence and presence of the cytosine DNA methylation inhibitor 5-azaC were compared. 5-azaC was found to be an effective DNA demethylating agent in E. coli. 5-azaC induced expression changes in 63 E. coli genes. The majority of the changes occurred at early stationary phase and were up-regulations. There are at least two mechanisms by which 5-azaC can induce gene expression changes. The first pathway is a DNA methylation-independent pathway that likely involves a DNA damage response. The second pathway is a DNA methylation-dependent pathway which occurs in the absence of changes in rpoS expression. The precise mechanism by which cytosine DNA methylation influences gene expression will be the focus of future studies.

DNA methylation analysis via PspGI digestion
DNA was isolated from 2 mL of bacterial cells grown to logarithmic phase and early stationary phase using the Qiagen DNeasy kit. DNA quality and quantity was measured using a NanoDrop 1000 spectrophotometer. To assess methylation at 5′CCWGG3′ sites, DNA (0.5-1 μg) was incubated with 10 units of the restriction enzymes BstNI or PspGI for 2 h at 60°C [12]. Reactions were monitored by electrophoresis on 1 % agarose gels and stained with 0.5 μg/mL ethidium bromide or 1x GelRed (VWR).

Sodium bisulfite sequencing
1 μg of DNA was treated with sodium bisulfite according to the instructions in the EpiTect Kit (Qiagen). Bisulfite-treated DNA was used as a template for PCR. The conditions were 94°C for 2 min, 94°C for 1 min, 55°C for 1 min, and 72°C for 1 min for 30 cycles, followed by one cycle of 72°C for 10 min. Primer sequences are osmE-F1-D 5′GAAAAGATAAAATTTT TTTAAAGTTAATT3′ and osmE-R1-D5′ACACT-CAAAATTCCTACCATATTCTTATT3′. PCR products were inserted into the pGEM-T Easy plasmid (Promega), introduced into E. coli JM109, isolated using the alkaline lysis procedure (Wizard SV kit, Promega), and analyzed by Sanger sequencing (GeneWiz, New Jersey).

Total RNA Isolation
At logarithmic phase, early stationary phase, and stationary phase, RNA was harvested from 4 mL of cells using the MasterPure RNA isolation kit according to the manufacturer's instructions (EpiCentre). The optional DNase step was included in all preparations. RNA was analyzed by spectrophotometry on a NanoDrop 1000 spectrophotometer. For microarray experiments, RNA quality was assessed by bioanalysis on a RNA 6000 Nano chip at the University of Rochester Genomics Center. RIN values were between 8 and 10. For RT-qPCR experiments, RNA was analyzed by bioanalysis or traditional agarose gel electrophoresis after heating the RNA for 2 min at 65°C in a solution containing 47.5 % formamide and 0.01 % sodium dodecyl sulfate.

DNA Microarray analysis
Microarray experiments were performed with RNA from untreated cells and cells treated with 5 μg/mL 5-azaC. Total RNA was treated with the reagents in the mRNA Only Prokaryotic mRNA Isolation Kit to reduce rRNA levels and polyadenylate the mRNA according to the manufacturer's instructions (EpiCentre). Modified RNA was used as input for array experiments. cRNA was made in the presence of either Cy3-CTP (RNA from untreated cells) or Cy5-CTP (RNA from 5-azaC-treated cells) using the Quick Amp 2color Labeling Kit (Agilent) according to the manufacturer's instructions. Equal amounts of dye-labeled cRNA were hybridized to E. coli 8*15 K DNA microarrays arrays (Agilent technologies). Microarrays were scanned on an Agilent Microarray Scanner. There were five biological replicates for both the logarithmic phase and early stationary phase samples. Data for non E. coli K-12 genomes on the array were removed prior to analysis. P values were determined with one-sample t-tests. The data were considered significant if the fold-change was > or < 2 fold and the Benjamini-Hochberg false discovery rate (FDR) was < 0.1. Pathway analysis was performed using DAVID with the high stringency setting and the same FDR cutoff [44].

RT-qPCR analysis
RT-qPCR was performed as described previously [11]. Fold-changes were determined using the ΔΔCT method [45]. Three or more independent samples were used for each experiment for all targets except fimA (n = 2). The data were considered significant if the fold change was > or < 2-fold and the p value was <0.05 in one sample ttests. Primer sequences are listed below.

Statistical analysis
All statistical tests were performed in R [46].

Additional files
Additional file 1: Microarray fold-changes, p values, and false discovery rates (FDR) for all E. coli genes from logarithmic phase 5-azaC treated cells. Fold-change is expression in the presence of 5-azaC divided by expression in the absence of 5-azaC and is log 2 transformed. P values are from one sample t-tests, and the FDR is the Benjamini-Hochberg false discovery rate. (CSV 153 kb) Additional file 2: Microarray fold-changes, p values, and false discovery rates (FDR) for all E. coli genes from early stationary phase 5-azaC treated cells. Fold-change is expression in the presence of 5-azaC divided by expression in the absence of 5-azaC and is log 2 transformed. P values are from one sample t-tests, and the FDR is the Benjamini-Hochberg false discovery rate. (CSV 153 kb) Additional file 3: Sanger DNA sequencing analysis of the BW25113 rpoS promoter and transcription start sites. (PDF 509 kb) Abbreviations 5-azaC, 5-azacytidine; RT-qPCR, reverse transcription quantitative PCR