Contribution of DNA adenine methylation to gene expression heterogeneity in Salmonella enterica

Expression of Salmonella enterica loci harboring undermethylated GATC sites at promoters or regulatory regions was monitored by single cell analysis. Cell-to-cell differences in expression were detected in ten such loci (carA, dgoR, holA, nanA, ssaN, STM1290, STM3276, STM5308, gtr and opvAB), with concomitant formation of ON and OFF subpopulations. The ON and OFF subpopulation sizes varied depending on the growth conditions, suggesting that the population structure can be modulated by environmental control. All the loci under study except STM5308 displayed altered patterns of expression in strains lacking or overproducing Dam methylase, thereby confirming control by Dam methylation. Bioinformatic analysis identified potential binding sites for transcription factors OxyR, CRP and Fur, and analysis of expression in mutant backgrounds confirmed transcriptional control by one or more of such factors. Surveys of gene expression in pairwise combinations of Dam methylation-dependent loci revealed independent switching, thus predicting the formation of a high number of cell variants. This study expands the list of S. enterica loci under transcriptional control by Dam methylation, and underscores the relevance of the DNA adenine methylome as a source of phenotypic heterogeneity.