Genomic characterization of Listeria monocytogenes isolates reveals that their persistence in a pig slaughterhouse is linked to the presence of benzalkonium chloride resistance genes
BMC Microbiologyvolume 18, Article number: 220 (2018)
The aim of this study was to characterize the genomes of 30 Listeria monocytogenes isolates collected at a pig slaughterhouse to determine the molecular basis for their persistence.
Comparison of the 30 L. monocytogenes genomes showed that successive isolates (i.e., persistent types) recovered from thew sampling site could be linked on the basis of single nucleotide variants confined to prophage regions. In addition, our study revealed the presence among these strains of the bcrABC cassette which is known to produce efflux pump-mediated benzalkonium chloride resistance, and which may account for the persistence of these isolates in the slaughterhouse environment. The presence of the bcrABC cassette was confirmed by WGS and PCR and the resistance phenotype was determined by measuring minimum inhibitory concentrations. Furthermore, the BC-resistant strains were found to produce lower amounts of biofilm in the presence of sublethal concentrations of BC.
High resolution SNP-based typing and determination of the bcrABC cassette may provide a means of distinguishing between resident and sporadic L. monocytogenes isolates, and this in turn will support better management of this pathogen in the food industry.
Listeria monocytogenes is a bacterial pathogen that causes listeriosis, the foodborne illness with the highest fatality rate [1, 2]. The ability of this organism to colonize and persist in food processing plants increases the risk of food contamination . Persistent strains have been involved in some listeriosis outbreaks [4,5,6], making persistent contamination of L. monocytogenes an important consideration for public health. The importance of the persistence of L. monocytogenes in food processing plants has been previously discussed [7,8,9]. Whereas some L. monocytogenes strains can persist within the food processing environment for years , some strains are merely transitory in these environments, and the approaches for their mitigation will be very different (i.e., different sanitation regimens). Many genetically distinct L. monocytogenes strains can be isolated sporadically during different sampling periods while other persistent strains can be isolated over the same periods [10,11,12]. Although several explanations have been proposed for the persistence of specific L. monocytogenes strains , the reasons are not fully understood. It has been postulated that L. monocytogenes strains persist due to opportunistic colonization of harborage sites in food processing facility premises and on equipment . Alternatively, it has been suggested that the persistence of L. monocytogenes is linked to its ability to produce biofilms [14, 15]. Other studies using phenotypic methods have observed resistance to quaternary ammonium compounds (QAC) in persistent strains [16,17,18,19]. Recently, genetic factors associated with the QAC-resistance phenotype and resistance to benzalkonium chloride (BC), a QAC largely used in the food industry, were identified [17, 20,21,22,23]. These include the benzalkonium chloride resistance cassette bcrABC , in which bcrAB code for the small multidrug resistance protein family (SMR) transporters and bcrC codes for a transcriptional factor . In addition, the qacH gene which codes for a SMR transporter associated with the export of benzalkonium chloride and acquired by the Tn6188 transposon [19, 20, 24], and other QAC determinants originally observed in Staphylococci .
Although there is no consensus on the definition of a persistent strain, a persistent status should be considered when the same molecular type is isolated many times at different visits within the same processing plant for a defined period of time . Therefore, identification of recurring highly genetically related isolates is the first step necessary in studying the association between resistance to BC and persistence. The persistence designation used in previous studies was based on pulsed-field gel electrophoresis (PFGE) typing [19, 24], an international standardized method for L. monocytogenes characterization. However, recent studies using whole-genome sequencing (WGS) enabling the core genome/whole genome MLST (cgMLST/wgMLST) analyses and single nucleotide variant (SNV) analyses have demonstrated that a better discriminatory power can be achieved [26, 27]. WGS-based typing methods may allow for better clustering of highly-related persistent isolates.
In this study, we set to characterize three groups of L. monocytogenes strains – two persistent and one sporadic – isolated from a pig slaughterhouse, with a particular emphasis on the BC resistance genotype/phenotype of these strains. We also performed complementary analyses to identify high quality single nucleotide variants (hqSNVs) and predict affected functions on the basis of these SNVs.
Materials and methods
Sampling, bacterial isolation, and pulsed field gel electrophoresis characterization
Sampling procedures, bacterial isolation and pulsed-field gel electrophoresis (PFGE) typing were performed according to Lariviere-Gauthier et al., . Briefly, bacterial strains were isolated from repeated environmental sampling at a pork slaughterhouse and cutting facility. The samples used in this study were collected during three separate visits over a one-month period, with an interval of 2 weeks between visits 1 and 2 and 1 week between visits 2 and 3. Samples were collected using pre-moistened swabs, after cleaning and sanitation procedures, on environmental surfaces in contact or not with the products, a list of surface types is detailed in Table 1. The disinfectant concentration used in this slaughterhouse, and specifically in cutting facility, was from 150 to 200 ppm.
Another strain, Lis8316 belonging to the group B was isolated 3 years after in the same slaughterhouse.
Detection of Listeria monocytogenes was performed using Health Canada’s MFHPB-30 based-method as described previously , briefly, two enrichment broths were used. The first enrichment was conducted in University of Vermont medium 1 broth (UVM-1; Lab M, Heywood, United Kingdom) for 24 h at 30 C°, and the second in Fraser broth (Lab M, United Kingdom) for 48 h at 37 C°. Selective Listeria Ottavani and Agosti (ALOA; AES Chemunex, Bruz, France) agar was used to isolate L. monocytogenes strains from the two broths and the identification of isolates were based on the use of carbohydrates, Christie-Atkins-Munch-Petersen (CAMP test), hemolysis and motility.
The pulsed field gel electrophoresis characterization of isolates was based on the Center disease and control center (CDC) PulseNet protocol . Pattern analysis was performed using Bionumerics software (version 6.5, Applied Maths, Kortrijk, Belgium) and strains were clustered with the unweighted pair group method with arithmetic means.
For whole genome sequencing, strains were selected according to their recurrence during the sampling period (transient vs persistent). PFGE type 1 and 16 belonging to the sequence type (ST) 5 and 9 respectively, were considered persistent since they were isolated at each visit and after cleaning and sanitation procedures (respectively named persistent B and A). PFGE types 6, 10, 9, 18, 8 respectively with a ST of 6, 5, 5, 9 and 1 were considered to be transient or sporadic strains since they were isolated only at one visit. The persistent group A and B and sporadic strains are phylogenetically unrelated, even if the sporadic strain p961 (pulsotype 18) showed a polymorphism of three different bands with the group A considering Apa1 and AscI restriction profiles. A PFGE profiles of a representative strains from the group A and B and p961 strain (pulsotype 18) were given in a Additional file 1: Figure S1. The numbers and characteristic selection of isolates are detailed in Table 1.
DNA extraction and sequencing
One mL of Brain Heart Infusion (BHI) broth (Becton and Dickinson, Franklin Lakes, NJ, USA) was inoculated with a single colony of L. monocytogenes and incubated for 4 to 6 h at 37 °C. Bacterial cells were collected from 400 μL of culture by centrifuging at 15,000 x g for 2 min. Bacterial pellets were resuspended in 200 μL of lysis buffer containing 1% Triton X-100 (Sigma-Aldrich, ON, Canada) and 10 mg/mL lysozyme (Fisher scientific, ON, Canada) and incubated for 20 min at 37 °C. Genomic DNA (gDNA) was extracted using the Maxwell® 16 Cell LEV DNA Purification kit (Promega, Madison, WI) as recommended by the manufacturer. DNA was then quantified using the Quant-iT™ High-Sensitivity DNA Assay Kit (Life Technologies Inc., Burlington, ON). Sequencing libraries were constructed from 1 ng of gDNA using Nextera XT DNA Sample Preparation and Index Kits (Illumina, Inc., San Diego, CA) according to manufacturers’ instructions. Genomic sequencing was performed on the Illumina MiSeq Platform (Illumina, Inc.) using a 600-cycle MiSeq Reagent kit v3 (Illumina, Inc.).
Genome and plasmid assembly and annotation
Reads were quality checked and de novo assembled using SPAdes 3.12.0  with default settings. The resulting draft genomes were analyzed for assembly quality using Quast 4.1 . For plasmid assembly, plasmidSPAdes  pipeline was used to assemble only the plasmid from the WGS data using the default settings. To identify the presence of large plasmid pLM80 harboring the bcrABC cassette, the resulting plasmid assemblies were blasted online using the basic local alignment search tool (Blastn) from the National Center for Biotechnology Information’s (NCBI) website (https://blast.ncbi.nlm.nih.gov/Blast.cgi) with default settings using the nucleotide collection (nr/nt) database.
Core genome MLST (cgMLST) characterization
cgMLST was conducted using a genome based characterization method developed by Moura et al. . The assembled 30 genomes were submitted to the Institut Pasteur’s website (http://bigsdb.pasteur.fr/listeria/listeria.html) for cgMLST typing using 1748 loci and a cgMLST type (CT) was defined as a group of cgMLST profiles that differ by up to seven allelic mismatches out of 1748 loci.
High quality core genome single nucleotide variant (hqSNV) analysis
The paired end reads of closely related isolates were quality trimmed and analyzed for hqSNV detection by using the SNVPhyl core genome hqSNV pipeline (https://github.com/apetkau/core-phylogenomics) , available as a tool shed in Galaxy (https://toolshed.g2.bx.psu.edu/). Briefly, the reference genomes were analyzed with PHASTER [32, 33] for identifying phage regions. The paired end reads were mapped against the reference genomes using SMALT v0.7.5 (http://www.sanger.ac.uk/science/tools/smalt-0). The FreeBayes 0.9.20 was used to call the high-quality variants with a minimum depth coverage of 15 and a minimum mapping quality of 30. The minimum SNV abundance ratio was adjusted at 0.75. The uncalled positions were analyzed for sufficient coverage with the BCFtools that came in the SAMtools package. Reference genomes EGD-e (NCBI reference number NC003210.1) and 1/2b str 10–0810 (NCBI reference number NZ_CP007168.1) were used for groups A and B respectively, according to their phylogenetically relatedness throughout all analyses.
The resulting variants (hqSNVs) were annotated and filtered to the reference genome Listeria monocytogenes EGD-e (NC003210.1) with snpEff and snpSift program respectively  for the prediction of affected functions based on the identification of genes in which the mutations were found. In this study, only the high impact hqSNVs were taken into consideration (i.e. SNVs leading to a stop gained, start loss, or a frameshift).
Screening of benzalkonium chloride resistance genes in L. monocytogenes strains
Identification of genetic determinants for resistance was conducted at the Institut Pasteur, according to their established protocol . Briefly, a set of resistance genes was identified based on previous studies [35,36,37] and PATRIC, a public resistance gene database . A BLASTn algorithm was used to identify the resistant genes with a minimum nucleotide identity of 70% and alignment length coverage of 70%. The results were represented in a presence-absence heatmap using the phandango visualization tool .
Minimal inhibitory concentration of BC and PCR screening for bcrABC
To determine if the persistent groups were resistant to BC, Minimal inhibitory concentrations (MICs) were determined according to the clinical and laboratory standards institute (CLSI) guidelines with minor modifications. Briefly, tryptic soy broth with yeast extract (TSBYE) supplied with BC (Sigma-Aldrich, ON, Canada) at a final concentration covering a range from 0 to 200 ppm was inoculated with 10− 5 CFU from an overnight culture. The plates were incubated at 30 °C for 48 h. The MIC was defined as the lowest concentration at which a growth inhibition was observed. The assay was conducted in triplicate.
PCR was used to confirm the presence of bcrABC genes. DNA was extracted from sporadic and persistent strains using Chelex (Biorad, Mississauga, ON, CA) as follows: isolates were plated on blood sheep agar and colonies were picked and resuspended in 50 μL of Chelex at 0.6% (w/v). The suspensions were heated first at 55 °C for 30 min and then at 98 °C for 15 min and were vortexed between steps. The solutions were centrifuged at 12000 rpm and the upper phase containing DNA was transferred to a clean tube and stored at 4 °C for quick use. The PCR parameters consisted of an initial denaturation of 98 °C for 30 s followed by 30 cycles of denaturation at 98 °C for 30 s, annealing at 60 °C for 40 s, elongation at 72 °C for 1.30 min, and final elongation at 72 °C for 5 min. The primers used here were designed in a previous study [21, 40].
Biofilm formation ability
Isolates from all groups were tested for biofilm formation ability using the microtiter plate assay in the presence and absence of BC (Sigma-Aldrich, ON, Canada), briefly, an overnight culture in tryptic soy broth supplemented with 0.6% (w/v) of yeast extract was inoculated to BHI (Becton and Dickinsen, Franklin Lakes, NJ, USA) supplemented or not with 3.125 and 0,78 ppm of BC (at final concentration) at 1/1000. 100 μL of the inoculated BHI was distributed into each well of 96 well microplates (Costar® 3370; Corning, NY, USA) and incubated at 30 °C for 48 h. The supernatant was discarded from the mature biofilm and this latter washed to eliminate planktonic cells before stained with Cristal Violet (CV) 0.1% (w/v) for 20 min and washed two times. Ethanol 95% was used to solubilize CV within biofilm and Optical density (Od) was taken at 595 nm. Experiment was repeated in triplicate at different days.
To determine if there were differences in biofilm formation capacity between strains from the three groups in the presence and absence of BC, linear mixed model was performed on the data of groups (A, B, C) and treatments (control, 3.125 and 0.78 ppm) as fixed factors and isolate as random factor. A series of a priori contrasts were conducted to compare the mean in the control groups with the mean in the other two treatments in each group and to compare the mean in each group for each treatment. Alpha level was adjusted with Benjamini-Hochberg procedure. Results were significantly different if p < 0.05.
A total of 30 L. monocytogenes strains was isolated: 25 were persistent, present during all visits, and 5 were sporadic, present during a single visit. The 25 persistent isolates were sub-divided into 2 Groups, A, 6 isolates, and B, 19 isolates, based on their pulsotypes. The 5 sporadic isolates were assigned to Group C (Table 1).
The cgMLST analysis showed eight cgMLST types (CT) corresponding to the six Pulsotypes. It is noteworthy that all of the CTs observed in our study were new in the BIGSdb database of the Institut Pasteur (at time of analysis). All of the strains from the group B (pulsotype 1) belonged to the new CT2802 whereas the strains from group A (pulsotype 16) showed two different CTs, with one strain belonging to CT4145 and the rest of the group to CT2825 (Fig. 1). Interestingly, the sporadic strain p961 (pulsotype 18) from sporadic strains group was closely related to persistent group A, and sporadic strains p959 and p976 were closely related to group B (Fig. 1).
The SNVphyl analysis showed some differences in hqSNV’s numbers within the persistent strains: with 0 to 25 hqSNVs in group A and 0 to 13 hqSNVs in group B. Thus, in group A, only the p984 strain exhibited a high number of hqSNVSs (25) and was more distant from the remaining strains (Fig. 2a). Almost all the strains from group B harbored 4 hqSNV differences or fewer, except P978, p965 and Lis8316 strains which had more hqSNV differences compared to the rest of group B. This was reflected by the close positions in the phylogenetic tree (Fig. 2b).
SNVs confined mostly to prophage regions
To study different functions affected by mutations, hqSNVs were analyzed with the snpEff program. The inferred function of impacted proteins by the hqSNVs were mostly related to bacteriophage proteins, particularly in group A (Fig. 3). Most of the hqSNVs were found within phage A118 (Table 2A). Furthermore, the position of these phages within the genome was the same between the unrelated persistent groups A and B (Table 2A and B). There were 8 hqSNVs that were detected in genes coding for known functions of the L. monocytogenes core genome in group B, but five out of these eight were present in only one to four strains (Table 2B). Otherwise, in group B, a SNV in the comk gene was identified at a position that resulted in a disruption of this gene due to loss of a stop codon (Table 2B). A stop codon in the inlA gene was detected as well as the SNVs, and this was confirmed by Sanger sequencing of the internalin A gene (inlA) in group B strains in a previous study . These SNVs lead to a premature stop codon (PMSC). Interestingly, another SNV was detected in an internalin B gene (inlB) confirmed by BLASTn. The rest of the SNVs were detected in genes coding for hypothetical proteins. Among them, many shared similarities with the A118 phage proteins (Additional file 2: Table S1).
Presence of bcrABC in a persistent strains group B
In order to study the presence or absence of genes associated with the persistence phenotype, identification of accessory genes was conducted on the persistent groups as well as on the five different sporadic strains isolated from the same slaughterhouse. Our results showed the presence of bcrABC, a resistance gene cassette to BC, a commonly used agri-food sanitizers [17, 20,21,22,23]. More importantly, this genetic determinant was present in only persistent group B and in its closely related sporadic strain p959, while the stress response island genes (ssi1) were present in all strains except p957 and p983 (Fig. 4). The bcrABC cassette was identified in a plasmid with 80 kbp in all the persistent strain group B with 100% similarity with plasmid pLM80-cont 507, previously described  with accession number of AADR01000010 (Fig. 5). In addition to bcrABC, this plasmid was identified as a vector for cadmium and arsenic resistance as well . To confirm the presence of the bcrABC genes in only the group B, a PCR was conducted in all studied isolates. The PCR confirmed the presence of the bcrABC cassette in the group B and its absence in the group A and in sporadic strains (Table 3). Furthermore, using the CMI test, we observed growth in all isolates from the group B at 6.25 ppm of BC, indicating decrease of sensitivity to this compound, while the rest of isolates from the group A and C did not grow at this concentration (Table 3).
BC effect on biofilm ability is strain and concentration dependent
In order to investigate the ability of biofilm formation by strains from the persistent group A, B and the sporadic strains in the presence and absence of BC, microtiter plate assay for biofilm growth was performed on all isolates. Since only strains from group B grew at 3.125 ppm and strains from groups A and C grew in the presence of 0.78 ppm of BC, comparisons were done at the latter concentration. Results showed no significant difference between isolates from the three groups in absence or presence of BC (Fig. 6a and b, respectively), in fact all the strains were able to form biofilm and seem to have an equal ability for biofilm formation. To study if BC affected the biofilm formation of strains in each group, biofilms formed in the presence of BC at 3.125 ppm (for group B strains) and 0,78 ppm were compared to the biofilms formed in the absence of BC. Results showed that biofilms formed in presence of BC decreased for all groups at all concentrations of BC, but the differences were not significant for group A and C strains grown in the presence of a sublethal concentration of 0.78 ppm compared to the biofilms grown in the absence of BC (Fig. 6c). In contrast, for group B strains, biofilm formation in the presence of BC at a sublethal concentration of 3.125 ppm was significantly lower than growth in the absence of BC (Fig. 6c).
The aim of this work was to identify genetic determinants that could be associated with persistence of L. monocytogenes strains within the slaughterhouse environment. We were especially interested in the potential role of genetic determinants related to resistance to quaternary ammonium, an antimicrobial compound commonly used as a disinfectant in slaughterhouses.
The persistence phenotype of some L. monocytogenes strains has been previously reported [3, 9, 13], however, the reason why these strains can persist in the food processing environment is yet to be fully understood . Here, we describe genomic characteristics of two unrelated persistent groups of strains isolated from slaughterhouse environments at three different times. Although our sampling interval was relatively short, we believe that the isolation of the same CTs at the same location after daily cleaning and sanitation procedures, is a sufficient parameter to attribute the persistence character to our isolates. Furthermore, a sampling completed 3 years after our sampling for this study at the same location allowed us to isolate the same CT from the group B (Lis8316).
The cgMLST and the SNVphyl analysis confirmed the close genetic relatedness between isolates of pulsotype 1 from group B belonging to ST5 and clonal complex (CC) CC5, a clonal group involved in a variety of food products and processing environments contamination as well as in at least two outbreak listeriosis in the United states . Two different CTs, CT2825 and CT4146 were distinguished in group A strains (pulsotype 16) leading to improved discriminatory power, as reported in previous studies [26, 43]. These CTs belonged to ST9 which was described previously as highly prevalent in the food associated L. monocytogenes strains [44,45,46]. The SNVphyl analysis showed the presence of 0 to 13 hqSNV differences in the persistent groups, except for strain p984, in group A, where the number of different hqSNVs reached 25 relatives to the rest of the strains. This result was in accordance with cgMLST typing as this strain had a different CT from the rest of group A. The number of hqSNVs did not vary significantly between each period and between strains within the two persistent groups, indicating that the strains were closely related. Indeed, there were no differences between them based on the number of SNVs. Such results are in accordance with other studies that observed limited accumulation of SNVs between two closely related strains isolated from a food processing plant at a 12 year interval . In the current study, a great proportion of SNVs were confined to prophage regions in the genomes and there were very few biologically significant mutations affecting the core genome. Thus, while previous studies showed that prophage insertion into the L. monocytogenes genome was the most important element contributing to the diversification of the strains , here we demonstrated in silico that the genes most affected by non-synonymous SNVs were also related to the prophage genes. These findings support the largely recognized view of the high stability of the L. monocytogenes genome [26, 37, 47, 48]. The presence of the comk prophage in the group B strains is in accordance with the persistence character of these strains and the association of the ST5 strains to food processing environment as described in previous study . This study allowed the identification of bcrABC in persistent strains group B. These results complement the Ortiz et al. study  where they found the tn6188 transposon, which is known to serve as a vector for the qacC resistance gene transfer in persistent strains. Efflux pump bcrABC was described in L. monocytogenes for the first time by Elhanafi et al. in a large Listeria monocytogenes plasmid, pLM80 . This plasmid was also previously identified as the vector for cadmium and arsenic resistance determinants in L. monocytogenes .
Since then, many studies found the bcrABC in strains isolated from retail foods [23, 50, 51]. The presence of bcrABC in strain Lis8216, isolated within a 3-years interval, may indicate the stability of the plasmid harboring this gene in environments where BC is commonly used. Moreover, persistent group B was isolated from conveyor surfaces, which have the greatest exposure to BC. However, in our analysis, if these strains were confirmed as less sensitive to BC treatments, they should still be considered as sensitive to BC concentrations used in the slaughterhouses where the sampling was done. This apparent paradox was previously reported [17, 25, 52] and was attributed to a dilution effect during sanitation procedures on wet surfaces, leading to sub-inhibitory concentrations of the disinfectants as reported previously . Persistent group A did not harbor the bcrABC cassette, and strains from this group were more sensitive to BC as compared to the ones in group B. Therefore, group A isolates were susceptible to BC. It is noteworthy that these strains were collected from equipment (e.g., saws), places known to be hard to clean, one may suppose that strains residing in such areas could not be exposed frequently to disinfectant. It is therefore easy to hypothesize that strains from group A did not submit pressure selection for bcrABC acquisition. Results for biofilm formation ability showed no significant differences between strains from persistent and non-persistent groups and all isolates were able to form a 48 h mature biofilm. Similar results were shown in previous studies where persistent and non-persistent L. monocytogenes strains isolated from different sources showed an equal biofilm formation capacity . Same results were obtained from environmental and food product strains in other recent studies where authors suggest that biofilm formation could not be a determinant for persistence  or the persistence of studied L. monocytogenes strains was not linked to biofilm formation ability and other determinants should be investigated . The presence of BC at 0.78 ppm during biofilm formation seems not to affect biofilm amount between the three groups, and when strains were compared for their ability to form biofilm in the presence and absence of BC, only isolates from group B seemed to be affected by the sublethal concentration of BC at 3.125 ppm, which may indicate the strain-specific variable responsiveness to the effect of BC on biofilm formation, and it appears that BC-resistant strains produced significantly less biofilm than susceptible strains as reported previously . It is not understood why BC-resistant L. monocytogenes strains had less ability to form biofilm in the presence of BC, though previous investigators suggested that modification of cellular morphology, phospholipids and fatty acids could be involved in biofilm formation .
Further studies using a greater number of unrelated persistent and non persistent strains should be conducted in different conditions, particularly those which mimic the slaughterhouse environment.
Our findings support the hypothesis that persistence, in the case of the strains from the group A, was related to biofilm formation and/or ecological niches formation in hard-to-clean places. Indeed, some authors suggested that persistence could be related to colonization of ecological niches less accessible during sanitation procedures as reported in previous reviews by Carpentier and Cerf (2011) and Ferreira et al. (2014) where authors concluded that there was no association between persistence and any phenotypic characteristics like biofilm formation and resistance to QAC [7, 13]. While we cannot exclude that some environmental sites were less efficiently cleaned, our findings support that there is not always only one determinant explaining the persistence of L. monocytogenes strains and that this phenomenon could be due to multiple factors facilitating the survival of this bacterium. This could explain why some studies focusing on only one determinant, like resistance to QAC, could not find a causal effect with persistence [9, 19]. We propose that the localized environmental conditions in which strains were isolated in the food processing plant could give clues to identify which determinant(s) enhances the strain’s survival. Thus, isolates collected from a conveyor known to be the most exposed surfaces to QAC disinfectants may benefit from QAC resistance while isolates recovered from difficult to clean surfaces may benefit from better biofilm formation phenotypes. We believe that persistent strains possess at least one, but most likely many persistence-associated determinants leading to survival in their specific environment conditions.
In this study, we described the distribution of non-synonymous SNVs in persistent L. monocytogenes genomes from a food processing environment and found evidence that the presence of genetic determinants associated with resistance to BC could be linked to the persistence of L. monocytogenes strains. The ability to differentiate between sporadic and persistent environmental contamination events using the methods described will be a key element in guiding risk management actions, such as determining the type of sanitary interventions (e.g., aggressive versus standard) required to mitigate identified hazards and adapt cleaning and sanitation procedures to avoid recurrent food contamination.
Brain heart infusion
Center disease and control center
Clinical and laboratory standards institute
High quality single nucleotide variant
Minimal inhibitory concentrations
National Center for Biotechnology Information
Pulsed field gel electrophoresis
Premature stop codon
Quaternary ammonium compound
Small multidrug resistance
Single nucleotide variant
Triptic soy broth yeast extract
University of Vermont Medium
Whole genome sequencing
Charlier C, Perrodeau E, Leclercq A, Cazenave B, Pilmis B, Henry B, Lopes A, Maury MM, Moura A, Goffinet F, et al. Clinical features and prognostic factors of listeriosis: the MONALISA national prospective cohort study. Lancet Infect Dis. 2017;17(5):510–9.
de Noordhout CM, Devleesschauwer B, Angulo FJ, Verbeke G, Haagsma J, Kirk M, Havelaar A, Speybroeck N. The global burden of listeriosis: a systematic review and meta-analysis. Lancet Infect Dis. 2014;14(11):1073–82.
Orsi RH, Borowsky ML, Lauer P, Young SK, Nusbaum C, Galagan JE, Birren BW, Ivy RA, Sun Q, Graves LM. Short-term genome evolution of Listeria monocytogenes in a non-controlled environment. BMC Genomics. 2008;9(1):539.
Ericsson H, Eklow A, Danielsson-Tham ML, Loncarevic S, Mentzing LO, Persson I, Unnerstad H, Tham W. An outbreak of listeriosis suspected to have been caused by rainbow trout. J Clin Microbiol. 1997;35(11):2904–7.
Nakari UM, Rantala L, Pihlajasaari A, Toikkanen S, Johansson T, Hellsten C, Raulo SM, Kuusi M, Siitonen A, Rimhanen-Finne R. Investigation of increased listeriosis revealed two fishery production plants with persistent Listeria contamination in Finland in 2010. Epidemiol Infect. 2014;142(11):2261–9.
Independent Investigator into the Listeriosis O, Weatherill S. Report of the independent investigator into the 2008 listeriosis outbreak. Canada; 2009.
Ferreira V, Wiedmann M, Teixeira P, Stasiewicz MJ. Listeria monocytogenes persistence in food-associated environments: epidemiology, strain characteristics, and implications for public health. J Food Prot. 2014;77(1):150–70.
Orsi RH, Borowsky ML, Lauer P, Young SK, Nusbaum C, Galagan JE, Birren BW, Ivy RA, Sun Q, Graves LM, et al. Short-term genome evolution of Listeria monocytogenes in a non-controlled environment. BMC Genomics. 2008;9:539.
Stasiewicz MJ, Oliver HF, Wiedmann M, den Bakker HC. Whole-genome sequencing allows for improved identification of persistent Listeria monocytogenes in food-associated environments. Appl Environ Microbiol. 2015;81(17):6024–37.
Simmons C, Stasiewicz MJ, Wright E, Warchocki S, Roof S, Kause JR, Bauer N, Ibrahim S, Wiedmann M, Oliver HF. Listeria monocytogenes and Listeria spp. contamination patterns in retail delicatessen establishments in three U.S. states. J Food Prot. 2014;77(11):1929–39.
Ortiz S, Lopez V, Villatoro D, Lopez P, Davila JC, Martinez-Suarez JN. A 3-year surveillance of the genetic diversity and persistence of Listeria monocytogenes in an Iberian pig slaughterhouse and processing plant. Foodborne Pathog Dis. 2010;7(10):1177–84.
Lariviere-Gauthier G, Letellier A, Kerouanton A, Bekal S, Quessy S, Fournaise S, Fravalo P. Analysis of Listeria monocytogenes strain distribution in a pork slaughter and cutting plant in the province of Quebec. J Food Prot. 2014;77(12):2121–8.
Carpentier B, Cerf O. Review--persistence of Listeria monocytogenes in food industry equipment and premises. Int J Food Microbiol. 2011;145(1):1–8.
Nowak J, Cruz CD, Tempelaars M, Abee T, van Vliet AHM, Fletcher GC, Hedderley D, Palmer J, Flint S. Persistent Listeria monocytogenes strains isolated from mussel production facilities form more biofilm but are not linked to specific genetic markers. Int J Food Microbiol. 2017;256:45–53.
Wang J, Ray AJ, Hammons SR, Oliver HF. Persistent and transient Listeria monocytogenes strains from retail deli environments vary in their ability to adhere and form biofilms and rarely have inlA premature stop codons. Foodborne Pathog Dis. 2015;12(2):151–8.
Lundén J. Persistent Listeria monocytogenes contamination in food processing plants; 2004.
Fox EM, Leonard N, Jordan K. Physiological and transcriptional characterization of persistent and nonpersistent Listeria monocytogenes isolates. Appl Environ Microbiol. 2011;77(18):6559–69.
Aase B, Sundheim G, Langsrud S, Rorvik LM. Occurrence of and a possible mechanism for resistance to a quaternary ammonium compound in Listeria monocytogenes. Int J Food Microbiol. 2000;62(1–2):57–63.
Ortiz S, López-Alonso V, Rodríguez P, Martínez-Suárez JV. The connection between persistent, disinfectant-resistant Listeria monocytogenes strains from two geographically separate Iberian pork processing plants: evidence from comparative genome analysis. Appl Environ Microbiol. 2016;82(1):308–17.
Müller A, Rychli K, Muhterem-Uyar M, Zaiser A, Stessl B, Guinane CM, Cotter PD, Wagner M, Schmitz-Esser S. Tn6188-a novel transposon in Listeria monocytogenes responsible for tolerance to benzalkonium chloride. PLoS One. 2013;8(10):e76835.
Elhanafi D, Dutta V, Kathariou S. Genetic characterization of plasmid-associated benzalkonium chloride resistance determinants in a Listeria monocytogenes strain from the 1998-1999 outbreak. Appl Environ Microbiol. 2010;76(24):8231–8.
Dutta V, Elhanafi D, Kathariou S. Conservation and distribution of the benzalkonium chloride resistance cassette bcrABC in Listeria monocytogenes. Appl Environ Microbiol. 2013;79(19):6067–74.
Rosa Ebner R, Stephan R, Althaus D, Brisse S, Maury M, Tasara T. Phenotypic and genotypic characteristics of Listeria monocytogenes strains isolated during 2011–2014 from different food matrices in Switzerland, vol. 57; 2015.
Ortiz S, López V, Martínez-Suárez JV. The influence of subminimal inhibitory concentrations of benzalkonium chloride on biofilm formation by Listeria monocytogenes. Int J Food Microbiol. 2014;189:106–12.
Xu D, Li Y, Shamim Hasan Zahid M, Yamasaki S, Shi L, Li J-r, Yan H. Benzalkonium chloride and heavy-metal tolerance in Listeria monocytogenes from retail foods. Int J Food Microbiol. 2014;190:24–30.
Moura A, Criscuolo A, Pouseele H, Maury MM, Leclercq A, Tarr C, Björkman JT, Dallman T, Reimer A, Enouf V. Whole genome-based population biology and epidemiological surveillance of Listeria monocytogenes. Nat Microbiol. 2016;2:16185.
Hingston P, Chen J, Dhillon BK, Laing C, Bertelli C, Gannon V, Tasara T, Allen K, Brinkman FSL, Truelstrup Hansen L, et al. Genotypes associated with Listeria monocytogenes isolates displaying impaired or enhanced tolerances to cold, salt, acid, or desiccation stress. Front Microbiol. 2017;8. https://doi.org/10.3389/fmicb.2017.00369.
Graves LM, Swaminathan B. PulseNet standardized protocol for subtyping Listeria monocytogenes by macrorestriction and pulsed-field gel electrophoresis. Int J Food Microbiol. 2001;65.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–5.
Petkau A. SNVPhyl: a single nucleotide variant phylogenomics pipeline for microbial genomic epidemiology. Microb Genom. 2017;3(6). https://doi.org/10.1099/mgen.0.000116.
Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: A Fast Phage Search Tool. Nucleic Acids Res. 2011;39(suppl_2):W347–52.
Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, Wishart DS. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44(W1):W16–21.
Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012;6(2):80–92.
Maury MM, Tsai Y-H, Charlier C, Touchon M, Chenal-Francisque V, Leclercq A, Criscuolo A, Gaultier C, Roussel S, Brisabois A, et al. Uncovering Listeria monocytogenes hypervirulence by harnessing its biodiversity. Nat Genet. 2016;48:308.
Gahan CGM, Hill C. Listeria monocytogenes: survival and adaptation in the gastrointestinal tract. Front Cell Infect Microbiol. 2014;4.
Kuenne C, Billion A, Mraheil MA, Strittmatter A, Daniel R, Goesmann A, Barbuddhe S, Hain T, Chakraborty T. Reassessment of the Listeria monocytogenespan-genome reveals dynamic integration hotspots and mobile genetic elements as major components of the accessory genome. BMC Genomics. 2013;14(1):47.
Wattam AR, Davis JJ, Assaf R, Boisvert S, Brettin T, Bun C, Conrad N, Dietrich EM, Disz T, Gabbard JL, et al. Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. Nucleic Acids Res. 2017;45(D1):D535–42.
Hadfield J, Croucher NJ, Goater RJ, Abudahab K, Aanensen DM, Harris SR. Phandango: an interactive viewer for bacterial population genomics. Bioinformatics. 2017:btx610.
Sidhu MS, Heir E, Leegaard T, Wiger K, Holck A. Frequency of disinfectant resistance genes and genetic linkage with beta-lactamase transposon Tn552 among clinical staphylococci. Antimicrob Agents Chemother. 2002;46(9):2797–803.
Nelson KE, Fouts DE, Mongodin EF, Ravel J, DeBoy RT, Kolonay JF, Rasko DA, Angiuoli SV, Gill SR, Paulsen IT, et al. Whole genome comparisons of serotype 4b and 1/2a strains of the food-borne pathogen Listeria monocytogenes reveal new insights into the core genome components of this species. Nucleic Acids Res. 2004;32(8):2386–95.
Chen Y, Luo Y, Curry P, Timme R, Melka D, Doyle M, Parish M, Hammack TS, Allard MW, Brown EW. Assessing the genome level diversity of Listeria monocytogenes from contaminated ice cream and environmental samples linked to a listeriosis outbreak in the United States. PLoS One. 2017;12(2):e0171389.
Maury MM, Tsai YH, Charlier C, Touchon M, Chenal-Francisque V, Leclercq A, Criscuolo A, Gaultier C, Roussel S, Brisabois A, et al. Uncovering Listeria monocytogenes hypervirulence by harnessing its biodiversity. Nat Genet. 2016;48(3):308–13.
Ebner R, Stephan R, Althaus D, Brisse S, Maury M, Tasara T. Phenotypic and genotypic characteristics of Listeria monocytogenes strains isolated during 2011–2014 from different food matrices in Switzerland. Food Control. 2015;57:321–6.
Chenal-Francisque V, Lopez J, Cantinelli T, Caro V, Tran C, Leclercq A, Lecuit M, Brisse S. Worldwide distribution of major clones of Listeria monocytogenes. Emerg Infect Dis. 2011;17(6):1110–2.
Martin B, Perich A, Gomez D, Yanguela J, Rodriguez A, Garriga M, Aymerich T. Diversity and distribution of Listeria monocytogenes in meat processing plants. Food Microbiol. 2014;44:119–27.
Orsi RH, Borowsky ML, Lauer P, Young SK, Nusbaum C, Galagan JE, Birren BW, Ivy RA, Sun Q, Graves LM, et al. Short-term genome evolution of Listeria monocytogenes in a non-controlled environment. BMC Genomics. 2008;9(1):539.
Holch A, Webb K, Lukjancenko O, Ussery D, Rosenthal BM, Gram L. Genome sequencing identifies two nearly unchanged strains of persistent Listeria monocytogenes isolated at two different fish processing plants sampled 6 years apart. Appl Environ Microbiol. 2013;79(9):2944–51.
Verghese B, Lok M, Wen J, Alessandria V, Chen Y, Kathariou S. Knabel S: comK prophage junction fragments as markers for Listeria monocytogenes genotypes unique to individual meat and poultry processing plants and a model for rapid niche-specific adaptation, biofilm formation, and persistence. Appl Environ Microbiol. 2011;77(10):3279–92.
Meier AB, Guldimann C, Markkula A, Pöntinen A, Korkeala H, Tasara T. Comparative phenotypic and genotypic analysis of Swiss and Finnish Listeria monocytogenes isolates with respect to benzalkonium chloride resistance. Front Microbiol. 2017;8:397.
Xu D, Nie Q, Wang W, Shi L, Yan H. Characterization of a transferable bcrABC and cadAC genes-harboring plasmid in Listeria monocytogenes strain isolated from food products of animal origin. Int J Food Microbiol. 2016;217:117–22.
Moretro T, Schirmer BCT, Heir E, Fagerlund A, Hjemli P, Langsrud S. Tolerance to quaternary ammonium compound disinfectants may enhance growth of Listeria monocytogenes in the food industry. Int J Food Microbiol. 2017;241:215–24.
Martínez-Suárez JV, Ortiz S, López-Alonso V. Potential impact of the resistance to quaternary ammonium disinfectants on the persistence of Listeria monocytogenes in food processing environments. Front Microbiol. 2016;7.
Nilsson RE, Ross T, Bowman JP. Variability in biofilm production by Listeria monocytogenes correlated to strain origin and growth conditions. Int J Food Microbiol. 2011;150(1):14–24.
Koreňová J, Oravcová K, Véghová A, Karpí Ková R, Kuchta T. Biofilm formation in various conditions is not a key factor of persistence potential of Listeria monocytogenes in food-processing environment. J Food Nutr Res. 2016;55(2).
Costa A, Bertolotti L, Brito L, Civera T. Biofilm formation and disinfectant susceptibility of persistent and nonpersistent Listeria monocytogenes isolates from Gorgonzola cheese processing plants. Foodborne Pathog Dis. 2016;13(11):602–9.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5.
Zhou Z, Alikhan N-F, Sergeant MJ, Luhmann N, Vaz C, Francisco AP, Carrico JA, Achtman M. GrapeTree: Visualization of core genomic relationships among 100,000 bacterial pathogens. bioRxiv. 2017.
We thank Dr. Alexandra Moura from Institut Pasteur for providing the cgMLST of the L. monocytogenes isolates, Dr. Guy Beauchamp from the University of Montreal for statistical analysis, Paul Manninger from the Canadian Food Inspection Agency for the whole-genome sequences of the isolates and Dr. Austin Markell from the Canadian Food Inspection Agency for critical reading of the manuscript.
This work was funded by the Canadian Food Inspection Agency (Food 1020 and 1231) and Natural Sciences and Engineering Research Council of Canada (NSERC) RDCPJ 520873–17. The funding organisms did not have any influence neither in the design of the study, in the analysis and the interpretation of data, nor in the writing of the manuscript.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
All genomes (except for the strain Lis8316 for which the genome is not published yet as it belongs to another current project) were submitted and published in the GenBank database on the NCBI web site (https://www.ncbi.nlm.nih.gov/assembly) and in the Institut Pasteur database (https://bigsdb.pasteur.fr/listeria/listeria.html) with the accession number of the bioproject PRJNA433177 and detailed in the Additional file 3: Table S2.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Pulsed field gel electrophoresis profile of AscI and Apa1 restriction enzymes in L. 20 monocytogenes strains isolated from a slaughterhouse over a four weeks period. (PDF 144 kb)
Table S1. complete annotation of the high impact SNV and function prediction for the group A (spreedsheet group A) and group B (spreed sheet group B). (XLSX 14 kb)
Table S2. PRJNA433177 Assembly Details and Institut Pasteur IDof the 29 Listeria monocytogenes strains isolated from an industrial slaughterhouse over a 4 weeks period. (XLSX 11 kb)