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

HAL
open science

# DNA Adenine Methyltransferase (Dam) Overexpression Impairs Photorhabdus luminescens Motility and Virulence

Amaury Payelleville, Anne Lanois-Nouri, Marie Gislard, Emeric Dubois, David Roche, Stéphane Cruveiller, Alain Givaudan, Julien Brillard

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

## To cite this version:

Amaury Payelleville, Anne Lanois-Nouri, Marie Gislard, Emeric Dubois, David Roche, et al., DNA Adenine Methyltransferase (Dam) Overexpression Impairs Photorhabdus luminescens Motility and Virulence. Frontiers in Microbiology, 2017, 8, 14 p. 10.3389/fmicb.2017.01671. [http://dx.doi.org/10.3389/fmicb.2017.01671](http://dx.doi.org/10.3389/fmicb.2017.01671)

HAL Id: hal-01628662

https://hal.science/hal-01628662

Submitted on 3 Nov 2017

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_3.png" /></figure>

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

Check for
updates

# DNA Adenine Methyltransferase (Dam) Overexpression Impairs *Photorhabdus luminescens* Motility and Virulence

Amaury Payelleville¹, Anne Lanois¹, Marie Gislard², Emeric Dubois², David Roche³, Stéphane Cruveiller³, Alain Givaudan¹ and Julien Brillard¹*

¹Diversité, Génomes Interactions Microorganismes Insectes (DGIMI), Institut National De La Recherche Agronomique, Université de Montpellier, Montpellier, France, ²MGX-Montpellier Genomix, Institut de Génomique Fonctionnelle, Montpellier, France, ³Le Commissariat à l'énergie atomique et aux énergies alternatives (CEA), Genoscope, Université d'Evry, Centre National De La Recherche Scientifique-UMR8030, Université Paris-Saclay, Evry, France

## OPEN ACCESS

**Edited by:**
Michael Sauer,
University of Natural Resources and
Life Sciences, Vienna, Austria

**Reviewed by:**
Louis S. Tisa,
University of New Hampshire,
United States
Pablo Ivan Nikel,
The Novo Nordisk Foundation Center
for Biosustainability (DTU Biosustain),
Denmark
Marjan Van Der Woude,
University of York, United Kingdom

*Correspondence:
Julien Brillard
julien.brillard@umontpellier.fr

**Specialty section:**
This article was submitted to
Microbial Physiology and Metabolism,
a section of the journal
Frontiers in Microbiology

**Received:** 24 May 2017
**Accepted:** 17 August 2017
**Published:** 01 September 2017

**Citation:**
Payelleville A, Lanois A, Gislard M, Dubois E, Roche D, Cruveiller S, Givaudan A and Brillard J (2017) DNA Adenine Methyltransferase (Dam) Overexpression Impairs *Photorhabdus luminescens* Motility and Virulence. *Front. Microbiol.* 8:1671.
doi: 10.3389/fmicb.2017.01671

Dam, the most described bacterial DNA-methyltransferase, is widespread in gamma-proteobacteria. Dam DNA methylation can play a role in various genes expression and is involved in pathogenicity of several bacterial species. The purpose of this study was to determine the role played by the *dam* ortholog identified in the entomopathogenic bacterium *Photorhabdus luminescens*. Complementation assays of an *Escherichia coli dam* mutant showed the restoration of the DNA methylation state of the parental strain. Overexpression of *dam* in *P. luminescens* did not impair growth ability *in vitro*. In contrast, compared to a control strain harboring an empty plasmid, a significant decrease in motility was observed in the *dam*-overexpressing strain. A transcriptome analysis revealed the differential expression of 208 genes between the two strains. In particular, the downregulation of flagellar genes was observed in the *dam*-overexpressing strain. In the closely related bacterium *Xenorhabdus nematophila*, *dam* overexpression also impaired motility. In addition, the *dam*-overexpressing *P. luminescens* strain showed a delayed virulence compared to that of the control strain after injection in larvae of the lepidopteran *Spodoptera littoralis*. These results reveal that Dam plays a major role during *P. luminescens* insect infection.

**Keywords:** entomopathogenic bacterium, MTase, insect, pathogenicity, flagellar genes, RNA-seq

## INTRODUCTION

Enterobacteria of the genus *Photorhabdus* are highly pathogenic to insects and are also symbiotically associated to nematodes of the family *Heterorhabditidae*. After invasion of the insect host by the nematodes, bacteria are released into the hemocoel of the insect prey where they multiply until septicemia occurs, employing a broad range of virulence factors that kill the insect. The bacterial symbionts contribute to maintain suitable conditions for nematode reproduction. During the final stages of development, the bacteria and the nematode reassociate and subsequently leave the insect carcass in search of a new insect host (Nielsen-LeRoux et al., 2012). As described for many other microbial pathogens, which constantly alternate between their host and the compartment they disperse in Avery (2006), *Photorhabdus* displays phenotypic heterogeneity (Boemare and Akhurst, 1988; Somvanshi et al., 2012). Such phenomena provide
diverse phenotypes within a bacterial population, which increases the adaptive potential to rapidly changing environmental conditions (Avery, 2006; Grimbergen et al., 2015). For instance, the production of the *Photorhabdus luminescens* Mad (Maternal adhesion) pilus is controlled by a reversible DNA switch of the promoter of the *mad* operon (Somvanshi et al., 2012). Two forms can co-exist depending on the promoter orientation: the bacteria are in a M-form (for mutualist) when the Mad pilus is produced, because it allows the colonization of the nematode. When the Mad pilus is not produced, the bacteria are in P-form (for pathogenic). In addition, it has been recently demonstrated that the virulence strategy of *P. luminescens* involves the generation of a bacterial subpopulation which causes septicemia in insects by displaying resistance to cationic antimicrobial peptides (Mouammine et al., 2017). Other mechanisms of phenotypic heterogeneity within clonal bacterial cultures exist, such as bistability (Veening et al., 2008). Phenotypic heterogeneity involving a positive-feedback loop regulation of a transcriptional regulator has been reported in the closely related genus *Xenorhabdus*, where it controls the expression of motility and virulence determinants (Jubelin et al., 2013). DNA methylation is another process that has also been described as responsible for bacterial phenotypic heterogeneity (Casadesus and Low, 2013), but this has not been investigated in *Photorhabdus*.

Numerous DNA methyltransferases (MTases) are a component of restriction-modification (RE) systems that protect the bacterial cell from invasion by bacteriophages DNA (Marinus and Løbner-Olesen, 2014). Other MTases are not linked to restriction endonuclease and are classified as “orphan” MTases (Casadesus, 2016). The best characterized orphan MTase in bacteria is Dam (DNA adenine methyltransferase), originally identified in *Escherichia coli* and widespread among gamma-proteobacteria (Lobner-Olesen et al., 2005). Dam MTase transfers a methyl group to an adenosine localized in sites 5'-GATC-3' of the DNA. After DNA replication, the newly synthetized DNA strand is unmethylated and the GATC sites are therefore transiently hemimethylated (Wion and Casadesus, 2006; Marinus and Casadesus, 2009) until the action of a Dam protein. The Dam enzyme has key roles in bacterial genome maintenance, and *E. coli* *dam* mutants have defects in many important physiological processes such as DNA replication initiation, chromosome partitioning, nucleoid structure, and mismatch repair (Lobner-Olesen et al., 2005).

Changes in DNA methylation can also alter the affinity of regulatory proteins to their target DNA binding sequences (Hale et al., 1994; Tavazoie and Church, 1998). Conversely, DNA-binding proteins can represent a hindrance to MTase for reaching specific DNA sequences, and therefore can inhibit DNA methylation. Both mechanisms may lead to alterations in gene expression. Thus, DNA-methylation can play major roles in transcriptional regulation, including genes involved in bacterial virulence (Heusipp et al., 2007). An illustrative example is the Pap (pyelonephritis-associated pili) pilus phase variation of uropathogenic *E. coli*. The *pap* operon transcription depends on the binding of the transcriptional repressor Lrp. Such binding depends on the methylation state (on several

GATC sites) of the *pap* promoter because of a better Lrp affinity for non-methylated DNA (Blyn et al., 1990; Braaten et al., 1994). The regulation of gene that mediates adhesion to uroepithelial cells is therefore governed by DNA methylation state. In *E. coli*, Dam regulates transcription of several other pili operons (van der Woude and Low, 1994), and expression of a major outer-membrane protein (Ag43) (Henderson and Owen, 1999).

In several bacterial species that possess a *dam* gene, Dam has been described as an important virulence gene regulator. Mutants of *Salmonella* Typhimurium lacking the Dam enzyme are avirulent in mice (Garcia-Del Portillo et al., 1999; Heithoff et al., 1999). The impact of Dam inactivation on bacterial virulence has also been reported in *Haemophilus influenzae* (Watson et al., 2004), *Klebsiella* (Mehling et al., 2007), *Actinobacillus* (Wu et al., 2006) and *Yersinia pestis* (Robinson et al., 2005). In *Yersinia pseudotuberculosis*, *Vibrio cholera*, and *Aeromonas hydrophila*, inactivation of the *dam* gene was shown to be a lethal mutation (Julio et al., 2001; Erova et al., 2006b; Demarre and Chattoraj, 2010). However, plasmid-mediated overexpression of the *dam* gene in *Y. pseudotuberculosis* resulted in a virulence decrease in mice compared to wild-type (Julio et al., 2002) and in a defect in colonization of *V. cholerae* in a suckling mouse model compared to wild-type (Julio et al., 2001). Similarly, *dam*-overexpressing strains of *Salmonella*, *Pasteurella multocida*, or *A. hydrophila* were also highly attenuated in mice (Heithoff et al., 1999; Chen et al., 2003; Erova et al., 2006b).

Phenotypes associated with an alteration of the DNA methylation state have mostly been described in mammalian pathogens, but have not yet been reported in an insect-pathogenic bacterium. The gene plu0087 of the *P. luminescens* TT01 genome is annotated “DNA adenine methylase (Deoxyadenosyl-methyltransferase)” (Duchaud et al., 2003) and displays 70.7% of identity and 86.3% of similarity with the *E. coli* *dam* gene. The purpose of this study was to determine the role of deregulating the *dam* orthologous gene found in *P. luminescens* TT01. In this work, we confirmed the N<sup>6</sup>-Adenine methyltransferase function of the *dam* gene product. Investigation of several phenotypes of a *P. luminescens* dam-overexpressing strain revealed a major role for Dam in motility as well as during virulence in the insect.

# MATERIALS AND METHODS

## Strains and Growth Conditions

The bacterial strains and plasmids used in this study are listed in Table 1. *E. coli*, *P. luminescens* and *Xenorhabdus nematophila* cells were routinely grown in Luria broth (LB) medium with a 180 rpm agitation at 37 and 28°C, respectively. As required, antibiotic concentrations used for bacterial selection were gentamycin at 15 µg mL<sup>-1</sup>, rifampicin at 100 µg mL<sup>-1</sup> and kanamycin at 20 µg mL<sup>-1</sup>. IPTG was added at 0.2 mM when required.

## *In silico* Analysis

Primer sequence was designed using the Primer3 software (Untergasser et al., 2012). The REBASE database (Roberts et al., 2015) was used to identify a putative Dam methyltransferase,
<table><caption>**TABLE 1** | Strains and plasmids used in this work.</caption><thead><tr><th>Strain or plasmid</th><th>Relevant genotype and characteristics<sup>a</sup></th><th>Reference or source</th></tr></thead><tbody><tr><td colspan="3"><strong>STRAINS</strong></td></tr><tr><td>Photorhabdus luminescens TT01</td><td>Wild type</td><td>Duchaud et al., 2003</td></tr><tr><td>Escherichia coli XL1 blue MRF<sup>f</sup></td><td>Δ(mcrA)183 Δ(mcrCB-hsdSMR-mrr)173 endA1 supE44 thi-1 recA1 gyrA96 relA1 lac [F<sup>r</sup> proAB lacIqZΔM15 Tn10 (Tetr)]</td><td>Agilent Technologies</td></tr><tr><td>E. coli WM3064</td><td>thrB1004 pro thi rpsL hsdS lacZΔM15 RP4-1360Δ(araBAD)567 ΔdapA1341::[ferm pir (wt)]</td><td>Paulick et al., 2009</td></tr><tr><td>E. coli MG1655</td><td>Wild type</td><td>Lobner-Olesen and von Freiesleben, 1996</td></tr><tr><td>E. coli Dam::16KM</td><td>MG1655 dam16::Km<sup>r</sup></td><td>Lobner-Olesen and von Freiesleben, 1996</td></tr><tr><td>Xenorhabdus nematophila F1</td><td>Wild type</td><td>Lanois et al., 2013</td></tr><tr><td>Micrococcus luteus</td><td>Wild type</td><td>Pasteur Institute Culture collection, Paris, France</td></tr><tr><td colspan="3"><strong>PLASMIDS</strong></td></tr><tr><td>pBBR1MCS-5</td><td>Cloning vector, Gm<sup>r</sup></td><td>Kovach et al., 1995</td></tr><tr><td>P<sub>tet</sub>-MCS</td><td>Cloning vector, Km<sup>r</sup></td><td>Jubelin et al., 2013</td></tr><tr><td>P<sub>tet</sub>-MCS-Dam</td><td>853 pb PCR fragment (dam gene) inserted between KpnI and BamHI site of P<sub>tet</sub>-MCS plasmid</td><td>This study</td></tr><tr><td>pBB-Dam</td><td>864 pb fragment (RBS and dam gene) isolated from P<sub>tet</sub>-MCS-Dam and inserted between EcoRI and BamHI site of pBBR1MCS-5 plasmid</td><td>This study</td></tr><tr><td>pJQ200KS</td><td>Mobilizable vector, Gm<sup>r</sup></td><td>Quandt and Hynes, 1993</td></tr><tr><td>pJQ-Δdam</td><td>Region overlapping the dam gene disrupted by a Cm<sup>r</sup> cassette and inserted between PstI and XbaI site of pJQ200KS plasmid</td><td>This study</td></tr></tbody></table>

M.PluTDamP (REBASE Enzyme Number 7410), in the *P. luminescens* TT01 genome. Alignments between Dam from various organism were performed using the Multalin tool (Corpet et al., 1998).

## Nucleic Acid Manipulations

The extraction of plasmid DNA from *E. coli* was performed using the GenElute™HP Plasmid® miniprep purification kit as recommended by the manufacturer (Sigma). Chromosomal DNA was extracted from bacterial cells using the QIAamp DNA Mini kit (Qiagen). Restriction enzymes and T4 DNA ligase were used as recommended by the manufacturer (New England Biolabs and Promega, respectively). Oligonucleotide primers were synthesized by Eurogentec (Seraing, Belgium) and are listed in Table S1. PCR was performed in a T100 thermal cycler (Biorad) using the iProof high-fidelity DNA polymerase (Biorad). Amplified DNA fragments were purified using a PCR purification kit (Roche) and separated on 0.7% agarose gels after digestion as previously described (Brillard and Lereclus, 2007). Digested DNA fragments were extracted from agarose gels with a centrifugal filter device (DNA gel extraction kit, Millipore, Molsheim, France). All constructions were confirmed by DNA sequencing (Eurofins Genomics).

Cloning the *P. luminescens* dam gene was performed as follows. The plu0087 gene was PCR amplified using two primers mapping immediately upstream and downstream (Cp-plu0087-F and Cp-plu0087-R, respectively, Table S1) the 804 bp ORF (open reading frame), using the following cycling conditions: 98°C, 10 s; 56°C, 30 s; 72°C, 30 s for 35 cycles. The 853 bp-long amplified DNA fragment was then digested according to

the endonuclease sites introduced in the primers (KpnI and BamHI). Because the *P. luminescens* dam gene has no clear ribosome binding site (RBS), the generated DNA fragment was first inserted immediately downstream of the RBS already present in the plasmid P<sub>tet</sub>-MCS (Jubelin et al., 2013), to create the P<sub>tet</sub>-MCS-Dam (Table 1). A 864 bp-long fragment corresponding to this insert together with the RBS was then isolated from the P<sub>tet</sub>-MCS-Dam plasmid by EcoRI and BamHI endonucleases, and was inserted between the corresponding sites of the low-copy plasmid pBBR1MCS-5 (Kovach et al., 1995) downstream of the P<sub>lac</sub> promoter. The recombinant plasmid (pBB-dam) was introduced in *E. coli* strains by electroporation, or transferred in *P. luminescens* and *X. nematophila* by conjugative mating as previously described (Givaudan and Lanois, 2000). Transconjugants harboring the pBBR1MCS-5 empty plasmid were used as a control.

Attempts to construct a *dam* mutant were performed as follows. Briefly, DNA fragments of the plu0087 upstream (540 pb) and downstream (569 bp) regions were PCR-amplified using the primer pairs upF-plu0087/upR-plu0087 and dnF-plu0087/dnR-plu0087, respectively (Table S1). PCR products were digested with PstI/BamHI and BamHI/XbaI using the primer-incorporated restriction sites (Table S1). In parallel, the Ω interposon harboring a Cm<sup>r</sup> cassette was digested with BamHI, as previously described (Brillard et al., 2002). The three digested DNA fragments were purified, ligated in PstI/XbaI-digested pJQ200KS (Table 1), and introduced by electroporation in *E. coli* XL1. The resulting pJQ-Δdam plasmid was transferred in *P. luminescens* by conjugative mating. Four independent transconjugants clones were then subjected to allelic exchange

<footnote><sup>a</sup>Km, kanamycin; Gm, gentamicin; Cm, chloramphenicol.</footnote>
in LB at 28°C, following the protocol routinely used in the laboratory (Brillard et al., 2002; Derzelle et al., 2004b; Brugirard-Ricaud et al., 2005). Because several attempts were unsuccessful, the transconjugants were additionally subjected to allelic exchange in a M9 minimal medium instead of LB, or incubated at room temperature or at 15°C instead of 28°C.

## DNA Methylation State Analysis in *E. coli* *dam* Mutant Expressing the *P. luminescens* *dam* Gene

The methylation state of GATC sites was assessed in an *E. coli* *dam* mutant strain Dam 16::KM (Lobner-Olesen and von Freiesleben, 1996) harboring pBB-Dam or the pBBR1MCS-5 empty plasmid, in order to determine the *P. luminescens* *dam* functionality in *E. coli*. These strains were grown in LB supplemented with gentamycin (to maintain the plasmid) and IPTG (Isopropyl β-D-1-thiogalactopyranoside) to allow the activation of the *P*<sub>lac</sub> promoter controlling the expression of the *dam* gene. Fifty nanogram of plasmid DNA extracted from these strains were then digested during 2 h at 37°C by MboI (which digests only non-methylated GATC sites) or DpnI (which digests only methylated GATC sites) and DNA digestions were analyzed after electrophoresis on a 1% agarose gel.

## Phenotype Analysis of *P. luminescens*

Bromothymol blue adsorption was determined after growth on NBTA (nutrient agar supplemented with 25 mg of bromothymol blue and 40 mg of triphenyltetrazolium chloride per liter). It allows the identification of variant forms (Boemare and Akhurst, 1988). Antibiotic production was assessed by measuring antibacterial activity against *Micrococcus luteus* (Table 1). Hemolysis was determined by the observation of a clearing surrounding bacteria grown on standard sheep blood agar plates as previously described (Brillard et al., 2001). Bioluminescence production, lipase activity on Tween 20, 40, 60, 80, and 85 were also assessed as previously described (Boemare and Akhurst, 1988).

For motility assays, agar plates were prepared with LB broth supplemented with 0.35% agar and inoculated using 5 µL of cells grown in exponential phase (OD<sub>540 nm</sub> = 0.8), as previously described (Givaudan et al., 1995). The diameter of the halo size of swimming motility was measured 24 h and 30 h after incubation. Data from 3 independent experiments (with 10 plates used in each condition) were analyzed using Wilcoxon test.

Growth of *P. luminescens* was monitored with a TECAN automated turbidimetric system (Infinite M200 TECAN®). Estimation of maximum specific growth rate (µ<sub>max</sub>) was performed on 4 independent biological replicates for each strain, using serial dilution of the inoculum as previously described (Augustin et al., 1999).

The *P. luminescens* biofilm formation was determined as follows. Five milliliters of LB medium in glass tubes were inoculated at 10% with an overnight culture and incubated for 12 days at 28°C in static conditions. The tubes were then rinsed with PBS before the addition of 7 mL of Crystal violet solution at 0.01% (in PBS) to stain the biofilms during 15 min. Biofilms were rinsed with PBS and then dissolved 3 h in 7 mL ethanol. The

OD<sub>570 nm</sub> measurement allowed the quantification of the biofilm-associated crystal violet. Data from 3 independent experiments with replicates (totalizing 17 tubes for each strain) were analyzed using Wilcoxon test.

The *P. luminescens* spontaneous mutation rate was assessed by quantifying the emergence of rifampicin-resistant CFUs as follows. *P. luminescens* was grown overnight in 100 mL of LB medium supplemented with gentamycin (for plasmid maintenance) before plating on LB and LB with rifampicin. The mutation rate was calculated as the rifampicin-resistant population divided by the total population. Data from 3 independent experiments were compared using the Student t-test.

## Plasmid Curing of *P. luminescens* Strains

Plasmid curing was performed as follows. For each *P. luminescens* transconjugant strain, a fresh colony was used to inoculate 5 mL LB and incubated overnight with shaking in the absence of antibiotic pressure. These cultures were used to inoculate 100 mL fresh LB at an OD<sub>540 nm</sub> = 0.05, and incubated with shaking until OD<sub>540 nm</sub> = 0.8 was reached. These cultures were then diluted and spread on LB agar-plates prior incubation until CFU were visible. For each strain, 50 CFU were then streaked on LB Gm and LB without antibiotic in parallel to check for the plasmid stability. The loss of the pBB-Dam plasmid, as well as that of the pBBR1MCS-5, was observed for all the 50 tested CFU. Four Gm<sup>S</sup> clones from each strain were then tested for their motility ability after inoculation on low-agar plates as described above. Finally, one clone from each strain was tested for insect virulence.

## Insect Virulence Assay

The virulence-related properties of *dam*-overexpression were assessed by comparing the killing effect of *P. luminescens* transconjugants harboring either the pBB-Dam or the pBBR1MCS-5 empty plasmid during infection in the common cutworm *Spodoptera littoralis* as previously described (Brillard et al., 2002). Four to five independent pathogenicity assays were performed for each bacterial strain. Briefly, 20 µL of exponentially growing bacteria (OD<sub>540 nm</sub> = 0.3) diluted in LB corresponding to about 1 × 10<sup>4</sup> CFU (1.4 × 10<sup>4</sup>, mean value of 4 experiments or 1.2 × 10<sup>4</sup>, mean value of 5 experiments, for *P. luminescens* harboring pBBR1MCS-5 or pBB-Dam, respectively), were injected into the hemolymph of 20 fifth-instar larvae of *S. littoralis* reared on an artificial diet. Insect larvae were then individually incubated at 23°C. Altogether, the survival rate of 80–100 larvae for each bacterial strain were analyzed. The CFU of bacteria were determined by plating dilutions on LB agar. Insect death was monitored over time for up to 60 h. The time for killing 50% of the insect larvae (LT<sub>50</sub>) was calculated. Statistical analysis (Wilcoxon test) was performed as previously described (Givaudan and Lanois, 2000; Brillard et al., 2002) using SPSS version 14.0 (SPSS, Inc., Chicago, IL) to compare the mortality state.

## RNA Preparation

Total RNA extraction was performed on cells harvested at OD<sub>540 nm</sub> = 0.5, from nine independent cultures for each strain, using RNeasy miniprep Kit (Qiagen), according to the
manufacturer's instructions. An additional incubation step with DNase I (Qiagen) was performed. The quantity and quality of RNA were assessed with an Agilent 2100 Bioanalyzer with the RNA 6000 Nano LabChip kit. Lack of DNA contamination was controlled by carrying out a PCR on each RNA preparation.

## RNA Sequencing

The RNA-sequencing was performed as previously described (Mouammine et al., 2017) with the following changes. Equal amounts of total RNA from three independent samples per strain were pooled together to generate one final biological RNA sample per strain. Thus, from the initial nine independent RNA samples per strain, three final RNA samples were generated for each strain and subsequently treated as follows prior sequencing. Ribo-Zero rRNA Removal Kit Bacteria (illumina, San Diego, CA) was used to remove ribosomal RNA from 4 µg of total RNA. For each sample, 100 ng of rRNA-depleted RNA was used to construct sequencing libraries using Illumina's TruSeq Stranded mRNA Sample Prep Kit (Low throughput). The mRNA was chemically fragmented. The first cDNA strand was generated by reverse transcription with random hexamer primers, SuperScript IV Reverse Transcriptase (Life Technologies), Actinomycine D. The Second strand cDNA was synthesized by replacing dTTP with dUTP. A single "A" nucleotide was added to the 3' end and ligation was carried out with Illumina's indexed adapters. After 15 cycles of PCR, libraries were validated on a Fragment Analyzer (AATI, Ankeny, IA) and quantified with a KAPA qPCR kit. On a sequencing lane of a flowcell V4, nine libraries were pooled in equal proportions, denatured with NaOH and diluted to 8 pM before clustering. Clustering and 50 nt single read sequencing were performed according to the manufacturer's instructions. Cluster formation, primer hybridisation and single end-read 50 cycles sequencing were performed on cBot and HiSeq2500 (Illumina, San Diego, CA), respectively. Image analyses and baseCALLing were performed using the Illumina HiSeq Control Software and Real-Time Analysis component. Demultiplexing was performed using Illumina's conversion software (bcl2fastq 2.17). The quality of the data was assessed using FastQC from the Babraham Institute and the Illumina software SAV (Sequencing Analysis Viewer). Potential contaminants were investigated with the FastQ Screen software from the Babraham Institute.

## RNA-Seq Analysis

High-throughput transcriptomic sequencing data were processed with a bioinformatic pipeline implemented at the Microscope platform (Vallenet et al., 2013). The reads were mapped onto the *P. luminescens* subsp. *laumondi* TT01 genome sequence (EMBL accession number: BX470251) with BWA software (v. 0.7.4) (Li and Durbin, 2009). We then used SAMtools (v.0.1.12) (Lister et al., 2009) to lower the false-positive discovery rate and to extract reliable alignments from BAM-formatted files. The number of reads matching each genomic object harbored by the reference genome was then calculated with the Bioconductor-GenomicFeatures package (Lawrence et al., 2013). For reads matching several genomic objects, the count number was weighted so as to keep the total number of reads constant.

Finally, we used the Bioconductor-DESeq package (Anders and Huber, 2010) with default parameters to analyze raw count data, to normalize the samples to the reliable reads and to evaluate differential expression between conditions, as previously described (Jubelin et al., 2013). Between 14 and 19 million Illumina sequences (50-base reads) were obtained for each sample and between 80 and 93% of high-quality sequences mapped to at least one site in the reference genome. The complete dataset from this study has been deposited in NCBI's Gene Expression Omnibus (GEO) database, under accession number GSE100650.

## RT-qPCR Analysis

For the validation of RNA-seq data, quantitative reverse transcription-PCR (RT-qPCR) were carried out as previously described (Mouammine et al., 2017). Briefly, RNA samples from 9 biological replicates for each strain were used for cDNA synthesis. The SuperScript II reverse transcriptase (Invitrogen) was used on 1 µg of total RNA with random hexamers (100 ng/µl; Roche Diagnostics). qPCR analyses were performed using SYBR green Master kit (Roche Diagnostics) with 1 µl of cDNA and specific gene primers at 1 µM (Table S1). The reactions were performed in triplicate at 95°C for 10clones from each strain min, followed by 45 cycles at 95°C for 5 s, 61°C for 10 s, and 72°C for 15 s and monitored in the LightCycler 480 system (Roche). Melting curves were analyzed and always contained a single peak. The data analyzed with the REST software 2009 (Pfaffl et al., 2002) using the pairwise fixed randomization test with 2,000 permutations are presented as a ratio with respect to the reference housekeeping gene *gyrB*, as previously described (Jubelin et al., 2013).

# RESULTS

## Conservation of Major Amino Acids in *P. luminescens* Dam Protein

Functional characterization of Dam proteins has been described in several organisms, and allowed the identification of several amino-acids essential for the Dam function (Yang et al., 2003; Erova et al., 2006a; Horton et al., 2006). Comparing these protein sequences with that of M.PluTDamP by a multiple alignment revealed that all major amino acids described in these organisms were conserved in *P. luminescens* Dam protein (Figure S1). They were also conserved in a protein encoded by a *dam* ortholog (XNC3v2_1950011) (83.3% of identity and 92.6% of similarity with M.PluTDamP) found in the closely related bacterium X. nematophila (Lanois et al., 2013). These findings strongly suggest that M.PluTDamP plays a role in DNA adenine methylation, as described for other Dam proteins.

## Complementation by M.PluTDamp of an *E. coli* dam Mutant

In order to confirm the *P. luminescens* Dam function, the plu0087 gene was cloned and introduced in an *E. coli dam* mutant or in its *E. coli* parental strain. The pBBR1MCS-5 empty plasmid was used as a control. The recombinant strains harboring pBB-Dam or pBBR1MCS-5 were assessed for their ability to methylate GATC
sites. Plasmid DNA extracted from these strains were digested by enzymes sensitive to DNA methylation. Results presented in **Figure 1** show that the DNA extracted from the *E. coli* MG1655 strain harboring either the pBB-Dam or the empty plasmid was digested by DpnI, but not by MboI, as expected, indicating that the GATC sites of the DNA are methylated in this strain, whatever the plasmid introduced. However, the DNA extracted from the *E. coli dam* mutant harboring the control empty plasmid was digested by MboI, but not by DpnI, confirming that the GATC sites are not methylated in this mutant strain. In contrast, the DNA extracted from the *E. coli dam* mutant harboring the pBB-Dam was digested by DpnI, but not by MboI, revealing a methylation on GATC sites. This indicates that the *P. luminescens dam* gene was able to complement the *E. coli dam* mutant, and therefore confirms that it is a genuine *dam* ortholog, with M.PluTDamP being able to methylate adenine on GATC sites of DNA.

## dam Overexpression in *P. luminescens* Does Not Alter Growth Nor Several Major Phenotypes

Construction of a *P. luminescens dam*-mutant failed despite several attempts. The role of *P. luminescens dam* gene was therefore investigated by using a strain overexpressing *dam*. In *P. luminescens*, genes placed under the control of the *P*<sub>lac</sub> promoter are constitutively expressed (Lanois et al., 2011; Mouammine et al., 2014). Therefore, the additional copy of the *dam* gene caused by the presence of pBB-Dam plasmid, together with constitutive expression of the strong *P*<sub>lac</sub> promoter are supposed to induce a *dam* overexpression in *P. luminescens*. This postulate was confirmed by quantification of mRNA

corresponding to the *dam* gene. RT-qPCR experiments showed an average of 23.1-fold induction of expression of *dam* in *P. luminescens* harboring pBB-Dam when compared to the control strain (ie, harboring a pBBR1MCS-5 empty plasmid).

Considering that this *dam* overexpression may modify *P. luminescens* physiology, growth of both strains was monitored with an automated turbidimetric system, and the maximum specific growth rate (μ<sub>max</sub>) was estimated (**Figure 2**). The growth curves of both strains overlapped with the same shape: their slope were similar during the exponential phase and they reached the same maximum OD during stationary phase. No lag-phase was observed. Moreover, the calculated μ<sub>max</sub> were not different: 0.647 h<sup>-1</sup> for *P. luminescens* harboring pBB-Dam vs. 0.636 h<sup>-1</sup> for the control strain (p = 0.92, Student t-test). Several phenotypes were also assessed to compare the *dam*-overexpressing strain to the control strain. No significant difference was observed between *P. luminescens* harboring pBB-Dam when compared to the control *P. luminescens* strain for bromothymol blue adsorption on NBTA, bioluminescence production, antibiotic production, hemolysis, lipase activities and mutation rate (**Table 2**). Altogether, these findings indicate that many *P. luminescens* phenotypes are not altered by *dam* overexpression.

## *dam* Overexpression Impairs the Bacterial Motility

The halo size of motility, assessed on low agar LB medium, was much smaller in *P. luminescens* harboring pBB-Dam when compared to that of the control strain (**Figure 3**). Data revealed that for *P. luminescens* harboring pBB-Dam, the median halo size after 30 h of incubation was 59.3% that of the control strain

<figure><img src="image_5.png" /><figcaption>**FIGURE 1** | Differential plasmid DNA digestion from *E. coli*. An *E. coli dam* mutant (Dam::16KM) or its parental strain (MG1655) were complemented with either a plasmid harboring the *P. luminescens dam* gene (Dam) or the pBBR1-MCS5 empty vector (MCS5). Plasmid DNA was extracted and digested by Mbol (active on unmethylated 5'-GATC-3' sites) or DpnI (active on 5'-GmeATC-3' sites). pBBR1-MCS5 (4,768 bp) has 17 GATC sites and pBB-Dam (5,618 bp) has 19 GATC sites.</figcaption></figure>
<figure><img src="image_6.png" /></figure>

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

**FIGURE 2** | Growth comparison between *P. luminescens* overexpressing *dam* gene (red) and control harboring an empty plasmid (blue). (A) Growth curves in LB of the *P. luminescens* Dam-overexpressing strain and control strain (inoculated at 1.61 × 10⁵ and 1.86 × 10⁵ CFU/ml, respectively). Mean values ± SDs of at least 3 independent biological replicates for each strain are shown. (B) Two-fold serial dilutions of cultures containing 1.25 × 10⁴ CFU/ml for each strain have been performed in LB. Growth was quantified by absorbance at 600 nm and measured every 30 min. Four independent biological replicates for each strain were used. The time of growth detection was defined as an increase of 0.1 unit of absorbance at 600 nm, and recorded for each dilution. The Ln of the dilutions as a function of the time for growth detection is indicated. Similar slopes indicate similar growth rates (see Methods section for details). Differences were not significant (*p* = 0.92, Student *t*-test).

(17.5 vs. 29.5 mm, respectively). This difference was significant as early as 24 h after inoculation (*p* < 0.01, Wilcoxon test), and highly significant after 30 h of incubation (*p* < 0.001, Wilcoxon test). This indicates that motility is significantly reduced but not abolished in the *dam*-overexpressing *P. luminescens* strain.

Because introducing the pBB-Dam plasmid led to an impaired motility, we checked if the loss of the plasmid would restore the WT phenotype. After growth in a liquid medium followed by growth on agar plates, both in the absence of antibiotic pressure, the loss of the pBB-Dam plasmid, as well as that of the pBBR1MCS-5, was observed for all the 50 tested CFU of *P. luminescens* (not shown). Four Gm-sensitive CFU originating from each strain were then tested for their motility ability after inoculation on low-agar plates, as described above. No significant difference in the halo size of motility was observed between the cured strains originating either form the *dam*-overexpressing *P. luminescens* strain or from the control strain, after 30 h of incubation (*p* = 0.44, Wilcoxon test). These results revealed that motility was fully restored in pBB-Dam cured cells, confirming that *dam*-overexpression causes an impaired motility in *P. luminescens*.

We wondered if *dam* overexpression in another species closely related to *P. luminescens* would also cause an impaired motility. We therefore introduced the pBB-Dam plasmid, or the control empty plasmid, in *X. nematophila* and checked the motility on a low agar medium. Results from 5 independent experiments showed that the median halo size after 30 h of incubation was 13 mm for *X. nematophila* harboring pBB-Dam vs. 19 mm for the control strain. This difference was significant (*p* = 0.001, Wilcoxon test) indicating that swimming motility is reduced, but not abolished, in the *dam*-overexpressing *X. nematophila* strain, similarly as it was observed for *P. luminescens*.

## *dam* Overexpression Increases the Biofilm Formation Ability

Determination of the biofilms formed in glass tubes by the *P. luminescens* strain overexpressing *dam* and the control strain

was analyzed by a crystal violet staining method (Figure 4). The results revealed a significant increase of biofilm-associated crystal violet measured for the *P. luminescens* overexpressing *dam* strain (*p* = 0.003, Wilcoxon test), suggesting an increase in adhesion properties for this strain.

## Effect of *dam* Overexpression on *P. luminescens* Insect Virulence

The insect virulence of the *P. luminescens* strain overexpressing *dam* was compared to that of the control strain. It was assessed by injection of bacterial cells in *S. frugiperda* (Figure 5). Both strains were pathogenic, being able to cause death of all injected larvae in <60 h. However, while the time needed to kill 50% of infected larvae (LT₅₀) was about 38 h for the control strain, it was significantly increased in *P. luminescens* pBB-dam, reaching 44.5 h (*p* < 0.001, Wilcoxon test). This delay in killing insect larvae indicates a reduction in virulence properties of the *P. luminescens* strain overexpressing *dam*.

After loss of the *dam* overexpressing plasmid as described above, cured strains displayed a significantly lower LT₅₀ compared to that of the *dam* overexpressing *P. luminescens* strain (data not shown), confirming that *dam*-overexpression causes an impaired virulence in *P. luminescens*.

## Flagellar Genes Are Downregulated in the *P. luminescens* *dam*-Overexpressing Strain

We wondered if the reduced motility and the delayed virulence in insects observed for the *P. luminescens* strain overexpressing *dam* were associated to changes in gene expression. RNA sequencing was therefore performed on the *dam* overexpressing strain and the control strain during exponential phase of growth. The transcriptome analysis revealed significant differences (log2 fold change ≥1; adjusted *p* ≤ 0.005) in expression for 208 genes between the two strains, with 121 down-regulated and 87 up-regulated genes in the *P. luminescens* strain overexpressing *dam* (Table S2). Both up- and down-regulated genes were
<table><caption>TABLE 2 | Phenotypes of *P. luminescens* TT01 transconjugants, overexpressing *dam* gene (pBB-dam) and control (pBBR1MCS-5).</caption><thead><tr><th rowspan="3">Strain</th><th rowspan="3">Btb<br>adsorption<sup>b</sup></th><th rowspan="3">Bioluminescence<sup>c</sup></th><th rowspan="3">Antibiotic<br>production<sup>d</sup></th><th rowspan="3">Sheep blood<br>hemolysis<sup>e</sup></th><th rowspan="3">Motility<sup>f</sup></th><th colspan="5">Tested phenotypes<sup>a</sup></th></tr><tr><th colspan="2">Tween 20</th><th colspan="2">Tween 40</th><th colspan="2">Lipolysis of<sup>g</sup></th><th rowspan="2">Growth rate<br>(h<sup>-1</sup>)<sup>h</sup></th><th rowspan="2">Mutation<br>rate<sup>i</sup></th></tr><tr><th>Tween 20</th><th>Tween 40</th><th>Tween 60</th><th>Tween 80</th><th>Tween 85</th></tr></thead><tbody><tr><td>TT01 WT</td><td>G</td><td>+</td><td>+</td><td>—</td><td>++</td><td>+</td><td>+</td><td>+</td><td>—</td><td>ND</td><td>ND</td></tr><tr><td>TT01+pBBR1MCS-5</td><td>G</td><td>+</td><td>+</td><td>—</td><td>++</td><td>+</td><td>+</td><td>+</td><td>—</td><td>0.636</td><td>4.80 × 10<sup>-8</sup></td></tr><tr><td>TT01+pBB-dam</td><td>G</td><td>+</td><td>+</td><td>—</td><td>+</td><td>+</td><td>+</td><td>+</td><td>—</td><td>0.647</td><td>3.84 × 10<sup>-8</sup></td></tr></tbody></table>

<sup>a</sup>All plates were incubated for 2 days at 28°C before assays were interpreted, unless otherwise indicated. Routinely tested phenotypes on the WT strain are indicated for comparison.

<sup>b</sup>Btb, bromothymol blue; G, green-blue colonies on NBTA medium.

<sup>c</sup>+, Luminescence detected by visual observation in a dark room.

<sup>d</sup>+, Halo size (>25 mm) of growth inhibition of Micrococcus luteus.

<sup>e</sup>, No halo of hemolysis detected

<sup>f</sup>+, Large spreading area (halo size >20 mm); +, reduced spreading area (halo size < 20 mm) after 30 h of incubation.

<sup>g</sup>+, Halo of precipitation; —, no halo of precipitation

<sup>h</sup>Growth rate was estimated on 4 independent biological replicates, using serial dilution of the inoculum as described in the method section. Differences were not significant (p = 0.92, Student t-test). ND, not done.
(Spontaneous mutation rate per CFU was assessed by quantifying the number of rifampicin-resistant CFU arising in the total population after overnight growth in liquid medium. Presented values are the mean of 3 independent experiments. Differences were not significant (p = 0.73, Student t-test).

distributed all over the chromosome (Table S2). The percentages by COG class in TT01 genome revealed that the 208 differentially expressed genes belong to various COG categories (Figure 6). Interestingly, we observed an enrichment in genes encoding proteins belonging to the N category, and putatively involved in “cell motility” (Figure 6 and Table S2). Most of them are either flagellar genes or encode putative pili/fimbrial proteins (Table S2). When focusing on the 49 flagellar genes that are found in the *P. luminescens* TT01 genome, 48 of them displayed a lower expression in the strain overexpressing *dam* when compared to the control strain, although not always reaching a significant adjusted *p*-value to be considered as differentially expressed (Table S3). Interestingly, 16 genes encoding putative pili or fimbrial proteins were upregulated in the *dam*-overexpressing strain (Table S2), including *pilL* (plu1049), the first gene of an operon encoding a type IV pilus. This result may explain the observed increased biofilm formation in this strain.

However, while an impaired virulence was observed in the *P. luminescens* strain overexpressing *dam*, only a limited number of genes (*n* = 5) encoding effectors influencing the infection process were differentially regulated (Table S2). In order to confirm these observation, a few genes related to the two major impaired phenotypes (i.e., motility and insect virulence) were therefore selected and their level of expression was quantified by RT-qPCR in the *P. luminescens* strain overexpressing *dam* in comparison to the control strain. As mentioned above, the *dam* gene was up-regulated about 23-fold. Results also showed that 7 of the 10 tested flagellar genes were significantly downregulated in the *dam*-overexpressing strain (Figure S2), in agreement with the results of the RNA-seq analysis. This result confirmed that a reduced expression at the transcriptional level of flagellar genes was responsible for the impaired motility. In contrast, 5 additional genes (*manA*, *sodA*, *luxS*, *tcaZ*, *lopT*) likely involved in insect virulence displayed similar level of expression between the two strains (Figure S2). Finally, RT-qPCR confirmed the overexpression of a fimbrial gene (*madA*) in the *dam*-overexpressing strain (Figure S2).

DISCUSSION

Given the high degree of conservation of Dam methyltransferase among enterobacteria and several other Gram-negative bacteria (Lobner-Olesen et al., 2005; Casadesus and Low, 2013; Marinus and Løbner-Olesen, 2014), the identification of a *dam* ortholog in *P. luminescens* was expected. We first showed that the amino-acids previously described as essential for the Dam function in several organisms (Erova et al., 2006a; Horton et al., 2006) were conserved in the predicted amino-acid sequence of the plu0087 *dam* orthologous gene. This result strongly suggested that the predicted encoded enzyme, M.PluTDamP, has the same DNA methylation function in *P. luminescens* as that described for Dam proteins from other Gram-negative bacteria. The finding that the *P. luminescens dam* gene was able to complement the *E. coli dam* mutant, confirmed that M.PluTDamP was able to induce DNA-adenine methylation of GATC sites. Since our results confirmed the putative function of M.PluTDamP, the enzyme should now be named M.PluTDam (or M.PluTII) (Roberts et al., 2003).
<figure><img src="image_8.png" /></figure>

**FIGURE 3** | Swimming motility of *P. luminescens* overexpressing *dam* gene (Dam) and control (MCS5). (A) Swimming halos were observed on low agar LB medium inoculated by 5 µL of exponentially growing cells. (B) Boxplots of the diameter of the halo size of motility of each strain measured after 30 h of incubation (see Methods section for details). Difference between the two strains is significant (*p* < 0.001, Wilcoxon test).

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

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

**FIGURE 4** | Biofilm formation ability of *P. luminescens* overexpressing *dam* gene (Dam) and control (MCS5). (A) Biofilms formed in a glass tube after 12 days of incubation in LB medium were stained with crystal-violet. (B) Boxplots of the biofilm-associated crystal violet measured at OD<sub>570 nm</sub> (see Methods section for details). Difference between the two strains is significant (*p* < 0.01, Wilcoxon test).

The pleiotropic role of DNA-methylation by Dam has been illustrated in many bacterial species (Marinus and Casadesus, 2009). In some of them, it was proposed an essential function of *dam* for cell viability (Julio et al., 2001; Erova et al., 2006b). In particular, the role of Dam during bacterial-host interactions has been reported in several bacterial species (Heusipp et al., 2007), and mostly in the mammalian pathogen *Salmonella* (García-Del Portillo et al., 1999; Heithoff et al., 1999), but also in *E. coli*, *H. influenzae*, *Mycobacterium tuberculosis*, *Campylobacter*, *Klebsiella*, *Actinobacillus*, *Yersinia pseudotuberculosis*, and *Y. pestis* (Julio et al., 2001; Watson et al., 2004; Robinson et al., 2005; Wu et al., 2006; Mehling et al., 2007; Kim et al., 2008; Shell et al.,

2013). However, the role of DNA methylation in host-pathogen interaction remains unexplored in bacterial insect pathogens. Construction of a *P. luminescens* *dam*-mutant failed despite several attempts, suggesting that the presence of Dam itself is possibly required for some essential mechanisms, as described in other bacterial species (Julio et al., 2001; Erova et al., 2006b; Demarre and Chattoraj, 2010). We therefore investigated the role of Dam in *P. luminescens* by using a strain overexpressing *dam*. Growth in standard conditions together with the observation of most major phenotypes were found unmodified by the 23-fold increase of *dam* expression. Remarkably, two major phenotypes were impaired when compared to the control strain: motility and
virulence properties were significantly reduced in *P. luminescens* overexpressing *dam*. *P. luminescens* harbors a large repertoire of factors involved in bacterial-host interaction (Clarke, 2008),

and several genes have been described to contribute to insect virulence (Nielsen-LeRoux et al., 2012). For instance, isogenic mutants of genes encoding the SodA or the LuxS proteins display reduced virulence (Krin et al., 2006; Chalabaev et al., 2007). The Tc toxins and the type III secretion system (T3SS) are also considered to play key roles during *P. luminescens*-insect interactions (Bowen et al., 1998; Brugirard-Ricaud et al., 2005; Gatsogiannis et al., 2013). In addition, a *manA*-mutant is impaired in motility, insect virulence, but also in biofilm formation (Amos et al., 2011). Despite a significantly impaired virulence in *P. luminescens* overexpressing *dam*, no change in gene expression was detected for *manA*, *sodA*, *luxS*, *tcaZ*, *lopT* (encoding an effector of the T3SS), neither by RNA-seq nor RT-qPCR analysis. Thus, the precise factors involved in the observed impaired virulence of the *P. luminescens* strain overexpressing *dam* remains to be elucidated.

Our findings suggest that genome-wide alterations of methylation states may significantly impact some major phenotypes, although the specific mechanisms by which DNA methylation regulates the expression of the genes involved in these phenotypes remains unknown in *P. luminescens*. Strikingly, there is no GATC site in the promoter region of *flhD*, the gene encoding the master regulator of flagellar cascade. Therefore, the *flhD* downregulation (and consequently, of other flagellar genes) in the Dam overexpressing strain cannot be directly related to a difference in GATC methylation states in this locus. However, the observed downregulation of flagellar genes in *P. luminescens dam*-overexpressing strain opens new fields of investigation. The impaired motility coupled with an impaired ability to kill insects identified in this study may illustrate the occurrence of a direct mechanism, yet undescribed, involving these two

<figure><img src="image_11.png" /><figcaption>**FIGURE 5** | Infection of *Spodoptera littoralis* larvae by *P. luminescens* overexpressing *dam* and the control strain. Proportion of survival of *S. littoralis* after injection of 10⁴ CFU of *P. luminescens* overexpressing *dam* (TT01 + pBB-Dam, green) or carrying the vector control (TT01 + pBBR1-MCS5, blue). Graph represents the results from 4 to 5 independent experiments (with 20 insect larvae per experiment). The survival of 50% of the infested larvae is represented by the dotted line. The time needed to kill 50% of infected larvae (LT₅₀) was significantly different between the two strains (*p* < 0.001, Wilcoxon test).</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**FIGURE 6** | Classification by COG (cluster of orthologous group) annotation of the 208 genes differentially expressed between *P. luminescens* overexpressing *dam* and the control strain. Results show the percentage of genes from each COG class differentially expressed between the *P. luminescens* *dam*-overexpressing strain and the control strain, according to the 2014 update (ftp://ftp.ncbi.nih.gov/pub/COG/COG2014/static/lists/homeCOGs.html).</figcaption></figure>
phenotypes. Indeed, a complex interplay between the regulation
of flagellar motility and the expression of virulence factors exists
in several bacterial pathogens (Josenhans and Suerbaum, 2002).
In particular, these two phenotypes were previously shown as
being linked in the closely related genus *Xenorhabdus* (Givaudan
and Lanois, 2000, 2017). In this bacterium, the global regulators
encoded by the class I operon *flhDC*, controls the expression of
class II genes, including most of the structural genes for the
flagellar hook-basal body, but also the *fliAZ* operon encoding
the alternative sigma factor FliA and another flagellar regulator,
FliZ (Givaudan and Lanois, 2000, 2017; Park and Forst, 2006).
FliZ was found to upregulate many genes, including genes
encoding 2 hemolysins, and a toxin complex (Tc) protein, all of
them being considered as virulence factors (Lanois et al., 2008).
In *P. luminescens*, no mutant in flagellar regulators has been
described (Givaudan and Lanois, 2017). However, the mutation
of a response regulator of a two-component system (i.e., AstR)
in *P. luminescens* TT01, which regulates the *flhDC* transcription
level, causes an impaired motility but has no impact on virulence
in insects (Derzelle et al., 2004a). In addition, it was shown
that two distinct flagellar genes (*flgG* and *motAB*) deletion
mutants, that were consequently non-motile, were as efficient
as their parental strain in killing insects (Easom and Clarke,
2008). Our result revealed that *dam* overexpression causes
both a reduced motility linked to a downregulation of flagellar
genes (i.e., *flhDC*, *motAB...*) and an impaired virulence. Thus,
further investigations are required in *P. luminescens* in order to
determine if a common regulator (such as FlhD, FlhC, FliA, or
FliZ) is involved in both motility and the ability to kill insects, as
already described in *X. nematophila*. Then, it will be important to
determine if the changes in the global methylation state caused
by *dam*-overexpression involve such common regulator, or in
contrast if it impairs distinct mechanisms that are involved in
each of these two phenotypes.

A differential expression of several genes encoding virulence factors such as *pap*, *agn43*, *sci1* in *E.coli* (Blyn et al., 1990; Henderson and Owen, 1999; Brunet et al., 2011), but also *gtr*, or *opvAB* in *Salmonella* (Broadbent et al., 2010; Cota et al., 2016), is related to a differential methylation state of the GATC sites found in their respective promoters. This is caused mostly by differences in DNA affinity, depending on the DNA methylation state, of various transcriptional repressors (OxyR, Lrp, or Fur). In *X. nematophila*, virulence attenuation in insects was shown to be associated with an *lrp* mutation and Lrp positively regulates the expression of the *flhD* gene encoding the master flagellar regulator (Cowles et al., 2007; Lanois et al., 2008). Thus, because of the presence of *lrp* (plu1600) but also of *oxyR* (plu4740) and *fur* (plu1327) orthologs in the *P. luminescens* genome (Duchaud et al., 2003), it remains to be determined if similar mechanisms exist in *P. luminescens*.

It was shown that GATC methylation itself but also the
level of Dam, have multiple functions in the cell (Low
and Casadesus, 2008). These functions are correlated with
three DNA transactions: DNA mismatch repair, initiation of
chromosome replication, and regulation of gene expression
(Marinus and Løbner-Olesen, 2014). Our results confirm that
modifying the *dam* level of expression in *P. luminescens*
causes major phenotypes, presumably in the regulation of gene

expression. However, *dam* overexpression did not significantly
cause difference in the growth rate (Figure 2). In addition,
this overexpression did not lead to the identification of drug-
induced mutants in the tested conditions (Table 2). This is in
contrast to what was described in other species (Chen et al.,
2003) for which the general mutation rate is however similar
to what is observed in *P. luminescens*. Our result suggests that
after DNA replication, the Mismatch Repair (MMR) apparatus is
highly efficient and therefore still able to discriminate between
the nascent (error-containing) DNA strand and the mother
strand. Many genes putatively involved in regulation of the
chromosome replication, in MMR (including *mutH*, the gene
encoding the nuclease which relies on the DNA methylation
state to identify the correct strand), in double strand breaks
(DSBs) repair and/or required for viability of *dam* mutants in
other bacterial species are conserved in *P. luminescens* TT01
(Table S4). Such mechanisms may require the presence of Dam
for cell viability in *P. luminescens* during the growth conditions
tested. In contrast, our results suggest that an elevated level
of *dam* expression does not drastically affect the mechanisms
involved in chromosomal replication or in DNA mismatch
repair.

This study revealed that DNA methylation state on GATC sites
seems to be critical in *P. luminescens* for phenotypes involving
interactions with the insect host. Eleven additional MTases are
found in the genome and probably contribute to the global DNA
methylation state which may also account for some of the *P.
luminescens* phenotypes during bacterial-insect interaction.

ETHICS STATEMENT

According to the EU directive 2010/63, this study is exempt from
the above requirements because experiments were performed on
invertebrates animals (insects).

AUTHOR CONTRIBUTIONS

AP, AL, AG, and JB analyzed the data and designed the
experiments; AP performed the experiments; MG and ED
performed the RNA-seq; DR and SC analyzed the RNA-seq data;
JB drafted the manuscript. All authors revised the manuscript
and have approved its final version.

ACKNOWLEDGMENTS

The authors thank N. Ginibre, S. Pagès, and W. Levrat for
technical assistance and Philippe Clair from the Montpellier
GenomiX facility for assistance with real-time PCR. They are
grateful to Anders Løbner-Olesen for gift of E. coli strains.
The authors thank INRA Plant Health and Environment (SPE)
division for financial support (SPE2015-65359) and GAIA
doctoral school #584 for supporting AP's PhD.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found
online at: http://journal.frontiersin.org/article/10.3389/fmicb.
2017.01671/full#supplementary-material
REFERENCES

Amos, M. R., Sanchez-Contreras, M., Jackson, R. W., Munoz-Berbel, X., Ciche, T. A., Yang, G., et al. (2011). Influence of the *Photorhabdus luminescens* phosphomannose isomerase gene, manA, on mannose utilization, exopolysaccharide structure, and biofilm formation. *Appl. Environ. Microbiol.* 77, 776-785. doi: 10.1128/AEM.02326-10

Anders, S., and Huber, W. (2010). Differential expression analysis for sequence
count data. *Genome Biol.* 11:R106. doi: 10.1186/gb-2010-11-10-r106

Augustin, J. C., Rosso, L., and Carlier, V. (1999). Estimation of temperature dependent growth rate and lag time of *Listeria monocytogenes* by optical density measurements. *J. Microbiol. Methods* 38, 137-146. doi: 10.1016/S0167-7012(99)00089-5

Avery, S. V. (2006). Microbial cell individuality and the underlying sources of heterogeneity. *Nat. Rev. Microbiol.* 4, 577-587. doi: 10.1038/nrmicro1460

Blyn, L. B., Braaten, B. A., and Low, D. A. (1990). Regulation of *pap* pilin phase variation by a mechanism involving differential dam methylation states. *EMBO J.* 9, 4045-4054.

Boemare, N. E., and Akhurst, R. J. (1988). Biochemical and physiological characterization of colony form variants in *Xenorhabdus* spp. (*Enterobacteriaceae*). *J. Gen. Microbiol.* 134, 751-761. doi: 10.1099/00221287-134-3-751

Bowen, D., Rocheleau, T. A., Blackburn, M., Andreev, O., Golubeva, E., Bhartia, R., et al. (1998). Insecticidal toxins from the bacterium *Photorhabdus luminescens*. *Science* 280, 2129-2132. doi: 10.1126/science.280.5372.2129

Braaten, B. A., Nou, X., Kaltenbach, L. S., and Low, D. A. (1994). Methylation patterns in *pap* regulatory DNA control pyelonephritis-associated pili phase variation in *E. coli*. *Cell* 76, 577-588. doi: 10.1016/0092-8674(94)90120-1

Brillard, J., Duchaud, E., Boemare, N., Kunst, F., and Givaudan, A. (2002). The Phla hemolysin from the entomopathogenic bacterium *Photorhabdus luminescens* belongs to the two-partner secretion family of hemolysins. *J. Bacteriol.* 184, 3871-3878. doi: 10.1128/JB.184.14.3871-3878.2002

Brillard, J., and Lereclus, D. (2007). Characterization of a small PlcR-regulated gene co-expressed with cereolysin O. *BMC Microbiol.* 7:52. doi: 10.1186/1471-2180-7-52

Brillard, J., Ribeiro, C., Boemare, N., Brehelin, M., and Givaudan, A. (2001). Two distinct hemolytic activities in *Xenorhabdus nematophila* are active against immunocompetent insect cells. *Appl. Environ. Microbiol.* 67, 2515-2525. doi: 10.1128/AEM.67.6.2515-2525.2001

Broadbent, S. E., Davies, M. R., and van der Woude, M. W. (2010). Phase variation controls expression of *Salmonella* lipopolysaccharide modification genes by a DNA methylation-dependent mechanism. *Mol. Microbiol.* 77, 337-353. doi: 10.1111/j.1365-2958.2010.07203.x

Brugirard-Ricaud, K., Duchaud, E., Givaudan, A., Girard, P. A., Kunst, F., Boemare, N., et al. (2005). Site-specific antiphagocytic function of the *Photorhabdus luminescens* type III secretion system during insect colonization. *Cell. Microbiol.* 7, 363-371. doi: 10.1111/j.1462-5822.2004.00466.x

Brunet, Y. R., Bernard, C. S., Gavioli, M., Lloubes, 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

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

Casadesus, J., and Low, D. A. (2013). Programmed heterogeneity: epigenetic mechanisms in bacteria. *J. Biol. Chem.* 288, 13929-13935. doi: 10.1074/jbc.R113.472274

Chalabaev, S., Turlin, E., Charles, J. F., Namane, A., Pages, S., Givaudan, A., et al. (2007). The HcaR regulatory protein of *Photorhabdus luminescens* affects the production of proteins involved in oxidative stress and toxemia. *Proteomics* 7, 4499-4510. doi: 10.1002/pmic.200700609

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* 149, 2283-2290. doi: 10.1099/mic.0.26251-0

Clarke, D. J. (2008). *Photorhabdus:* a model for the analysis of pathogenicity and mutualism. *Cell. Microbiol.* 10, 2159-2167. doi: 10.1111/j.1462-5822.2008.01209.x

Corpet, F., Gouzy, J., and Kahn, D. (1998). The ProDom database of protein domain families. *Nucleic Acids Res.* 26, 323-326. doi: 10.1093/nar/26.1.323

Cota, I., Bunk, B., Sproer, C., Overmann, J., Konig, C., and Casadesus, J. (2016). OxyR-dependent formation of DNA methylation patterns in *OprABOFF* and *OprABON* cell lineages of *Salmonella enterica*. *Nucleic Acids Res.* 44, 3595-3609. doi: 10.1093/nar/gkv1483

Cowles, K. N., Cowles, C. E., Richards, G. R., Martens, E. C., and Goodrich-Blair, H. (2007). The global regulator Lrp contributes to mutualism, pathogenesis and phenotypic variation in the bacterium *Xenorhabdus nematophila*. *Cell. Microbiol.* 9, 1311-1323. doi: 10.1111/j.1462-5822.2006.00873.x

Demarre, G., and Chattoraj, D. K. (2010). DNA adenine methylation is required to replicate both *Vibrio cholerae* chromosomes once per cell cycle. *PLoS Genet.* 6:e1000939. doi: 10.1371/journal.pgen.1000939

Derzelle, S., Ngo, S., Turlin, E., Duchaud, E., Namane, A., Kunst, F., et al. (2004a). AstR-AstS, a new two-component signal transduction system, mediates swarming, adaptation to stationary phase and phenotypic variation in *Photorhabdus luminescens*. *Microbiology* 150, 897-910. doi: 10.1099/mic.0.26563-0

Derzelle, S., Turlin, E., Duchaud, E., Pages, S., Kunst, F., Givaudan, A., et al. (2004b). The PhoP-PhoQ two-component regulatory system of *Photorhabdus luminescens* is essential for virulence in insects. *J. Bacteriol.* 186, 1270-1279. doi: 10.1128/JB.186.5.1270-1279.2004

Duchaud, E., Rusniok, C., Frangeul, L., Buchrieser, C., Givaudan, A., Taourit, S., et al. (2003). The genome sequence of the entomopathogenic bacterium *Photorhabdus luminescens*. *Nat. Biotechnol.* 21, 1307-1313. doi: 10.1038/nbt886

Easom, C. A., and Clarke, D. J. (2008). Motility is required for the competitive fitness of entomopathogenic *Photorhabdus luminescens* during insect infection. *BMC Microbiol.* 8:168. doi: 10.1186/1471-2180-8-168

Erova, T. E., Fadl, A. A., Sha, J., Khajanchi, B. K., Pillai, L. L., Kozlova, E. V., et al. (2006a). Mutations within the catalytic motif of DNA adenine methyltransferase (Dam) of *Aeromonas hydrophila* cause the virulence of the Dam-overproducing strain to revert to that of the wild-type phenotype. *Infect. Immun.* 74, 5763-5772. doi: 10.1128/IAI.00994-06

Erova, T. E., Pillai, L., Fadl, A. A., Sha, J., Wang, S., Galindo, C. L., et al. (2006b). DNA adenine methyltransferase influences the virulence of *Aeromonas hydrophila*. *Infect. Immun.* 74, 410-424. doi: 10.1128/IAI.74.1.410-424.2006

Garcia-Del Portillo, F., Pucciarelli, M. G., and Casadesus, 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

Gatsogiannis, C., Lang, A. E., Meusch, D., Pfaumann, V., Hofnagel, O., Benz, R., et al. (2013). A syringe-like injection mechanism in *Photorhabdus luminescens* toxins. *Nature* 495, 520-523. doi: 10.1038/nature11987

Givaudan, A., Baghdiguian, S., Lanois, A., and Boemare, N. (1995). Swarming and swimming changes concomitant with phase variation in *Xenorhabdus nematophilus*. *Appl. Environ. Microbiol.* 61, 1408-1413.

Givaudan, A., and Lanois, A. (2000). FlhDC, the flagellar master operon of *Xenorhabdus nematophilus*: requirement for motility, lipolysis, extracellular hemolysis, and full virulence in insects. *J. Bacteriol.* 182, 107-115. doi: 10.1128/JB.182.1.107-115.2000

Givaudan, A., and Lanois, A. (2017). Flagellar Regulation and Virulence in the Entomopathogenic Bacteria-*Xenorhabdus nematophila* and *Photorhabdus luminescens*. *Curr. Top. Microbiol. Immunol.* 402, 39-51. doi: 10.1007/82_2016_53

Grimbergen, A. J., Siebring, J., Solopova, A., and Kuipers, O. P. (2015). Microbial bet-hedging: the power of being different. *Curr. Opin. Microbiol.* 25, 67-72. doi: 10.1016/j.mib.2015.04.008

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* 284, 967-970. doi: 10.1126/science.284.5416.967

Henderson, I. R., and Owen, P. (1999). The major phase-variable outer membrane protein of *Escherichia coli* structurally resembles the immunoglobulin A1
protease class of exported protein and is regulated by a novel mechanism involving Dam and OxyR. *J. Bacteriol.* 181, 2132–2141.

Heusipp, G., Falker, S., and Schmidt, M. A. (2007). DNA adenine methylation and bacterial pathogenesis. *Int. J. Med. Microbiol.* 297, 1–7. doi: 10.1016/j.ijmm.2006.10.002

Horton, J. R., Liebert, K., Bekes, M., Jeltsch, A., and Cheng, X. (2006). Structure and substrate recognition of the *Escherichia coli* DNA adenine methyltransferase. *J. Mol. Biol.* 358, 559–570. doi: 10.1016/j.jmb.2006.02.028

Josenhans, C., and Suerbaum, S. (2002). The role of motility as a virulence factor in bacteria. *Int. J. Med. Microbiol.* 291, 605–614. doi: 10.1078/1438-4221-00173

Jubelin, G., Lanois, A., Severac, D., Rialle, S., Longin, C., Gaudriault, S., et al. (2013). FliZ is a global regulatory protein affecting the expression of flagellar and virulence genes in individual *Xenorhabdus nematophila* bacterial cells. *PLoS Genet.* 9:e1003915. doi: 10.1371/journal.pgen.1003915

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*. *Infect. Immun.* 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. *Infect. Immun.* 70, 1006–1009. doi: 10.1128/IAI.70.2.1006-1009.2002

Kim, J. S., Li, J., Barnes, I. H., Baltzegar, D. A., Pajaniappan, M., Cullen, T. W., et al. (2008). Role of the *Campylobacter jejuni* Cj1461 DNA methyltransferase in regulating virulence characteristics. *J. Bacteriol.* 190, 6524–6529. doi: 10.1128/JB.00765-08

Kovach, M. E., Elzer, P. H., Hill, D. S., Robertson, G. T., Farris, M. A., Roop, R. M. 2nd, et al. (1995). Four new derivatives of the broad-host-range cloning vector pBBR1MCS, carrying different antibiotic-resistance cassettes. *Gene* 166, 175–176. doi: 10.1016/0378-1119(95)00584-1

Krin, E., Chakroun, N., Turlin, E., Givaudan, A., Gaboriau, F., Bonne, I., et al. (2006). Pleiotropic role of quorum-sensing autoinducer 2 in *Photorhabdus luminescens*. *Appl. Environ. Microbiol.* 72, 6439–6451. doi: 10.1128/AEM.00398-06

Lanois, A., Jubelin, G., and Givaudan, A. (2008). FliZ, a flagellar regulator, is at the crossroads between motility, haemolysin expression and virulence in the insect pathogenic bacterium *Xenorhabdus*. *Mol. Microbiol.* 68, 516–533. doi: 10.1111/j.1365-2958.2008.06168.x

Lanois, A., Ogier, J. C., Gouzy, J., Laroui, C., Rouy, Z., Givaudan, A., et al. (2013). Draft genome sequence and annotation of the entomopathogenic bacterium *Xenorhabdus nematophila* strain F1. *Genome Announc.* 1:e00342–13. doi: 10.1128/genomeA.00342-13

Lanois, A., Pages, S., Bourot, S., Canoy, A. S., Givaudan, A., and Gaudriault, S. (2011). *Transcriptional analysis of a Photorhabdus sp.* variant reveals transcriptional control of phenotypic variation and multifactorial pathogenicity in insects. *Appl. Environ. Microbiol.* 77, 1009–1020. doi: 10.1128/AEM.01696-10

Lawrence, M., Huber, W., Pages, H., Aboyoup, P., Carlson, M., Gentleman, R., et al. (2013). Software for computing and annotating genomic ranges. *PLoS Comput. Biol.* 9:e1003118. doi: 10.1371/journal.pcbi.1003118

Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. *Bioinformatics* 25, 1754–1760. doi: 10.1093/bioinformatics/btp324

Lister, R., Pelizzola, M., Dowen, R. H., Hawkins, R. D., Hon, G., Tonti-Filippini, J., et al. (2009). Human DNA methylomes at base resolution show widespread epigenetic differences. *Nature* 462, 315–322. doi: 10.1038/nature08514

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

Lobner-Olesen, A., and von Freiesleben, U. (1996). Chromosomal replication incompatibility in Dam methyltransferase deficient *Escherichia coli* cells. *EMBO J.* 15, 5999–6008.

Low, D. A., and Casadesus, J. (2008). Clocks and switches: bacterial gene regulation by DNA adenine methylation. *Curr. Opin. Microbiol.* 11, 106–112. doi: 10.1016/j.jmb.2008.02.012

Marinus, M. G., and Casadesus, J. (2009). Roles of DNA adenine methylation in host-pathogen interactions: mismatch repair, transcriptional regulation, and more. *FEMS Microbiol. Rev.* 33, 488–503. doi: 10.1111/j.1574-6976.2008.00159.x

Marinus, M. G., and Løbner-Olesen, A. (2014). DNA methylation. *EcoSal Plus* 6. doi: 10.1128/ecosalplus.ESP-0003-2013

Mehling, J. S., Lavender, H., and Clegg, S. (2007). A Dam methylation mutant of *Klebsiella pneumoniae* is partially attenuated. *FEMS Microbiol. Lett.* 268, 187–193. doi: 10.1111/j.1574-6968.2006.00581.x

Mouammine, A., Lanois, A., Pages, S., Lafay, B., Molle, V., Canova, M., et al. (2014). Ail and PagC-related proteins in the entomopathogenic bacteria of *Photorhabdus* genus. *PLoS ONE* 9:e110060. doi: 10.1371/journal.pone.0110060

Mouammine, A., Pages, S., Lanois, A., Gaudriault, S., Jubelin, G., Bonabaud, M., et al. (2017). An antimicrobial peptide-resistant minor subpopulation of *Photorhabdus luminescens* is responsible for virulence. *Sci. Rep.* 7:43670. doi: 10.1038/srep43670

Nielsen-LeRoux, C., Gaudriault, S., Ramarao, N., Lereclus, D., and Givaudan, A. (2012). How the insect pathogen bacteria *Bacillus thuringiensis* and *Xenorhabdus/Photorhabdus* occupy their hosts. *Curr. Opin. Microbiol.* 15, 220–231. doi: 10.1016/j.mib.2012.04.006

Park, D., and Forst, S. (2006). Co-regulation of motility, exoenzyme and antibiotic production by the EnvZ-OmpR-FlhDC-FliA pathway in *Xenorhabdus nematophila*. *Mol. Microbiol.* 61, 1397–1412. doi: 10.1111/j.1365-2958.2006.05320.x

Paulick, A., Koerdt, A., Lassak, J., Huntley, S., Wilms, I., Narberhaus, F., et al. (2009). Two different stator systems drive a single polar flagellum in *Shewanella oneidensis* MR-1. *Mol. Microbiol.* 71, 836–850. doi: 10.1111/j.1365-2958.2008.06570.x

Pfaffl, M. W., Horgan, G. W., and Dempfle, L. (2002). Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. *Nucleic Acids Res.* 30:e36. doi: 10.1093/nar/30.9.e36

Quandt, J., and Hynes, M. F. (1993). Versatile suicide vectors which allow direct selection for gene replacement in gram-negative bacteria. *Gene* 127, 15–21. doi: 10.1016/0378-1119(93)90611-6

Roberts, R. J., Belfort, M., Bestor, T., Bhagwat, A. S., Bickle, T. A., Bitinaite, J., et al. (2003). A nomenclature for restriction enzymes, DNA methyltransferases, homing endonucleases and their genes. *Nucleic Acids Res.* 31, 1805–1812. doi: 10.1093/nar/gk7274

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. doi: 10.1093/nar/gku1046

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

Shell, S. S., Prestwich, E. G., Baek, S. H., Shah, R. R., Sassetti, C. M., Dedon, P. C., et al. (2013). DNA methylation impacts gene expression and ensures hypoxic survival of *Mycobacterium tuberculosis*. *PLoS Pathog.* 9:e1003419. doi: 10.1371/journal.ppat.1003419

Somvanshi, V. S., Sloup, R. E., Crawford, J. M., Martin, A. R., Heidt, A. J., Kim, K. S., et al. (2012). A single promoter inversion switches *Photorhabdus* between pathogenic and mutualistic states. *Science* 337, 88–93. doi: 10.1126/science.1216641

Tavazoe, S., and Church, G. M. (1998). Quantitative whole-genome analysis of DNA-protein interactions by *in vivo* methylase protection in *E. coli*. *Nat Biotechnol* 16, 566–571. doi: 10.1038/nbt0698-566

Untergasser, A., Cutcutache, I., Koressar, T., Ye, J., Faircloth, B. C., Remm, M., et al. (2012). Primer3-new capabilities and interfaces. *Nucleic Acids Res.* 40:e115. doi: 10.1093/nar/gks596

Vallenet, D., Belda, E., Calteau, A., Cruveiller, S., Engelen, S., Lajus, A., et al. (2013). MicroScope—an integrated microbial resource for the curation and comparative analysis of genomic and metabolic data. *Nucleic Acids Res.* 41, D636–D647. doi: 10.1093/nar/gks1194

van der Woude, M. W., and Low, D. A. (1994). Leucine-responsive regulatory protein and deoxyadenosine methylase control the phase variation and expression of the *sfa* and *daa* pili operons in *Escherichia coli*. *Mol. Microbiol.* 11, 605–618. doi: 10.1111/j.1365-2958.1994.tb00340.x

Veening, J. W., Smits, W. K., and Kuipers, O. P. (2008). Bistability, epigenetics, and bet-hedging in bacteria. *Annu. Rev. Microbiol.* 62, 193–210. doi: 10.1146/annurev.micro.62.081307.163002

Watson, M. E. Jr., 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
Wion, D., and Casadesus, J. (2006). N6-methyl-adenine: an epigenetic signal for DNA-protein interactions. *Nat. Rev. Microbiol.* 4, 183–192. doi: 10.1038/nrmicro1350

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

Yang, Z., Horton, J. R., Zhou, L., Zhang, X. J., Dong, A., Zhang, X., et al. (2003). Structure of the bacteriophage T4 DNA adenine methyltransferase. *Nat. Struct. Biol.* 10, 849–855. doi: 10.1038/nsb973

**Conflict of Interest Statement:** 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 © 2017 Payelleville, Lanois, Gislard, Dubois, Roche, Cruveiller, Givaudan 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) or licensor 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.