Ligilactobacillus murinus Strains Isolated from Mice Intestinal Tract: Molecular Characterization and Antagonistic Activity against Food-Borne Pathogens

Considering the objectives of “One Health” and the Sustainable development Goals “Good health and well-being” for the development of effective strategies to apply against bacterial resistance, food safety dangers, and zoonosis risks, this project explored the isolation and identification of Lactobacillus strains from the intestinal tract of recently weaned mice; as well as the assessment of antibacterial activity against clinical and zoonotic pathogens. For molecular identification, 16S rRNA gene-specific primers were used and, via BLAST-NCBI, 16 Ligilactobacillus murinus, one Ligilactobacillus animalis, and one Streptococcus salivarius strains were identified and registered in GenBank after the confirmation of their identity percentage and the phylogenetic analysis of the 16 Ligilactobacillus murinus strains and their association with Ligilactobacillus animalis. The 18 isolated strains showed antibacterial activity during agar diffusion tests against Listeria monocytogenes ATCC 15313, enteropathogenic Escherichia coli O103, and Campylobacter jejuni ATCC 49943. Electrophoretic and zymographic techniques confirmed the presence of bacteriolytic bands with a relative molecular mass of 107 kDa and another of 24 kDa in Ligilactobacillus murinus strains. UPLC-MS analysis allowed the identification of a 107 kDa lytic protein as an N-acetylmuramoyl-L-amidase involved in cytolysis and considered a bacteriolytic enzyme with antimicrobial activity. The 24 kDa band displayed similarity with a portion of protein with aminopeptidase function. It is expected that these findings will impact the search for new strains and their metabolites with antibacterial activity as an alternative strategy to inhibit pathogens associated with major health risks that help your solution.


Introduction
Zoonotic diseases, food safety dangers, and bacterial resistance are related and affect human, animal, and environmental health. Hence, the "One Health" comprehensive approach acknowledges the significant need to develop effective strategies to solve these global health issues [1]. In that sense, food-borne diseases, especially zoonotic diseases, are the cause of morbidity and mortality globally. In accordance with the European Food Safety Authority, Campylobacter spp. is the main zoonotic agent that causes diarrheic disease, with 120,946 cases, 8605 hospital admissions, and 45 deaths, which represents over 60% of the cases reported as of 2020 in the European Union. Enteropathogenic E. coli is acknowledged as the most important cause of acute diarrheal disease in children under four years old, and Shiga toxin-producing strains of E. coli (STEC) are highly virulent and associated with these pathogens of clinical and zoonotic interest. Furthermore, the presence of enzymes with broad-spectrum bacteriolytic activity can be safely applied as promoters of animal health growth as an alternative to the use of antibiotics for the control of antibiotic-resistant pathogenic bacteria.

Isolation and Selection of LAB
LAB isolation was performed from just-weaned CD1 mice collected from the Animal Facility, Multidisciplinary Research Unit in FES Cuautitlan, UNAM. The project was evaluated and approved by the CICUAE-Institutional Committee for the Care and Use of Experimental Animals with protocol number C18_20, and the euthanasia was performed in accordance with NOM-062-ZOO-1999 [22]. Under sterile conditions, the duodenum, jejunum, and colon were extracted from each mouse in order to obtain the intestinal contents by washing through the lumen of each intestine portion with sterile peptone water (NaCl 8.5 mg·mL −1 , proteose peptone 1 mg·mL −1 ). The collected intestinal matter was inoculated in Lactobacilli MRS broth (De Man Rogosa and Sharp, BD Difco TM , USA), in proportion (1:10) and incubated at 30 ± 2 • C for 24 h in a jar anaerobic system (BD GasPak TM , USA). Subsequently, streak seeding was performed in MRS agar incubation at 30 ± 2 • C for 48 h per intestine portion. The round, white colonies of 1-2 mm that, when stained, were Gram-positive were selected.

Strain Handling
Each of the isolated Gram-positive Bacillus strains was cultivated in MRS broth at 37 ± 2 • C for 24 h, with a pass of 2% (v/v) in the same medium incubated for 12 h. The cells were collected by centrifugation at 5000× g for 10 min at 4 • C (Centurion Scientific K2015 centrifuge, UK), resuspended in MRS broth and skim milk at 10% (w/v), ultrafrozen at −50 • C (Thermo Scientific Revco Ultima II, USA), and freeze-dried (Labconco ® FreeZone 4.5L, KC, USA) for their preservation at the long term and in MRS-glycerol broth (50-50% v/v) for their use in subsequent tests.

Genomic DNA Extraction
Twelve-hour cultures from each one of the LAB strains were carried out for DNA extraction through a phenol-chloroform method [24], verifying its purity and concentration in a spectrophotometer (NanoDrop ® ND-1000, USA) (A 260-280 nm ), in addition to its integrity in agarose gels at 1.5%. The DNA of L. lactis ATCC 0205V and L. rhamnosus ATCC 53103 strains equally cultivated were used as control.

Primer Sequences Design
Species-specific primers for BAL strains were designed (Table 1) from the strains reported in the literature, usually found in mouse intestinal microbiota. The partial or total sequences of the 16S rRNA gene were obtained from the GenBank from each of these strains. These sequences were used for the design of primers in the Primer-BLAST software by setting the following parameters: 18 to 25 nucleotides length, percentage of guanine-cytokine (%GC) between 40% and 50%, and difference of fusion temperature less than or equal to 5 • C. The designed primers were synthesized by Sigma Aldrich Química S de RL de CV.

Sequence Analysis and Identity Percentage Matrix
The 16S rRNA sequences of each bacterial isolation obtained from the chromatograms were edited with BioEdit 7.2 software, and the sense and antisense sequences of each strain were aligned with ClustalW to create their consensus sequence. Every consensus sequence was aligned with BLAST from NCBI and compared with the sequences from the GenBank database. The similarity percentage between the sequences was obtained with multiple alignments in ClustalW and the option of an identity matrix.

Phylogenetic Tree Building
The sequences of the 18 isolated LAB strains were used for the building of phylogenetic trees, in addition to 12 reference sequences obtained from the alignment in BLAST-NCBI software. The phylogenetic analysis was performed in MEGA11 software, with the Maximum-Likelihood method of the suggested Jukes-Cantor model [25]. Confidence levels were determined by bootstrap analysis with 1000 repetitions.

Antagonistic Activity by Agar Diffusion Tests
Inhibitory activity was detected by diffusion tests in agar plates [26], with some modifications. MRS-based agar plates (BD, Difco TM ) at 1.5% (w/v) were used, upon which 5µL aliquots of a logarithmic phase culture of each isolated LAB were placed, corresponding to 7-9 log CFU mL −1 . In parallel, for each test, a 10 mL soft agar overlay (0.6% w/v), BHI agar for L. monocytogenes, LB agar for E. coli, and Brucella agar for C. jejuni was prepared, added with 100 µL from a 12 h culture of each test pathogen. Each overlay was conducted at 45 • C, and it was poured over every MRS-based agar plate. Solidified plates were incubated at 37 ± 2 • C for 24 h, and the test was conducted with C. jejuni. The plates were incubated with the use of the Anaerocult ® C system (5-7% O 2 , 8-10% CO 2 ) at 40 ± 2 • C for 48 h.

Bacterial Growth Parameters
The growth kinetics of every isolated strain was created from new cultivation in MRS broth at 37 ± 2 • C, and OD and CFU mL −1 were obtained at 0, 4, 8, 12, 16, 20, and 24 h. Growth parameters were calculated with the CFU mL −1 from the logarithmic phase and the following formulas [27]: Specific growth rate (µ): means the growth rate per unit of CFU or biomass  Cultures of every LAB strain in logarithmic phase (4 h), in 1 mL MRS broth, in Eppendorf tubes, were centrifuged at 5000× g, for 10 min at 4 • C (Centurion Scientific K2015, UK), neutralized at pH 7 (NaOH 6N) and sterilized by membrane filtration 0.22 µm (Merck Millipore TM , Ireland) [28].

Protein Quantification
Proteins from Sterile-filtered, neutralized supernatants were precipitated with the Methanol/Chloroform method [29]. The protein concentration of each sample was determined by the Bradford method [30], with Bradford reactive (Bio-Rad, cat. no. 5000205, USA). In every test, a standard curve was built via Bovine Seric Albumin 1 mg·mL −1 (Bio-Rad cat no. 5000206).

Electrophoretic Profile of Supernatants Proteins (SDS-PAGE)
Electrophoresis was performed on denaturing and reducing polyacrylamide gels at 12% (30%T,2.67%C) [31]. Fifteen µL of the sample containing the precipitated proteins (1 µg/µL) and 5 µL of Laemmli sample buffer 4× (Bio-Rad, cat no. 1610747) were heated to 95 • C for 5 min and placed in each well of the gel. Precision Plus Protein TM All Blue molecular Weight marker (Bio-Rad, cat. no. 1610373) was used. The gels were run in an electrophoresis chamber (MiniPROTEAN III, Bio-Rad) at 90 volts for 180 min and 20 mA (Bio-Rad, PowerPac Basic). After electrophoresis, the gels were washed with sterile distilled water, stained with Coomassie brilliant blue (Bio-Rad cat. no. 1610400) for 1 h, and destained with a solution (20% methanol (J.T. Baker ® , USA), 15% acetic acid (J.T. Baker ® , USA) and 65% deionized water) for 12 h. Gels were analyzed in the Gel Logic 100 imaging system (Kodak, USA) to establish the electrophoretic profile of protein bands.

Zymograms
The zymography tests were carried out with the cells obtained from 60 mL of a culture of L. monocytogenes, E. coli REPEC, or C. jejuni, washed twice with sterile saline solution (0.9% w/v), which were used as a substrate when adding them in 12% polyacrylamide mixtures. Micrococcus lysodeikticus ATCC 4698 was used as a positive control at 0.2% (w/v) for the identification of PGH activity [32]. Precipitated proteins of every LAB strain were placed in the wells, and the SDS-PAGE was carried out under the same conditions previously described. When the electrophoresis was finished, the gels were disassembled and put in a container with deionized water for 30 min, proceeded by incubation in a renaturing solution (25 mM Tris-HCl pH 8 with 1% v/v Triton) for 6 h with the purpose of identifying lytic bands. In order to obtain a better definition, the gels were stained with methylene blue at 1% (w/v) in KOH solution at 0.1% (w/v) for 30 min and were subjected to destaining with distilled water [33].

UPLC-MS Fingerprinting of L. murinus B1 Lytic Proteins
The protein bands that demonstrated lytic activity separated in SDS-PAGE gels, stained with Coomassie brilliant blue G-250 0.02% (w/v), were cut out of the gel with a sterile scalpel and transferred to an Eppendorf tube. The preparation of the samples was carried out according to the protocol developed at the Research and Industry Support Services Unit (USAII), Faculty of Chemistry, UNAM [34]. The gel fragments were incubated in 100 µL of destaining solution (80 mg ammonium bicarbonate, 20 mL analytically pure water (J.T. Baker ® ), and 20 mL acetonitrile) for 15 min in an orbital shaker at 50 rpm.
To reduce and alkylate the proteins in the samples, the liquid phase was discarded, and 5. Enzymatic digestion was performed with 1 µL of trypsin (Promega cat #V528A) (1µg/µL) and 100 µL of 25 mM ammonium bicarbonate, incubated for 18 h at 37 • C. Peptides were extracted with 50% (v/v) acetonitrile and 5% formic acid (J.T. Baker ® ), and the sample placed in a new tube was dried in a vacuum concentrator for 2 h. Peptides were resuspended in 10 µL of 0.1% (v/v) formic acid, and sample salts were removed with a protein concentrator (Bond Elut OMIX, C18, No. A5700310K) and eluted in 25 µL of mobile phase (3% acetonitrile, 97% analytically pure water and 0.1% formic acid).
Peptide analysis was performed on an SYNAPT G2-S integrated nano-LC-ESI-/MS system (Waters Corporation, UK) equipped with a NanoLockSpray ion source and coupled online to a nanoACQUITY Ultra Performance Liquid Chromatograph (UPLC; Waters Corporation). The binary solvent system was composed of 2% acetonitrile with 0.1% formic acid (mobile phase A) and 98% acetonitrile with 0.1% formic acid (mobile phase B). The sample was injected onto a nanoACQUITY UPLC 2G Trap column, 5 µm, 180 µm × 20 mm (S/N 03073720616583) and subsequently washed with mobile phase A at a flow rate of 0.4 µL/min. Peptides were separated on a nanoACQUITY UPLC HSS T3 C18 column, 1.8 µm, 75 µm × 150 mm (S/N 01973716416209) using a linear gradient to 40% mobile phase B with a flow rate of 0.4 µL/min.
Data were acquired and processed using Protein Lynxs Global Selver 2.5.1 TM software, Waters (PLGS). For the identification of the proteins, the databases extracted from the website were used: https://www.uniprot.org, accessed on 5 January 2023, for Ligilactobacillus.

Statistical Analysis
Three independent replicates of each assay were performed in duplicate for each one. A variance analysis ANOVA was applied to the obtained data, expressed as means and standard deviation. The Tukey test was conducted for the comparison of the media with a significance level (p < 0.05) with the Minitab ® statistical software.

LAB Isolation and Microbiological Tests
Eighteen LAB strains were isolated under anaerobiosis conditions from the intestine of just-weaned CD1 mice per intestine section with a distribution of 28% duodenum (5 isolations), 22% jejunum (4 isolations), and 50% colon (9 isolations). In biochemical tests, the 18 isolated strains were negative for catalase and oxidase. In oxidation-fermentation tests, all strains fermented glucose, while only E4 and H2 strains fermented sorbitol. Moreover, A2, C1, and D1 strains were negative to xylose; E1, F3, and H1 strains were negative to lactose, and D1, E1, F2, and H1 strains were negative to saccharose.

Identity Percentage Matrix
The identity matrix allowed the authors to establish the similarity of the isolated BAL sequences. As a result, 100% similarity was obtained between L. murinus C1 and E1 and between L. murinus E2 and D1. The similarity for the rest of the L. murinus strains was between 92% and 99% ( Table 2).

Phylogenetic Tree
The phylogenetic tree built by the Maximum-Likelihood based on the Jukes-Cantor model shows that strains with similar sequences are clustered in the same clade level [35]. In particular, the isolated L. murinus E3 and I1 strains in this study were clustered in the same clade as the L. murinus DSM 100193, TCD6, NM28, and V10 (isolated from mice), L. animalis P38 and FR12 (isolated from chicken) and L. animalis JCM8692 (isolated from pork) strains found in GenBank. Likewise, the remaining L. murinus and L. animalis H2 strains isolated in this study were located in different clades (Figure 1).

Bacterial Growth Parameters
The growth parameters for every strain are listed in Table 3. L. m salivarius E4 displayed the highest µ values with 1.25 h −1 and 1.13 h −1, re other hand, L. murinus D1, E5, H1, and H3 were clustered with intervals 0.51 h −1 . Regarding τ, L. murinus H3 had the highest value at 2.12 h, and were observed for L. murinus A1, A2, E3, and S. salivarius E4 with inter and 0.84 h.

Bacterial Growth Parameters
The growth parameters for every strain are listed in Table 3. L. murinus E3 and S. salivarius E4 displayed the highest µ values with 1.25 h −1 and 1.13 h −1, respectively; on the other hand, L. murinus D1, E5, H1, and H3 were clustered with intervals between 0.39 and 0.51 h −1 . Regarding τ, L. murinus H3 had the highest value at 2.12 h, and the lowest values were observed for L. murinus A1, A2, E3, and S. salivarius E4 with intervals between 0.59 and 0.84 h. Table 3. Kinetic parameters of the isolated LAB strains. Data are expressed as a mean and standard deviation. abcdef shows a significant difference per column (Tukey).

Agar Diffusion Tests
The inhibitory activity against L. monocytogenes was observed in the 18 isolated strains, highlighting the activity of L. animalis H2 with inhibition diameters of 17 mm (p < 0.5) and L. murinus A1, A2, B1, C1, E2, F1 and F3 with inhibition diameters of 15 mm. Additionally, the isolated LAB strains showed an inhibitory effect against E. coli REPEC, with larger inhibition halos for L. animalis H2 with 25 mm and L. murinus G1, A1, and F3 with 23 mm. During the tests performed with isolated LAB strains against C. jejuni, a significant inhibitory effect was observed, highlighting the activity of L. animalis H2 with inhibition diameters of 44 mm (p < 0.5) ( Table 4).

SDS-PAGE and Zymograms
As a result, in protein profiles, some differences were observed in the number and intensity of protein bands. Zymograms against the different testing microorganisms evidenced two bands with an apparent molecular mass of 108 and 25 kDa, corresponding to two lytic bands of 107 and 24 kDa. The 107 kDa lytic bands were identified in the 16 L. murinus strains against M. lysodeikticus ATCC4698, E. coli REPEC O103, and C. jejuni. Moreover, the lytic band was identified in renatured zymograms after two hours against C. jejuni ATCC 49943, and after renaturation, for six hours, the lytic activity was observed as diffuse. On the other hand, the 24 kDa lytic bands were detected with different intensities in the L. murinus strains against L. monocytogenes, E. coli REPEC O103, and C. jejuni. It is important to highlight that the 107 kDa lytic band found in L. murinus B1 supernatants displayed a higher density against E. coli REPEC, and 24 kDa lytic bands in L. murinus H3 and I1 supernatants against L. monocytogenes (Figure 2). L. animalis and S. salivarius strains did not show lytic activity against the testing microorganisms.

UPLC-MS Analysis
The UPLC-MS analysis confirmed the presence of several specific peptides from the protein bands with lytic activity (Figure 3). Through the detection of 20 (Table 5), the lytic band with a molecular mass of 107 kDa was identified as a P N-acetylmuramoyl-L-alanine amidase activity. This PGH has a molecular structu of a catalytic domain formed by two 1,4-beta-N-acetylmuramoylhydrolase reg one domain bound to the peptidoglycan formed by nine repeated sequences

UPLC-MS Analysis
The UPLC-MS analysis confirmed the presence of several specific peptides obtained from the protein bands with lytic activity (Figure 3). Through the detection of 20 peptides (Table 5), the lytic band with a molecular mass of 107 kDa was identified as a PGH with N-acetylmuramoyl-L-alanine amidase activity. This PGH has a molecular structure made of a catalytic domain formed by two 1,4-beta-N-acetylmuramoylhydrolase regions and one domain bound to the peptidoglycan formed by nine repeated sequences of amino acids called LysM (lysin-like motif) domains (UniProt accession A0A4S2EQ32_9LACO).   Eight peptides were identified in the 24 kDa lytic band via the UPLC-MS analysis. Such peptides showed similarities with 12% of coverage with a protein fragment containing the GA (protein G-related Albumin-binding) domain within the genomic sequence of L. salivarius with aminopeptidase molecular function (UniProt accession A0A1V9RCV3_9LACO).

Discussion
Several niches were explored to obtain LAB, including human and animal gastrointestinal tract [16,17]. In this study, LAB strains were isolated under anaerobiosis conditions from recently weaned CD1 mice intestines. The 18 catalase and oxidase-negative isolated strains showed differences in their fermentation profile. All the strains fermented the glucose (hexose) used by LAB as the main source of energy. Meanwhile, most strains fermented pentose (xylose) and disaccharides (lactose and saccharose). On the other hand, the strains identified as L. murinus did not ferment sorbitol (polyol sugar). These results match those reported for the L. murinus 313 reference strain in terms of glucose, lactose, and saccharose fermentation, which does not ferment xylose and sorbitol [36]. Based on the fermentation type, lactobacilli are classified into three groups: obligately homofermentative, fermenting via the glycolytic pathway; facultative heterofermentative, fermenting via the glycolytic/pentose phosphoketolase pathway and obligately heterofermentative, fermenting via the pentose phosphoketolase pathway. Therefore, obligately homofermentative lactobacilli cannot assimilate pentoses [37]. However, there are exceptions for some strains, such as the homolactic fermentation of a pentose, for instance, phosphoketolase is induced in the presence of xylose, and the metabolism may shift between these two pathways according to the xylose concentration in the medium [38].
Molecular identification allowed the detection of a high proportion of L. murinus and one L. animalis. Specifically, L. murinus strains have been isolated from different niches, such as rodent gastrointestinal tract [39], dog feces [40], fermented dough [41], and marine sediments [19], while L. animalis has been found in the gastrointestinal tract of several animal species [42]. Furthermore, the identity matrix and the phylogenetic analysis determined the similarity between the sequences using the V2 variable region of the 16S rRNA gene from the 16 L. murinus strains and confirmed its close relationship with L. animalis, which has been described by several authors, identifying these two species at the same phylogenetic level with similarity percentages of up to 99% between 16S rRNA gene regions [41,43].
The isolated LAB inhibitory activity against the testing pathogens was demonstrated by the agar diffusion method. L. animalis H2 and L. murinus A1, A2, B1, C1, E2, F1 and F3 strains displayed more inhibition (p < 0.5) against L. monocytogenes, while L. animalis H2 and L. murinus G1, A1 and F3 showed larger halos against E. coli REPEC. In these studies, the inhibitory activity of the 18 isolated strains was highlighted against C. jejuni. According to reports, several LABs have displayed antagonistic activity against Campylobacter strains. For instance, Pediococcus pentosaceus CWBI B73, Lactiplantibacillus pentosus CWBI B78, and Enterococcus faecium THT inhibited C. jejuni and C. coli species [44]. On the other hand, L. salivarius SMXD51, MMS122, and MMS151 showed inhibition capacity against C. jejuni and C. coli [45]. In these studies, besides evaluating the antagonistic activity of LAB through the agar-diffusion method, cell supernatants were assessed, confirming that the inhibition capacity of the isolated strains was related to the production of protein compounds.
In general, antibacterial compounds with a protein nature are produced by LAB in low quantities, which makes them difficult to detect and analyze. Therefore, in this study, SDS-PAGE electrophoretic and zymography techniques were used, which are highly sensitive methods for the detection of enzymatic activity [46], which showed the presence of two proteins with lytic activity existing in the supernatants of the 16 strains of L. murinus isolated.
The 107 kDa lytic band, found in copolymerized gels with M. lysodeikticus, E. coli REPEC O103, and C. jejuni, was identified via the UPLC-MS fingerprinting method as an extracellular PGH with N-acetylmuramoyl-L-alanine amidase activity, considered as a bacteriolytic enzyme with antibacterial activity. According to the literature, the PGH are enzymes that hydrolyze peptidoglycan (PG) links, known as the main component of the cellular wall for Gram-positive and Gram-negative bacteria. The PG is formed by chains of alternate N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc) units, bound by B-1,4 links and cross-linked by peptides formed out of L-and D-amino acids. Specifically, the bacteriolytic activity from N-acetylmuramoyl-L-alanine amidase originated due to the specific site of hydrolysis in the PG, which cleaves the amide link between MurNAc and the L-alanine of the peptide [47]. Furthermore, the efficiency of the N-acetylmuramoyl-L-alanine amidase activity identified can be attributed to the number of binding LysM modules, which recognize the GlcNAc scraps of the PG and contribute to their additive junction [48].
The 24 kDa band with lytic activity against L. monocytogenes, E. coli REPEC O103, and C. jejuni displayed 12% similarity with a GA domain-containing protein fragment, reported in the L. salivarius genome. Particularly, the GA domain identified in this lytic band via the IDQMLELTVDQKDNFNK peptide is a domain localized in a variety of proteins of the cellular surface in different bacteria [49].
Similar to this research, several studies have reported the lytic activity from PGH against Gram-positive and Gram-negative pathogens [28,50]. One of the most researched PGH is lysostaphin, a 25 kDa peptidase produced by Staphylococcus simulans, that disrupts the peptide link between the third and fourth glycine residues of the pentaglycine crosslink in the S. aureus PG, which was demonstrated to be efficient in infections caused by Staphylococcus spp. [51].
The findings of this study can be applied from the perspective of the sustainable health and well-being objective [7] and the multidisciplinary effort with the One Health framework to reduce health risks such as sanitary problems, zoonotic illnesses, and antimicrobial resistance because these issues represent a growing threat for human and animal wellbeing. In the search for alternative solutions, there is the use of new probiotic strains and their bioactive components to reduce the use of antimicrobials because they provide health benefits and do not generate risks of resistance mechanisms. In this context, our research is important because it explores a new niche in the pursuit of Lactobacillus strains, which displayed inhibitory capacity against test pathogen microorganisms, C. jejuni, which is considered the main bacterial cause of gastroenteritis in humans around the world and the zoonotic pathogen responsible for gastroenteritis in humans associated to chicken meat consumption [52]. E. coli EPEC was identified as one of the main causes of diarrhea in humans and animals [2] and acknowledged as the primary bacteria that causes diarrhea in children under five years old [53], and L. monocytogenes considered as one emergent Food-borne pathogen, provokes severe health issues that appear sporadically or as an epidemic [54]. Due to the inhibitory capacity of the 16 Ligilactobacillus murinus isolated strains presented in this study, as well as the identification of their antibacterial metabolites, a feasible alternative is suggested for the prevention and control of these pathogens in order to promote animal health and well-being and, subsequently, human health as well.

Conclusions
The molecular identification of 16 strains of Ligilactobacillus murinus isolated from mice intestinal tract, its wide spectrum antibacterial activity against clinical and zoonotic pathogens, and the presence of bacteriolytic enzymes, one of them being PGH are important finding that contributes to the antagonist potential of these strains. In addition, these enzymes have a potential application as alternative antibacterial agents and treatment of infections caused by antibiotic resistant bacteria because the molecular mechanisms of resistance to antibiotics are not related to the hydrolysis of the PG.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/microorganisms11040942/s1, Figure S1: Products obtained through PCR with Ligilactobacillus animalis primers for the 18 isolated LAB strains. Table S1: Molecular identification (BLAST-NCBI) of isolated LAB strains from mice intestinal tract.