# FepR as a Central Genetic Target in the Adaptation to Quaternary Ammonium Compounds and Cross-Resistance to Ciprofloxacin in Listeria monocytogenes

Pierre Emmanuel Douarre, Yann Sévellec, Patricia Le Grandois, Christophe Soumet, Arnaud Bridier, Sophie Roussel

► To cite this version:

Pierre Emmanuel Douarre, Yann Sévellec, Patricia Le Grandois, Christophe Soumet, Arnaud Bridier, et al.. FepR as a Central Genetic Target in the Adaptation to Quaternary Ammonium Compounds and Cross-Resistance to Ciprofloxacin in Listeria monocytogenes. Frontiers in Microbiology, 2022, 13, pp.864576. 10.3389/fmicb.2022.864576. anses-03707762

HAL Id: anses-03707762

https://anses.hal.science/anses-03707762

Submitted on 28 Jun 2022

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_1.png" /></figure>

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

Check for
updates

# FepR as a Central Genetic Target in the Adaptation to Quaternary Ammonium Compounds and Cross-Resistance to Ciprofloxacin in *Listeria monocytogenes*

Pierre-Emmanuel Douarre<fnref n="1" />†, Yann Sévellec<fnref n="1" />†, Patricia Le Grandois², Christophe Soumet², Arnaud Bridier²† and Sophie Roussel<fnref n="1" />*†

**OPEN ACCESS**

**Edited by:**

Catherine D. Carrillo,
Canadian Food Inspection Agency,
Canada

**Reviewed by:**

Ulas Tezel,
Boğaziçi University, Turkey
Xiaobing Jiang,
Henan Normal University, China

*Correspondence:

Sophie Roussel
sophie.roussel@anses.fr

**Specialty section:**
This article was submitted to *Food Microbiology*, a section of the journal *Frontiers in Microbiology*

**Received:** 28 January 2022

**Accepted:** 27 April 2022

**Published:** 18 May 2022

**Citation:** Douarre P-E, Sévellec Y, Le Grandois P, Soumet C, Bridier A and Roussel S (2022) FepR as a Central Genetic Target in the Adaptation to Quaternary Ammonium Compounds and Cross-Resistance to Ciprofloxacin in *Listeria monocytogenes*. *Front. Microbiol.* 13:864576.
doi: 10.3389/fmicb.2022.864576

The foodborne pathogen, *Listeria monocytogenes*, (*Lm*), frequently undergoes selection pressure associated with the extensive use of disinfectants, such as quaternary ammonium compounds, which are widely used in food processing plants. The repeated exposure to sub-inhibitory biocide concentrations can induce increased tolerance to these compounds, but can also trigger the development of antibiotic resistance, and both increase the risk of food contamination and persistence in food production environments. Although the acquisition of genes can explain biocide tolerance, the genetic mechanisms underlying the adaptive cross-resistance to antibiotics remain unclear. We previously showed that repeated exposure to benzalkonium chloride (BC) and didecyldimethyl ammonium chloride (DDAC) led to reduced susceptibility to ciprofloxacin in *Lm* strains from diverse sources. Here, we compared the genomes of 16 biocide-adapted and 10 parental strains to identify the molecular mechanisms of fluoroquinolone cross-resistance. A core genome SNP analysis identified various mutations in the transcriptional regulator *fepR* (Imo2088) for 94% of the adapted strains and mutations in other effectors at a lower frequency. FepR is a local repressor of the MATE fluoroquinolone efflux pump FepA. The impact of the mutations on the structure and function of the protein was assessed by performing *in silico* prediction and protein homology modeling. Our results show that 75% of the missense mutations observed in *fepR* are located in the HTH domain of the protein, within the DNA interaction site. These mutations are predicted to reduce the activity of the regulator, leading to the overexpression of the efflux pump responsible for the ciprofloxacin-enhanced resistance.

**Keywords:** *Listeria monocytogenes*, biocide adaptation, ciprofloxacin, cross-resistance, *fepR* gene, benzalkonium chloride, didecyldimethylammonium chloride

<footnote n="1">¹ Maisons-Alfort Laboratory for Food Safety, Salmonella and Listeria Unit, University of Paris-Est, French Agency for Food, Environmental and Occupational Health and Safety (ANSES), Maisons-Alfort, France, ² Antibiotics, Biocides, Residues and Resistance Unit, French Agency for Food, Environmental and Occupational Health and Safety (ANSES), Fougères, France</footnote>

<footnote>† These authors have contributed equally to this work</footnote>
<table><caption>**TABLE 1** | Description of the 10 *Listeria monocytogenes* parental strains and the corresponding strains adapted to the two biocides, benzalkonium chloride (BC), and didecyldimethylammonium chloride (DDAC), or exposed to water (referred to as Control).</caption><thead><tr><th rowspan="2">Strain</th><th rowspan="2">Listadapt id<sup>a</sup></th><th rowspan="2">Sourcecategory</th><th rowspan="2">Lineage</th><th rowspan="2">CC</th><th rowspan="2">CIP MIC<sup>c</sup>(mg/L)</th><th colspan="2">Biocide adaptation</th></tr><tr><th>DDAC</th><th>BC</th></tr></thead><tbody><tr><td>Parental</td><td>FR-DA-B-CH-295</td><td>F</td><td>II</td><td>CC21</td><td>1</td><td></td><td></td></tr><tr><td>Control</td><td>FR-BIO-T-500</td><td></td><td></td><td></td><td>1</td><td></td><td></td></tr><tr><td>Adapted</td><td>FR-BIO-DC-501</td><td></td><td></td><td></td><td>8</td><td>X</td><td></td></tr><tr><td>Parental</td><td>FR-ME-U-UN-465</td><td>F</td><td>I</td><td>CC2</td><td>2</td><td></td><td></td></tr><tr><td>Control</td><td>FR-BIO-T-502</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>FR-BIO-BC-503</td><td></td><td></td><td></td><td>8</td><td></td><td></td></tr><tr><td>Adapted</td><td>FR-BIO-DC-504</td><td></td><td></td><td></td><td>16</td><td>X</td><td>X</td></tr><tr><td>Parental</td><td>FR-ME-P-UN-410</td><td>F</td><td>II</td><td>CC37</td><td>2</td><td></td><td></td></tr><tr><td>Control</td><td>FR-BIO-T-507</td><td></td><td></td><td></td><td>1</td><td></td><td></td></tr><tr><td>Adapted</td><td>FR-BIO-DC-508</td><td></td><td></td><td></td><td>8</td><td>X</td><td></td></tr><tr><td>Parental</td><td>FR-VE-U-UN-470</td><td>F</td><td>II</td><td>CC37</td><td>2</td><td></td><td></td></tr><tr><td>Control</td><td>FR-BIO-T-505</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>FR-BIO-DC-506</td><td></td><td></td><td></td><td>8</td><td>X</td><td></td></tr><tr><td>Parental</td><td>CZ-NAT-SO-22</td><td>NE</td><td>I</td><td>CC1</td><td>2</td><td></td><td></td></tr><tr><td>Control</td><td>CZ-BIO-T-276</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>CZ-BIO-BC-277</td><td></td><td></td><td></td><td>8</td><td></td><td></td></tr><tr><td>Adapted</td><td>CZ-BIO-CD-278</td><td></td><td></td><td></td><td>8</td><td>X</td><td>X</td></tr><tr><td>Parental</td><td>NO-DEE-FE-3</td><td>A</td><td>II</td><td>CC11</td><td>1</td><td></td><td></td></tr><tr><td>Control</td><td>NO-BIO-T-40</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>NO-BIO-BC-41</td><td></td><td></td><td></td><td>4</td><td></td><td>X</td></tr><tr><td>Parental</td><td>CZ-OTH-UN-8</td><td>A</td><td>II</td><td>CC7</td><td>2</td><td></td><td></td></tr><tr><td>Control</td><td>CZ-BIO-T-282</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>CZ-BIO-BC-283</td><td></td><td></td><td></td><td>8</td><td></td><td>X</td></tr><tr><td>Adapted</td><td>CZ-BIO-CD-284</td><td></td><td></td><td></td><td>16</td><td>X</td><td></td></tr><tr><td>Parental</td><td>CZ-NAT-SO-15</td><td>NE</td><td>II</td><td>CC20</td><td>2</td><td></td><td></td></tr><tr><td>Control</td><td>CZ-BIO-T-286</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>CZ-BIO-BC-287</td><td></td><td></td><td></td><td>8</td><td></td><td>X</td></tr><tr><td>Adapted</td><td>CZ-BIO-CD-288</td><td></td><td></td><td></td><td>8</td><td>X</td><td></td></tr><tr><td>Parental</td><td>IT-FOX-FE-63</td><td>A</td><td>I</td><td>CC1</td><td>4</td><td></td><td></td></tr><tr><td>Control</td><td>IT-BIO-T-137</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>IT-BIO-BC-138</td><td></td><td></td><td></td><td>8</td><td></td><td>X</td></tr><tr><td>Adapted</td><td>IT-BIO-CD-139</td><td></td><td></td><td></td><td>8</td><td>X</td><td></td></tr><tr><td>Parental</td><td>NL-GOA-UN-2</td><td>A</td><td>II</td><td>CC26</td><td>2</td><td></td><td></td></tr><tr><td>Control</td><td>NL-BIO-T-65</td><td></td><td></td><td></td><td>2</td><td></td><td></td></tr><tr><td>Adapted</td><td>NL-BIO-BC-66</td><td></td><td></td><td></td><td>8</td><td></td><td>X</td></tr><tr><td>Adapted</td><td>NL-BIO-CD-67</td><td></td><td></td><td></td><td>8</td><td>X</td><td></td></tr></tbody></table>

were filtered because these mutations were not associated with the cross-resistance phenotype.

The annotated genome feature format (GFF) files produced by Prokka v1.13.3 (Seemann, 2014) were used to identify the pangenome of the 36 genomes through the Panaroo workflow (Tonkin-Hill et al., 2020) in strict mode for maximal contamination removal. Differences observed in the gene content between parental and control or adapted strains were characterized with a BLAST analysis.

## Prediction of *fepR* Protein Structure

FepR is a transcriptional regulator of the fluoroquinolone efflux pump (FepA). Based on the different FepR amino acid sequences, the protein structure of the regulator FepR was predicted by

homology using the SWISS-MODEL application (Waterhouse et al., 2018)<fnref n="4" /> based on the highest scoring template (Lmo2088 from EGD-e accession: 5ZTC on https://www.rcsb.org/structure/ 5ZTC) (Samad et al., 2019). Model quality was estimated as recommended in Studer et al. (2019).

Domains were predicted using the InterProScan (Jones et al., 2014)<fnref n="5" /> and PROSITE (De Castro et al., 2006)<fnref n="6" /> webtools. Ligand-binding pockets were predicted using P2rank software (Krivák and Hoksza, 2018). Ligand docking was evaluated using the COACH-D online tool

<footnote><sup>a</sup>The strains were isolated as part of the European project "LISTADAPT" (Adaptive traits of *Listeria monocytogenes* to its diverse ecological niches; https://onehealthejp.eu/jrp-listadapt/).</footnote>

<footnote>CZ, Czech Republic; NO, Norway; NL, the Netherlands; FR, France; IT, Italy.</footnote>

<footnote><sup>b</sup>F, Food; A, Animal; NE, Natural environment.</footnote>

<footnote><sup>c</sup>Ciprofloxacin minimum inhibitory concentration (CIP MIC).</footnote>

<footnote n="4">⁴https://swissmodel.expasy.org/interactive</footnote>

<footnote n="5">⁵https://www.ebi.ac.uk/interpro/search/sequence/</footnote>

<footnote n="6">⁶https://prosite.expasy.org/</footnote>
<table><caption>**TABLE 2** | Genes presenting mutations in the 16 adapted strains [7 strains adapted to benzalkonium chloride (BC) and 9 strains adapted to didecyldimethylammonium chloride (DDAC)].</caption><thead><tr><th>Gene</th><th>Product</th><th>EGD-e locus</th><th>Frequency<sup>a</sup></th><th>BC</th><th>DDAC</th></tr></thead><tbody><tr><td>-</td><td>Hypothetical protein</td><td>Absent in EGD-e</td><td>2</td><td>1</td><td>1</td></tr><tr><td>inJ</td><td>Internalin-J precursor</td><td>Absent in EGD-e</td><td>1</td><td>-</td><td>1</td></tr><tr><td>-</td><td>Hypothetical protein</td><td>Absent in EGD-e</td><td>1</td><td>-</td><td>1</td></tr><tr><td>-</td><td>Hypothetical protein</td><td>Absent in EGD-e</td><td>1</td><td>-</td><td>1</td></tr><tr><td>-</td><td>Hypothetical protein</td><td>Absent in EGD-e</td><td>2</td><td>1</td><td>1</td></tr><tr><td>-</td><td>Hypothetical protein</td><td>Absent in EGD-e</td><td>1</td><td>1</td><td>-</td></tr><tr><td>nagC</td><td>N-acetylglucosamine repressor</td><td>Imo0178</td><td>1</td><td>-</td><td>1</td></tr><tr><td>ftsH</td><td>ATP-dependent zinc metalloprotease FtsH</td><td>Imo0220</td><td>1</td><td>-</td><td>1</td></tr><tr><td>walR</td><td>Transcriptional regulatory protein WalR</td><td>Imo0287</td><td>3</td><td>2</td><td>1</td></tr><tr><td>walK</td><td>Sensor histidine kinase WalK</td><td>Imo0288</td><td>1</td><td>1</td><td>-</td></tr><tr><td>manZ</td><td>PTS system mannose-specific EIID component</td><td>Imo0781</td><td>1</td><td>1</td><td>-</td></tr><tr><td>rsbS</td><td>RsbT antagonist protein RsbS</td><td>Imo0890</td><td>2</td><td>2</td><td>-</td></tr><tr><td>rsbU</td><td>Phosphoserine phosphatase RsbU</td><td>Imo0892</td><td>3</td><td>2</td><td>1</td></tr><tr><td>mntH</td><td>Divalent metal cation transporter MntH</td><td>Imo1424</td><td>1</td><td>-</td><td>1</td></tr><tr><td>murG</td><td>UDP-N-acetylglucosamine-N-acetylmuramyl- (pentapeptide) pyrophosphoryl-undecaprenol N-acetylglucosamine transferase</td><td>Imo2035</td><td>1</td><td>1</td><td>-</td></tr><tr><td>fepR</td><td>HTH-type transcriptional regulator AcrR</td><td>Imo2088</td><td>15</td><td>6</td><td>9</td></tr><tr><td>-</td><td>Internalin</td><td>Imo2396</td><td>1</td><td>-</td><td>1</td></tr><tr><td>kdpD</td><td>Sensor protein KdpD</td><td>Imo2679</td><td>1</td><td>-</td><td>1</td></tr></tbody></table>

<table><caption>**TABLE 3** | Amino acid changes in the FepR protein.</caption><thead><tr><th>Listadapt id</th><th colspan="16">Position/amino acid change in FepR</th><th>Biocide adaptation</th></tr><tr><th>EGD-e (Imo2088)</th><th>19</th><th>22</th><th>27</th><th>37</th><th>44</th><th>46</th><th>100</th><th>103</th><th>108</th><th>116</th><th>121</th><th>130</th><th>157</th><th>168</th><th></th></tr><tr><th></th><th>Phe</th><th>Thr</th><th>Ile</th><th>Ser</th><th>Asp</th><th>Asp</th><th>Trp</th><th>Tyr</th><th>Pro</th><th>Gln</th><th>Glu</th><th>Phe</th><th>Gly</th><th>Tyr</th><th></th></tr></thead><tbody><tr><td>FR-BIO-DC-501</td><td></td><td>Ile</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>FR-BIO-BC-503</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Stop</td><td></td><td>BC</td></tr><tr><td>FR-BIO-DC-504</td><td></td><td></td><td></td><td></td><td></td><td></td><td>FS</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>FR-BIO-DC-508</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Stop</td><td></td><td></td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>FR-BIO-DC-506</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Stop</td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>CZ-BIO-BC-277</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Tyr</td><td></td><td></td><td></td><td></td><td></td><td></td><td>BC</td></tr><tr><td>CZ-BIO-CD-278</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Stop</td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>NO-BIO-BC-41</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Leu</td><td></td><td></td><td></td><td>BC</td></tr><tr><td>CZ-BIO-BC-283</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>BC</td></tr><tr><td>CZ-BIO-CD-284</td><td>Leu</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>CZ-BIO-BC-287</td><td></td><td></td><td></td><td></td><td></td><td>Gly</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>BC</td></tr><tr><td>CZ-BIO-CD-288</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Stop</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>IT-BIO-CD-139</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Tyr</td><td></td><td></td><td></td><td></td><td></td><td></td><td>DDAC</td></tr><tr><td>NL-BIO-BC-66</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Leu</td><td></td><td></td><td></td><td></td><td>BC</td></tr><tr><td>NL-BIO-CD-67</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Stop</td><td></td><td></td><td>DDAC</td></tr></tbody></table>

Positions in orange correspond to the HTH domain and position in green locates in the vicinity of the ligand-binding pocket. Missense mutations are indicated in blue and the frameshift (FS) and premature stop codon (Stop) are indicated in red.

(Wu et al., 2018). To further explore the docking site, the secondary structure of the fepA promoter region (76 bp) extracted from EGD-e sequence (NC_003210.1) was generated using the DNA Sequence to Structure tool (available at http://www.scfbio-iitd.res.in/software/drugdesign/bdna.jsp) and

docked on the helix-turn-helix (HTH) domain of FepR using the HADDOCK web server (V. 2.4) (Van Zundert et al., 2016). The protein surface was visualized using pymol.<fnref n="7" />

<footnote><sup>a</sup>The frequency is defined here as the number of strains showing a mutation in this gene.</footnote>

<footnote n="7">⁷https://pymol.org/2/</footnote>
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Structure of the FepR protein. The carbon of the side chain of the mutated residues inside the HTH domain are shown in pink. Oxygen atoms are shown in red and nitrogen in blue. Hydrogen atoms are hidden.</figcaption></figure>

# RESULTS

## Single Nucleotide Polymorphism Identification

The comparison of the whole genome sequences between the 16 adapted strains and their parental strains identified a total of 37 SNPs and 3 deletion sequences in 18 genes as well as 10 SNPs and two insertion sequences in 11 intergenic regions (Supplementary Table 2). Overall, 18 and 21 mutations were observed in the strains adapted to BC and DDAC, respectively, and the number of mutations in each strain varied from one (CZ-BIO-DC-288) to seven mutations (FR-BIO-DC-501).

Among the 18 mutated genes, 13 encoded proteins of known functions and five encoded hypothetical proteins (Table 2). Our results revealed that mutations in the *fepR* gene were detected in

15 out of the 16 adapted strains, in contrast to the other genes that were identified at much lower frequencies. The *fepR* gene encodes a *tetR*-like transcriptional regulator that represses the activity of the multidrug efflux pump FepA. The 15 mutations were spread out in 14 different positions in *fepR*. Six of these were non-sense mutations resulting in the insertion of a premature stop codon (PMSC) (six strains), one was a frameshift mutation (one strain), while missense mutations were identified for the eight remaining strains (Table 3).

The BC-adapted strain IT-BIO-BC-138 that did not display mutation in the *fepR* gene contained a single missense mutation (319A319T/Asn107Tyr) in the transcriptional regulator *walR*. This gene encodes a two-component response regulator involved in membrane metabolism and antibiotic resistance. Mutations in this gene were identified in two other strains (NL-BIO-BC-66 and NL-BIO-DC-67) (Table 2).

The other genes in which mutations were observed were mainly stress/antibiotic-resistance-related genes (*rbsS*, *rsbU* and *kdpD*) and genes involved in cell wall metabolism (*nagC*, *manZ*, *murG*). Some of these genes were only detected in BC- (*walK*, *manZ*, *rsbS*, and *murG*) or in DDAC-adapted (*nagC* and *kdpD*) strains while five mutated genes (*walR*, *rbsU*, *fepR*, and two genes encoding hypothetical proteins) were identified in both (Table 2 and Supplementary Table 2).

## Investigation of *fepR* Mutations

The impact of the eight missense mutations on the structure and function of the protein was analyzed by performing *in silico* prediction and protein homology modeling.

The quality parameters are reported in Supplementary Table 3. Most mutations (6 out of 8) were located in the HTH domain of the protein located in residues 19, 22, 27, 37, 44, and 46 and are predicted to significantly modified the protein surface in this area (Figure 1). These changes potentially reduce the capacity of FepR to bind DNA leading to an upregulation of the *fepA* gene.

Based on the main allele of *fepR*, DNA interaction was predicted to occur at two locations inside the HTH domain: between residues 1, 23, 24, 35, 39, 44, and 45 and between residues 2, 33–37, 40, and 41 (confidence score 1) (Figure 2). Docking was attempted with a dsDNA sequence extracted from the promoter region of EGD-e. The best binding pose predicted that the interacting residues were 33, 34, 36, 37, 39, 40, and 41 (Figure 2A). The mutations detected in residues 19, 22, 27, 37, 44, and 46 may have a significant impact on protein binding to DNA at the recognition sequence reducing the efficiency of this regulator (Figure 2B). The docking analysis showed a reduction of the HTH domain affinity for the *fepA* promoter sequence for these six mutations (affinity score from −95.2 to −82.9) in comparison to the affinity score of −104.2 for the native *fepR* structure (Supplementary Table 3). To investigate the affinity of FepR to QACs, pocket prediction was performed and a pocket was predicted in the vicinity of neighboring residues 60, 100, 101, 104, 105, 119, 123, 126, 156, 159, 160, and 163 with an affinity of −6.3 for BC and DDAC (Figure 3). The two remaining mutations did not seem to affect the binding of these ligands either because
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Impact of the Ser37Tyr mutation observed in the DDAC-adapted IT-BIO-DC-139 strain on DNA binding. Carbon atoms are shown in green for the parental strain (IT-FOX-FE-63) (A) and in blue for the DDAC-adapted strain IT-BIO-DC-139 (B). The calculated hydrogen bonds are indicated by a yellow dotted line. Only the scaffold is shown for the double-stranded DNA (dsDNA) molecule. The dsDNA is docked in the parental strain; the modification of the residue 37 alters the conformation of the protein at the DNA-binding site, leading to a reduction of affinity for the DNA strands (affinity score of −82.9 compared with −104.2 from the *fepR* model of the parental strain). (DDAC, didecyldimethylammonium chloride).

the mutation is located far from the pocket (Pro108Leu) or because it did not seem to affect the affinity for BC or DDAC (Phe130Leu).

## Pangenome

The pangenome of the control and adapted strains was compared with that of the parental strains to detect potential variations in gene content. The different presence-absence matrices and summary statistics generated by Panaroo are available in **Supplementary Table 4**.

A deletion of 62 genes was observed in the BC-adapted strain FR-BIO-DC-506. This 41 kb region shares 93.5% identity (72% coverage) with the *L. monocytogenes* prophage LP-101 (NC_024387.1) (Denes et al., 2014). The prophage identified in this study is integrated in the parental and control genomes at the tRNA-Arg locus. The analysis also revealed that 83 genes were missing in the adapted strain CZ-BIO-BC-287. The element carrying these genes displayed a strong homology (>99% identity and coverage) with the pLmN1546 plasmid isolated from the *Listeria monocytogenes* strain Lm N1546 (CP013725.1) (Tasara et al., 2016). This 86 kb plasmid encoding the cadmium-resistance transcriptional regulatory protein CadC was present in all the other strains derived from the parental strain. A similar plasmid was also lost by the control and BC-adapted strains NL-BIO-T-65 and NL-BIO-BC-66, respectively. The loss of this plasmid in the control strain indicates it is not linked with the resistance phenotype. A few additional missing genes encoding hypothetical protein, transposase and internalin are indicated in **Supplementary Table 4**.

## DISCUSSION

CIP belongs to the fluoroquinolone family, a family of antibiotics extensively used in human medicine. Different types of selective pressure (such as antibiotics, biocides, or heavy metals) may play a role in the prevalence of antibiotic resistance of *Lm* in the food chain. Exposure to disinfectant biocides such as the QACs BC and DDAC, extensively used in the food processing industry, has the potential to promote cross-resistance to CIP, as indicated by an increasing number of studies (Rakic-Martinez et al., 2011; Yu et al., 2018; Guérin et al., 2021; Kode et al., 2021; Bland et al., 2022). These biocides may co-select for strains resistant to CIP, due to shared mechanisms of resistance (cross-resistance) (Kampf, 2019). However, very limited information is currently available on the genetic mechanisms underlying CIP resistance following exposure to QACs. Here, we studied this issue by analyzing genetic traits linked to this cross-resistance.

The most striking finding of this study was the recurrence of the mutations observed in *fepR*, the multi-antimicrobial extrusion (MATE) transporter (*fepA*) regulator for all the CIP resistant DDAC-adapted strains and for six out of the seven CIP resistant BC-adapted strains. This is consistent with very recent findings showing that mutations in *fepR* play a role in the reduction in antibiotic susceptibility including CIP following low-level adaptation to a commercial sanitizer containing QAC (Bland et al., 2022). Interestingly, overexpression of the MATE efflux pump FepA through inactivation of the *fepR* repressor has been described to be responsible for fluoroquinolone resistance in *Lm* (Guérin et al., 2014). The involvement of efflux pump systems in adaptation to QAC and CIP resistance found here is also consistent with previous studies reporting their involvement
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Docking of benzalkonium chloride (pink) in the predicted pocket of FepR using the P2Rsankweb application. The predicted interacting residues are shown in violet. Oxygen atoms are shown in red, nitrogen in blue, sulfur in yellow.</figcaption></figure>

in cross-resistance between various biocides (especially QAC) and antibiotics in various pathogenic Gram-positive and Gram-negative species including *Lm*, *Staphylococcus aureus*, *Escherichia coli*, and *Salmonella enterica* (Buffet-Bataillon et al., 2016). Nevertheless, this study is the first to report the *fepR* gene as a central target in the adaptation to DDAC and cross-resistance development to CIP in *Lm*.

The FepR protein contains an HTH domain in aa residues 1–60 (predicted through InterProScan and PROSITE) and a tetracyclin repressor-like C-terminal domain in residues 99–192. The analysis of the protein structure revealed that 75% of the missense mutations were located in the DNA interaction site within the HTH domain thus reducing the affinity of FepR for the *fepA* promoter region and leading to the overexpression of the efflux pump FepA responsible for the reduced susceptibility to CIP.

The mutation Glu21PMSC observed by Guérin et al. (2014) in a clinical *Lm* isolate was also detected inside the DNA interaction site. qRT-PCR analysis demonstrated the overexpression (64-fold higher) of the MATE efflux pump *fepA* in the mutant strain and a clean *fepR* deletion mutant,

confirming the role of FepR as a local repressor of FepA (Guérin et al., 2014). Another mutation has been reported by Wilson et al. (2018) in one CIP-resistant CC1 strain, isolated from food production supplies in Australia (Wilson et al., 2018). This mutation was located at nucleotide position 181 and resulted in a PMSC in place of aspartic acid (Asp61PSMC). More recently, Bland et al. (2022) identified five mutations (His41Tyr/Leu96PMSC/Pro107Gln/Pro108Ser/Glu176PMSC) in five strains belonging to three hypervirulent CCs (Bland et al., 2022). All these mutations are different from those found in this study. Our results also demonstrate that *fepR* mutations occurred independently of the genetic background (CC) or the source, as the 15 strains developing cross-resistance were very diverse.

In addition to the recurrent mutations observed in *fepR*, other genes were mutated in QAC-adapted CIP-resistant strains. Mutations in the *walR* and *walK* genes were identified in four QAC-adapted CIP resistant strains. The two-component system is known to regulate cell wall metabolism, virulence regulation, biofilm production, oxidative stress resistance, and antibiotic resistance via direct or indirect regulation of autolysins (Takada and Yoshikawa, 2018). Interestingly, the BC-adapted strain IT-BIO-BC-138 was the only strain that did not display mutation in *fepR* but a missense mutation in *walR*. This result suggests that such mutation is responsible for the CIP resistance observed in IT-BIO-BC-138 and that WalR regulator could be another target for QAC adaptive CIP cross-resistance selection.

Some QAC-adapted CIP-resistant strains also displayed mutations in the *rsbS* or *rsbU* genes, both involved in the signal transduction responsible for the regulation of alternative sigma factor $\sigma^B$, which is known as a mediator of the general stress response in *Lm* (Guerreiro et al., 2020). These mutations may thereby participate in the CIP resistance observed in QAC-adapted strains through the upregulation of stress response factors acting on various cellular processes including cell envelope modification, pH homeostasis, osmoregulation, quorum sensing, or antibiotic resistance (Guerreiro et al., 2020).

## CONCLUSION

In conclusion, several genetic markers associated with CIP cross-resistance in QAC-resistant Listeria have been identified. In particular, the DNA binding region of the FepR regulator appears to be an important target in CIP cross-resistance. In the present study, the mutations observed in the FepR HTH domain are predicting to affect the promoter activity of *fepA* and resulting in the deregulation of the efflux pump FepA, responsible for the cross-resistance phenotype. Other stress-related genes such as the transcriptional regulator WalR may contribute to the CIP resistance mechanism and will need further attention.

## DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.
## AUTHOR CONTRIBUTIONS

P-ED contributed substantially to data analysis, quality control of the genomic dataset, writing, and editing the manuscript. YS contributed to study design, acquisition of strains and the corresponding genomes, data analysis as well as writing, and editing the manuscript. PL provided technical assistance. CS participated in study design and contributed to writing and editing of the manuscript. AB provided technical assistance and participated in designing the study and in writing and revising the manuscript. SR designed, coordinated the study, and contributed to writing and editing the manuscript. All authors read and approved the final manuscript.

## FUNDING

This work was supported by the One Health European Joint Programme- European Union's Horizon 2020 Research and Innovation Programme (Grant Agreement No 773830).

## REFERENCES

Bland, R., Waite-Cusic, J., Weisberg, A. J., Riutta, E. R., Chang, J. H., and Kovacevic, J. (2022). Adaptation to a Commercial Quaternary Ammonium Compound Sanitizer Leads to Cross-Resistance to Select Antibiotics in Listeria monocytogenes Isolated From Fresh Produce Environments. *Front. Microbiol.* 12:782920. doi: 10.3389/fmicb.2021.782920

Buffet-Bataillon, S., Tattevin, P., Maillard, J. Y., Bonnaure-Mallet, M., and Jolivet-Gougeon, A. (2016). Efflux pump induction by quaternary ammonium compounds and fluoroquinolone resistance in bacteria. *Fut. Microbiol.* 11, 81–92. doi: 10.2217/fmb.15.131

De Castro, E., Sigrist, C. J., Gattiker, A., Bulliard, V., Langendijk-Genevaux, P. S., Gasteiger, E., et al. (2006). ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. *Nucleic Acids Res.* 34, W362–W365. doi: 10.1093/nar/gkl124

Denes, T., Vongkamjan, K., Ackermann, H. W., Moreno Switt, A. I., Wiedmann, M., and Den Bakker, H. C. (2014). Comparative genomic and morphological analyses of Listeria phages isolated from farm environments. *Appl. Environ. Microbiol.* 80, 4616–4625. doi: 10.1128/AEM.00720-14

EFSA, and ECDC. (2021). The European Union One Health 2019 Zoonoses Report. *EFSA J.* 19:e06406. doi: 10.2903/j.efsa.2021.6406

Gerba, C. P. (2015). Quaternary ammonium biocides: efficacy in application. *Appl. Environ. Microbiol.* 81, 464–469. doi: 10.1128/AEM.02633-14

Guérin, A., Bridier, A., Le Grandois, P., Sévellec, Y., Palma, F., Félix, B., et al. (2021). Exposure to Quaternary Ammonium Compounds Selects Resistance to Ciprofloxacin in Listeria monocytogenes. *Pathogens* 10:220. doi: 10.3390/ pathogens10020220

Guérin, F., Galimand, M., Tuambilangana, F., Courvalin, P., and Cattoir, V. (2014). Overexpression of the novel MATE fluoroquinolone efflux pump FepA in Listeria monocytogenes is driven by inactivation of its local repressor FepR. *PLoS One* 9:e106340–e106340. doi: 10.1371/journal.pone.0106340

Guerreiro, D. N., Arcari, T., and O'byrne, C. P. (2020). The σ(B)-Mediated General Stress Response of Listeria monocytogenes: Life and Death Decision Making in a Pathogen. *Front. Microbiol.* 11:1505. doi: 10.3389/fmicb.2020.01505

Jones, P., Binns, D., Chang, H.-Y., Fraser, M., Li, W., and Mcanulla, C. (2014). InterProScan 5: genome-scale protein function classification. *Bioinformatics* 30, 1236–1240. doi: 10.1093/bioinformatics/btu031

Kampf, G. (2019). Antibiotic ResistanceCan Be Enhanced in Gram-Positive Species by Some Biocidal Agents Used for Disinfection. *Antibiotics* 8:13. doi: 10.3390/ antibiotics8010013

## ACKNOWLEDGMENTS

We would like to thank Renáta Karpíškova (Veterinary Research Institute, Brno, Czech Republic) for providing Czech strains from soil and ruminants, isolated in the Czech Republic and Francesco Pomilio (IZSAM, Istituto Zooprofilattico Sperimentale dell'Abruzzo e del Molise Teramo, Italy) for providing an Italian strain isolated from a wild animal. We also thank Taran Skjerdal (Norwegian Veterinary Institute, Food Safety and Animal Health Research Oslo, Norway) for providing a Norwegian strain, isolated from a wild animal and Bart Wullings (Food Safety Research, Wageningen University and Research, Netherlands) for providing a Dutch strain, isolated from a farm animal.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.864576/full#supplementary-material

Kode, D., Nannapaneni, R., and Chang, S. (2021). Low-Level Tolerance to Antibiotic Trimethoprim in QAC-Adapted Subpopulations of Listeria monocytogenes. *Foods* 10:1800. doi: 10.3390/foods10081800

Krivák, R., and Hoksza, D. (2018). P2Rank: machine learning based tool for rapid and accurate prediction of ligand binding sites from protein structure. *J. Cheminform.* 10:39. doi: 10.1186/s13321-018-0285-8

Lunden, J., Autio, T., Markkula, A., Hellstrom, S., and Korkeala, H. (2003). Adaptive and cross-adaptive responses of persistent and non-persistent Listeria monocytogenes strains to disinfectants. *Int. J. Food Microbiol.* 82, 265–272. doi: 10.1016/s0168-1605(02)00312-4

Moura, A., Criscuolo, A., Pouseele, H., Maury, M. M., Leclercq, A., and Tarr, C. (2016). Whole genome-based population biology and epidemiological surveillance of Listeria monocytogenes. *Nat. Microbiol.* 2:16185. doi: 10.1038/ nmicrobiol.2016.185

Olaimat, A. N., Al-Holy, M. A., Shahbaz, H. M., Al-Nabulsi, A. A., Abu Ghoush, M. H., Osaili, T. M., et al. (2018). Emergence of Antibiotic Resistance in Listeria monocytogenes Isolated from Food Products: A Comprehensive Review. *Compr. Rev. Food Sci. Food Saf.* 17, 1277–1292. doi: 10.1111/1541-4337.12387

Palma, F., Brauge, T., Radomski, N., Mallet, L., Felten, A., Mistou, M.-Y., et al. (2020). Dynamics of mobile genetic elements of Listeria monocytogenes persisting in ready-to-eat seafood processing plants in France. *BMC Genomics* 21:130. doi: 10.1186/s12864-020-6544-x

Rakic-Martinez, M., Drevets, D. A., Dutta, V., Katic, V., and Kathariou, S. (2011). Listeria monocytogenes strains selected on ciprofloxacin or the disinfectant benzalkonium chloride exhibit reduced susceptibility to ciprofloxacin, gentamicin, benzalkonium chloride, and other toxic compounds. *Appl. Environ. Microbiol.* 77, 8714–8721. doi: 10.1128/AEM.0598714

Samad, A., Li, Y., Zhang, C., Chen, F., Zeng, W., Fan, X., et al. (2019). X-ray crystal structure of putative transcription regulator Imo2088 from Listeria monocytogenes. *Biochem. Biophys. Res. Commun.* 520, 434–440. doi: 10.1016/j. bbrc.2019.10.033

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. *Bioinformatics* 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

Studer, G., Rempfer, C., Waterhouse, A. M., Gumienny, R., Haas, J., and Schwede, T. (2019). QMEANDisCo—distance constraints applied on model quality estimation. *Bioinformatics* 36, 1765–1771. doi: 10.1093/bioinformatics/btz828

Takada, H., and Yoshikawa, H. (2018). Essentiality and function of WalK/WalR two-component system: the past, present, and future of research. *Biosci. Biotechnol. Biochem.* 82, 741–751. doi: 10.1080/09168451.2018.1444466
Tasara, T., Klumpp, J., Bille, J., and Stephan, R. (2016). Genome Sequences of Listeria monocytogenes Strains Responsible for Cheese- and Cooked Ham Product-Associated Swiss Listeriosis Outbreaks in 2005 and 2011. Genome Announc. 4:e106-e116. doi: 10.1128/genomeA.00106-16

Tonkin-Hill, G., Macalasdair, N., Ruis, C., Weimann, A., Horesh, G., and Lees, J. A. (2020). Producing Polished Prokaryotic Pangenomes with the Panaroo Pipeline. Genome Biol. 21:180. doi: 10.1186/s13059-020-02090-4

Van Zundert, G. C. P., Rodrigues, J. P. G. L. M., Trellet, M., Schmitz, C., Kastritis, P. L., and Karaca, E. (2016). The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes. J. Mol. Biol. 428, 720-725. doi: 10.1016/j.jmb.2015.09.014

Vila Nova, M., Durimel, K., La, K., Felten, A., Bessières, P., Mistou, M.-Y., et al. (2019). Genetic and metabolic signatures of Salmonella enterica subsp. enterica associated with animal sources at the pangenomic scale. BMC Genomics 20:814. doi: 10.1186/s12864-019-6188-x

Waterhouse, A., Bertoni, M., Bienert, S., Studer, G., Tauriello, G., and Gumienny, R. (2018). SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 46, W296-W303. doi: 10.1093/nar/gky427

Wiktorczyk-Kapischke, N., Skowron, K., Grudlewska-Buda, K., Wałecka-Zacharska, E., Korkus, J., and Gospodarek-Komkowska, E. (2021). Adaptive Response of Listeria monocytogenes to the Stress Factors in the Food Processing Environment. Front. Microbiol. 12:710085. doi: 10.3389/fmicb.2021.710085

Wilson, A., Gray, J., Chandry, P. S., and Fox, E. M. (2018). Phenotypic and Genotypic Analysis of Antimicrobial Resistance among Listeria monocytogenes Isolated from Australian Food Production Chains. Genes 9:80. doi: 10.3390/genes9020080

Wu, Q., Peng, Z., Zhang, Y., and Yang, J. (2018). COACH-D: improved protein-ligand binding sites prediction with refined ligand-binding poses through molecular docking. Nucleic Acids Res. 46, W438-W442. doi: 10.1093/nar/gky439

Yu, T., Jiang, X., Zhang, Y., Ji, S., Gao, W., and Shi, L. (2018). Effect of Benzalkonium Chloride Adaptation on Sensitivity to Antimicrobial Agents and Tolerance to Environmental Stresses in Listeria monocytogenes. Front. Microbiol. 9:2906. doi: 10.3389/fmicb.2018.02906

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

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Douarre, Sévellec, Le Grandois, Soumet, Bridier and Roussel. 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.