Epidemiological study on Listeria monocytogenes in Egyptian dairy cattle farms’ insights into genetic diversity of multi-antibiotic-resistant strains by ERIC-PCR

Listeria monocytogenes (L. monocytogenes) is frequently detected in ruminants, especially dairy cattle, and associated with the sporadic and epidemic outbreak of listeriosis in farms. In this epidemiological study, the prevalence, virulence, antibiotic resistance profiles, and genetic diversity of L. monocytogenes in three Egyptian dairy cattle farms were investigated. The risk factors associated with the fecal shedding of L. monocytogenes were analyzed. The L. monocytogenes strains from the three farms were categorized into distinct genotypes based on sampling site and sample type through enterobacterial repetitive intergenic consensus polymerase chain reaction (ERIC-PCR). A total of 1896 samples were collected from animals, environments, and milking equipment in the three farms. Results revealed that 137 (7.23%) of these samples were L. monocytogenes positive. The prevalence of L. monocytogenes in the animal samples was high (32.1%), and the main environmental source of prevalent genotypes in the three farms was silage. For all sample types, L. monocytogenes was more prevalent in farm I than in farms II and III. Risk factor analysis showed seasonal variation in production hygiene. For all sample types, L. monocytogenes was significantly more prevalent in winter than in spring and summer. The level of L. monocytogenes fecal shedding was high likely because of increasing age, number of parities, and milk yield in dairy cattle. Two virulence genes, namely, hlyA & prfA, were also detected in 93 strains, whereas only one of these genes was found in 44 residual strains. Conversely, iap was completely absent in all strains. The strains exhibited phenotypic resistance to most of the tested antibiotics, but none of them was resistant to netilmicin or vancomycin. According to sample type, the strains from the animal samples were extremely resistant to amoxicillin (95.2%, 80/84) and cloxacillin (92.9%, 78/84). By comparison, the strains from the environmental samples were highly resistant to cefotaxime (86.95%, 20/23). Furthermore, 25 multi-antibiotic resistance (MAR) patterns were observed in L. monocytogenes strains. All strains had a MAR index of 0.22–0.78 and harbored antibiotic resistance genes, including extended-spectrum β-lactamase (blaCTX-M [92.7%] and blaDHA-1 [66.4%]), quinolones (qnrS [91.2%], qnrA [58.4%], parC [58.4%], and qnrB [51%]), macrolides (erm[B] [76.6%], erm(C) [1.5%], and msr(A) [27%]), trimethoprim (dfrD [65.7%]), and tetracyclines (tet(M) [41.6%], tet(S) [8%], and int-Tn [26.3%]). ERIC-PCR confirmed that the strains were genetically diverse and heterogeneous. A total of 137 isolated L. monocytogenes strains were classified into 22 distinct ERIC-PCR groups (A–V). Among them, ERIC E (10.2%) was the most prevalent group. These results indicated that environment and milking equipment served as reservoirs and potential transmission ways of virulent and multidrug-resistant L. monocytogenes to dairy animals, consequently posing threats to public health. Silage is the main environmental source of prevalent genotypes on all three farms. Therefore, hygienic measures at the farm level should be developed and implemented to reduce L. monocytogenes transmission inside dairy cattle farms. Supplementary Information The online version contains supplementary material available at 10.1007/s11356-022-19495-2.


Introduction
Listeria monocytogenes is a facultative intracellular Grampositive bacterium, which has been widely explored because of its association with numerous outbreaks of listeriosis worldwide (Vazquez-Boland et al. 2001). It is included in the World Health Organization list of foodborne pathogens and mostly related to raw milk, unpasteurized milk, and other dairy products (Shamloo et al. 2019). L. monocytogenes infections in humans and animals include eye infections, uveitis, keratitis (Nightingale et al. 2004), septicemia, encephalitis, uterine infections (causing abortion and still birth) (Papić et al. 2019), and subclinical mastitis (Constable et al. 2016).
Listeria monocytogenes has been extensively characterized in many animal species although farm livestock have the greatest pretentious (Chow et al. 2021). It is shed by most diseased ruminants, which are asymptomatic carriers, into their environment via their feces (Nightingale et al. 2004). It is a global environmental bacterium present in diverse farm environments, including water, silage, and feces, which are the main resources and potential reservoirs of L. monocytogenes in dairy farms. Consequently, the dairy cattle farms have different genotypes of L. monocytogenes associated with human listeriosis outbreaks (Castro et al. 2018).
The ecology of L. monocytogenes in the farm environment is complicated and inadequately appreciated. Furthermore, factors influencing the persistence of genotypes in dairy farms are unknown. Therefore, our knowledge about the persistence patterns and contamination methods of L. monocytogenes in dairy cattle farms should be improved to elucidate the contaminant source and provide crucial data for reducing transmission and developing intervention strategies against L. monocytogenes at the level of farms and from animals to humans (Walland et al. 2015;Castro et al. 2018;Chow et al. 2021).
Fecal shedding of L. monocytogenes has many risk factors, including inadequate hygiene, sanitation, housing conditions, silage, antibiotic therapy (approved by veterinarians), and atmospheric season (Bandelj et al. 2018). In addition, bulk tank milk (BTM), milk filter, milking machine, milk handler, fecal contamination, poor on-farm hygiene during milking, storage, and transportation are considered contamination sources (Pantoja et al. 2012).
The virulence potential of L. monocytogenes is determined by several molecular determinants (Camejo et al. 2011), which are classified as follows: listeriolysin O (encoded by hlyA), internalins (coded by inlA, inlC, and inlJ), virulence regulator (represented by prfA), actin assembly (coded by actA), phosphatidylinositol-phospholipase C (encoded by plcA), and invasion-associated protein (represented by iap) (Liu et al. 2007). Listeriolysin O (LLO) is a pore-forming toxin produced by hlyA that contributes to the lysis of bacterium-containing phagocytic vacuoles; as a result, bacterial cells are released into the host cytoplasm.
The persistent increase in antibiotic resistance among Listeria spp., notably L. monocytogenes, has been linked to ongoing selective pressure from widespread antibiotic usage in agriculture, animals, and humans during the last few decades (Baquero et al. 2020). L. monocytogenes is resistant to β-lactam, fosfomycin, third-generation cephalosporins, quinolones, erythromycin, tetracycline-minocycline, and trimethoprim (Iwu and Okoh 2020). Multiantibiotic-resistant (MAR) L. monocytogenes strains are determined in food, environmental, and clinical samples (Baquero et al. 2020;Iwu and Okoh 2020;Swetha et al. 2021). Antibiotic resistance develops in these pathogens through gene mutations or acquisition of mobile genetic elements. Antibiotic resistance genes are rapidly being acquired by L. monocytogenes, which may have originated from commensal organisms prevalent in foods and food-processing areas. Some antibiotic resistance genes identified in L. monocytogenes are β-lactamase genes (blaSHV, blaCTX-M, blaOXA, blaIMP, blaCMY, and blaTEM), quinolone-resistant genes (qnrA, qnrB, qnrS, gyrA, and parC), macrolide-resistant genes (erm (A), erm (B), erm(C), erm (TR), mef (A), and msr (A)), and tetracycline-resistant genes (tetA, tetK, tetL, tetM, tetS, and int-Tn) (Baquero et al. 2020). However, additional information about the mechanisms of antibiotic resistance among L. monocytogenes strains should be obtained to design strategies that can prevent the emergence and spread of resistance and to develop innovative therapeutic approaches against multidrug-resistant organisms.
One of the critical points to avoid the distribution of healthcare-associated infections and develop infection control is distinguishing the genetic relatedness among different pathogenic strains. Repetitive element sequence-based polymerase chain reaction (PCR) is a molecular approach used to study repetitive nucleotide sequences within a bacterial genome and to cluster bacterial strains. Enterobacterial repetitive intergenic consensus (ERIC) is one of the repetitive basics characterized by the diverse patterns and numbers of bacterial genomes. L. monocytogenes strains possess short ERIC sequences (Jersek et al. 1999). Enterobacterial repetitive intergenic consensus polymerase chain reaction (ERIC-PCR) is a rapid, dependable, and cost-effective technique for molecular typing that distinguishes the genetic diversity among strains (Parsaie Mehr et al. 2017).
Few epidemiological surveys have been performed on L. monocytogenes in dairy cattle farms in Egypt and its resistance to various antibiotics. Hence, this research was designed to (a) systematically survey the prevalence and distribution of L. monocytogenes in three dairy cattle farms, (b) categorize the risk factors associated with L. monocytogenes fecal shedding in the three dairy farms, (c) determine the virulence and phenotypic and genotypic antibiotic resistance profiles of distinct types of L. monocytogenes strains and (d) establish the genetic diversity of L. monocytogenes strains through Enterobacterial repetitive intergenic consensus polymerase chain reaction (ERIC-PCR).

Study areas
A cross-sectional study was performed to investigate the prevalence of L. monocytogenes and the associated risk factors in the environment of dairy cattle farms in Dakahlia Governorate, Egypt. Three dairy cattle farms were chosen on the basis of their owners' willingness to permit recurrent sample collection. The map of Dakahlia Governorate was constructed to highlight the location of the three selected dairy cattle farms in relation to the rest of Dakahlia (Supplementary Fig. 1).

Questionnaire preparation and data collection
For basic purposes, a structured questionnaire was prepared. The schedule was pretested with necessary adjustment to confirm the importance of the questions and the nature of the sample producers. The dairy farms were selected on the basis of owners' willingness to allow frequent sample collection. The following data were obtained using the questionnaire: general farm information, animal movement, purchase of animals, visitors and staff, health status, and milking process, as shown in Supplementary Table 1. The farmers were asked questions, and their responses were recorded in the questionnaire. Other details, such as the sanitary condition of the farms, were collected through direct observation.

Assessment of production hygiene
The production hygiene of all farms was assessed depending on the premises hygiene, which was evaluated at each visit during the study period. Aspects assessed included milk room, milking station, waiting area, manure passage, resting area, cow cleanliness, feed troughs, and water troughs. Full scores of 1 to 3 were assigned for each of the evaluated sites in the three farms, where "1" referred to a major deficit in production hygiene, "2" denoted a minor deficit in production hygiene, and "3" indicated no notable deficit in production hygiene according to Castro et al. (2018). During farm visits, the production hygiene of each farm was evaluated on the basis of the cleanliness of the premises, particularly overall cleanliness, hygienic design and condition of materials, drainage, air quality, lighting, and insect control.

Collection and processing of samples
A total of 1896 samples were collected from the three examined dairy cattle farms. Among them, 858, 432, and 606 samples were obtained from farms I, II, and III, respectively. Each farm was visited three times (in February, April, and June 2021) so that seasonal variation could be examined. At each farm level, the samples were collected from animals (milk and feces), environment (water, silage, manure, and soil), and milking equipment (teat cup swabs, milk filters, BTM samples, and floor swabs in storage areas). In farm I, 300 animal samples were obtained for each type; 15 environmental samples were acquired for each type; and 150, 30, 12, and 6 samples were obtained from teat cup swabs, milk filters, BTM samples, and floor swabs in storage areas, respectively. In farm II, 150 animal samples were collected for each type; 9 environmental samples were obtained for each type; and 75, 12, 6, and 3 samples were gathered from teat cup swabs, milk filters, BTM samples, and floor swabs in storage areas, respectively. In farm III, 210 animal samples were obtained for each type; 12 environmental samples were collected for each type; and 105, 18, 9, and 6 samples were acquired from teat cup swabs, milk filters, BTM samples, and floor swabs in storage areas, respectively.

Animal samples
A. Milk samples Milk samples were collected before milking early in the morning in accordance with previously described methods (Scha et al. 1971). Briefly, before the samples were collected, udders, particularly teats, were cleansed and dried. Each teat end was cleaned with a pledget of cotton wet with 70% ethyl alcohol. A separate pledget of cotton was used for each teat. The first few streams of milk were discarded. About 10 ml of milk was collected into 15-ml sterile glass vials and labeled as RF (right front), and RR (right rear), LF (left front), LR (left rear).

B. Fecal samples
Fecal samples were obtained directly from the rectum of each cow by using a separate clean plastic sleeve for each sample. The plastic sleeves were inverted, and the content was aseptically transferred into sterile plastic vials (conical 50-ml propylene screw top; VWR International, Inc., West Chester, PA).

Environmental samples
A. Water Water samples were collected from common water troughs in a 50-cm sterile glass bottle (American Public Health Association 1971).

B. Soil
Soil samples were obtained from various sites of yards, especially from the wetted region with high moisture 1 3 and organic matter load at a depth of 5 cm in a sterile glass bottle fitted with a sterile glass stopper (Clegg et al. 1983).

C. Silage
Silage samples were gathered from silage bunkers of the three dairy farms and placed in a sterile Whirl-Pak bag (Vongkamjan et al. 2012).

D. Manure
Manure composite samples were aseptically obtained from various sites in each pen by using a sterile plastic sleeve. The plastic sleeve content was homogenized, and an aliquot was aseptically transferred into 50-ml plastic vials.

Milking equipment
Bulk tank milk samples (100 ml) were aseptically collected. Milk filter samples were obtained and aseptically transported into a sterile sealable plastic bag. For the milking equipment, a sponge was utilized to wipe the inner side of a certain area. Aseptic cotton swabs were used to collect teat cups and floor swabs in the storage area when the routine washing cycle was complete. The sponges were deposited in sterilized bags containing neutralizing buffer, and the cotton swabs were placed in aseptic tubes with 3 ml of neutralizing buffer. Each sample was loaded in coolers with ice packs and transferred to the laboratory for bacterial examination.

Sample analysis for L. monocytogenes
The soil, silage, fecal, and manure composite samples were weighed, diluted with 1% buffered peptone water (BPW), and pummeled (BagMixer; Weymouth, MA, Interscience Laboratories, Inc.) for bacterial analysis as formerly pronounced by Latorre et al. (2010). All samples were pre-enriched by adding 5 ml of BPW or 5 ml of milk to 5 ml of half Fraser broth (Oxoid, Basingstoke, UK) and incubated for 24 h at 37 °C. For the enrichment, 1 ml of the pre-enriched half Fraser broth was added to 9 ml of Fraser broth (Oxoid, Basingstoke, UK) for each sample and incubated at 37 °C for 24 h. Each enriched Fraser broth culture was cultured onto Palcam agar (Oxoid, Basingstoke, UK) and incubated for 48 h at 37 °C. Presumptive L. monocytogenes colonies (gray with black center) were biochemically identified using the following tests: catalase test, oxidase test, evaluation of hemolysis type, motility at 25 °C and 37 °C, in addition to sugar fermentation test (Van Kessel et al. 2004). The strains expressing these standard features were further tested using the API Listeria test (BioMerieux) for confirmation.

Molecular characterization of L. monocytogenes
The genomic DNA was extracted from overnight culture of brain heart infusion broth using the boiled lysate method (Agersborg et al. 1997). Amplification of 16S RNA gene, virulence genes (listeriolysin O (hlyA), positive regulatory factor (prfA)), and adherence (iap) genes for the L. monocytogenes species identification was done by the standard PCR assays in Applied Biosystem, 2720 Thermal Cycler (USA), in a total volume of 25 μL consisted of 12.5 μL of 2 × PCR master mix (Promega, Madison, USA), 1 μL of individual primer (Metabion, Germany), 4.5 μL PCRgrade water, and 6 μL DNA template. The primers and PCR conditions were used as earlier described (Wang et al. 1992;Bohnert et al. 1992;Germini et al. 2009;Soni et al. 2014). The amplified PCR products were organized on a 1.5% agarose gel that was tainted by 1% ethidium bromide and photo-documented under UV illumination. L. monocytogenes ATCC 35,152 strain was utilized as a positive control.

Genetic diversity analysis using ERIC-PCR
The L. monocytogenes strains were genotyped via the enterobacterial repetitive intergenic consensus polymerase chain reaction (ERIC-PCR) fingerprinting assay as described in previous study (Bilung et al. 2018). Genomic DNA was extracted using a QIAamp DNA Mini Kit (Qiagen, Germany). The oligonucleotide forward primer sequence was ERIC1 primer 5′ATG TAA GCT CCT GGG GAT TCAC-3′, and the reverse primer sequence was ERIC1 primer 5′ AAG TAA GTG ACT GGG GTG AGCG-3′ (Versalovic et al. 1991). Each 25 μL of the PCR mixture was composed of 12.5 μL of 2 × PCR master mix (Promega, Madison, USA), 1 μL of individual primer (Metabion, Germany), 4.5 μL of PCR-grade water, and 6 μL of DNA template. PCR amplification was performed with the following thermal cycles (Biometra): primary denaturation at 94 °C for 5 min; 35 cycles of denaturation at 94 °C for 30 s, annealing at 52 °C for 1 min, and extension at 72 °C for 1 min; and a final cycle of 72 °C for 12 min. The amplified PCR products were resolved by electrophoresis on 1.5% agarose gels and photo-documented under UV illumination (Alpha Innotech). ERIC fingerprinting data were converted into a binary code based on the presence or deficiency of each band. Dendrogram was created by the unweighted pair group approach with arithmetic average (UPGMA) and Ward's hierarchical clustering routine. Cluster analysis and dendrogram construction were presented with SPSS, version 22 (IBM 2013) (Hunter 1990). Similarity index (Jaccard/Tanimoto coefficient and number of intersecting elements) among all samples was analyzed by the online tool (https:// plane tcalc. com/ 1664/).

Statistical analysis
Data was recorded using Microsoft Excel spreadsheet (version 15.0), and the investigation was conducted using SPSS (Statistical Set for Social Science) software version 22. Accordingly, descriptive statistics such as percentages and frequency distribution were utilized to determine the prevalence. The antibiotic sensitivity patterns were presented in percentages. Probability level (P) was calculated by using non-parametric test (Chi-square test).

Prevalence and ecology of L. monocytogenes in the dairy cattle farms
Of 1896 samples that were collected, 7.23% (farm I, 3.4%; farm II, 15.2%; and farm III, 6.9%) were positive for L. monocytogenes (Table 1). L. monocytogenes (7.23%) was more predominant in animal samples (32.1%) than in milking equipment (21.9%) and environmental samples (16.8%). L. monocytogenes was determined in all the sampled areas of the three investigated farms except water and floor swabs from the storage area of farm I which were L. monocytogenes negative. Overall, L. monocytogenes was more prevalent in farm I (for all types of samples) than in farms II and III. The prevalence of L. monocytogenes in fecal samples was slightly higher than that in milk samples and varied among the three farms (I, 3% and 2.7%; II, 14.7% and 12%; III, 6.7% and 6.2%, respectively

Associated risk factors of the fecal shedding of L. monocytogenes
The prevalence of L. monocytogenes showed seasonal variation in different samples collected from the animals, environment, and milking equipment of the three dairy cattle farms (Table 2). Analyses revealed that L. monocytogenes was significantly (P = 0.000) more prevalent in all samples obtained in winter (farm I, 6.6%; farm II, 25.7%; and farm III, 13.7%) than in spring (farm I, 3%; farm II, 13.9%; and farm III, 6.8%) and summer (farm I, 0.35%; farm II, 6.3%; farm III, 0%). The seasonal prevalence of L. monocytogenes in all sample categories for each farm was significantly higher (P = 0.043-0.000) in winter than in spring and summer. In winter, the prevalence of L. monocytogenes was the highest in the environmental samples followed by the milking equipment. Conversely, the lowest prevalence was observed in the animal samples, on the level of three farms. The three farms had different hygiene levels ( Table 3). The overall hygiene scores (2.8 to 2.5) of farms I and III were higher than that of farm II (1.4). The rank of each farm according to hygiene score had a strong inverse connection with the prevalence of L. monocytogenes. Differences in production hygiene stated in the questionnaire might be due to the considerably greater incidence of L. monocytogenes in farm II (15.2%) than in farms I (6.9%) and III (3.4%). The fecal shedding of L. monocytogenes was positively associated (P < 0.05) with age, parity, and
All L. monocytogenes strains (100%) had MAR to at least 4 of the 18 antibiotics examined. Their MAR index ranged from 0.22 to 0.78. The highest MAR index of 0.78 was recorded in 24.1% (33/137) isolates from which 27 strains were found in farm II. Furthermore, a MAR index of 0.72 was detected in 16.8% (23/137) of isolates.
The ERIC-PCR dendrogram (Fig. 1) showed that L. monocytogenes obtained from the three farms were genetically diverse and heterogeneous, as indicated by their categorization into specific genotype by sampling site and sample type. The ecological distribution of the 22 ERIC types of L. monocytogenes in the examined farms is shown in Table 5. Furthermore,13, 16, and 12 ERIC types were observed in farm I (A, E-G, I-N, M, Q, S, and T), farm II (A, C, E-N, Q, and S-U), and farm III (B-E, I, L, O, P, R-T, and V). Only ERIC E, I, L, S, and T were common in the three farms. ERIC A, E, F, G, I-N, Q, and S-T were present in farms I and II. ERIC C, E, I, L, S, and T were found in farms II and III. In farm I, the predominant type 1 3 was ERIC F (24%, 7/29), followed by ERIC G (13.8%, 4/29) at which ERIC F was common in milk, feces, silage, and BTM isolates, whereas ERIC G was common in the isolates of feces and silage only. The predominant type in farm II was ERIC M (13.6%, 9/66), followed by ERIC A (12.1%, 8/66), ERIC J (10.6%, 7/66), ERIC Q (10.6%, 7/66), ERIC F (9%, 6/66), ERIC I (7.6%, 5/66), and ERIC C (6%, 4/66), which accounted for 69.7% (46/66) of the isolates in farm II (Table 5). ERIC A was common in isolates of milk, feces, and silage; ERIC F and C were common in feces and silage isolates; and ERIC M was common in milk, manure, and BTM isolates. The five predominant types in farm III were ERIC E (23.8%, 10/42), ERIC C (11.9%, 5/42), ERIC V (11.9%, 5/42), ERIC D (9.5%, 4/42), and ERIC L (9.5%, 5/42). ERIC E and C were evident in the isolates of faces and silage. ERIC V was found in the isolates of milk, manure, and floor swabs of the milking room. It was also common in fecal and water isolates.

Discussion
The epidemiology of L. monocytogenes in clinical human, animal, and food specimens has been widely explored (Castro et al. 2018;Swetha et al. 2021). To our knowledge, this study was the first to elucidate the ecology of L. monocytogenes in dairy cattle farms and identify the sources of the pathogen and risk factors in the environment of dairy cattle farms. Various samples were collected from clinically normal cattle, the environment, and the milking system. The results showed that the prevalence of L. monocytogenes was 7.23% (farm I, 3.4%; farm II, 15.2%; farm III, 6.9%) in all the three studied farms. Mohammed and Abdel Aziz (2017) found a greater prevalence (28.1%) in Egypt possibly because of differences in farm size, management practices, and hygiene ranking. Conversely, the prevalence of L. monocytogenes was low in the USA (4.48%) (Van Kessel et al. 2011), Iran (2.02%) (Sohrabi et al. 2013), and Italy (1.6%) (Bianchi et al. 2013).
The prevalence of L. monocytogenes was higher in the environment than in dairy cattle. Moreover, silage (27.8%) and manure (19.4%) followed by water and soil (8.3%, each) were considered the principal source of L. monocytogenes in the dairy farm environment. These values confirmed the importance of environmental sources of L. monocytogenes in dairy farms, especially silage (Mohammed and Abdel Aziz 2017). The prevalence of L. monocytogenes in milk filters was higher (6.7%) than that in the BTM samples (5.2%) partly because of the concentrations of the bacteria in the filter and dilution in the tanks (Vilar et al. 2007;Bandelj et al. 2018).
Data analysis showed that five risk factors were significantly (P < 0.05) associated with the increased shedding of L. monocytogenes in dairy cattle farms. These risk factors included two environmental factors (season and hygienic condition in the farms) and three animal factors (age of cows, parity, and milk yield). The prevalence of L. monocytogenes had significant seasonal variation in all types of samples (animal, environmental, and milking equipment samples) from the three examined farms. In particular, the prevalence was higher in winter than in spring and summer. Although the analysis of seasonal variation was limited by the short study period, this finding was consistent with those of other investigations (Dalzini et al. 2016; Bandelj et al. 2018) which identified significant seasonal differences and revealed the higher prevalence of L. monocytogenes in cold seasons (winter and early spring) than in other seasons. This high prevalence in cold months could be due to several factors, including the crowding of cattle in indoor facilities, difficulty in maintaining excellent hygiene practices under such circumstances, and the ingestion of spoiled silage (Ryser and Marth, 2007). Other studies have identified insignificant seasonal differences (Hassan et al. 2001;Mohammed et al. 2010). This disparity in seasonal prevalence among different studies might be attributed to variations in study design, climate, and management practices in various geographical regions.
The environmental prevalence of L. monocytogenes in farm II was higher than that in farms I and III. As such, it could be due to poor milking hygiene and contamination pressure from the environment of farms to dairy cattle in our study (Castro et al. 2018).
In the three farms studied, the fecal shedding of L. monocytogenes had a positive significant relationship with the age, parity, and milk yield of dairy cattle. Dairy cattle are affected by several stressors, such as pregnancy, parturition, and lactation, which may suppress host immunity, increase L. monocytogenes fecal shedding in healthy ruminants, and increase their susceptibility to pathogen infection (Roberts and Wiedmann 2003).
The presence of key virulence factors, hlyA and prfA, confers pathogenicity to L. monocytogenes strains (Poimenidou et al. 2018). The hlyA, which encodes listeriolysin O (LLO), is a basic L. monocytogenes pathogenicity gene that helps release bacterial cells through the host cell vacuole (Roberts et al. 2005). PrfA is a protein required for the transcription of the prfA-regulated virulence gene cluster and prfA itself. The iap is a surface protein that acts as a murein hydrolase. In our study, the majority of L. monocytogenes strains (93) were positive for hlyA and prfA although 44 tested strains were positive for only one virulence gene (hlyA or prfA). Conversely, all the tested bacterial strains were negative for iap. These virulence genes are linked to L. monocytogenes strains from clinical and food samples in Ireland (Poimenidou et al. 2018), animals in Egypt (Elbar et al. 2020), and the environment in South Africa (Iwu and Okoh 2020). Consistently, the rates of prfA and hly in L. monocytogenes strains obtained from environmental water in South Africa are high (Kayode et al. 2021). Therefore, the presence of the virulence genes strongly suggested that L. monocytogenes strains from the studied dairy cattle farms could cause listeriosis in humans. MAR L. monocytogenes from different sources, such as human, food, and environmental samples, has been widely described (Castro et al. 2018;Swetha et al. 202). In this work, all L. monocytogenes strains from dairy farms were resistant to penicillin, neomycin, cefoxitin, and nalidixic acid. They were also resistant to amoxicillin, cloxacillin, cefotaxime, amikacin, erythromycin, norfloxacin, tetracyclines, and gentamicin, which are frequently useful in the treatment of human listeriosis. However, they were highly susceptible to chloramphenicol and ciprofloxacin. These findings confirmed previous observations, which demonstrated that L. monocytogenes strains have variable resistance to commonly used antibiotics in the medication of clinical and veterinary infections (Su et al. 2016;Swetha et al. 2021). Tahoun et al. (2017) reported high tetracycline, clindamycin, and rifampicin resistance in L. monocytogenes isolated from an Egyptian dairy farm. However, a previous investigation showed that L. monocytogenes isolated from environmental water has high resistance rates against sulfamethoxazole, oxytetracycline, and amoxicillin but not against ampicillin (Kayode et al. 2021). These differences in the susceptibility patterns of L. monocytogenes strains could be dependent on geographical variations and antibiotic use for humans and animals.
All L. monocytogenes strains exhibited multiple resistances to four classes of antibiotics, namely, β-lactam (particularly second-generation cephalosporins), aminoglycosides, quinolones, and macrolides, which pose risks to public health because of challenges in the treatment of listeriosis. Furthermore, 25 antimicrobial resistance patterns were observed among the MAR L. monocytogenes strains isolated from dairy farms, viewing resistance to antibiotics ranging from 4 to 14. This observation was consistent with the findings of Iwu and Okoh (2020) on higher multiple resistance than single resistance. The observed resistance might be attributed to medication use or feed additives in the livestock industry (Zeitoun et al. 2015). It is significant in the context of the incidence of temporal and spatial changes in antibiotic resistance (Yan et al. 2010). Thus, the emergence of antibiotic resistance should be continuously monitored, and other treatment methods should be developed. In the current work, all L. monocytogenes had a MAR index of > 0.20, indicating that the strains isolated from the three dairy farms originated from high-risk sources in which they were constantly exposed to antibiotics and had a high-risk potential (Bilung et al. 2018).
The mechanisms conferring resistance to different classes of antibiotics were explored to determine whether each MAR L. monocytogenes strain in this investigation had at least one antimicrobial resistance gene. The results showed more prevalence of blaCTX-M gene responsible to produce CTX-M β-lactamases among cefotaxime-resistant strains in comparison to the gene encoding DHA-type β-lactamases (Iwu and Okoh 2020). The first L. monocytogenes strains isolated from dairy farms harboring the plasmid-mediated AmpC β-lactamase DHA-1 gene were designed a few years after the statement of gram-negative bacteria harboring the blaDHA-1 gene. The widespread use of β-lactam antibiotics has resulted in a surge in the occurrence rate of ESBLs because of their low toxicity and effectiveness in the treatment of various infectious diseases, thereby posing a severe threat to global health (Livermore 1996).
Plasmid-mediated quinolone resistance genes (qnrS, qnrA, and qnrB) and parC were observed in quinoloneresistant strains. Though quinolones are not recommended treatment options for Listeria infections, they can indirectly disseminate the emergence of resistant L. monocytogenes strains because of their massive usage for the medication of multiple infections (Godreuil et al. 2003). The resistance of Gram-positive bacteria to quinolones is due to adjustments in the quinolone resistance-determining regions of the intracellular targets of quinolones, DNA gyrase encoded by gyrA and gyrB, and topoisomerase IV encoded by parC and parE (Hooper and Jacoby 2016). The erm (erythromycin ribosome methylase) encodes a 23S rRNA methyltransferase responsible for the modification of the macrolide-lincosamide-streptogramin B (MLSB) antibiotic binding site (Leclercq 2002).
Three MLSB resistance genes, namely, ermB, ermC, and msrA, were determined in the erythromycin-resistant strains. Among them, ermB is the most predominant (Morvan et al. 2010). However, erm (A), erm (TR), and mef (A) encoding the recorded efflux pumps in Gram-positive bacteria were not observed in the tested strains (Leclercq 2002;Granier et al. 2011). The high occurrence of dfrD, encoding a resistant dihydrofolate reductase, has been reported in L. monocytogenes strains from the environment and humans (Morvan et al. 2010). In our study, the prevalence of tet (M) was higher than that of tet (S) that was also noticed by Escolar et al. (2017). The MAR strains exhibited resistance to tetracycline, indicating ribosome protection because of tet (M) or tet(S) (Charpentier and Courvalin 1999). As a determinant of resistance to tetracyclines, tet (M) is highly prevalent in Gram-positive bacteria resistant to tetracyclines and commonly related to conjugative elements of the Tn916 family (Leclercq et al. 2005). In our study, the existence of int-Tn gene in the strains harboring tet (M) confirmed that L. monocytogenes was partly resistant to tetracycline because of the acquisition of conjugative transposons (Poyart-Salmeron et al., 1989).
ERIC-PCR analysis revealed that the L. monocytogenes strains from the three dairy farms were genetically diverse and heterogeneous. Heterogeneity was indicated by the various origins of the strains (animal, environment, and milking equipment) and sampling locations. The potential contamination routes of L. monocytogenes in the dairy cattle farm environment were detected by investigating the prevalence of genotypes among the three farms and sampling areas inside each farm through ERIC-PCR. In the three farms, the sources of predominant genotypes were milk, feces, silage, manure, floor swabs of the milking room, and water, suggesting that various locations in the farm environment could be considered ecological niches for the persistence of L. monocytogenes inside dairy cattle farms (Morvan et al. 2010). This study suggested that the consumption of contaminated silage was the potential source of L. monocytogenes in the investigated farms (Nightingale et al. 2005;Ho et al. 2007). L. monocytogenes contamination in BTM was likely attributed to the milking system (teat cups and milking filters), as observed in other studies (Terentjeva et al. 2021). Our study also demonstrated low similarity in the L. monocytogenes genotypes among the three farms possibly because of the geographical location of the three farms (Dakahlia Governorate). Moreover, strains with the same genotype had diverse antibiotic susceptibility. For example, the strains (n = 7) with a similar ERIC A type had different antibiotic resistance patterns, antibiotic resistance genes, and virulence genes.

Conclusion
For the first time, this work investigated the prevalence of virulent and multi-antibiotic-resistant L. monocytogenes strains in dairy cattle farms in Egypt. Genotyping analysis through ERIC-PCR revealed that all L. monocytogenes strains were diverse and heterogenous, as they were categorized into specific genotypes in terms of sampling sites and sample type. Our study confirms the importance of hygienic practices with respect to silage production and milking hygiene should be developed and implemented to prevent not only the introduction and spread of L. monocytogenes into herds but also its entry in milk. These findings elucidated the epidemiology of L. monocytogenes in dairy cattle farms and served as a basis for implementing control strategies to reduce the risks of L. monocytogenes dissemination in dairy cattle and the environment.
Availability of data and materials All data generated or analyzed during this study are included in this published article and its supplementary information files.

Declarations
Ethics approval and consent to participate This protocol was performed by following the animal ethics guidelines and approved by Medical Research Ethics Committee of Mansoura University with code number (R/75).

Consent for publication Not applicable.
Competing interests The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.