Abstract
Listeriosis is a serious foodborne disease of a global concern, and can effectively be controlled by a continuous surveillance of the virulent and multidrug-resistant strains of Listeria monocytogenes. This study was planned to investigate prevalence of L. monocytogenes in bovine raw milk samples. A total of 457 raw milk samples collected from 15 major cities in Rajasthan, India, were analyzed for the presence of L. monocytogenes by using standard microbiological and molecular methods. Five of the 457 samples screen tested positive for L. monocytogenes. Multiplex serotyping showed that 3/5 strains belonged to serotype 4b followed by one strain each to 1/2a and to 1/2c. Further virulence potential assessment indicated that all strains possessed inlA and inlC internalins, and, in addition, two strains also possessed the gene for inlB. All strains were positive for Listeriolysin O (LLO) and showed phosphatidylinositol-specific phospholipase C (PI-PLC) activity on an in vitro agar medium with variations in production levels among the strains. A good correlation between the in vitro pathogenicity test and the chick embryo test was observed, as the strains showing higher LLO and PI-PLC activity were found to be lethal to fertilized chick embryos. All strains were resistant to the majority of antibiotics and were designated as multidrug-resistant strains. However, these strains were susceptible to 9 of the 22 tested antibiotics. The maximum zone of inhibition (mm) and acceptable minimum inhibitory concentration were observed with azithromycin, and thus it could be the first choice of a treatment. Overall, the presence of multidrug-resistant L. monocytogenes strains in the raw milk of Rajasthan region is an indicator of public health hazard and highlighting the need of consumer awareness in place and implementation of stricter food safety regulations at all levels of milk production.
Introduction
L
Although pasteurization destroys the bacterium in raw milk, this process does not eliminate the risk of contamination in the dairy chain. In a study from Finland, it was noticed that the dairy products actually prepared from pasteurized milk might become contaminated with the cells of L. monocytogenes at various stages of production (Lyytikäinen et al., 2000). Besides, the ability of L. monocytogenes to survive under harsh conditions and to form biofilm further intensifies the threat to dairy industry. This in turn poses extra economical burden on the industry in terms of analysis cost and possible product recalls.
As on today, 13 serotypes of L. monocytogenes are prevalent, and majority (>95%) of those obtained from contaminated food and infected humans include the 1/2a, 1/2b, and 4b serotypes (Orsi et al., 2011); and of these, 4b is most virulent (Martins and Germano, 2011). Thus, the characterization of isolates at the serogroup level is crucial in the prevention of listeriosis spread.
Of the various available sero-grouping methods, polymerase chain reaction (PCR)-based method is convenient in operation and economical. However, nowadays more superior, accurate, fast, and better discriminatory methods based on whole genome sequencing (WGS) are becoming the method of choice in national reference laboratories to compare patient, food, and environmental isolates than conventional methods (Hyden et al., 2016). Likewise, the characterization of virulence genes (inlA, inlB, inlC and inlJ) responsible for establishing infectious cycle (Bierne et al., 2007) and assessing the production of toxic enzymes (Listeriolysin O, LLO) under in vitro and in vivo conditions in suitable infectious models help understand its virulence and pathogenic potential.
Generally, Listeria strains are susceptible to most of the antibiotics, but indiscriminate antibiotic use has led to the emergence of antibiotic-resistant Listeria spp., which has now become a public health risk. Resistant strains of L. monocytogenes have been obtained from raw milk samples (Jamali et al., 2013) and from dairy farms (Srinivasan et al., 2005). The monitoring of multidrug-resistant pathogens isolated from milk samples could be useful to update the present antibiotic resistance patterns.
India is the largest producer of milk in the world and Rajasthan holds second position among all the Indian states. Furthermore, only 15% of the total produced milk in the country is processed by dairy sector, rest 85% through informal channels (FAO, 2013). Thus, in such cases there are ample chances that the milk might be consumed as raw at homes or traditionally processed into different dairy foods (Doijad et al., 2011). Besides this, the prevalent hygiene levels in dairies and the awareness rate in milk handlers about the safe milk production in India, including Rajasthan are still below the satisfactory level (Lingathurai and Vellathurai, 2010; Sharma et al., 2015). There are no prerequisite standards for quality milk productions prescribed to the farmers in the country.
Earlier many reports concerning prevalence of L. monocytogenes have been published from other parts of India (Barbuddhe et al., 2016), but there is no such report from the state of Rajasthan. Therefore, the present study was designed to investigate (1) the presence of L. monocytogenes in raw bovine milk samples in Rajasthan, India, and (2) to characterize them into serogroups, their genetic variability and virulotypes, and (3) assessing their pathogenicity potential and antimicrobial resistance.
Materials and Methods
Sampling
A total of 457 raw milk samples were collected from different dairy farms situated in 15 major cities in Rajasthan, India, from May 2014 to November 2015. The samples were quickly transported to the laboratory at 4°C and analyzed immediately.
Isolation and identification
A slightly modified method of McClain and Lee (1987) was used for the isolation of L. monocytogenes strains. Briefly, 25 mL of homogenized raw milk sample was added to 225 mL of Listeria enrichment broth (HiMedia, Mumbai India), mixed for 2 min and incubated at 37°C for 24 h. After that, 1 mL of sample was transferred to 9 mL of Fraser broth and incubated for 24 h at 37°C. At the end of incubation, 10-fold serial dilutions (10–1 to 10–5) were prepared in 0.1% sterile peptone water and a 100 μL aliquot from each dilution was poured on the Listeria chromogenic agar supplemented with Listeria selective supplement (HiMedia, Mumbai India). After 24–48 h of incubation at 37°C, typical bluish-green colonies with yellow halo were picked and transferred to brain heart infusion tubes to proceed for Gram staining and catalase test. Presumptive positive results were also tested by the API Listeria Kit (bioMérieux, Marcy-l'Étoile France) for confirmative purposes.
Molecular characterization
Only the presumptively identified L. monocytogenes were processed for PCR. For that, genomic DNA was extracted as described by Pospiech and Neumann (1995). U1/LI1 (938 bp) primers were used for the identification of Listeria spp. (Table 1), whereas primers LM1/LM2 (702 bp, Table 1) corresponding to Listeriolysin O were used for L. monocytogenes identification. The PCR cycling conditions were those as described by the respective authors. A T100™ Thermal cycler (Bio-Rad Laboratories, Gurgaon, Haryana, India) was used for all PCR applications during the whole study. Amplified PCR products were stained with ethidium bromide and visualized using a gel documentation system (Bio-Rad Laboratories).
16S rRNA sequencing
For further confirmation, all PCR identified L. monocytogenes isolates were subjected to 16S rRNA region sequencing by using the bacterial universal primers (Lane, 1991; Table 1) and PCR protocol described by (Rao et al., 2015). The sequencing facility was outsourced from Eurofins Genomics, Bengaluru, India. As per supplied information the amplified PCR products were purified by sodium acetate–ethanol precipitation and a 2 μL of product was sequenced for both primers using the Big Dye Terminator Cycle Sequencing Kit (Applied Biosystem, Foster City, CA, USA). The reactions were performed on an Applied Biosystem ABI 3730 (Applied Biosystem) automated genetic analyzer. The generated partial sequences were then searched for percent homology on NCBI by doing nBLAST and finally submitted to NCBI GenBank for accession number.
Serogrouping
Serogroups were categorized based on the multiplex PCR method (Doumith et al., 2004). The primer pairs used are mentioned in Table 1.
Pulsed-field gel electrophoresis
Subtyping of L. monocytogenes isolates was performed by pulsed-field gel electrophoresis (PFGE) as per CDC PulseNet procedure (Graves and Swaminathan, 2001). The prepared agarose plugs were digested with 25 U of AscI (Thermo Scientific, Waltham, MA, USA) at 37°C for 3 h and 50 U of ApaI (Thermo Scientific). The plugs were then loaded on 1% PFGE certified agarose gel and run in 0.5 × TBE buffer on a CHEF-Mapper System (Bio-Rad Laboratories). A lambda DNA (λ) ladder (Bio-Rad Laboratories) was used as molecular DNA marker. The gels were stained with ethidium bromide and images were captured by Syngene G-Box imagining system (Syngene Bioimaging Private Ltd., Gurgaon, Haryana, India). The obtained PFGE patterns were analyzed using BioNumerics (version 7.6) software package (Applied Maths, Sint-Martens-Latem, Belgium). Dendrogram was prepared by using the Dice band-based similarity coefficient (position tolerance 1.0%) and the unweighted pair group method with arithmetic averages as cluster analysis method.
Virulence determinants
The standardized multiplex PCR method (Liu et al., 2007) was used to determine the virulence genes. The primers used for inlA, inlB, inlC, and inlJ virulence genes are listed in Table 1.
Hemolytic activity
The hemolytic activity in the supernatant of isolates was quantified as described by Sampathkumar et al. (1998) and Upadhyay et al. (2012) with minor modification. The concentration of culture was adjusted to McFarland standard 1.0. Percent hemolysis was calculated as per the formula provided by Bhakdi et al. (1984).
%Hemolysis = (1 − ODs/ODt) × 100, where ODs and ODt represents the differences in optical density at 600 nm between the sample and positive control (100% hemolysis) and between negative control (nonhemolyzed) and positive control, respectively.
Phosphatidylinositol-specific phospholipase C activity
A slightly modified method of Leclercq (2004) was used to ascertain the phosphatidylinositol-specific phospholipase C (PI-PLC) activity of L. monocytogenes. A single colony was streaked on L. mono differential agar (HiMedia, Mumbai, India), and analyzed for the presence of light greenish centered colonies with yellow to transparent halo after 24–48 h (37°C) of incubation (Negi et al., 2015).
Chick embryo toxicity test
The in vivo pathogenicity of the strains was assessed by performing chick embryo toxicity test with minor modifications (Olier et al., 2002). Ten-day-old fertilized eggs were purchased from Khatipura Poultry Farm, Jaipur, India and five eggs were used for one strain. Eggs inoculated only with phosphate-buffered saline (pH 7.2) were taken as control. Strains were considered pathogenic, if they caused mortality within 5 days postinoculation.
Antibiotic susceptibility
Following antibiotic E-strips (bioMérieux, Marcy-l'Etoile, France); penicillin G, vancomycin, ciprofloxacin, rifampicin, gatifloxacin, chloramphenicol, teicoplanin, amikacin, ampicillin, azithromycin, amoxicillin–clavulanic acid, nalidixic acid, gentamicin, piperacillin, oxacillin, kanamycin, tetracycline, ceftazidime, ceftriaxone, norfloxacin, linezolid, and streptomycin were used to determine the minimum inhibitory concentration (MIC) by following the manufacturer's instructions. Briefly, cultures were grown to 0.5 McFarland turbidity; thereafter, a sterile swab soaked in the bacterial suspension was streaked on the Mueller-Hinton agar plates (HiMedia, Mumbai, India). The E-strips were gently placed on the plates and incubated at 37°C for 18–24 h. The strains were categorized resistant (R), intermediately resistant (I) and sensitive (S) on the basis of zone of inhibition(s) (mm) (Clinical and Laboratory Standards Institute (CLSI), 2015). As CLSI has given breakpoints for fewer antimicrobials such as ampicillin and penicillin; therefore guidelines given for Staphylococci and Enterococcus were employed (Conter et al., 2009; Li et al., 2007). The MIC values were interpreted based on the manufacturer's instructions.
Results
Isolation and identification of L. monocytogenes
Of the 457 raw milk samples, 5 samples (1.1%) were found presumably positive for L. monocytogenes on the basis of colony color characteristic and API identification (Table 2). The PCR confirmation of five isolates revealed that all were L. monocytogenes as they give band at the expected size in independent PCR assays.
16S rRNA sequencing
The nBLAST results showed a sequence similarity of 92–97% and all isolates were found closely related to L. monocytogenes (GeneBank Accession Numbers: KX911246–KX911250).
Multiplex PCR serogrouping
Of the five L. monocytogenes, three strains belonged to serogroup “4b, 4d, 4e” and subsequently one each to “1/2c, 3c” and” 1/2a 3a” (Table 3).
Chick embryo test in vivo.
Strains caused mortality in fertilized chick embryos.
Does mortality after 3 days of inoculation.
Does mortality after 4 days of inoculation.
BP, Bharatpur; JPC, Jaipur; LLO, Listeriolysin O; LM, Listeria monocytogenes; PI-PLC, phosphatidylinositol-specific phospholipase C; RM, Raw milk; SH, Sirohi.
PFGE analysis
PFGE analysis by AscI and ApaI enzymes revealed five distinct pulsotypes (Fig. 1). The isolates were considered to have similar pulsotypes when the numbers and position of bands were indistinguishable (Gianfranceschi et al., 2009).

Dendrogram obtained from PFGE profile of AscI and ApaI macrorestriction among Listeria monocytogenes isolates obtained from raw milk samples. PFGE, pulsed-field gel electrophoresis.
Determination of virulence genes
The results demonstrated that all five L. monocytogenes isolates were harboring inlA and inlC genes. In two isolates, besides the aforementioned genes, inlB gene was also present. No isolate was found to have inlJ gene (Table 3).
Hemolysis and PI-PLC activity
Of the five strains, one strain was found to be highly hemolytic (>66% hemolysis), whereas three strains showed moderate hemolysis (>33–66% hemolysis). The remaining one strain showed low hemolysis (below 33% hemolysis) (Table 3). Likewise, strong PI-PLC activity was also demonstrated by one strain only and by others mild to weak activity was shown (Table 3).
Chick embryo pathogenicity test
Out of the five isolates, three were considered as pathogenic as they caused mortality in the chick embryos within 5 days after inoculation (Table 3).
Antimicrobial resistance
Overall, all strains showed 100% resistance to four antibiotics (Table 4) Four, three, and one of the strains showed resistance to three, one, and five antibiotics, respectively. One out of five strains exhibited intermediate level of resistance to three antibiotics. All strains were susceptible to other antibiotics.
Antibiotic concentration impregnated in E-strips: Pencillin G to Gatifloxacin 0.002–32; Chloramphenicol to Linezolid 0.016–256; Streptomycin 0.064–1024.
I, intermediate resistant; R, resistant; S, susceptible.
Discussion
Milk and milk products get inadvertently contaminated with L. monocytogenes at various stages of production; and consumption of such products might always be a risk factor for listeriosis (Jamali et al., 2013). In the present study, L. monocytogenes were isolated, identified, and confirmed from bovine raw milk samples using microbiological media plating, PCR, and 16S rRNA approaches. Earlier, many authors (Cetinkaya et al., 2014; Soni and Dubey, 2014: Jamali et al., 2013) have used these ideal approaches for the isolation, genotypic identification, and validation of L. monocytogenes from different sources and coincide with our findings.
Herein, 1.1% of the tested samples were found to be contaminated with L. monocytogenes. This concords with earlier findings from India on raw milk samples that reported presence of L. monocytogenes (Bhilegaonkar et al., 1997; Barbuddhe et al., 2002; Kalorey et al., 2008; Doijad et al., 2011; Karthikeyan et al., 2015). However, the noticed prevalence is lower than the mentioned studies. Several factors, such as soil, milking equipments, milking environment, sick animals, manure, and poor forage quality at dairy farms, have been recognized as the source of Listeria contamination in raw milk (Oliver et al., 2005).
Serogrouping is important for two reasons: first, it helps epidemiologists screen the involvement of a serogroup in an outbreak; and second in linking the evidences with cases. In the present study, most strains belonged to the serogroup “4b, 4d, 4e” and few to other serogroups (“1/2a, 3a” and “1/2c, 3c”). Similar variations in serogroups among L. monocytogenes strains have been observed in earlier studies (Aurora et al., 2009; Cetinkaya et al., 2014). The very high percentage reported in the present study concerning serogroup “4b, 4d, 4e” concords well with Hamdi et al. (2007), where all of the analyzed raw milk L. monocytogenes isolates belonged to serogroup “4b, 4d, 4e.” Similarly, Karthikeyan et al. (2015) observed that out of 15 L. monocytogenes isolates, eight belonged to serogroup “4b, 4d, 4e.” However, Jamali et al. (2013) reported that majority of L. monocytogenes strains isolated in their case were of serogroup “1/2a, 3a” type.
It has been reported that at a global level, serogroups 1/2a, 1/2b, and 4b are involved in >98% of human foodborne listeriosis (Liu, 2006). Despite the fact that serotype 1/2a is largely retrieved from food samples; the majority of human listeriosis cases have a preponderance of serotype 4b. The disparity noticed between the occurrence of serotype 4b in food and human listeriosis cases suggests that serotype 4b is more virulent. Other serotypes (3a, 3b, 4d, and 4e) have also often been reported from food samples (Cetinkaya et al., 2014). Here, one strain was also found to belong to serogroup 1/2c, which is not frequently detected in early milk epidemiological studies.
To have deeper epidemiological insight among the isolates, all strains were subjected for their genetic diversity analysis using PFGE. From long ago, worldwide PFGE is used by several researchers and diagnostic laboratories for subtyping and comparison of L. monocytogenes strains obtained from diverse regions and sources (Negi et al., 2015). Because of the high discriminatory power and reproducibility of PFGE it is considered as gold standard method for the subtyping of L. monocytogenes (Gerner-Smidt et al., 2006). In this study, PFGE investigation provides discrimination among the isolates of the same serogroup as well as among different serogroups isolated from different regions. Furthermore, the dendrogram analysis revealed that two isolates (LMRM28 and LMRMSH66) belonging to serogroup 4b, 4d, and 4e isolated were very closer in genetic similarity provided the fact that they were isolated from different places. These observations pointed a possible epidemiological association between these strains.
Similar findings were reported by pervious investigators (Aurora et al., 2009; Fox et al., 2012; Negi et al., 2015) who found close similarity among the L. monocytogenes strains of the same serogroup. On the other hand, two isolates (LMRMJPC185 and LMRM79) having closer PFGE profile belonging to different serogroups (4b, 4d, 4e and 1/2c, 3c) were isolated from different areas. Similar findings were reported from others (Mazurier and Wernars, 1992; Aurora et al., 2009), where strains belonged to different sergroups and were isolated from different areas having identical profiles. The remaining one isolate recovered from different area has a totally different PFGE profile than other four isolates. Thus, our results clearly demonstrate the contamination by genetically different strains to the raw milk from surveyed region and that may cause disease in animals and humans. In fact such strains may be liable for posing illness in humans and animals, and their adjoining environment. In future, analyzing of such strains by current modern techniques based on WGS will surely help us in getting a clear-cut resolution that enhance outbreak detection and investigation due to better interpretation, as PFGE only provides limited genetic resolution (Jackson et al., 2016).
Generally, L. monocytogenes strains are considered pathogenic, but their virulence potential varies from one strain to another and depends on certain genetic determinants (Neves et al., 2008). Several methods have been developed to analyze the virulence potential, but among them, PCR-based identification of virulent genes has emerged as a preferable one due to some added advantages. Genes (inlA and inlB) help the bacterium internalize mammalian cells, whereas inlC promotes the intracellular spread of infection. The exact function of inlJ is still unknown. Liu et al. (2007) suggested that inlJ helps in discriminating virulent L. monocytogenes from avirulent ones. However, others have reported some strains in which no inlJ gene is present, but still they cause mortality in mouse through intraperitoneal route (Liu et al., 2006; Roberts et al., 2006). The presence of these genes in strains, however, does not certainly indicate that these strains will cause infections in humans through conventional oral ingestion.
In our study, genes for invasiveness and infection process were present in all five isolates (Table 3). Indeed, in two isolates an additional internalin gene inlB was also present that further increases their affinity for liver and spleen tissues (Bierne et al., 2007). Surprisingly, gene inlJ was not detected among the tested strains; this might be due to a defect or mutation present at primer binding sites resulting in the failure of the primers to amplify the concerned gene (Lomonaco et al., 2013; Negi et al., 2015). Similar findings have also been observed in other studies dealing with different virulence genes (Lomonaco et al., 2013). The present results indicate that these strains have the potential to establish human listeriosis and thereby may cause mortality, if enter through the intraperitoneal route.
Earlier we have noted the presence of hemolysin/LLO-producing gene (hly) by means of PCR, but it might be possible that the protein formed by this gene is truncated. Functional LLO helps the bacterium escape from the phagosomes into the cytosol of the infected cell and protects it from the extracellular immune response. However, this vacuolar disruption is facilitated by PI-PLC group of proteins: PLC-A and PLC-B (Poussin et al., 2009). In this study, all analyzed isolates showed the ability to produce LLO and PI-PLC proteins, which further affirms them being pathogenic. We observed variations in the LLO- and PI-PLC-producing abilities of strains; this disparity in activities might be due to noncoordination of genes responsible for these functions (Shakuntala et al., 2006; Kaur et al., 2010). In addition, hemolysin activity is also interfered during successive subculturing through culture media (Dudhe et al., 2014). The obtained results were in full agreement with those of earlier studies (Momtaz and Yadollahi, 2013; Dudhe et al., 2014; Osman et al., 2014).
We observed that some of the strains tested presented lethal activity in chicken embryos after 5 days of inoculation. Thus these strains indirectly reflected their potential to cause human infections. Similar findings have also been reported by others (Shoukat et al., 2014; Negi et al., 2014).
A successful treatment of listeriosis lies in the early administration of antibiotics with a good bactericidal activity (Hof, 2004). The present observations revealed that the isolates were having high resistance to penicillin G, azithromycin, piperacillin, oxacillin, ceftriaxone, ampicillin, amoxicillin–clavulanic acid, and nalidixic acid and to a lesser extent to ceftazidime, linezolid, tetracycline, teicoplanin, chloramphenicol, and rifampicin (Table 4). These antibiotics are routinely used in the clinical management of listeriosis (Nwachukwu et al., 2010). It can be inferred that this high resistance profile might be due to frequent use of these antibiotics in veterinary practices in controlling infections and to cure human diseases in Rajasthan. The acquisition of mobile genetic elements like self-transferable plasmids and conjugative transposons can be considered as most probable cause of resistance development among these strains.
In addition, environmental factors that bring about mutations in the chromosomal genes might also aid in conferring resistance to L. monocytogenes. Results concerning the resistance of these strains to the aforementioned antibiotics are in line with the earlier findings (Rahimi et al., 2010; Soni et al., 2013; Dudhe et al., 2014). However, all isolates were sensitive to most of the tested antibiotics, similar to other reports (Aureli et al., 2003; Arslan and Özdemir, 2008; Filiousis et al., 2009; Cetinkaya et al., 2014). Interestingly, some authors have also reported vancomycin-resistant strains of L. monocytogenes (Harakeh et al., 2009; Yan et al., 2010). The importance of these antibiotics in the treatment of infections in humans and animals adds to our findings a great value. Furthermore, multidrug-resistant strains of L. monocytogenes have been reported in many previous studies (Jamali et al., 2013; Soni et al., 2013; Dudhe et al., 2014) and again our results in this respect are consistent with those. Importantly, the presence of multidrug-resistant strains in this region puts a question mark about listeriosis treatment and warns a potential threat to public health.
Conclusion
The presence of virulent strains of L. monocytogenes in raw milk samples indicates that unpasteurized milk and traditional dairy products could serve as a vehicle for foodborne listeriosis especially in high-risk groups. The PFGE results obtained in this study may be helpful in reinterpreting the epidemiological prevalence of L. monocytogenes in the studied region and India. Besides, multiple drug resistance noticed against commonly used clinical antimicrobials for listeriosis treatment poses a serious threat to public health. The information generated in this study is highly useful for the veterinarians engaged in controlling the risks of L. monocytogenes and tackling antibiotic resistance among strains.
Footnotes
Acknowledgment
This study was part of a project on milk quality and safety, funded by Rashtriya Krishi Vikas Yojana (RKVY), Government of Rajasthan, India.
Disclosure Statement
No competing financial interests exist.
