<figure><img src="image_1.png" /></figure>

HAL
open science

# Novel Identification of Bacterial Epigenetic Regulations Would Benefit From a Better Exploitation of Methylomic Data

Amaury Payelleville, Julien Brillard

► To cite this version:

Amaury Payelleville, Julien Brillard. Novel Identification of Bacterial Epigenetic Regulations Would Benefit From a Better Exploitation of Methylomic Data. Frontiers in Microbiology, 2021, 12, 10.3389/fmicb.2021.685670. [https://doi.org/10.3389/fmicb.2021.685670](https://doi.org/10.3389/fmicb.2021.685670). [hal-03227954](https:// HAL Id: hal-03227954)

https://hal.inrae.fr/hal-03227954

Submitted on 10 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Novel Identification of Bacterial Epigenetic Regulations Would Benefit From a Better Exploitation of Methylomic Data

Amaury Payelleville¹,² and Julien Brillard¹*

¹DGIMI, INRAE, Univ. Montpellier, Montpellier, France, ²Cellular and Molecular Microbiology, Faculté des Sciences, Université Libre de Bruxelles, Gosselies, Belgium

DNA methylation can be part of epigenetic mechanisms, leading to cellular subpopulations with heterogeneous phenotypes. While prokaryotic phenotypic heterogeneity is of critical importance for a successful infection by several major pathogens, the exact mechanisms involved in this phenomenon remain unknown in many cases. Powerful sequencing tools have been developed to allow the detection of the DNA methylated bases at the genome level, and they have recently been extensively applied on numerous bacterial species. Some of these tools are increasingly used for metagenomics analysis but only a limited amount of the available methylomic data is currently being exploited. Because newly developed tools now allow the detection of subpopulations differing in their genome methylation patterns, it is time to emphasize future strategies based on a more extensive use of methylomic data. This will ultimately help to discover new epigenetic gene regulations involved in bacterial phenotypic heterogeneity, including during host-pathogen interactions.

**Keywords:** bacterial epigenetics, DNA methylation, methylome, transcriptome, phenotypic heterogeneity

## OPEN ACCESS

**Edited by:**

Xin Deng,
City University of Hong Kong,
Hong Kong

**Reviewed by:**

Josep Casadesús,
Seville University, Spain

Chao Chen,
Wuhan University, China

*Correspondence:

Julien Brillard
julien.brillard@umontpellier.fr

**Specialty section:**

This article was submitted to
Evolutionary and Genomic
Microbiology,
a section of the journal
Frontiers in Microbiology

Received: 25 March 2021

Accepted: 22 April 2021

Published: 14 May 2021

**Citation:**

Payelleville A and Brillard J (2021)
Novel Identification of Bacterial
Epigenetic Regulations Would Benefit
From a Better Exploitation
of Methylomic Data.
Front. Microbiol. 12:685670.
doi: 10.3389/fmicb.2021.685670

## INTRODUCTION

Epigenetic regulations have been studied mainly in eukaryotes where they are involved in cell differentiation or disease occurrence, through diverse mechanisms such as histone modifications or DNA methylation. However, evidences for the existence of epigenetic regulation in prokaryotes are increasingly reported (Sánchez-Romero and Casadesús, 2020). Such gene regulations can occur by feedback loops (positive or negative) but most examples involve DNA methylation (Adhikari and Curtis, 2016). DNA methylation is a base modification system that acts by the addition of a methyl group from an S-adenosyl-methionine molecule to an Adenine or a Cytosine in the DNA. In the growing cell, this process usually occurs shortly after the DNA replication on the newly synthesized strand. Enzymes responsible for DNA methylation are called DNA-methyltransferases (MTases) and catalyze three types of DNA methylation modifications: N<sup>6</sup>-methyl-adenine (m6A), C<sup>5</sup>-methyl-cytosine (m5C) and N<sup>4</sup>-methyl-cytosine (m4C). While all three types are described in archaea and bacteria, m4C modification is not reported in eukaryotes. Genes encoding putative MTases are found in almost all bacterial species, and most bacterial genomes analyzed so far display a DNA methylation pattern
that is species or strain-specific (Blow et al., 2016). Many bacterial MTases belong to restriction-modification systems, as they are genetically and functionally associated to a restriction endonuclease (REase) that protects the bacterial cell from exogenous DNA (Marinus and Løbner-Olesen, 2014; Roberts et al., 2015). In addition, “solitary” or “orphan” MTases are frequently found in the genomes of bacteria (Blow et al., 2016) and some of them carry key roles in genome maintenance (Løbner-Olesen et al., 2005). DNA methylation can also affect the interaction of DNA-binding proteins with their cognate sites, either directly (e.g., steric hindrance) or by changes in DNA topology (Casadesús, 2016), resulting in epigenetic regulations (Casadesús and Low, 2013). This mechanism, among others, is responsible for prokaryotic phenotypic heterogeneity (Sánchez-Romero et al., 2020), a phenomenon of critical importance for successful infection by several major pathogens (Balaban et al., 2004; Weigel and Dersch, 2018).

Thanks to recent advances in long-read sequencing technologies, data of the methylated bases over the complete genome (methylome) of bacteria is progressively reported (Rand et al., 2017; Beaulaurier et al., 2019; García-Pastor

et al., 2019). However, most of these technologies are exclusively used to generate new genomics data, while the methylomic data are often set aside. Thus, there is a gap between the scarce reports of bacterial gene regulation associated to DNA methylation and the increasing availability of unexploited methylomic data. Given the pervasiveness of DNA methylation in prokaryotes, we believe that a deeper analysis of methylomic data could lead to identifying new candidates of epigenetically regulated genes. This review illustrates the importance of DNA methylation associated to epigenetic regulations in bacteria and aims to raise awareness on the available yet underused tools, in the field of bacterial epigenetics.

## BACTERIAL TRANSCRIPTIONAL FACTORS INVOLVED IN EPIGENETIC REGULATION

While DNA methylation occurs on motifs located anywhere on the DNA, typical bacterial epigenetic regulations are located in

<figure><img src="image_4.png" /></figure>
promoter regions (Figure 1). The main examples are linked to m6A modifications and either involve (i) CcrM, an MTase found in several α-proteobacteria, where it plays a critical role in controlling the cell cycle (Mouammine and Collier, 2018); or (ii) Dam, an MTase conserved in Enterobacteriaceae, which is associated to the formation of subpopulations with distinct phenotypes (van der Woude and Baumler, 2004; Marinus and Løbner-Olesen, 2014; Adhikari and Curtis, 2016; Sánchez-Romero and Casadesús, 2020).

The affinity of some DNA-binding proteins, including transcriptional regulators, can be affected by the DNA methylation state within or in the vicinity of their binding sites. After DNA replication (Figure 1A), the future methylation state of the locus will depend on whether the regulator or the MTase binds the newly synthesized (unmethylated) DNA strand first (Figure 1B). Gene expression will subsequently differ based on if (or where) the regulator is bound in these promoter regions. Thereby, this mechanism can therefore give rise to two bacterial subpopulations with distinct transcription patterns and consequently distinct phenotypes (Figure 1C).

Examples of epigenetic regulations associated to Dam often involve genes playing critical roles during bacterial-host interaction. Several transcriptional regulators were shown to be sensitive to the methylation state of the promoter region they control (for a recent review, see Sánchez-Romero and Casadesús, 2020). A canonic example is the regulation of the pilus-encoding pap operon by Lrp, extensively reviewed elsewhere (van der Woude et al., 1992, 1996; Peterson and Reich, 2008). Nevertheless, other examples of epigenetic gene regulation involving important transcriptional regulators deserve to be mentioned here (Figure 1D). A competition between the OxyR transcriptional regulator and Dam has been described for several promoter regions, as *agn43* (encoding an outer membrane protein) in *E. coli* or *gtr* or *opvAB* (encoding proteins involved in LPS modification) in *Salmonella* (Waldron et al., 2002; Wallecha et al., 2002; Cota et al., 2016; Sánchez-Romero et al., 2020). Such epigenetic regulations direct the formation of subpopulations that have a distinct fitness depending on the environment. For instance in *Salmonella*, control of the *opvAB* operon by OxyR produces one subpopulation that is resistant to infection by many bacteriophages but avirulent, and another that is phage sensitive but able to infect animal hosts (Cota et al., 2015). HdfR is a transcriptional regulator involved in the epigenetic regulation of the *std* fimbrial operon (García-Pastor et al., 2019). Here, the regulation mechanism is complexified by a positive feedback loop involving two additional regulators (StdE and StdF) encoded by the *std* operon. The Fur transcriptional regulator is critical for the iron stress response (Escolar et al., 1999) and is also involved in the epigenetic regulation of the *sci1* operon, encoding a Type VI Secretion System (T6SS) in *E. coli* (Brunet et al., 2011; Brunet et al., 2020). In this case, an intraspecific bacterial competition occurs between the two subpopulations dependent on the active expression of the T6SS. In addition to these few selected examples, a list of other regulators sensitive to DNA methylation can be found elsewhere (Sánchez-Romero and Casadesús, 2020). Altogether,

they illustrate that the DNA methylation pattern is important for several bacterial phenotypes.

## BACTERIAL PHENOTYPES ASSOCIATED TO CHANGES IN DNA-METHYLATION PATTERN

The DNA methylation pattern can be significantly impacted by environmental conditions, as observed in some α or γ-proteobacteria (Hale et al., 1994; Ichida et al., 2007; Doberenz et al., 2017). Yet, the methylome of several γ-proteobacteria appeared very stable despite various growth conditions tested (Cohen et al., 2016; Westphal et al., 2016; Payelleville et al., 2018). In many examples, major phenotypic modifications are driven by the modification of the DNA methylation pattern through mutation or overexpression of an MTase gene. For instance, *dam* mutation caused impaired virulence in *Salmonella*, *Klebsiella*, *Haemophilus*, *Yersinia*, and *Actinobacillus* (García-Del Portillo et al., 1999; Heithoff et al., 1999; Watson et al., 2004; Robinson et al., 2005; Taylor et al., 2005; Wu et al., 2006; Mehling et al., 2007). While *dam* deletion can sometimes be lethal for the bacterium (as in *V. cholerae* or *A. hydrophila*) (Julio et al., 2001; Erova et al., 2006) its overexpression often leads to a decreased virulence phenotype (Heithoff et al., 1999; Julio et al., 2002; Chen et al., 2003; Erova et al., 2012; Payelleville et al., 2017).

Other examples below illustrate the major phenotypes found in association with DNA-methylation pattern deregulation, whatever the type of methylated based involved. For most of these examples, the underlying mechanisms responsible for the modified phenotype has not been described yet. The CamA orphan m6A MTase conserved among all *C. difficile* is involved in sporulation and persistence and many genes in the MTase mutant are shown to be differentially regulated (Oliveira et al., 2020). Dcm, a broadly distributed m5C MTase in *Enterobacteriaceae*, plays a role during the stationary growth phase as its deletion leads to an increase in the RpoS sigma factor expression (Kahramanoglou et al., 2012; Militello et al., 2012, 2016, 2020) and it may also be linked to antimicrobial compound resistance (Militello et al., 2014). Furthermore, the deletion of m4C MTases encoding genes can impair virulence and cause broad transcriptional changes as shown in *H. pylori* (Kumar et al., 2012, 2018). In another human pathogen, *Leptospira interrogans*, the impaired virulence caused by the deletion of an MTase is associated to the dysregulation of an extracytoplasmic function (ECF) sigma factor (Gaultney et al., 2020).

The so-called “Phasevarion” (for phase-variable regulon) allows broad changes in DNA-methylation patterns. It occurs in bacteria that express MTases with phase variation character (i.e., a reversible switch in expression leading to single-cell phenotypic heterogeneity). This phenomenon provides a way to modify the expression of numerous genes simultaneously through epigenetic regulation by a single event of phase variation. This way, phasevarions confer the ability for the bacterium to adapt to broader environmental conditions. The phenomenon has been described in several major pathogens [*H. influenzae*
(De Bolle et al., 2000), *H. pylori* (de Vries et al., 2002; Srikhanta et al., 2017b), *Streptococcus suis* (Atack et al., 2018b), *S. pneumoniae* (Manso et al., 2014), *Neisseria meningitidis* (Jen et al., 2014; Seib et al., 2015), *N. gonorrheae* (Srikhanta et al., 2009; Seib et al., 2017), *Moraxella catarrhalis* (Blakeway et al., 2014; Blakeway et al., 2019), *Kingella kingae* (Srikhanta et al., 2017a)] and is reviewed elsewhere (Atack et al., 2018a; Seib et al., 2020). Moreover, new bioinformatic tools allowing for the detection of genetic signatures of those events have been recently developed, giving hope for a further increase in the description of the underlying mechanism (Atack et al., 2020a,b). All the phasevarions described until now involve MTases that are part of RM systems. Given the significant proportion of type II solitary MTases in prokaryotes [more than 50%, (Blow et al., 2016)], and considering that they are involved in most described epigenetic regulations, we can hypothesize that phasevarion involving solitary MTases might be discovered in the future. Again, even though the expression switch of the MTase is well described, the exact epigenetic regulation mechanisms of the various regulon members remain to be elucidated. Modern technologies enabling the analysis of the cells' methylome may help to answer these questions.

## SEQUENCING METHODS TO ANALYZE BACTERIAL METHYLOMES

Three main methods are currently used for the base resolution sequencing of DNA methylation (reviewed in Beaulaurier et al., 2019). They all require powerful bioinformatic tools to have an accurate and complete view of the modified bases in the whole genome.

(I) The Whole Genome Bisulfite Sequencing (WGBS) has been used for years, mostly in eukaryotes. In this approach, unmethylated cytosines are converted into uracils by the bisulfite treatment. After sequencing (usually short-read Illumina sequencing) and alignment, m5C can then be identified (Carless, 2009). This tool which is specific to the detection of m5C has found scarce application for the bacterial methylome analysis so far (Kahramanoglou et al., 2012). (II) The possibility to use the long-range method SMRT sequencing (for Single Molecule Real Time) to detect DNA methylation at a genomic scale was described in 2010 (Flusberg et al., 2010; Cloney, 2016). This technology makes use of a DNA polymerase sensitive to base modifications. Upon recognition of a modified base on the ssDNA matrix, a delay in polymerization is generated during the recording of the DNA sequence (InterPulse Duration or IPD). While it can easily detect m4C and m6A, this technology has strong limitations to detect m5C (requires a substantial coverage rate). Therefore, performing both WGBS and SMRT sequencing can allow to determine an exhaustive methylome, something rarely done (Payelleville et al., 2018; Vandenbussche et al., 2021). (III) More recently, another long-range sequencing method, the Oxford Nanopore technology (ONT), was shown to efficiently detect modified bases. While ssDNA crosses nanopores embedded in a lipid membrane,

a voltage potential is applied. Analysis of the electrolytic current signals, which are sensitive to base modifications, reveals both the DNA sequence and the methylation state of the DNA matrix (Rand et al., 2017; Simpson et al., 2017). The earliest studies using ONT were focused on m5C detection of CpG islands, found in some eukaryotes (Laszlo et al., 2013), but various bioinformatic models have been developed since to increase the accuracy of other DNA methylation predictions. In certain DNA motifs, such models are now able to detect m5C and m6A using a low read coverage (as low as twofold) (Ni et al., 2019) with a significant precision on *E. coli* data (Liu et al., 2019). Recently, the methylomes of eight microbial reference species have been validated using various methods, including ONT sequencing for m6A detection (McIntyre et al., 2019).

The increasing number of bacterial metagenomic studies which use such sequencing technologies (Petersen et al., 2019; Tourancheau et al., 2021) is about to expand the set of data that could be used in parallel for analyzing DNA methylation. This, in turn, could stimulate bacterial epigenetic research. In addition, because the DNA methylation pattern is often strain-dependent, incorporation of methylomic information into shotgun metagenomic analyses was proposed as a new tool for distinguishing genomes of closely related strains, hence providing a much more accurate clustering of metagenomic sequences (Beaulaurier et al., 2018). If such strategy became widespread in the future, it would reinforce the importance of exploiting methylomic data.

## COMBINING METHYLOME ANALYSIS WITH OTHER APPROACHES TO IDENTIFY PUTATIVE EPIGENETIC REGULATION

Numerous examples, as illustrated above, have demonstrated that the DNA methylation pattern at a given site can impact gene expression (Sánchez-Romero and Casadesús, 2020). Therefore, in a bacterial genome, each subset of unmethylated recognition motifs could be considered as a putative epigenetic regulatory site (Figure 1). Such approach was recently proved to be strikingly efficient in *Salmonella* and led to the identification of several new genes displaying expression heterogeneity controlled by Dam-methylation (Sánchez-Romero et al., 2020). A broad conservation among prokaryotes of unmethylated sites that are usually recognized by conserved MTases (Dam in γ-proteobacteria, or CcrM in α-proteobacteria) (Blow et al., 2016) strengthens the hypothesis of a widespread occurrence of this regulatory mechanism in bacteria.

Given the large amount of transcriptomic data released in the databases, together with the increasing acquisition of genomic data by the help of technologies that also allow to identify the DNA methylation pattern (SMRT or ONT), we propose the coupling of transcriptome analysis with extensive methylome analysis. It may lead to the identification
<figure><img src="image_5.png" /></figure>

of putative epigenetic regulation networks. This would be particularly true if unmethylated motifs located in gene regulatory regions (i.e., promoters) correlate with differential regulation of the genes in an altered environmental condition (or in a mutant strain, **Figure 2**).

Up to now, this combined methylome/transcriptome strategy has been rarely employed and led to the identification of only a limited number of candidate genes under a putative epigenetic regulation, which still need to be confirmed by additional mechanistic studies (Blow et al., 2016; Cohen et al., 2016; Doberenz et al., 2017; Payelleville et al., 2018; Nye et al., 2019). This low occurrence of candidates has to be balanced by the fact that a complex gene regulation may require multiple factors for a fine tuning of expression, as exemplified by the epigenetic regulation of the Std fimbriae in *Salmonella* described above (García-Pastor et al., 2019). Furthermore, the contribution of nucleoid associated proteins, such as H-NS, in epigenetic mechanisms involving DNA-methylation has also been reported in various cases and therefore also contribute to the formation of bacterial subpopulations (Nicholson and Low, 2000; Camacho et al., 2005; Cota et al., 2016).

## BACTERIAL PHENOTYPIC HETEROGENEITY: METHYLOME AND TRANSCRIPTOME ANALYSES OF SUBPOPULATIONS

The methylome analysis methods described above determine the DNA methylation to the nucleotide resolution at the population level. However, distinct DNA methylation patterns can drive the emergence of different subpopulations with different expression profiles (Figure 2). In 2015, an improvement in SMRT sequencing, SMALR for Single Molecule modification Analysis of Long Reads, was proposed (Beaulaurier et al., 2015). The enhancement resides on the ability of SMALR to identify epigenetic heterogeneity, where a subpopulation displays a distinct methylation pattern compared to the rest of the population. Despite its potential to identify subpopulations with different methylation patterns, currently few studies reporting the use of SMALR can be found in the literature (Modlin et al., 2020). This may be due to an ongoing need for improved or analysis-specific tools based
on modern sequencing technologies to decipher the mechanisms that give rise to the formation of subpopulations in an isogenic bacterial culture.

Importantly, two tools were recently designed to study complete gene expression at the subpopulation level: (i) Record-seq (Schmidt et al., 2018) can report a change in gene expression during a bottleneck situation where the amount of events is too low for a classic global transcriptomic analysis (Schmidt et al., 2018; Tanna et al., 2020); (ii) PETRI-seq (Blattman et al., 2020) allows to detect subpopulations with a different transcriptomic profile where less mRNA is needed, compared to classic transcriptional studies. Combining the SMALR and one of those transcriptomic tools should drastically improve the detection of candidate genes subject to DNA methylation regulation and heterogeneously expressed among a population.

In parallel, as prokaryotic phenotypic heterogeneity is of critical importance for a successful infection by various major pathogens, it is crucial that, in the near future, evolution of methylome analysis (supported by both improved sequencing coverage rates, and development of appropriate computational tools) will be more sensitive to allow for the distinction of differential DNA methylation patterns among a single DNA sample. This will increase the possibility to identify heterogeneity in epigenetic marks between bacterial subpopulations, including *in vivo* during bacterial infections.

## CONCLUSION

It is now time to consider the large amount of available data that could be thoroughly exploited in order to identify new

## REFERENCES

Adhikari, S., and Curtis, P. D. (2016). DNA methyltransferases and epigenetic regulation in bacteria. *FEMS Microbiol. Rev.* 40, 575–591. doi: 10.1093/femsre/fuw023

Atack, J. M., Guo, C., Litfin, T., Yang, L., Blackall, P. J., Zhou, Y., et al. (2020a). Systematic analysis of REBASE identifies numerous Type I restriction-modification systems with duplicated, distinct *hsdS* specificity genes that can switch system specificity by recombination. *mSystems* 5, e00497–e00520.

Atack, J. M., Guo, C., Yang, L., Zhou, Y., and Jennings, M. P. (2020b). DNA sequence repeats identify numerous Type I restriction-modification systems that are potential epigenetic regulators controlling phase-variable regulons; phasevarions. *FASEB J.: Off. Publ. Federation Am. Societies Exp. Biol.* 34, 1038–1051. doi: 10.1096/fj.201901536rr

Atack, J. M., Tan, A., Bakaletz, L. O., Jennings, M. P., and Seib, K. L. (2018a). Phasevarions of bacterial pathogens: methylomics sheds new light on old enemies. *Trends Microbiol.* 26, 715–726. doi: 10.1016/j.tim.2018.01.008

Atack, J. M., Weinert, L. A., Tucker, A. W., Husna, A. U., Wileman, T. M., Hadjirin, N., et al. (2018b). *Streptococcus suis* contains multiple phase-variable methyltransferases that show a discrete lineage distribution. *Nucleic Acids Res.* 46, 11466–11476.

Balaban, N. Q., Merrin, J., Chait, R., Kowalik, L., and Leibler, S. (2004). Bacterial persistence as a phenotypic switch. *Science (New York, N.Y.)* 305, 1622–1625. doi: 10.1126/science.1099390

candidates of putative epigenetic regulation. Although it might often be challenging to confirm such mechanisms of regulation, the candidates detected would most likely unveil major roles in the life cycle of the bacteria. This assumption is after all exemplified by the mechanisms of epigenetic regulation which have been deciphered up to now (van der Woude et al., 1996; Wallecha et al., 2002; Camacho et al., 2005; van der Woude and Henderson, 2008; Brunet et al., 2011; Cota et al., 2016; García-Pastor et al., 2019).

## AUTHOR CONTRIBUTIONS

AP and JB designed and wrote the manuscript and designed the figures using Biorender.com. Both the authors contributed to the article and approved the submitted version.

## FUNDING

This project has received funding from the European Union's Horizon 2020 Research and Innovation Program under the Marie Skłodowska-Curie Grant Agreement No. 801505. AP also received a Grant from GAIA doctoral school #584. JB received funding from INRA Plant Health and Environment (SPE) division (SPE-IB17-DiscriMet) and from the French National Research Agency (EPI-PATH project, ANR-17-CE20-0005).

## ACKNOWLEDGMENTS

We extremely grateful to Tatjana Schlechtweg for her careful reading of the manuscript.

Beaulaurier, J., Schadt, E. E., and Fang, G. (2019). Deciphering bacterial epigenomes using modern sequencing technologies. *Nat. Rev. Genet.* 20, 157–172. doi: 10.1038/s41576-018-0081-3

Beaulaurier, J., Zhang, X.-S., Zhu, S., Sebra, R., Rosenbluh, C., Deikus, G., et al. (2015). Single molecule-level detection and long read-based phasing of epigenetic variations in bacterial methylomes. *Nat. Commun.* 6:7438.

Beaulaurier, J., Zhu, S., Deikus, G., Mogno, I., Zhang, X.-S., Davis-Richardson, A., et al. (2018). Metagenomic binning and association of plasmids with bacterial host genomes using DNA methylation. *Nat. Biotechnol.* 36, 61–69. doi: 10.1038/nbt.4037

Blakeway, L. V., Power, P. M., Jen, F. E. C., Worboys, S. R., Boitano, M., Clark, T. A., et al. (2014). ModM DNA methyltransferase methylome analysis reveals a potential role for *Moraxella catarrhalis* phasevarions in otitis media. *FASEB J.: Off. Publ. Federation Am. Societies Exp. Biol.* 28, 5197–5207. doi: 10.1096/fj.14-256578

Blakeway, L. V., Tan, A., Jurcisek, J. A., Bakaletz, L. O., Atack, J. M., Peak, I. R., et al. (2019). The *Moraxella catarrhalis* phase-variable DNA methyltransferase ModM3 is an epigenetic regulator that affects bacterial survival in an in vivo model of otitis media. *BMC Microbiol.* 19:276. doi: 10.1186/s12866-019-1660-y

Blattman, S. B., Jiang, W., Oikonomou, P., and Tavazoie, S. (2020). Prokaryotic single-cell RNA sequencing by in situ combinatorial indexing. *Nat. Microbiol.* 5, 1192–1201. doi: 10.1038/s41564-020-0729-6

Blow, M. J., Clark, T. A., Daum, C. G., Deutschbauer, A. M., Fomenkov, A., Fries, R., et al. (2016). The epigenomic landscape of prokaryotes. *PLoS Genet.* 12:e1005854. doi: 10.1371/journal.pgen.1005854
Brunet, Y. R., Bernard, C. S., and Cascales, E. (2020). Fur-dam regulatory interplay at an internal promoter of the enteroaggregative *Escherichia coli* Type VI secretion *sci1* gene cluster. *J. Bacteriol.* **202**, e00075-20

Brunet, Y. R., Bernard, C. S., Gavioli, M., Lloubès, R., and Cascales, E. (2011). An epigenetic switch involving overlapping *fur* and DNA methylation optimizes expression of a type VI secretion gene cluster. *PLoS Genet.* **7:e1002205**. doi: 10.1371/journal.pgen.1002205

Camacho, E. M., Serna, A., Madrid, C., Marqués, S., Fernández, R., De La Cruz, F., et al. (2005). Regulation of *finP* transcription by DNA adenine methylation in the virulence plasmid of *Salmonella enterica*. *J. Bacteriol.* **187**, 5691–5699. doi: 10.1128/jb.187.16.5691-5699.2005

Carless, M. (2009). Investigation of genomic methylation status using methylation-specific and bisulfite sequencing polymerase chain reaction. *Methods Mol. Biol.* (*Clifton, N.J.*) **523**, 217–234. doi: 10.1007/978-1-59745-190-1_15

Casadesús, J. (2016). Bacterial DNA methylation and methylomes. *Adv. Exp. Med. Biol.* **945**, 35–61. doi: 10.1007/978-3-319-43624-1_3

Casadesús, J., and Low, D. A. (2013). Programmed heterogeneity: epigenetic mechanisms in bacteria. *J. Biol. Chem.* **288**, 13929–13935. doi: 10.1074/jbc.r113.472274

Chen, L., Paulsen, D. B., Scruggs, D. W., Banes, M. M., Reeks, B. Y., and Lawrence, M. L. (2003). Alteration of DNA adenine methylase (Dam) activity in *Pasteurella multocida* causes increased spontaneous mutation frequency and attenuation in mice. *Microbiology (Reading, England)* **149**, 2283–2290. doi: 10.1099/mic.0.26251-0

Cloney, R. (2016). DNA methylation: a SMRT analysis of prokaryotic epigenomes. *Nat. Rev. Genet.* **17**:195. doi: 10.1038/nrg.2016.22

Cohen, N. R., Ross, C. A., Jain, S., Shapiro, R. S., Gutierrez, A., Belenky, P., et al. (2016). A role for the bacterial GATC methylome in antibiotic stress survival. *Nat. Genet.* **48**, 581–586. doi: 10.1038/ng.3530

Cota, I., Bunk, B., Spröer, C., Overmann, J., König, C., and Casadesús, J. (2016). OxyR-dependent formation of DNA methylation patterns in *OpvAB*<sup>OFF</sup> and *OpvAB*<sup>ON</sup> cell lineages of *Salmonella enterica*. *Nucleic Acids Res.* **44**, 3595–3609. doi: 10.1093/nar/gkv1483

Cota, I., Sánchez-Romero, M. A., Hernández, S. B., Pucciarelli, M. G., García-Del Portillo, F., and Casadesús, J. (2015). Epigenetic control of *Salmonella enterica* O-Antigen chain length: a tradeoff between virulence and bacteriophage resistance. *PLoS Genet.* **11:e1005667**. doi: 10.1371/journal.pgen.1005667

De Bolle, X., Bayliss, C. D., Field, D., Van De Ven, T., Saunders, N. J., Hood, D. W., et al. (2000). The length of a tetranucleotide repeat tract in *Haemophilus influenzae* determines the phase variation rate of a gene with homology to type III DNA methyltransferases. *Mol. Microbiol.* **35**, 211–222. doi: 10.1046/j.1365-2958.2000.01701.x

de Vries, N., Duinsbergen, D., Kuipers, E. J., Pot, R. G. J., Wiesenekker, P., Penn, C. W., et al. (2002). Transcriptional phase variation of a type III restriction-modification system in *Helicobacter pylori*. *J. Bacteriol.* **184**, 6615–6623. doi: 10.1128/jb.184.23.6615-6624.2002

Doberenz, S., Eckweiler, D., Reichert, O., Jensen, V., Bunk, B., Spröer, C., et al. (2017). Identification of a *Pseudomonas aeruginosa* PAO1 DNA methyltransferase, its targets, and physiological roles. *mBio* **8**, e02312–e02316.

Erova, T. E., Kosykh, V. G., Sha, J., and Chopra, A. K. (2012). DNA adenine methyltransferase (Dam) controls the expression of the cytotoxic enterotoxin (*act*) gene of *Aeromonas hydrophila* via tRNA modifying enzyme-glucose-inhibited division protein (GidA). *Gene* **498**, 280–287. doi: 10.1016/j.gene.2012.02.024

Escolar, L., Pérez-Martín, J., and De Lorenzo, V. (1999). Opening the iron box: transcriptional metalloregulation by the Fur protein. *J. Bacteriol.* **181**, 6223–6229. doi: 10.1128/jb.181.20.6223-6229.1999

Flusberg, B. A., Webster, D. R., Lee, J. H., Travers, K. J., Olivares, E. C., Clark, T. A., et al. (2010). Direct detection of DNA methylation during single-molecule, real-time sequencing. *Nat. Methods* **7**, 461–465. doi: 10.1038/nmeth.1459

García-Del Portillo, F., Pucciarelli, M. G., and Casadesús, J. (1999). DNA adenine methylase mutants of *Salmonella typhimurium* show defects in protein secretion, cell invasion, and M cell cytotoxicity. *Proc. Natl. Acad. Sci. U.S.A.* **96**, 11578–11583. doi: 10.1073/pnas.96.20.11578

García-Pastor, L., Sánchez-Romero, M. A., Jakomin, M., Puerta-Fernández, E., and Casadesús, J. (2019). Regulation of bistability in the std fimbrial operon of *Salmonella enterica* by DNA adenine methylation and transcription factors *HdfR*, *StdE* and *StdF*. *Nucleic Acids Res.* **47**, 7929–7941. doi: 10.1093/nar/gkz530

Gaultney, R. A., Vincent, A. T., Lorioux, C., Coppée, J.-Y., Sismeiro, O., Varet, H., et al. (2020). 4-Methylcytosine DNA modification is critical for global epigenetic regulation and virulence in the human pathogen *Leptospira interrogans*. *Nucleic Acids Res.* **48**, 12102–12115. doi: 10.1093/nar/gkaa966

Hale, W. B., Van Der Woude, M. W., and Low, D. A. (1994). Analysis of nonmethylated GATC sites in the *Escherichia coli* chromosome and identification of sites that are differentially methylated in response to environmental stimuli. *J. Bacteriol.* **176**, 3438–3441. doi: 10.1128/jb.176.11.3438-3441.1994

Heithoff, D. M., Sinsheimer, R. L., Low, D. A., and Mahan, M. J. (1999). An essential role for DNA adenine methylation in bacterial virulence. *Science* (New York, N.Y.) **284**, 967–970. doi: 10.1126/science.284.5416.967

Ichida, H., Matsuyama, T., Abe, T., and Koba, T. (2007). DNA adenine methylation changes dramatically during establishment of symbiosis. *FEBS J.* **274**, 951–962. doi: 10.1111/j.1742-4658.2007.05643.x

Jen, F. E. C., Seib, K. L., and Jennings, M. P. (2014). Phasevarious mediate epigenetic regulation of antimicrobial susceptibility in *Neisseria meningitidis*. *Antimicrob. Agents Chemother.* **58**, 4219–4221. doi: 10.1128/aac.00004-14

Julio, S. M., Heithoff, D. M., Provenzano, D., Klose, K. E., Sinsheimer, R. L., Low, D. A., et al. (2001). DNA adenine methylase is essential for viability and plays a role in the pathogenesis of *Yersinia pseudotuberculosis* and *Vibrio cholerae*. *Infection Immunity* **69**, 7610–7615. doi: 10.1128/iai.69.12.7610-7615.2001

Julio, S. M., Heithoff, D. M., Sinsheimer, R. L., Low, D. A., and Mahan, M. J. (2002). DNA adenine methylase overproduction in *Yersinia pseudotuberculosis* alters YopE expression and secretion and host immune responses to infection. *Infection Immunity* **70**, 1006–1009. doi: 10.1128/iai.70.2.1006-1009.2002

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

Kumar, R., Mukhopadhyay, A. K., Ghosh, P., and Rao, D. N. (2012). Comparative transcriptomics of *H. pylori* strains AM5, SS1 and their *hpyAVIBM* deletion mutants: possible roles of cytosine methylation. *PLoS One* **7:e42303**. doi: 10.1371/journal.pone.0042303

Kumar, S., Karmakar, B. C., Nagarajan, D., Mukhopadhyay, A. K., Morgan, R. D., and Rao, D. N. (2018). N4-cytosine DNA methylation regulates transcription and pathogenesis in *Helicobacter pylori*. *Nucleic Acids Res.* **46**, 3429–3445. doi: 10.1093/nar/gky126

Laszlo, A. H., Derrington, I. M., Brinkerhoff, H., Langford, K. W., Nova, I. C., Samson, J. M., et al. (2013). Detection and mapping of 5-methylcytosine and 5-hydroxymethylcytosine with nanopore MspA. *Proc. Natl. Acad. Sci. U.S.A.* **110**, 18904–18909. doi: 10.1073/pnas.1310240110

Liu, Q., Fang, L., Yu, G., Wang, D., Xiao, C.-L., and Wang, K. (2019). Detection of DNA base modifications by deep recurrent neural network on Oxford Nanopore sequencing data. *Nat. Commun.* **10**:2449.

Løbner-Olesen, A., Skovgaard, O., and Marinus, M. G. (2005). Dam methylation: coordinating cellular processes. *Curr. Opin. Microbiol.* **8**, 154–160. doi: 10.1016/j.mib.2005.02.009

Manso, A. S., Chai, M. H., Atack, J. M., Furi, L., De Ste Croix, M., Haigh, R., et al. (2014). A random six-phase switch regulates pneumococcal virulence via global epigenetic changes. *Nat. Commun.* **5**:5055.

Marinus, M. G., and Løbner-Olesen, A. (2014). DNA methylation. *EcoSal Plus* **6**:10.

McIntyre, A. B. R., Alexander, N., Grigorev, K., Bezdan, D., Sichtig, H., Chiu, C. Y., et al. (2019). Single-molecule sequencing detection of N6-methyladenine in microbial reference materials. *Nat. Commun.* **10**:579.

Mehling, J. S., Lavender, H., and Clegg, S. (2007). A Dam methylation mutant of *Klebsiella pneumoniae* is partially attenuated. *FEMS Microbiol. letters* **268**, 187–193. doi: 10.1111/j.1574-6968.2006.00581.x
Militello, K. T., Finnerty-Haggerty, L., Kambhampati, O., Huss, R., and Knapp, R. (2020). DNA cytosine methyltransferase enhances viability during prolonged stationary phase in *Escherichia coli*. *FEMS Microbiol. Lett.* 367:fnaa166.

Militello, K. T., Mandarano, A. H., Varechtchouk, O., and Simon, R. D. (2014). Cytosine DNA methylation influences drug resistance in *Escherichia coli* through increased *sugE* expression. *FEMS Microbiol. Lett.* 350, 100–106. doi: 10.1111/1574-6968.12299

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

Militello, K. T., Simon, R. D., Qureshi, M., Maines, R., Vanhorne, M. L., Hennick, S. M., et al. (2012). Conservation of Dcm-mediated cytosine DNA methylation in *Escherichia coli*. *FEMS Microbiol. Lett.* 328, 78–85. doi: 10.1111/j.1574-6968.2011.02482.x

Modlin, S. J., Conkle-Gutierrez, D., Kim, C., Mitchell, S. N., Morrissey, C., Weinrick, B. C., et al. (2020). Drivers and sites of diversity in the DNA adenine methylomes of 93 *Mycobacterium tuberculosis* complex clinical isolates. *Elife* 9:e58542.

Mouammine, A., and Collier, J. (2018). The impact of DNA methylation in Alphaproteobacteria. *Mol. Microbiol.* 110, 1–10. doi: 10.1201/978020348 7013.ch1

Ni, P., Huang, N., Zhang, Z., Wang, D.-P., Liang, F., Miao, Y., et al. (2019). DeepSignal: detecting DNA methylation state from Nanopore sequencing reads using deep-learning. *Bioinformatics* (Oxford, England) 35, 4586–4595. doi: 10.1093/bioinformatics/btz276

Nicholson, B., and Low, D. (2000). DNA methylation-dependent regulation of *pel* expression in *Salmonella typhimurium*. *Mol. Microbiol.* 35, 728–742. doi: 10.1046/j.1365-2958.2000.01743.x

Nye, T. M., Jacob, K. M., Holley, E. K., Nevarez, J. M., Dawid, S., Simmons, L. A., et al. (2019). DNA methylation from a Type I restriction modification system influences gene expression and virulence in *Streptococcus pyogenes*. *PLoS Pathogens* 15:e1007841. doi: 10.1371/journal.ppat.1007841

Oliveira, P. H., Ribis, J. W., Garrett, E. M., Trzilova, D., Kim, A., Sekulovic, O., et al. (2020). Epigenomic characterization of *Clostridioides difficile* finds a conserved DNA methyltransferase that mediates sporulation and pathogenesis. *Nat. Microbiol.* 5, 166–180. doi: 10.1038/s41564-019- 0613-4

Payelleville, A., Lanois, A., Gislard, M., Dubois, E., Roche, D., Cruveiller, S., et al. (2017). DNA Adenine Methyltransferase (Dam) overexpression impairs *Photorhabdus luminescens* motility and virulence. *Front. Microbiol.* 8:1671. doi: 10.3389/fmicb.2017.01671

Payelleville, A., Legrand, L., Ogier, J.-C., Roques, C., Roulet, A., Bouchez, O., et al. (2018). The complete methylome of an entomopathogenic bacterium reveals the existence of loci with unmethylated Adenines. *Sci. Rep.* 8:12091.

Petersen, L. M., Martin, I. W., Moschetti, W. E., Kershaw, C. M., and Tsongalis, G. J. (2019). Third-generation sequencing in the clinical laboratory: exploring the advantages and challenges of nanopore sequencing. *J. Clin. Microbiol.* 58, e01315–e01319.

Peterson, S. N., and Reich, N. O. (2008). Competitive Lrp and Dam assembly at the *pap* regulatory region: implications for mechanisms of epigenetic regulation. *J. Mol. Biol.* 383, 92–105. doi: 10.1016/j.jmb.2008.07.086

Rand, A. C., Jain, M., Eizenga, J. M., Musselman-Brown, A., Olsen, H. E., Akeson, M., et al. (2017). Mapping DNA methylation with high-throughput nanopore sequencing. *Nat. Methods* 14, 411–413. doi: 10.1038/nmeth. 4189

Roberts, R. J., Vincze, T., Posfai, J., and Macelis, D. (2015). REBASE—a database for DNA restriction and modification: enzymes, genes and genomes. *Nucleic Acids Res.* 43, D298–D299.

Robinson, V. L., Oyston, P. C. F., and Titball, R. W. (2005). A *dam* mutant of *Yersinia pestis* is attenuated and induces protection against plague. *FEMS Microbiol. Lett.* 252, 251–256.

Sánchez-Romero, M. A., and Casadesús, J. (2020). The bacterial epigenome. *Nat. Rev. Microbiol.* 18, 7–20. doi: 10.1038/s41579-019-0286-2

Sánchez-Romero, M. A., Olivenza, D. R., Gutiérrez, G., and Casadesús, J. (2020). Contribution of DNA adenine methylation to gene expression heterogeneity

in *Salmonella enterica*. *Nucleic Acids Res.* 48, 11857–11867. doi: 10.1093/nar/ gkaa730

Schmidt, F., Cherepkova, M. Y., and Platt, R. J. (2018). Transcriptional recording by CRISPR spacer acquisition from RNA. *Nature* 562, 380–385. doi: 10.1038/ s41586-018-0569-1

Seib, K. L., Jen, F. E. C., Scott, A. L., Tan, A., and Jennings, M. P. (2017). Phase variation of DNA methyltransferases and the regulation of virulence and immune evasion in the pathogenic *Neisseria*. *Pathogens Dis.* 75:ftx080. doi: 10.1093/femspd/ftx080

Seib, K. L., Jen, F. E. C., Tan, A., Scott, A. L., Kumar, R., Power, P. M., et al. (2015). Specificity of the ModA11, ModA12 and ModD1 epigenetic regulator N(6)-adenine DNA methyltransferases of *Neisseria meningitidis*. *Nucleic Acids Res.* 43, 4150–4162. doi: 10.1093/nar/gkv219

Seib, K. L., Srikhanta, Y. N., Atack, J. M., and Jennings, M. P. (2020). Epigenetic regulation of virulence and immunoevasion by Phase-Variable restriction-modification systems in bacterial pathogens. *Ann. Rev. Microbiol.* 74, 655–671. doi: 10.1146/annurev-micro-090817-062346

Simpson, J. T., Workman, R. E., Zuzarte, P. C., David, M., Dursi, L. J., and Timp, W. (2017). Detecting DNA cytosine methylation using nanopore sequencing. *Nat. Methods* 14, 407–410. doi: 10.1038/nmeth.4184

Srikhanta, Y. N., Dowideit, S. J., Edwards, J. L., Falsetta, M. L., Wu, H.-J., Harrison, O. B., et al. (2009). Phasevarions mediate random switching of gene expression in pathogenic *Neisseria*. *PLoS Pathogens* 5:e1000400. doi: 10.1371/journal.ppat. 1000400

Srikhanta, Y. N., Fung, K. Y., Pollock, G. L., Bennett-Wood, V., Howden, B. P., and Hartland, E. L. (2017a). Phasevarion-regulated virulence in the emerging pediatric pathogen *Kingella kingae*. *Infection Immunity* 85, e00319–e00417.

Srikhanta, Y. N., Gorrell, R. J., Power, P. M., Tsyganov, K., Boitano, M., Clark, T. A., et al. (2017b). Methylomic and phenotypic analysis of the ModH5 phasevarion of *Helicobacter pylori*. *Sci. Rep.* 7:16140.

Tanna, T., Schmidt, F., Cherepkova, M. Y., Okoniewski, M., and Platt, R. J. (2020). Recording transcriptional histories using Record-seq. *Nat. Protoc.* 15, 513–539. doi: 10.1038/s41596-019-0253-4

Taylor, V. L., Titball, R. W., and Oyston, P. C. F. (2005). Oral immunization with a *dam* mutant of *Yersinia pseudotuberculosis* protects against plague. *Microbiology (Reading, England)* 151, 1919–1926. doi: 10.1099/mic.0.27959-0

Tourancheau, A., Mead, E. A., Zhang, X. S., and Fang, G. (2021). Discovering multiple types of DNA methylation from bacteria and microbiome using nanopore sequencing. *Nat. Methods* 18, 491–498. doi: 10.1038/s41592-021- 01109-3

van der Woude, M., Braaten, B., and Low, D. (1996). Epigenetic phase variation of the *pap* operon in *Escherichia coli*. *Trends Microbiol.* 4, 5–9. doi: 10.1016/0966- 842x(96)81498-3

van der Woude, M. W., and Baumler, A. J. (2004). Phase and antigenic variation in bacteria. *Clin. Microbiol. Rev.* 17, 581–611. table of contents, doi: 10.1128/cmr. 17.3.581-611.2004

van der Woude, M. W., Braaten, B. A., and Low, D. A. (1992). Evidence for global regulatory control of pilus expression in *Escherichia coli* by Lrp and DNA methylation: model building based on analysis of *pap*. *Mol. Microbiol.* 6, 2429–2435. doi: 10.1111/j.1365-2958.1992.tb01418.x

van der Woude, M. W., and Henderson, I. R. (2008). Regulation and function of *Ag43* (*flu*). *Ann. Rev. Microbiol.* 62, 153–169. doi: 10.1146/annurev.micro.62. 081307.162938

Vandenbussche, I., Sass, A., Van Nieuwerburgh, F., Pinto-Carbo, M., Mannweiler, O., Eberl, L., et al. (2021). Detection of cytosine methylation in *Burkholderia cenocepacia* by single-molecule real-time sequencing and whole-genome bisulfite sequencing. *Microbiology* 167. doi: 10.1099/mic.0.001027

Waldron, D. E., Owen, P., and Dorman, C. J. (2002). Competitive interaction of the OxyR DNA-binding protein and the Dam methylase at the antigen 43 gene regulatory region in *Escherichia coli*. *Mol. Microbiol.* 44, 509–520. doi: 10.1046/j.1365-2958.2002.02905.x

Wallecha, A., Munster, V., Correnti, J., Chan, T., and Van Der Woude, M. (2002). Dam- and OxyR-dependent phase variation of *agn43*: essential elements and evidence for a new role of DNA methylation. *J. Bacteriol.* 184, 3338–3347. doi: 10.1128/jb.184.12.3338-3347.2002

Watson, M. E., Jarisch, J., and Smith, A. L. (2004). Inactivation of deoxyadenosine methyltransferase (*dam*) attenuates *Haemophilus influenzae* virulence. *Mol. Microbiol.* 53, 651–664. doi: 10.1111/j.1365-2958.2004.04140.x
Weigel, W. A., and Dersch, P. (2018). Phenotypic heterogeneity: a bacterial virulence strategy. *Microbes Infection* 20, 570–577. doi: 10.1016/j.micinf.2018.01.008

Westphal, L. L., Sauvey, P., Champion, M. M., Ehrenreich, I. M., and Finkel, S. E. (2016). Genomewide dam methylation in *Escherichia coli* during long-term stationary phase. *mSystems* 1, e00130–e00216.

Wu, H., Lippmann, J. E., Oza, J. P., Zeng, M., Fives-Taylor, P., and Reich, N. O. (2006). Inactivation of DNA adenine methyltransferase alters virulence factors in *Actinobacillus actinomycetemcomitans*. *Oral Microbiol. Immunol.* 21, 238–244. doi: 10.1111/j.1399-302x.2006.00284.x

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Payelleville and Brillard. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.