1. Introduction
Listeria monocytogenes is the primary cause of human cases and listeriosis outbreaks and has a considerable negative economic impact on society and the food industry [
1]. Although
L. monocytogenes is the only recognized human pathogen among
Listeria species, it is also pathogenic for animals [
2,
3].
L. ivanovii is the only other pathogen responsible primarily for listeriosis in animals [
4], but it has been reported to cause listeriosis in humans [
5].
Human listeriosis outbreaks have been documented globally, including the world’s largest outbreak reported in South Africa [
3,
6].
L. monocytogenes causes sporadic cases, protracted outbreaks, and even multi-country outbreaks, and the specific source may not be known [
7]. The European Food Safety Authority [
8] reported 2,183 confirmed invasive human cases of
L. monocytogenes in 2021. In Europe, the case fatality rate is high (13.7%), similar to 2020 [
8], confirming listeriosis as one of the most severe foodborne diseases.
L. monocytogenes is an important foodborne zoonotic agent, and it has been demonstrated to be present in several food types and, therefore, poses a food safety risk [
9]. Meat and meat products constitute a daily human diet due to the high nutritional value of their components, such as proteins, important amino acids, vitamins, and minerals [
10]. However, the nutritional components in meat function as ‘natural media’ for microorganisms such as
L. monocytogenes [
11]. The consumption of ready-to-eat (RTE) meat products has been described as a vehicle for approximately 30% of human listeriosis outbreaks between 2008 and 2015 [
12]. Contaminated RTE meat products are the main concern for public health [
13]. The ability of
L. monocytogenes to survive common food processing conditions like low pH levels, a high salt concentration, low water activity, and refrigeration temperatures facilitates its proliferation in the food environment [
14]. Due to the pathogen’s ubiquity, contamination of meat and meat products occurs at various processes, including RTE products [
15] and distribution stages [
16].
For decades, the traditional serotyping of
L. monocytogenes has been used to characterize isolates recovered from several sources for investigative purposes [
17]. However, researchers and diagnosticians now rely on more sensitive, specific, and accurate molecular methods to diagnose, confirm, and characterize
L. monocytogenes isolates. Some of these methods include the polymerase chain reaction (PCR), multi-virulence-locus sequence typing (MVLST), multi-locus sequence typing (MLST), multilocus variable number tandem repeat analysis (MLVA), pulse-field gel electrophoresis (PFGE) and whole-genome sequencing (WGS), which are now being used [
3,
18,
19].
The sequence types (STs) and the clonal complexes (CC) of
L. monocytogenes have been used to characterize the pathogen [
20,
21], and numerous STs have been identified in
L. monocytogenes isolates worldwide [
22]. Of significance is the frequent association of some STs and CCs with isolates of
L. monocytogenes that are implicated with human listeriosis, thus making the detection of these STs and CCs critical in epidemiological investigations [
21,
23,
24,
25,
26].
The pathogenicity of
L. monocytogenes has been associated with the possession of virulence factors, especially those present in the
Listeria Pathogenicity Islands (LIPIs) [
3,
27,
28]. The virulence factors in the LIPIs play vital roles in the pathogenicity of
L. monocytogenes. For example, the LIPI-1 and LIPI-3 clusters contain genes related to the infectious life cycle and survival in the food processing environment [
28]. The presence of several virulence factors, such as surface-associated internalins, listeriolysin O, and listeriolysin S (LLS) in
L. monocytogenes significantly regulate its pathogenicity [
29,
30].
Antimicrobial resistance (AMR) genes have been documented in
L. monocytogenes isolates are produced to facilitate the development of phenotypic resistance to antimicrobial agents [
31]. Variable frequencies of AMR genes have been reported for
L. monocytogenes isolated from cattle farms, abattoirs, and retail outlets [
32,
33,
34]. The abuse and overuse of antimicrobial agents in human and animal populations result in developing resistance to antimicrobial agents, which is facilitated by the production of appropriate resistance genes as an adaptive response by the pathogen [
35,
36]. It has also been documented that the leading cause of resistance of
L. monocytogenes to antimicrobials is horizontal gene transfer (HGT) of mobile genetic elements such as plasmids and transposons carrying resistant genes and the activation of efflux pump systems [
36].
Plasmids are found in several bacterial pathogens, including
Listeria spp. [
37], and of significance is their ability to carry genetic materials with the potential to encode AMR [
38]. In addition, some plasmids provide other benefits to the host cells with potential contribution to stress survival [
39]. Some of the plasmids detected in
L. monocytogenes include plasmid profiles (N1-011A, J1776, and pLM5578), which were detected in
L. monocytogenes isolates recovered from food processing environments in South Africa [
40].
Proviruses, prophages in bacterial organisms like
L. monocytogenes commonly found in the
Listeria genome, have been reported to play an essential role in bacterial evolution, survival, and persistence [
41]. Prophages/proviruses are known to mediate defense against phage infection through diverse mechanisms in bacteria [
42]. Several frequencies and types of prophages have been reported in
L. monocytogenes from many sources [
40,
43,
44,
45].
The Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) CRISPR-cas system exists in several bacteria, including
Listeria spp., which acts as an adaptive immune system of bacteria, is known to help invade the host immune system [
46]. Several types of CRISPR-Cas have been reported in
Listeria spp., which include Cas-type IA, Cas-type IB, and Cas-type IIA [
47]. In
L. monocytogenes, it has been found that 41.4% of some isolates contain putative
cas genes [
48,
49]. Various CRISPR-Cas systems in
L. monocytogenes isolates recovered from cattle farms, abattoirs, foods, food processing environments, and retail outlets have been found [
46,
50,
51,
52].
South Africa experienced a large outbreak of human listeriosis in 2017-2018 [
6] caused by a strain of
L. monocytogenes, ST6, due to the consumption of ‘polony,’ an RTE pork product [
53]. Earlier reports in the country have documented the occurrence of listeriosis in livestock [
54]. Reports exist using WGS to characterize
L. monocytogenes recovered from the large human listeriosis outbreak [
24], isolates of
Listeria spp. obtained from beef processing environments [
40] and the red meat and poultry value chain [
55]. Most recently, Gana et al. [
56] used WGS to characterize
L. innocua isolates from cattle farms, abattoirs, and retail outlets. To date, there is a dearth of comprehensive information on the WGS analysis of
L. monocytogenes circulating in the beef production chain's three levels (production, processing, and retailing) in Gauteng province, South Africa.
The specific objectives of the current study were, therefore, to apply WGS and bioinformatics analyses to characterize isolates of L. monocytogenes recovered from cattle farms, cattle abattoirs, and retail outlets in Gauteng province to unravel the diversity in the profiles of their sequence types, virulence factors, resistance genes, plasmids, CRISPR-Cas systems, and proviruses. We also investigated the potential effects of the origin of L. monocytogenes isolates (sources and sample/food types) on their profiles.
4. Discussion
In the most recent outbreak of
L. monocytogenes in South Africa, considered the largest in the world,
L. monocytogenes ST6 was determined to be responsible. It was due to consuming contaminated ‘polony,’ an RTE product [
6]. The epidemiology, WGS analysis, and the comparison of South Africa’s outbreak with reports from other countries have been documented [
24,
80,
81]. Beef and beef-based products have been reported to be responsible for listeriosis in other countries [
3,
82]. As a result of the outbreak in the country, WGS analyses have been used to investigate the population structure of
L. monocytogenes isolated in the meat value chain in South Africa [
40,
55]. However, the current study is the first to document the use of WGS and bioinformatics analyses to characterize
L. monocytogenes recovered from the three levels or industries (cattle farms, cattle abattoirs, and retail outlets) of the country's beef production chain.
In our study,
L. monocytogenes was detected at an overall frequency of 6.1% (60/990), comprising significantly different detection frequencies of 3.4%, 4.6%, and 9.3% for cattle farms, cattle abattoirs, and retail outlets, respectively, across Gauteng province. Variable frequencies of
L. monocytogenes have been reported for samples collected elsewhere from cattle farms, cattle abattoirs, and retail outlets [
12,
33,
83,
84,
85]. The differences in the frequencies of
L. monocytogenes across countries may be due to different management practices in the three industries and the prevalence of the the pathogen in these countries.
Of the
seven STs (ST1, ST2, ST14, ST31, ST204, ST224, and ST876) detected in our study, three were found at comparatively high frequencies for ST876 (11.7%), ST2 (21.7%), and ST204 (46.7%). Our findings are different from those reported for 217
L. monocytogenes isolates recovered from red meat and poultry value chain in South Africa [
55], where a total of 20 STs were detected, comprising ST204 (14.7%), ST2 (13.8%), ST1 (11.5%), ST9 (11.1%), and ST321 (9.7%). It is pertinent to mention that the current study and two other studies, one conducted on the food chain [
40] and the other on meat and meat products [
55], all of which were after the large human listeriosis outbreak of 2017-2018 [
6] failed to detect ST6 of
L. monocytogenes which was responsible for the outbreak. However, it cannot be over-emphasized that ST204, the most frequently detected ST in the three studies, may pose a potential food safety concern regarding human listeriosis in the country since it has been associated with cases of human listeriosis elsewhere [
21,
86,
87]. It has been documented that ST204 is the most common ST in meat products in Australia and France [
88,
89] and in food processing plants [
90]. Furthermore, other STs observed in our study have been detected in meat and meat products and other foods implicated in cases and outbreaks of listeriosis by others [
25,
89,
91].
It has been documented that
L. monocytogenes strains are delineated into sequence types (STs) based on conventional multilocus sequence typing (MLST), which utilizes seven alleles. STs are then grouped into clonal complexes (CCs) with strains in the same sharing at least six of the seven MLST alleles [
21]. The clonal complexes (CC) and lineages to which isolates of
L. monocytogenes have been determined to predict the vir-ulence or pathogenicity potential of the microorganism recovered from human cases or foods [
26,
92,
93]. This association of CCs and lineages of
L. monocytogenes with virulence is linked to the type and number of virulence factors they carry. It is, therefore, it is significant that
L. monocytogenes isolates allocated to CC1 and CC2in lineage, I were detected in 24 (40%) of the isolates, which have been reported to be frequently associated with human listeriosis [
26,
94,
95,
96], and CC204 in lineage II constituted 46.7% of our isolates and has been documented to be predominantly documented in foods [
40,
55] in South Africa and elsewhere [
97,
98]. It cannot be over-emphasized that in our study, all 60 isolates of
L. monocytogenes were carriers of 5 LIPI-1 virulence factors (
prfA,
plcA,
hly,
mpl, and
plcB) and LIPI-3 cluster, which are known to play a significant role in the virulence/pathogenicity of CC1 and CC2
L. monocytogenes [
26,
99,
100] were carriers of 8 virulence factors (
IlsA,
IlsB, I
lsD,
IlsG,
IlsH,
IlsP,
IlsX, and
IlsY) were detected in 13.3-20% of our isolates. This is relevant to food safety, considering that 31.7% (19/60) of the CC1 and CC2
L. monocytogenes in the The current study was based on beef and beef products, including RTE products. Unsurprisingly, CC204 was detected at an overwhelmingly high frequency (40.5%-58.3%) across samples from the three industries (cattle farms, abattoirs, and retail outlets) and more importantly, 40.5% for beef and beef products.
The distribution of STs of
L. monocytogenes within and across the three industries was significantly different, demonstrating that the industries were significantly associated with the STs detected. However, it is pertinent to mention that our findings that the industries were significantly associated with the STs detected may be limited to the current sampling scope, including the locations and sampling span. This is because the STs/CCs of L. monocytogenes are known to be frequently introduced and transmitted; therefore, some are expected to be found in only one location in one sampling effort [
21]. A cross-sectional or ‘snapshot’ study, like ours, will be unable to infer the persistence of the CCs in that location over multiple years, thus limiting our study. It was also important to observe that some STs (ST31, ST204, and ST876) were distributed across all the industries. At the same time, ST224 was found exclusively in abattoir isolates, ST1 and ST14 were detected only in the isolate from the retail industry, and ST2 was uniquely shared between the
L. monocytogenes isolates obtained from the Abattoir and Retail industries. The differences in the number and types of STs recovered in the three industries (cattle farm: 3 types, abattoirs: 5 types, retail industry: 6 types) may be explained in part by the number of isolates tested per industry, 11, 12, and 37, respectively. Furthermore, the variation in the number of sources could have contributed to the recovery of isolates of
L. monocytogenes (cattle farms versus abattoirs versus retail outlets). It was also significant that ST204 was predominantly detected in three industries. Other reports have documented differences in the number and frequency of STs in
L. monocytogenes from these industries by others [
101,
102]. Just as found with the effects of industries on the distribution of
L. monocytogenes STs, the frequency of STs varied significantly according to the eight sample/food types tested. Interestingly, ST204 was detected at the highest frequency across all the sample/food types tested. Therefore, there is a possibility that ST204 is widespread among
L. monocytogenes isolates in Gauteng province, as earlier documented by others [
40,
55]. It is also noteworthy that some STs, such as ST2, were found in as many as five sample/food types (HT abattoirs and the four types of retail outlets), whereas ST224 was found only in HT abattoirs. It cannot be over-emphasized that the distribution of the STs of
L. monocytogenes, according to the sample/food types, can influence human exposure to some of these STs [
47,
103].
It is noteworthy that the LIPI-3 genes were detected at frequencies ranging from 13.3% (Ils
P) to 20% (
IlsA, IlsB, IlsG, IlsH, IlsX, and
IlsY) of our isolates. This is because the LIPI-3 gene cluster is known to be involved in the infectious life cycle and survival in the food processing environment [
27,
28]. It has been documented that these virulence factors perform different roles and functions, such as being responsible for surface protein anchoring, adherence, invasion, immune modulation, and intracellular survival, among others; some virulence factors have been implicated in human listeriosis [
3,
29,
103,
104,
105]. Equally relevant is the finding that among our 60 isolates of
L. monocytogenes, 26 (59.1%) shared unique virulence factors, including virulence factors belonging to LIPI-1 and LIPI-3 genes. Matle et al. [
106] similarly reported the presence of 47 similar virulence factors in six sequenced isolates of
L. monocytogenes from RTE products in the country.
The current study's detection of 44 virulence factors is considerably lower than the 68 putative virulence factors earlier reported in the country [
40]. The differences in the critical virulence factors detected in both studies conducted in South Africa may be accounted for partly by the origins and types of the samples from where the
L. monocytogenes were isolated, and different
L. monocytogenes populations resident in these locations. In our study, the
L. monocytogenes isolates originated from cattle farms (faeces, feeds: grain, grass, and silage), cattle abattoirs (pre-evisceration and post-evisceration carcass swabs, chilled carcass swabs, and effluent), and from retail outlets (raw beef, offal & organs, milled beef, and RTE) in Gauteng province. On the contrary, the
L. monocytogenes isolates analyzed by Mafuna et al. [
40] were initially recovered from raw meat, processed meat, RTE meat products, and environmental samples collected from a commercial pig farm environment during a listeriosis outbreak [
55]. Studies on
L. monocytogenes recovered from meat, meat products, and other foods by others have similarly reported variable types and frequencies of virulence factors[
93,
107].
The clustering of virulence factors within and across the three beef industries and sample/food type, as well as the MST based on their STs, is not surprising but indicates that these variables affect the consumer’s exposure potential to virulent isolates of
L. monocytogenes in agreement with published reports [
108,
109]. Our finding of food safety and therapeutic importance is that RTE beef products were carriers of virulent
L. monocytogenes STs were also resistant genes. RTE products have been documented to be implicated in most human listeriosis cases or outbreaks [
3,
12]. Relevant to South Africa is that one of the three RTE beef products (Vienna, ‘polony,’ ‘biltong’) are two popular delicacies in the country. ‘Polony’ is a beef RTE product that was implicated in the recent large outbreak of human listeriosis in the country [
6]. Biltong is essentially raw meat spiced and dried and has been reported to be contaminated with
Salmonella [
110] and Shiga toxin
Escherichia coli [
111] in the country. It is, therefore, a source of concern that ‘polony’ constituted 5 of the 7 RTE products where isolates of
L. monocytogenes were carriers of 32-39 virulence factors, of which ST204 occurred at the highest frequency (3/5), which has been implicated in human listeriosis [
86], the detection of CCs 1 and 2 (2/5), 6 LIPI-1 cluster genes (5/5), and were all positive for the
fosX and
vga(G) AMR genes. Our findings further confirm the potential food safety and therapeutic importance of the overwhelming detection of ST204 across all our samples, considering that the isolates found in the RTE products carried several virulence genes. Matle et al. [
106] detected 142 virulence genes across the sequences of six RTE isolates, which are considerably higher than found in the RTE products in the current study. Other reports have documented the contamination of RTE meat products with
L. monocytogenes and characterized the isolates regarding their virulence factors and resistance genes [
3,
93,
107,
112].
It is interesting that only two AMR genes,
fosX, and vga(G), which encode phenotypic resistance to fosfomycin and lincosamide, respectively, were detected in our study, having each been found in all the 60 isolates of
L. monocytogenes studied, regardless of the industries and sample/food types. This is not a surprise because, in South Africa, some antimicrobial agents, including fosfomycin, tetracycline, and sulphonamides, are legally allowed to be sold over the counter for use in the livestock industry. This has been made possible through the Fertilizers, Farm Feeds, Agricultural Remedies and Stock Remedies Act 36 of 1947. Therefore, these farmers use antimicrobial agents to treat livestock and as growth promoters in the [
113,
114,
115] without veterinary oversight. Thus, our study's exceptionally high detection frequency of the
fosX gene (100%) may have therapeutic implications. In agreement with the current study was the detection of the
fosX gene in all (100%)
L. monocytogenes isolates recovered RTE products of animal origin (100%) [
99] and from the food chain [
40] in South Africa. However, unlike our current study, where only two AMR genes,
fosX and
vga(G), were detected, Mafuna et al. [
40] found four resistance genes (
fosX,
lin,
norB, and
mprF) in all isolates of
L. monocytogenes and
tetS obtained from the country's meat food chain. The origin of the samples and isolate genotypes in both studies may account for the differences in their findings. Additionally, it is known that different genotypes can have different AMR gene profiles [
116]. The detection of the AMR resistance gene,
fosX, in all (100%)
L. monocytogenes isolates in the current study may be due to a report based on the analysis of 1,696 isolates of
L. monocytogenes revealed the
fosX gene to be part of the
Listeria core genome, where all isolates harbored this gene [
117]. Mota et al. [
118] have also indicated that
L.monocytogenes is currently considered to be intrinsically resistant to fosfomycin because of the lack of expression in the membrane transport systems and natural resistance to lincomycin due to the ribosomal protection of an ATP-binding cassette F (ABC-F) protein. It is, therefore, no surprise because Parra-Flores et al. [
52] detected the
fosX gene in 100% of
L. monocytogenes recovered from RTE foods, and Hanes et al. [
119] reported that 97.8% of
L. monocytogenes isolated from 2010-2021 in the USA were carriers of the gene. Regardless of whether the high frequency of the
fosX gene resulted from the overuse of antimicrobial agents in the cattle industry or the genes being part of the genome of
Listeria spp., it is important to note that the findings may pose therapeutic complications should the gene be expressed. However, not all genes, including resistance genes, are expressed because they may be lost, limiting their application to assessing their therapeutic implications and significance [
120]. The potential therapeutic effect of the high frequency of
vga(G), which encodes lincosamide resistance, cannot be assessed because the antimicrobial agent is not routinely used on cattle in South Africa.
In our study of 60 isolates of
L. monocytogenes, only one AMR plasmid, NF033156 was found at a low frequency of 5% (3/60), and all were from ST204 isolates recovered from an abattoir and two from retail outlets. Matle et al. [
106] found no plasmid in a study on six isolates of
L. monocytogenes recovered from RTE meat products. Of relevance is that Mafuna et al. [
40] detected plasmids in 71% of the 143 isolates of
L. monocytogenes studied, and their detection was ST-specific. Although there were differences in the types and frequencies of AMR plasmids identified in both studies, the over-representation of the plasmids in some STs is common. The differences in the types and frequencies of plasmids recovered from
L. monocytogenes in both studies may be due to the origins of the samples that yielded the pathogen and the losses of plasmids. Plasmids are essential in the carriage of AMR genes and other genetic materials in
L. monocytogenes and other bacteria [
121,
122]. Furthermore, different types of AMR plasmids have been identified in
L. monocytogenes in meat products in studies by others, and it has been demonstrated that
L monocytogenes may gain or lose plasmids [
123].
Three conjugative plasmids were detected in 80% of the 60 isolates of
L. monocytogenes tested in our study, with their detection at statistically significant different frequencies for MOBP2 (38.3%), MOBV (16.7%), and FA_orf13; FA_orf17b (5%). Notably, the occurrence of the plasmids was associated with the STs and the in- dustry, with our observation that both MOBV and MOBP2 were associated with the STs and the beef retail industries. In agreement with our findings, Mao et al. [
124] demonstrated that a conjugative plasmid, pLM1686, was associated with four STs (ST87, ST59, ST9 and ST120) in China. The authors also reported that the plasmid had the self-transmissible ability and existed in various
L. monocytogenes isolates, providing
L. monocytogenes advantages of surviving in adverse environments. Others have described the types and roles of conjugative plasmids in
L. monocytogenes [
124,
125].
In our study, all (100%) 60 isolates of
L. monocytogenes were carriers of proviruses/ prophages. This is higher than the 90.9% (30/143) found in
L. monocytogenes isolated from the food chain in the country [
40]. A considerably lower likelihood of prophage- carrying isolates of
L.
monocytogenes was 14.4%, as suggested by Vu et al. [
45], who also indicated that the prophages in
L. monocytogenes are highly diverse and could be at least 16. The detection that all
L. monocytogenes isolates in our study were positive for provirus/prophages is important, considering that they are known to play critical roles in
L. monocytogenes including mediating defense against phage infection, bacterial survival, and persistence in stressful environments [
41,
42]. Interestingly, proviruses (prophages) in the class Caudoviricetes were detected in all 60 isolates of
L. monocytogenes isolates in our study, their co-location with AMR gene (
fosX), and being ST-specific (ST31 and ST204) indicate that the provirus /prophages may serve as the vector for the
fosX gene. This finding requires further investigation. It is pertinent to mention that it has been suggested that the
fosX gene is harboured on the
Listeria core genome [
126].
We detected the CRISPR-Cas subsystem (Class1-Subtype-I-B_1) in 10% (6/60) of the
L. monocytogenes isolates, which were evenly distributed across the three industries but were highly represented in ST31, 83.3% (5/6). Parra-Flores et al. [
52] recovered the CRISPR-Cas system from 71% of RTE foods sampled in Chile. Although, for comparison, no prior reports have been published on the occurrence of the CRISPR- Cas system in
L. monocytogenes isolated from any source in the country, Mafuna et al. [
47] detected three CRISPR-Cas system types (CAS-Type IIA system, CAS-Type IB system, and CAS-Type IIC system) in 41 non-pathogenic
Listeria spp. (
L. innocua and
L. welshimeri) recovered from meat and food processing facilities (FPF). Regardless of the
Listeria spp., the CRISPR-Cas system is known to degrade foreign genetic elements and has been documented in
L. monocytogenes isolates [
127,
128,
129]. It has been suggested that some of the CRISPR-Cas systems detected in L. monocytogenes are functional with spacers matching sequences of known
Listeria phages and plasmids [
129]. The CRISPR-Cas system, which exists in
L. monocytogenes, acts as an adaptive immune system and has been documented to help invade the host immune system [
46].
The data from our study demonstrate the occurrence of virulence factors (n=44) and STs, which were widely distributed across the three industries and eight sample/food types investigated. Others have documented similar findings [
40,
55,
93,
107]. The significantly higher distribution of certain STs and virulence factors in some food types, particularly RTE products, increases the risk of listeriosis in humans. The AMR gene,
fosX, was present in all 60 isolates of
L. monocytogenes from the country's beef production chain may pose a therapeutic threat since Fosfomycin, encoded by the gene, is commonly used in the country’s livestock industry. Isolates of
L. monocytogenes may pose a therapeutic threat since Fosfomycin, encoded by the gene, is widely used in the country's livestock industry. The AMR gene, the CRISPR-Cas system, and proviruses/prophages are known to play some roles in the survival of
L. monocytogenes in the food system [
38,
39,
46].