The aim of this study was to investigate the level of strain variability amongst food and clinical Listeria monocytogenes isolates growing at low temperatures (4 and 7 °C) in both laboratory media and real food matrices. Isolates (n = 150) grown in laboratory media demonstrated a large variation in growth profiles measured using optical density. Overall, it was noted that clinical isolates exhibited a significantly higher growth rate (p ≤ 0.05) at 7 °C than the other isolates. Analysis of variance (ANOVA) tests of isolates grouped using Multi Locus Sequence Typing (MLST) revealed that clonal complex 18 (CC18) isolates were significantly (p ≤ 0.05) faster growing at 4 °C than other CC-type isolates while CC101, CC18, CC8, CC37 and CC14 were faster growing than other CC types at 7 °C. Euclidean distance and Ward method-based hierarchical clustering of mean growth rates classified 33.33% of isolates as faster growing. Fast and slow growing representative isolates were selected from the cluster analysis and growth rates were determined using plate count data in laboratory media and model food matrices. In agreement with the optical density experiments, CC18 isolates were faster and CC121 isolates were slower than other CC types in laboratory media, UHT milk and fish pie. The same trend was observed in chocolate milk but the differences were not statistically significant. Moreover, pan-genome analysis (Scoary) of isolate genome sequences only identified six genes of unknown function associated with increased cold tolerance while failing to identify any known cold tolerance genes. Overall, an association that was consistent in laboratory media and real food matrices was demonstrated between isolate CC type and increased cold tolerance.

Precise classification of foodborne pathogen Listeria monocytogenes is a necessity in efficient foodborne disease surveillance, outbreak detection, and source tracking throughout the food chain. In this study, a total of 150 L. monocytogenes isolates from various food products, food processing environments, and clinical sources were investigated for variations in virulence, biofilm formation, and the presence of antimicrobial resistance genes based on their Whole-Genome Sequences. Clonal complex (CC) determination based on Multi-Locus Sequence Typing (MLST) revealed twenty-eight CC-types including eight isolates representing novel CC-types. The eight isolates comprising the novel CC-types share the majority of the known (cold and acid) stress tolerance genes and are all genetic lineage II, serogroup 1/2a-3a. Pan-genome-wide association analysis by Scoary using Fisher's exact test identified eleven genes specifically associated with clinical isolates. Screening for the presence of antimicrobial and virulence genes using the ABRicate tool uncovered variations in the presence of Listeria Pathogenicity Islands (LIPIs) and other known virulence genes. Specifically, the distributions of actA, ecbA, inlF, inlJ, lapB, LIPI-3, and vip genes across isolates were found to be significantly CC-dependent while the presence of ami, inlF, inlJ, and LIPI-3 was associated with clinical isolates specifically. In addition, Roary-derived phylogenetic grouping based on Antimicrobial-Resistant Genes (AMRs) revealed that the thiol transferase (FosX) gene was present in all lineage I isolates, and the presence of the lincomycin resistance ABC-F-type ribosomal protection protein (lmo0919_fam) was also genetic-lineage-dependent. More importantly, the genes found to be specific to CC-type were consistent when a validation analysis was performed with fully assembled, high-quality complete L. monocytogenes genome sequences (n = 247) extracted from the National Centre for Biotechnology Information (NCBI) microbial genomes database. This work highlights the usefulness of MLST-based CC typing using the Whole-Genome Sequence as a tool in classifying isolates.