Skip to main content

Advertisement

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

Article metrics

Abstract

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.

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 [57]. 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 [1416].

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.

Results

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 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.

Fig. 1
figure1

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

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 A600 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 A600 readings (p = 9.9E-5). Post hoc analysis indicated that there was no difference in A600 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 [2830]. 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 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.

Fig. 2
figure2

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

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.

Fig. 3
figure3

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. a) logarithmic phase DNA in untreated cells, b) logarithmic phase DNA in 5-azaC treated cells, c) early stationary phase DNA in untreated cells, and d) early stationary phase DNA in 5-azaC treated cells. The numbers represent the position of the cytosines with respect to the transcription start site (+1). The cytosines in the 5′CCWGG3′ site are at positions −13 and −12. Open boxes represent Cs in the osmE sequence (Ts after bisulfite sequencing), and closed boxes represent 5-MeCs in the osmE sequence (Cs after bisulfite sequencing). Rows represent different clones. Xs represent ambiguous sequence information

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 up-regulated 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.

Table 1 Gene expression changes in logarithmic phase, 5-azacytidine-treated Escherichia coli cells
Table 2 Gene expression changes in early stationary phase, 5-azacytidine-treated Escherichia coli cells
Table 3 RT-qPCR analysis of gene expression changes in 5-azacytidine treated and dcm knockout Escherichia coli

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.

Table 4 Overrepresented functional categories in 5-azacytidine-treated E. coli

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 up-regulation 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, 3234]. 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 methylation-independent 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, 3234]. 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 [1113]. 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-azaC-dependent gene expression changes. Not all of the genes induced by 5-azaC at stationary phase are predicted to 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 [1113]. The mechanism of Dcm-dependent 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 RpoS-dependent 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.

Methods

Bacterial strains and growth

E. coli K-12 wild-type strain BW25113 and dcm knockout strain JW1944-2 have been previously described [12, 43]. The BW25113 (wild-type strain) genotype is F, Δ(araD-araB)567, ΔlacZ4787(::rrnB-3), λ-, rph-1, Δ(rhaD-rhaB)568, hsdR415. The JW1994-2 (dcm knockout strain) genotype is F, Δ(araD-araB)567, ΔlacZ4787(::rrnB-3), λ-, Δdcm-735::kan, rph-1, Δ(rhaD-rhaB)568, hsdR415. Cells were grown in LB at 37 °C with shaking at 250 RPM to logarithmic phase (~2 h, A600 of ~0.45), early stationary phase (~8 h, A600 of ~2.7), and stationary phase (~24 h, A600 of 3.8). 5-azaC (Sigma) was dissolved fresh in 1x PBS at 5 mg/ml, sterilized by filtration through a 0.22 μM filter, and added to experimental samples at 0.5–50 μg/mL. The same volume of 1x PBS was added to control flasks.

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′GAAAAGATAAAATTTTTTTAAAGTTAATT3′ and osmE-R1-D5′ACACTCAAAATTCCTACCATATTCTTATT3′. 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 2-color 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 t-tests. Primer sequences are listed below.

  • atpH-F, 5′TGCCGTCGAACACCAAAGTGTAGA3′

  • atpH-R, 5′CAGTTGCTCACCACAAACTGCGAT3′

  • cspB-F, 5′CGGCTTTATTTCTCCTGTTGATG3′

  • cspB-R, 5′GGACCTTTAGCACCACTCTC3′

  • fimA-F, 5′ACTCTGGCAATCGTTGTTCT3′

  • fimA-R, 5′TGAACGGTCCCACCATTAAC3′

  • mdh-F, 5′AAGTTGAAGTGCCGGTTATTGGCG3′

  • mdh-R, 5′ATCAGCCACTTCCTGCTCGGTAAA3′

  • osmE-F, 5′TGCGGGTAAACCTTCGTCTGAAGT3′

  • osmE-R, 5′TTACTTCGCAGCCTGTGGATCAGT3′

  • recA-F, 5′CGAAAGCGGAAATCGAAGGCGAAA3′

  • recA-R, 5′TAGTGGTTTCCGGGTTACCGAACA3′

  • rpoS-F, 5′CAGCCGTATGCTTCGTCTTA3′

  • rpoS-R, 5′CGTCATCTTGCGTGGTATCT3′

  • rsmI-F, 5′AGGGCCAGCTTTACATTGTACCGA3′

  • rsmI-R, 5′TATCCTCGGCGGCAATCAGATCAA3′

DNA sequencing analysis of the rpoS promoter

PCR was used to amplify the main stationary phase specific rpoS promoter [36, 37] and four rpoS transcription start sites [38]. E. coli BW25113 genomic DNA was used as a template and the Promega GoTaq master mix was used as a source of enzyme, dNTPs, and buffer. The conditions were 94 °C for 2 min, 94 °C for 1 min, 54 °C for 1 min, and 72 °C for 45 s for 30 cycles, followed by one cycle of 72 °C for 10 min. The PCR products were purified using the Qiagen QIAquick PCR purification kit and analyzed by Sanger DNA sequencing of both strands (GeneWiz, New Jersey). Three forward/reverse primer pairs were utilized; rpoSp-F1/rpoSp-R1, rpoSp-F2/rpoSp-R2, rpoSp-F3/rpoSp-R3. The primers sequences are listed below.

  • rpoSp-F1, 5′ACCTTGCAGGTGGGTAATG3′

  • rpoSp-R1, 5′GAGATAGGCGTACTGGTTGATG3′

  • rpoSp-F2, 5′GCTACTCAACAACCGCAAATC3′

  • rpoSp-R2, 5′CGTAGGCACTCAGGTAATCATC3′

  • rpoSp-F3, 5′CAGCAAAGGACAGGCAATTATC3′

  • rpoSp-R3, 5′CCTCAACTCCGTTCTCATCAA3′

Statistical analysis

All statistical tests were performed in R [46].

Abbreviations

5-azaC, 5-azacytidine; RT-qPCR, reverse transcription quantitative PCR

References

  1. 1.

    Zemach A, McDaniel IE, Silva P, Zilberman D. Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010;328:916–9.

  2. 2.

    Sanchez-Romero MA, Cota I, Casadesus J. DNA methylation in bacteria: from the methyl group to the methylome. Curr Opin Microbiol. 2015;25:9–16.

  3. 3.

    Marinus M, Lobner-Olesen A. DNA Methylation. In: Böck A, Curtiss III R, Kaper J, Karp P, Neidhardt F, Nyström T, et al., editors. EcoSal - Escherichia coli and Salmonella: Cellular and Molecular Biology. Washington: ASM Press; 2009. p. 1–88.

  4. 4.

    Marinus MG, Morris NR. Isolation of deoxyribonucleic acid methylase mutants of Escherichia coli K-12. J Bacteriol. 1973;114:1143–50.

  5. 5.

    Hennecke F, Kolmar H, Brundl K, Fritz HJ. The vsr gene product of E. coli K-12 is a strand- and sequence-specific DNA mismatch endonuclease. Nature. 1991;353:776–8.

  6. 6.

    Robertson AB, Matson SW. Reconstitution of the very short patch repair pathway from Escherichia coli. J Biol Chem. 2012;287:32953–66.

  7. 7.

    Sohail A, Lieb M, Dar M, Bhagwat AS. A gene required for very short patch repair in Escherichia coli is adjacent to the DNA cytosine methylase gene. J Bacteriol. 1990;172:4214–21.

  8. 8.

    Ohno S, Handa N, Watanabe-Matsui M, Takahashi N, Kobayashi I. Maintenance forced by a restriction-modification system can be modulated by a region in its modification enzyme not essential for methyltransferase activity. J Bacteriol. 2008;190:2039–49.

  9. 9.

    Takahashi N, Naito Y, Handa N, Kobayashi I. A DNA methyltransferase can protect the genome from postdisturbance attack by a restriction-modification gene complex. J Bacteriol. 2002;184:6100–8.

  10. 10.

    Hattman S, Schlagman S, Cousens L. Isolation of a mutant of Escherichia coli defective in cytosine-specific deoxyribonucleic acid methylase activity and in partial protection of bacteriophage lambda against restriction by cells containing the N-3 drug-resistance factor. J Bacteriol. 1973;115:1103–7.

  11. 11.

    Militello KT, Mandarano AH, Varechtchouk O, Simon RD. Cytosine DNA methylation influences drug resistance in Escherichia coli through increased sugE expression. FEMS Microbiol Lett. 2014;350:100–6.

  12. 12.

    Militello KT, Simon RD, Qureshi M, Maines R, VanHorne ML, Hennick SM, et al. Conservation of Dcm-mediated cytosine DNA methylation in Escherichia coli. FEMS Microbiol Lett. 2012;328:78–85.

  13. 13.

    Kahramanoglou C, Prieto AI, Khedkar S, Haase B, Gupta A, Benes V, et al. Genomics of DNA cytosine methylation in Escherichia coli reveals its role in stationary phase transcription. Nat Commun. 2012;3:886.

  14. 14.

    Goll MG, Kirpekar F, Maggert KA, Yoder JA, Hsieh CL, Zhang X, et al. Methylation of tRNAAsp by the DNA methyltransferase homolog Dnmt2. Science. 2006;311:395–8.

  15. 15.

    Schaefer M, Pollex T, Hanna K, Tuorto F, Meusburger M, Helm M, et al. RNA methylation by Dnmt2 protects transfer RNAs against stress-induced cleavage. Genes Dev. 2010;24:1590–5.

  16. 16.

    Som S, Friedman S. Regulation of EcoRII methyltransferase: effect of mutations on gene expression and in vitro binding to the promoter region. Nucleic Acids Res. 1994;22:5347–53.

  17. 17.

    Kaminskas E, Farrell AT, Wang YC, Sridhara R, Pazdur R. FDA drug approval summary: azacitidine (5-azacytidine, Vidaza) for injectable suspension. Oncologist. 2005;10:176–82.

  18. 18.

    Stresemann C, Lyko F. Modes of action of the DNA methyltransferase inhibitors azacytidine and decitabine. Int J Cancer. 2008;123:8–13.

  19. 19.

    Yang X, Lay F, Han H, Jones PA. Targeting DNA methylation for epigenetic therapy. Trends Pharmacol Sci. 2010;31:536–46.

  20. 20.

    Barbe J, Gibert I, Guerrero R. 5-Azacytidine: survival and induction of the SOS response in Escherichia coli K-12. Mutat Res. 1986;166:9–16.

  21. 21.

    Lal D, Som S, Friedman S. Survival and mutagenic effects of 5-azacytidine in Escherichia coli. Mutat Res. 1988;193:229–36.

  22. 22.

    Doiron KM, Lavigne-Nicolas J, Cupples CG. Effect of interaction between 5-azacytidine and DNA (cytosine-5) methyltransferase on C-to-G and C-to-T mutations in Escherichia coli. Mutat Res. 1999;429:37–44.

  23. 23.

    Doskocil J, Sorm F. The effects of 5-azacytidine and 5-azauridine on protein synthesis in Escherichia coli. Biochem Biophys Res Commun. 1970;38:569–74.

  24. 24.

    Schaefer M, Hagemann S, Hanna K, Lyko F. Azacytidine inhibits RNA methylation at DNMT2 target sites in human cancer cell lines. Cancer Res. 2009;69:8127–32.

  25. 25.

    Ali IK, Ehrenkaufer GM, Hackney JA, Singh U. Growth of the protozoan parasite Entamoeba histolytica in 5-azacytidine has limited effects on parasite gene expression. BMC Genomics. 2007;8:7.

  26. 26.

    Qiu X, Hother C, Ralfkiaer UM, Sogaard A, Lu Q, Workman CT, et al. Equitoxic doses of 5-azacytidine and 5-aza-2'deoxycytidine induce diverse immediate and overlapping heritable changes in the transcriptome. PLoS One. 2010;5:e12994.

  27. 27.

    Ringquist S, Smith CL. The Escherichia coli chromosome contains specific, unmethylated dam and dcm sites. Proc Natl Acad Sci USA. 1992;89:4539–43.

  28. 28.

    Bhagwat AS, Roberts RJ. Genetic analysis of the 5-azacytidine sensitivity of Escherichia coli K-12. J Bacteriol. 1987;169:1537–46.

  29. 29.

    Friedman S. Bactericidal effect of 5-azacytidine on Escherichia coli carrying EcoRII restriction-modification enzymes. J Bacteriol. 1982;151:262–8.

  30. 30.

    Krasich R, Wu SY, Kuo HK, Kreuzer KN. Functions that protect Escherichia coli from DNA-protein crosslinks. DNA Repair (Amst). 2015;28:48–59.

  31. 31.

    Frommer M, McDonald LE, Millar DS, Collis CM, Watt F, Grigg GW, et al. A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proc Natl Acad Sci U S A. 1992;89:1827–31.

  32. 32.

    Salem AM, Nakano T, Takuwa M, Matoba N, Tsuboi T, Terato H, et al. Genetic analysis of repair and damage tolerance mechanisms for DNA-protein cross-links in Escherichia coli. J Bacteriol. 2009;191:5657–68.

  33. 33.

    Nakano T, Morishita S, Katafuchi A, Matsubara M, Horikawa Y, Terato H, et al. Nucleotide excision repair and homologous recombination systems commit differentially to the repair of DNA-protein crosslinks. Mol Cell. 2007;28:147–58.

  34. 34.

    Becket E, Chen F, Tamae C, Miller JH. Determination of hypersensitivity to genotoxic agents among Escherichia coli single gene knockout mutants. DNA Repair (Amst). 2010;9:949–57.

  35. 35.

    Grenier F, Matteau D, Baby V, Rodrigue S. Complete Genome Sequence of Escherichia coli BW25113. Genome Announc. 2014;2:e01038-14.

  36. 36.

    Lange R, Fischer D, Hengge-Aronis R. Identification of transcriptional start sites and the role of ppGpp in the expression of rpoS, the structural gene for the sigma S subunit of RNA polymerase in Escherichia coli. J Bacteriol. 1995;177:4676–80.

  37. 37.

    Takayanagi Y, Tanaka K, Takahashi H. Structure of the 5' upstream region and the regulation of the rpoS gene of Escherichia coli. Mol Gen Genet. 1994;243:525–31.

  38. 38.

    Mendoza-Vargas A, Olvera L, Olvera M, Grande R, Vega-Alvarado L, Taboada B, et al. Genome-wide identification of transcription start sites, promoters and transcription factor binding sites in E. coli. PLoS One. 2009;4:e7526.

  39. 39.

    Argaman L, Hershberg R, Vogel J, Bejerano G, Wagner EG, Margalit H, et al. Novel small RNA-encoding genes in the intergenic regions of Escherichia coli. Curr Biol. 2001;11:941–50.

  40. 40.

    Matzke MA, Mosher RA. RNA-directed DNA methylation: an epigenetic pathway of increasing complexity. Nat Rev Genet. 2014;15:394–408.

  41. 41.

    Salgado H, Peralta-Gil M, Gama-Castro S, Santos-Zavaleta A, Muniz-Rascado L, Garcia-Sotelo JS, et al. RegulonDB v8.0: omics data sets, evolutionary conservation, regulatory phrases, cross-validated gold standards and more. Nucleic Acids Res. 2013;41:D203–13.

  42. 42.

    Finkel SE. Long-term survival during stationary phase: evolution and the GASP phenotype. Nat Rev Microbiol. 2006;4:113–20.

  43. 43.

    Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, et al. Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006;2:1–11.

  44. 44.

    Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:P3.

  45. 45.

    Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.

  46. 46.

    R Core Team. R. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015. https://r-project.org.

Download references

Acknowledgements

We thank Jie Wang for advice on microarray data analysis. We thank Robert O’Donnell, Betsy Hutchison, Rebecca Huss, and Samantha Gage for critical reading of the manuscript.

Funding

We thank the Geneseo Foundation for financial support of the project.

Availability of data and materials

Microarray data were deposited in the Gene Expression Omnibus repository (accession number GSE73707).

Authors’ contributions

KTM and RDS conceived and designed the experiments. KTM, AHM, AD, SMH, JCL, and SC performed the experiments. All authors were responsible for data analysis. KTM prepared the manuscript. All authors read and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Correspondence to Kevin T. Militello.

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 log2 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 log2 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)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Militello, K.T., Simon, R.D., Mandarano, A.H. et al. 5-azacytidine induces transcriptome changes in Escherichia coli via DNA methylation-dependent and DNA methylation-independent mechanisms. BMC Microbiol 16, 130 (2016) doi:10.1186/s12866-016-0741-4

Download citation

Keywords

  • Stationary phase
  • 5-azacytidine
  • DNA methylation inhibitor
  • DNA methylation
  • 5-methylcytosine
  • Dcm
  • RpoS
  • Escherichia coli
  • Sodium bisulfite sequencing
  • DNA microarray