Association of Virulence, Biofilm, and Antimicrobial Resistance Genes with Specific Clonal Complex Types of Listeria monocytogenes

Abstract

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.