2.1 Morphological and Phenotypic Characteristics
A total of 20 isolates were obtained from Nipponia nippon feces (Table 1), characterized by their typical small pinpointed, creamy white colonies, Gram-positive nature, catalase, hippuric acid, and gelatin liquefaction negativity, and both coccus and rod-shaped morphology.
Table 1
Morphological, physiological and biochemical characteristics of lactic acid bacterial strains isolated from the fecal samples of Nipponia nippon
Characteristic | 1 | 2 |
---|
isolates | D1、D2、D3、D4、D5、D6、D7、D8、D9、D11、D12、D13、D14、D15、D16、D17、D18、D19、D20 | E7 |
Morphology | coccus | rod |
Gram reaction | ་ | ་ |
catalase | - | - |
gelatin liquefaction | - | - |
Sucrose | 6/19 | ་ |
Xylose | 17/19 | ་ |
salicin | ་ | ་ |
glucose | ་ | ་ |
lactose | 9/19 | ་ |
Inulin | 16/19 | ་ |
raffinose | - | ་ |
cellobiose | ་ | ་ |
esculin | ་ | ་ |
maltose | 10/19 | ་ |
sorbitol | 5/19 | ་ |
mannitol | 7/19 | ་ |
hippuric acid | - | - |
+: positive or weakly positive reaction |
−: negative reaction |
Number: the number of positive reactions |
2.2 Acid tolerance |
The tolerance test results for twenty LAB isolates at different pH levels were shown in Table 2. At pH 3.0, most isolates demonstrated high tolerance. At pH 2.0, all isolates except strain D8 maintained viable bacterial counts above 104 CFU/mL. Under pH 1.0 conditions, viable bacterial counts were detected only in strains D1, D6, E7 and D20, exceeding 103 CFU/mL while no viable counts were detected in the remaining strains.
Table 2
Acid resistance capacity of LAB isolates under different acid conditions (log CFU/mL)
strains | Initial concentration | pH3 | pH 2 | pH 1 |
---|
D1 | 8.67 ± 0.04 | 6.61 ± 0.04 | 4.96 ± 0.24 | 3.31 ± 0.54 |
D2 | 8.74 ± 0.14 | 7.15 ± 0.01 | 4.66 ± 0.10 | - |
D3 | 8.57 ± 0.07 | 6.22 ± 0.11 | 4.49 ± 0.15 | - |
D4 | 8.38 ± 0.46 | 6.31 ± 0.26 | 4.08 ± 0.11 | - |
D5 | 8.56 ± 0.15 | 5.66 ± 0.04 | 4.55 ± 0.32 | - |
D6 | 8.63 ± 0.08 | 6.72 ± 0.04 | 5.92 ± 0.13 | 3.26 ± 0.36 |
E7 | 8.44 ± 0.10 | 6.32 ± 0.02 | 5.16 ± 0.40 | 3.07 ± 0.39 |
D8 | 8.73 ± 0.06 | 6.85 ± 0.03 | - | - |
D9 | 8.63 ± 0.10 | 6.86 ± 0.06 | 4.96 ± 0.25 | - |
D10 | 8.78 ± 0.13 | 5.67 ± 0.27 | 4.42 ± 0.12 | - |
D11 | 8.69 ± 0.04 | 6.21 ± 0.20 | 5.42 ± 0.23 | - |
D12 | 8.58 ± 0.25 | 5.70 ± 0.18 | 4.88 ± 0.09 | - |
D13 | 8.73 ± 0.16 | 6.76 ± 0.05 | 5.20 ± 0.44 | - |
D14 | 8.48 ± 0.06 | 6.76 ± 0.06 | 5.04 ± 0.35 | - |
D15 | 8.66 ± 0.08 | 7.29 ± 0.15 | 5.32 ± 0.37 | - |
D16 | 8.39 ± 0.23 | 6.26 ± 0.22 | 4.31 ± 0.06 | - |
D17 | 8.73 ± 0.18 | 6.58 ± 0.18 | 4.32 ± 0.26 | - |
D18 | 8.62 ± 0.34 | 6.56 ± 0.06 | 4.55 ± 0.24 | - |
D19 | 8.66 ± 0.31 | 5.52 ± 0.29 | 4.43 ± 0.11 | - |
D20 | 8.58 ± 0.03 | 6.91 ± 0.03 | 5.32 ± 0.37 | 3.02 ± 0.31 |
1-, No viable bacteria were detected |
2.3 Bile tolerance |
The tolerance results of twenty LAB isolates under different concentrations of bile salt were shown in Table 3. All strains showed high tolerance at 0.3% bile salt, with viable bacterial counts above 106 CFU/m L. At 1% bile salt, all strains except D3, D4, D12, D17, and D18 maintained counts above 104 CFU/mL. At 2% bile salt, only isolates D1, D2, D6, E7, D8, and D9 had viable counts above 102 CFU/mL, with no viable bacteria detected in other strains. Strain E7 had a particularly high viable count exceeding 105 CFU/mL. Among the twenty strains, only six demonstrated good tolerance to low pH and high bile salts, warranting further analysis.
Table 3
Survivability of the LAB isolates under different concentrations of bile salt (log CFU/mL)
strains | Initial concentration | 0.3% | 1% | 2% |
---|
D1 | 8.51 ± 0.06 | 6.18 ± 0.44 | 5.22 ± 0.09 | 3.56 ± 0.53 |
D2 | 8.76 ± 0.17 | 6.05 ± 0.32 | 5.08 ± 0.32 | 2.29 ± 0.78 |
D3 | 8.56 ± 0.13 | 6.70 ± 0.13 | - | - |
D4 | 8.32 ± 0.16 | 6.31 ± 0.09 | - | - |
D5 | 8.60 ± 0.10 | 6.64 ± 0.24 | 5.61 ± 0.26 | - |
D6 | 8.64 ± 0.11 | 6.92 ± 0.17 | 5.08 ± 0.43 | 2.47 ± 0.27 |
E7 | 8.79 ± 0.25 | 6.93 ± 0.14 | 5.31 ± 0.16 | 5.24 ± 0.22 |
D8 | 8.66 ± 0.11 | 6.43 ± 0.38 | 5.70 ± 0.06 | 4.04 ± 0.19 |
D9 | 8.86 ± 0.11 | 7.08 ± 0.17 | 5.78 ± 0.15 | 3.86 ± 0.05 |
D10 | 8.61 ± 0.26 | 6.57 ± 0.29 | 4.69 ± 0.01 | - |
D11 | 8.31 ± 0.24 | 6.31 ± 0.10 | 5.63 ± 0.28 | - |
D12 | 8.65 ± 0.08 | 6.41 ± 0.19 | - | - |
D13 | 8.64 ± 0.08 | 6.42 ± 0.38 | 4.36 ± 0.39 | - |
D14 | 8.23 ± 0.26 | 6.20 ± 0.46 | 4.79 ± 0.18 | - |
D15 | 8.36 ± 0.22 | 6.42 ± 0.16 | 4.59 ± 0.11 | - |
D16 | 8.37 ± 0.14 | 6.34 ± 0.32 | 5.66 ± 0.09 | - |
D17 | 8.78 ± 0.16 | 6.10 ± 0.14 | - | - |
D18 | 8.65 ± 0.11 | 6.51 ± 0.17 | - | - |
D19 | 8.57 ± 0.12 | 6.70 ± 0.29 | 5.42 ± 0.12 | - |
D20 | 8.86 ± 0.11 | 6.57 ± 0.29 | 4.69 ± 0.01 | - |
The mean of 3 values of each sample are presented with ± SD
2.4 Antimicrobial Activity
Six LAB isolates exhibited varying degrees of antibacterial activity against three pathogens (Escherichia coli, Staphylococcus aureus and Salmonella) as shown in Table 4. All isolates, except D8 (Table 4), showed significant inhibition zones (> 14mm) against the pathogens, with D1 exhibiting the strongest inhibitory effect on Salmonella (> 19mm).
Table 4
Antibacterial activities of LAB isolates
strains | Escherichia coli | Salmonella | Staphylococcus aureus |
---|
D1 | + | 14.31 ± 0.21 | ++ | 19.83 ± 3.10 | ++ | 16.59 ± 0.76 |
D2 | + | 14.43 ± 2.08 | + | 15.43 ± 3.68 | ++ | 18.13 ± 2.13 |
D6 | ++ | 16.17 ± 0.60 | ++ | 17.36 ± 3.92 | ++ | 16.58 ± 3.07 |
E7 | ++ | 16.53 ± 0.86 | ++ | 16.59 ± 1.57 | ++ | 16.65 ± 1.28 |
D8 | ++ | 16.35 ± 2.17 | + | 15.60 ± 2.01 | + | 12.04 ± 1.19 |
D9 | ++ | 16.56 ± 2.39 | ++ | 16.87 ± 2.72 | ++ | 16.51 ± 0.47 |
The different letters represent significant diversity and the same letters represent no significant diversity in one queue. The mean of 3 values of each sample are presented with ± SD
1±: 8 mm < < zone diameters ≤ 12 mm;+༚12 mm < < zone diameters ≤ 16 mm༛++༚16 mm < < zone diameters ≤ 20 mm༛+++20 mm < < zone diameters
2.5 Antibiotic Susceptibility
Selected isolates exhibited multi-drug resistance. All strains were sensitive to penicillin and chloramphenicol and resistant to vancomycin. Isolates D1 and D9 showed resistance to 8 and 6 different antibiotics, respectively, while strains D2 and D8 showed resistance to 4 different antibiotics. Strain E7 was resistant to cefradine, ciprofloxacin and vancomycin. The D6 strain was only resistant to Carbenicillin, Tetracycline, vancomycin. Although LAB are generally considered safe, the possibility of gene transfer in the gut of the crested ibis suggests selecting strains with less resistance. Strains D6 and E7 were safer than others.
Table 5
Susceptibility to antimicrobial agent of LAB isolates
| antibiotics |
---|
strains | PNC | CFP | CE | CB | CHL | TE | DC | ERY | CIP | GM | KAN | VAN |
ATCC25923 | S | S | S | S | S | S | S | S | S | S | S | S |
D1 | S | M | R | R | S | R | R | M | R | R | R | R |
D2 | S | M | R | S | S | M | S | M | R | S | R | R |
D6 | S | S | M | R | S | R | S | M | S | M | S | R |
E7 | S | S | R | S | S | M | S | M | R | S | M | R |
D8 | S | S | R | S | S | S | S | S | R | S | R | R |
D9 | S | S | R | S | S | S | R | S | R | R | R | R |
PNC Penicillin, CFP Cefoperazone, CE Cephradine, CB Carbenicillin, CHL Chloramphenicol, TE Tetracycline, DC Doxycycline, ERY Erythromycin, CIP Ciprofloxacin, GM Gentamicin, KAN Kanamycin, VAN vancomycin
1R-Resistance;M-medium sensitivity༛S-Sensitive
2.6 Cell surface hydrophobicity
The cell surface hydrophobicity of LAB, a measure of adhesion to intestinal epithelial cells, varied among the six isolates (Fig. 1). The E7 isolate exhibited significant moderate hydrophobicity (54.62%), while D2, D6, D8 isolates also showed moderate hydrophobicity (31.04%, 39.66%, 30.74%, respectively). The D1, D9 isolates exhibited low hydrophobicity (27.80%, 34.70%, respectively). The cell hydrophobicity of the strain E7 was the most significant among the six isolates (P < 0.05). Based on the above tests, strain E7 was selected for subsequent studies.
Values expressed as mean ± SD. Different letters represent significant difference, P < 0.05.
2.7 Hemolytic activity
On the blood agar plates, ETEC CVCC196 was γ hemolyzed (Fig. 2A). However, the strain E7 showed no hemolytic activity (Fig. 2B).
2.8 Molecular identification of LAB strains
Strain E7 was further identified by 16S rRNA sequencing and phylogenetic analysis, showing 99% sequence similarity to L. plantarum based on BLASTn. The sequence of the strain E7 was uploaded to the China Center for Type Culture Collection (CCTCC M 2024523), and the phylogenetic tree was constructed using MEGA7.0 (Fig. 3).
2.9 Growth performance of mice
During experimental period, all mice were healthy with no organ abnormalities. Analysis of body weight, final body weight, and ADG showed that L. plantarum E7 supplementation had no negative effect on growth performance (Fig. 4). No significant differences were found in organ indices (heart, liver, spleen, lungs, and kidneys) between the NC, BL, and BH groups (P > 0.05) (Table 6).
Table 6
The organ indices of each group mice (mg/g).
| Heart | Liver | Spleen | lungs | kidneys |
---|
NC group | 7.5 ± 0.16 | 45.2 ± 0.29 | 3.7 ± 0.13 | 7.2 ± 0.13 | 15.5 ± 0.16 |
BL group | 8.8 ± 0.20 | 44.1 ± 0.25 | 3.1 ± 0.12 | 7.2 ± 0.09 | 14.9 ± 0.2 |
BH group | 7.7 ± 0.16 | 46 ± 0.77 | 3.1 ± 0.14 | 7.4 ± 0.14 | 15.5 ± 0.13 |
2.10 Analysis of hematological parameters
Hematological parameters are shown in Table 7. Except for PLT, there were no significant differences between BL and BH groups compared to the NC group (P > 0.05). The PLT in BL group and BH group were significantly higher than the NC group (P<0.05).
Table 7
The hematological parameters of each group of mice.
Parameters | NC group | BL group | BH group | Parameters | NC group | BL group | BH group |
---|
WBC, 109/L | 4.18 ± 1.29 | 4.93 ± 1.00 | 5.05 ± 0.85 | HGB, g/L | 149.33 ± 15.16 | 150.83 ± 6.4 | 147.60 ± 8.14 |
Lymph, 109/L | 3.33 ± 1.26 | 4.08 ± 0.64 | 4.32 ± 0.96 | HCT, % | 49.21 ± 4.88 | 48.81 ± 2.15 | 48.18 ± 3.13 |
Mon, 109/L | 0.12 ± 0.04 | 0.15 ± 0.05 | 0.14 ± 0.1 | MCV, fL | 50.18 ± 1.27 | 48.75 ± 1.56 | 49.62 ± 0.59 |
Gran, 109/L | 0.73 ± 0.27 | 0.80 ± 0.20 | 1.08 ± 0.45 | MCH, pg | 15.15 ± 0.36 | 15.05 ± 0.50 | 15.14 ± 0.32 |
Lymph,% | 78.9 ± 4.62 | 81.65 ± 2.99 | 77.36 ± 3.37 | MCHC,g/L | 303.83 ± 4.70 | 305.17 ± 2.14 | 306 ± 4.36 |
Mon, % | 3.05 ± 0.80 | 3.02 ± 0.34 | 3.23 ± 0.77 | RDW, % | 14.42 ± 0.30 | 14.5 ± 0.81 | 14.9 ± 0.45 |
PDW, % | 16.25 ± 0.3 | 16.43 ± 0.2 | 16.58 ± 0.16 | PLT, 109/L | 1613 ± 126.18a | 1819.66 ± 149.37b | 1825.4 ± 144.63b |
RBC,1012/L | 9.81 ± 0.94 | 10.10 ± 0.27 | 9.72 ± 0.74 | MPV, fL | 5.6 ± 0.32 | 5.63 ± 0.34 | 5.88 ± 0.34 |
Values were shown as mean ± SD and a, b indicated a significant difference compared to NC group (a, bp < 0.05).
WBC, white blood cell counts; Lymph, lymphocytes; Mon, monocytes; PDW, platelet; RBC, red blood cell counts; HGB, hemoglobin; HCT, hematocrit; MCV, mean corpuscular volume; MCH, mean corpuscular hemoglobin; MCHC, mean corpuscular hemoglobin concentration; RDW, red blood cell distribution width; PLT, platelet counts; MPV, mean platelet volume;
2.11 Complete genome information of L. plantarum E7
The L. plantarum E7 genome is a circular chromosome with a total length of 3197533 bp and an average GC content of 46.50%. There were two circular plasmids, one with a total length of 56484 bp and an average GC content of 39.74%, and the other with a total length of 2503 bp and an average GC content of 55.97% (Table S1). The chromosome contains 3024 protein-coding sequences, accounting for 83.54% of the total length, as well as six 5srRNAs, five 16srRNAs, five 23 sRNAs, 72 tRNAs, 54 ncRNAs, and one CRISPRs. Plasmid1 contained 58 protein-coding sequences, accounting for 83.37% of the total length, six ncRNAs, and no 5S sRNAs, 16srRNAs, 23SrRNA, tRNA, ncRNA, or CRISPR. Plasmid 2 contained three protein-coding sequences, accounting for 55.97% of the total length, one ncRNA, and no 5SrRNA, 16SrRNA, 23SrRNA, tRNA, ncRNA, or CRISPR (Table S1), and its genome lacks CRISPRs. From inside to outside, the diagrammatic of the genome circle shows that the first circle represents scale, the second circle represents GC Skew, the third circle represents GC-content, the fourth and seventh circles represent COG to which each CDS belongs, the fifth and sixth circles represent the positions of CDS, tRNA, and rRNA on the genome (Fig. 5). The complete genome sequence was uploaded to NCBI (accession number CP158575.1).
The genome of L. plantarum E7 was annotated with 113 genes in the CAZy (Fig. 6), which showed that 32 genes, accounting for 0.0106 percent, were associated with glycosyltransferase (GTS). There was one gene related to polysaccharide lyase (PLs), accounting for 0.0003%. Genes associated with sugar esterase (CEs) were 13, accounting for 0.0043%, secondary activity (AAs) related genes were eight, accounting for 0.0026%, and carbohydrates related genes (CBMs) were 11, accounting for 0.0036%. There were 48 genes related to glycoside hydrolases (GHs), accounting for 0.0159%; however, the plasmid had no carbohydrate active enzymes. These findings suggest that the metabolic activity of L. plantarum E7 is mainly focused on the breakdown of carbohydrates to provide energy for the body.
In total, 2688 KEGG annotation genes were identified in the genome of L. plantarum E7 (Fig. 7A). These genes were assigned to 8 different signaling pathways. The distributing of these genes among various categories was as follows: 877 genes related to metabolism, 217 genes related to environmental information processing, 196 genes related to genetic information processing, 78 genes related to cellular processes, 78 genes related to human diseases, 991 genes related to Brite Hierarchies, 218 genes were related to Not Included in Pathway or Brite, and 33 genes were related to Organismal System. Among the metabolism-related genes, the main annotations were carbohydrate metabolism (270), amino acid metabolism (164), metabolism of cofactors and vitamins (87), energy metabolism (87), lipid metabolism (50), nucleotide metabolism (70), metabolism of other amino acids (38), biosynthesis of other secondary metabolites (34), xenobiotics biodegradation and metabolism (29), metabolism of terpenoids and polyketides (23) and glycine biosynthesis and metabolism (28).
In total, 29 KEGG annotation genes were identified in the Plasmid 1 genome of L. plantarum E7 (Fig. 7B). These genes were assigned to 4 different signaling pathways. The distribution of these genes among various categories was as follows: 12 genes related to metabolism, 1 gene related to environmental information processing, 11 genes related to Brite Hierarchies, 5 genes related to Not Included in Pathway or Brite. Among the metabolism-related genes, the main annotations included carbohydrate metabolism (4), amino acid metabolism (1), energy metabolism (3) and nucleotide metabolism (4). The plasmid 2 genomes had no functional annotation genes These findings suggest that L. plantarum E7 has strong abilities in carbohydrates and amino acid metabolism, which could be valuable for applications in various industries such as Food Engineering and biologics.
Based on gene annotation results, the genome of L. plantarum E7 contains numerous probiotic marker genes related to acid and bile tolerance, temperature tolerance, oxidative stress, riboflavin synthesis, exopolysaccharide secretion, and cell adhesion (Table S2). AntiSMASH analysis detected four secondary metabolite synthesis genes clusters namely RiPP-like, T3PKS, terpene and cyclic-lactone-autoinducer in the chromosome, with no similar gene clusters found in the plasmids (Table 8, Fig. 8). There were no potential probiotic genes in the plasmids.
Similar gene clusters were analyzed using MiBIG database and ClusterBlast algorithm (Fig. 8). The RiPP-like gene cluster of strain E7 and Lactiplantibacillus plantarum strain XJ25, Lactiplantibacillus plantarum strain SRCM10347 and Lactiplantibacillus plantarum strain LpYC41 have 100% homology. The T3PKS gene cluster has only 97% homology with Lactiplantibacillus plantarum strain RI-113. There are no gene clusters in the plasmids. These results indicate that L. plantarum E7 has the ability to produce new antibacterial substances.
Table 8
Secondary metabolite synthesis gene cluster of L. plantarum E7
Region | Type | From | to | Most similar known cluster | Similarity |
---|
Region1 | RiPP-like | 361955 | 374105 | 0 | 0 |
Region2 | T3PKS | 1814322 | 1855491 | 0 | 0 |
Region3 | terpene | 2874410 | 2895090 | 0 | 0 |
Region4 | cyclic-lactone-autoinducer | 3132211 | 3152813 | 0 | 0 |
Compared to the CARD database, 40 genes related to antibiotic resistance accounting for 0.628% were identified in the genome, with no resistance genes in the plasmids. There were 11 genes related to antibiotic targets, accounting for 0.364%. There was 1 gene related to antibiotic biosynthesis, accounting for 0.033%. There were 25 genes related to the total genes, accounting for 0.827%. There were no resistance genes in the plasmids (Table 9).
Table 9 Antimicrobial resistance genes detected in the genome of L. plantarum E7
Property
|
Number of Genes
|
Percentage (%)
|
Antibiotic Resistance
|
19
|
0.628
|
Antibiotic Target
|
11
|
0.364
|
Antibiotic Biosynthesis
|
1
|
0.033
|
Total genes
|
25
|
0.827
|
Comparing the chromosome genome sequence of L. plantarum E7 with the VFDB database (Table S3). Only a few virulence genes related to adherence, stress survival, exoenzyme activity, immune modulation, and regulation were detected in the genome data of L. plantarum E7, with none in the plasmids (Table S3).
2.11 Statistical Analysis
All data were expressed as means and standard deviations and analyzed using SPSS version 27.0. The difference was evaluated by one-way ANOVA and statistical significance was set at P < 0.05.