RESEARCH ARTICLE

Impact of high hydrostatic pressure on the single nucleotide polymorphism of stress-related dnaK, hrcA, and ctsR in the Lactobacillus strains

Joanna Bucka-Kolendo1*, Barbara Sokołowska2, 3

1Department of Microbiology, Culture Collection of Industrial Microorganisms-Microbiological Resource Center, Prof. Wacław Dąbrowski Institute of Agricultural and Food Biotechnology—State Research Institute, Warsaw, Poland;

2Department of Microbiology, Prof. Wacław Dąbrowski Institute of Agricultural and Food Biotechnology—State Research Institute, Warsaw, Poland;

3Institute of High Pressure Physics, Polish Academy of Sciences, Warsaw, Poland

Abstract

Lactic acid bacteria (LAB) are widespread in environments and can either have a positive impact because their ability to survive in harsh conditions and influence the product (probiotic properties, change of structure-EPS [exopolysaccharides], etc.), or a negative impact, (so not needed) because of their spoilage ability (beer, juices). High hydrostatic pressure (HHP), one of the non-thermal preservation methods used in the food industry, can force the LAB to activate the adaptative mechanisms. Under pressurization, the changes in the bacteria cells can occur at the transcriptional or translational level. This study evaluated the HHP on the single nucleotide polymorphism (SNP) changes in three genes, dnaK, ctsR, and hrcA, related to the stress-response mechanism in LAB. The correlation between the DNA polymorphism and the gene expression under HHP stress was assessed. The applied pressure of 300 MPa resulted in a low ratio of nonsynonymous substitutions to the synonymous substitutions (0 to 1.12), and a lower number of mutations was observed for pressurized strains (from 6 in hrcA to 11 in dnaK) than in controlled (from 3 in ctsR to 92 in hrcA). In all pressurized strains, the expression of genes was observed, whereas, in control strains, the gene expression was detected in three out of five strains. Although there was a noticeable change in stress-related gene expression after HHP, there was no correlation with SNPs. At the same time, with a high frequency of synonymous changes in nucleotide and high diversity for hrcA and dnaK, a very low diversity was found in ctsR sequences. The LAB strains stress response mechanisms are much more complex. The study requires information on the general mechanism and changes in the membranes’ composition, proteome changes, and gene expression patterns. The mutations in genes related to stress can have important implications for the strains’ fitness effect and adaptive ability of LAB strains, especially considering their food industry implication where the HHP techniques are used.

Key words: high hydrostatic pressure, lactic acid bacteria, nonsynonymous mutation, single nucleotide polymorphism, stress response, synonymous mutation

*Corresponding Author: Joanna Bucka-Kolendo, Department of Microbiology, Culture Collection of Industrial Microorganisms-Microbiological Resource Center, Prof. Wacław Dąbrowski Institute of Agricultural and Food Biotechnology—State Research Institute, Warsaw 02-532, Poland. Email: joanna.bucka@ibprs.pl

Received: 7 June 2022; Accepted: 18 August 2022; Published: 13 September 2022

DOI: 10.15586/qas.v14i4.1132

© 2022 Codon Publications
This is an Open Access article distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International (CC BY-NC-SA 4.0). License (http://creativecommons.org/licenses/by-nc-sa/4.0/)

Introduction

Several studies describe the lactic acid bacteria’s (LAB) ability to survive and respond to various environmental stresses (Tsuda et al., 2019; Bucka-Kolendo and Sokołowska, 2017) as they acquire preserving, probiotic, or spoilage properties (Bangar et al., 2022; Zapaśnik et al., 2022; Han et al., 2015). The most frequent phenotype described in the literature is the adapted cell (Papadimitriou et al., 2016), while adaptation refers to the effort of the cell to resist and persist under stress. However, the stress response to environmental factors can differ between species and depend on the applied stress (Mahmmodi et al., 2021; Van de Guchte et al., 2002). Extensive studies revealed the mechanisms involved in heat shock (De Angelis et al., 2004), bile (Bron et al., 2006), oxidative (Serrano et al., 2007), pH, and ethanol (Parente et al., 2010), where bacteria, through activation of the mechanisms involved in the stress response, adapt to the new conditions (Bucka-Kolendo and Sokołowska, 2017). Although the activated mechanisms may partially overlap, they are not identical (Papadimitriou et al., 2016), which can cause opposed results in LAB species, at species or even sub-species level.

Furthermore, the combination of stresses can trigger a cross-protection response (Yang et al., 2021b; Bucka-Kolendo and Sokołowska, 2017; Papadimitriou et al., 2016; Van de Guchte et al., 2002). The molecular mechanisms underlying the adaptation potential and response are based on the coordinated gene expression that can affect processes in cells, like cell division, transport, membrane composition, and DNA metabolism. Therefore, gene transcription, expression levels, and mechanisms engaged in bacteria growth under diverse stress conditions are greatly valued. Lactobacillus is a significant group of widespread organisms in different environments.

High hydrostatic pressure (HHP) is not a common stress factor for the LAB, as they are not generally exposed (Bucka-Kolendo and Sokołowska, 2017). However, HHP is a popular nonthermal preservation technique used in the food and beverage industry that reduces the number of microorganisms (Yaman et al., 2020; Chen et al., 2016) while preserving organoleptic molecules and providing “fresh” food. There is still limited knowledge about how LAB responds to the HHP, and its molecular mechanisms are not fully understood. HHP can negatively affect all molecular mechanisms in bacteria where DNA is involved, such as replication, transcription, and recombination (Salvador-Castell et al., 2020). The interaction between DNA and proteins may be disturbed due to the changes in the electrostatic and hydrophobic interactions. Pressure can dissociate ribosomal subunits and disturb the cytoskeletal proteins, resulting in reversible morphological changes. Since some effects of different factors can be similar, it is assumed that the ability to react to HHP comes from the cross-protection system HHP due to the fluidity of a complex response mechanism. Membrane fluidity among the critical factors is responsible for bacteria’s survival and growth under high-pressure conditions (Molina-Hoppner et al., 2003), where membrane lipids are stabilized by HHP and increase the melting points of lipids and transition the lipid bilayer to the gel state. Bacteria adapted to the HHP can adjust the phospholipid composition of the membrane by increasing the quantity of unsaturated fatty acids.

In bacteria, 90% of the genome represents genes; the rest contains small intergenic regions occupied by regulatory sites (Rocha, 2018). Genes are usually organized into operons and have only a few or no introns, and the insertions and deletions of genetic material tend to be determinants of gene expression (Rocha, 2018; Price et al., 2006). Many phenotypic variations among species are assigned to single nucleotide polymorphisms (SNPs) (Bailey et al., 2021; Hunt et al., 2009). Within the population, single base changes occur with a frequency greater than 1%. SNPs can be either synonymous (SS) when they do not cause changes in the amino acids (AA) or non-synonymous (NSS) when the AA structure is altered (Bailey et al., 2021; Lebeuf-Taylor et al., 2019; Hunt et al., 2009). It has been demonstrated that mechanisms altering the proteins’ structure, function, and expression level by affecting mRNA splicing, stability, structure, and protein folding are now better understood (López-González et al., 2018). NSS is more frequent and has a more substantial effect than SS mutations (Rocha, 2018). Since the effect is usually negative, those mutations are progressively removed from the population by the natural (purifying) selection, causing the low ratio of nonsynonymous (Ka) to synonymous (Ks) substitutions (ω - Ka/Ks). At the same time, synonymous substitutions can also be influenced by purifying selection, especially in fast-growing bacteria (Rocha, 2018). Many experimental studies prove that SS mutations can have positive solid fitness effects and drive adaptive evolution (Bailey et al., 2021; Liu et al., 2019).

It is crucial to provide insight into the genomic dynamics and polymorphism that characterize the physiological state of cells after exposure to stress, such as HHP, and understand the relationship between gene functions and phenotypic characteristics. Five LAB strains (two Loigolactobacillus backii, two Lactiplantibacillus plan-tarum, and one Lacticaseibacillus rhamnosus) were used. Strains were previously identified and analyzed (BuckaKolendo et al., 2020, 2021). This work described the proteomic and transcriptomic changes in selected Lactobacillus treated with HHP to determine the general adaptive response resulting from pressurization. As in previous studies, the HHP affects the proteome of treated LAB strains causing differences in the mass spectra profiles analyzed with MALDI-TOF MS (Bucka-Kolendo et al., 2020). Pressurization of the LAB strains can lead to changes in the expression patterns of stress-related genes (Bucka Kolendo et al., 2021). Among many functional genes associated with adaptation in LAB (Bucka-Kolendo et al., 2017), the three genes (dnaK, hrcA, and ctsR) previously described in the literature as stress-related (BuckaKolendo et al., 2017), were selected. The genes, like dnaK, hrcA, and ctsR, involved in the stress response (BuckaKolendo et al., 2021) can have different expressions under the stress factor, resulting in changes in the phenotype.

The aim was to relate the effect of the HHP on strains’ fitness through the changes in those gene expressions using SNP analysis of the partially sequenced genes. The hypothesis was that differences due to HHP might contribute to the bacteria in their gene mutations associated with the stress response and be elucidated with the genes’ phylogenetic clusters. The phenotypic and phylogenetic characterization of the LAB strains allows rising new insights into the adaptive abilities under the HHP. The fitness changes of the cellular response to stress factors can give an understanding of the individual strains’ responses to different factors (Douillard et al., 2016). As LAB are essential probiotics, starter, commensal, and pathogenic microorganisms, the in-depth research about the physiology of LAB stress has a significant meaning.

Materials and Methods

Schematic overview of the experiment

The schematic workflow of the study is presented in Figure 1. The diagram shows the overall process for the observation of the changes in the nucleotide sequences in the stress-related genes dnaK, ctsR, and hrcA under the HHP 300 MPa/5’. LAB strains were screened for possible adjustments to the changing environment. The first step was the isolation of strains from spoiled food products and the identification of the bacteria. The second step consists of applying HHP, molecular analysis with stress-related gene sequencing and expression, and growth analysis. The final step include the statistical and analytical analysis of SNP, phylogenesis, and PPI (protein-protein interaction).

Figure 1. A schematic flow chart of the experimental design used to establish and monitor the SNP in stress-related genes dnaK, ctsR, and hrcA. SNP, single nucleotide polymorphism.

Lactic acid bacteria

Five LAB strains (KKP 3565 Loigolactobacillus backii, KKP 3566 Loigolactobacillus backii, KKP 3568 Lactiplantibacillus plantarum, KKP 3569 Lactiplantibacillus plantarum, and KKP 3570 Lacticaseibacillus rhamnosus) were isolated from the food products, beer, tomato juice, and bread, and the process was performed according to the ISO 15214:2000 as described by Bucka-Kolendo et al. (2021). In addition, bacteria were identified using genetic and proteomic methods, as defined by Bucka-Kolendo et al. (2020, 2021). The isolated strains were deposited in the Culture Collection of Industrial Microorganisms-Microbiological Resource Center (IAFB, Warsaw, Poland), supported by the European Horizon 2020 research and innovation programme under grant agreement No 871129-IS_MIRRI21 Project. Strains were given the collection numbers, and the 16S rDNA sequences of each strain were deposited in the GenBank NCBI database (Table 1).

Table 1. Isolated strains of lactic acid bacteria.

Strain GenBank accessions Origin Identification based on 16S rDNA New nomenclature
KKP 3565 OK2913330 Beer Lactobacillus backii Loigolactobacillus backii
KKP 3566 OK287375 Beer Lactobacillus backii Loigolactobacillus backii
KKP 3568 OK291331 Bread Lactobacillus plantarum Lactiplantibacillus plantarum
KKP 3569 OK297672 Tomato juice Lactobacillus plantarum Lactiplantibacillus plantarum
KKP 3570 OK297673 Tomato juice Lactobacillus rhamnosus Lacticaseibacillus rhamnosus

Bacteria growth conditions and count number analysis (PCA)

LAB were grown and incubated under anaerobic conditions at 30°C for 48–72 h on MRS agar (Lactobacillus Agar DeMan, Rogosa, and Sharpe, Merc KGaA, Darmstadt, Germany), as described by Bucka-Kolendo et al. (2021). For counting, plates containing less than 300 CFU/mL were selected according to ISO 4833-1:2013. The non-treated LAB were considered control strains.

Immediately after HHP processing, the viability of the LAB strains was evaluated by counting colony-forming units on MRS Agar. The difference between control and treated strains was counted based on the number of surviving bacteria. The analysis was performed in two replicates.

Application of HHP

Using U 4000/65 apparatus (Unipress, Warsaw, Poland), the stationary phase LAB were threatened with HHP, as mentioned previously by Bucka-Kolendo et al. (2020). Therefore, exposition to 300 MPa for 5 min was chosen based on the former analyses (data not shown). At the HHP higher than 300 MPa, a significant reduction of live cells was observed, and at lower HHP, there were no significant differences from the controls. Hence, the used parameters were selected to observe possible cell changes but did not cause the bacteria’s death.

The HHP chamber details were as follows: the 0.96 L working volume, 600 MPa of maximum working pressure, from −10°C to +80°C working temperature, and the pressure-transmitting fluid was (1:1, v/v) water-polypropylene glycol. The time needed to generate the 300 MPa pressure was 70–80 s, with a release time of 2–4 s. The pressurization times did not include the come-up and come-down times. The assays were performed under two independent processes, and unpressurized strains were used as a control.

DNA extraction

According to the instructions, DNA from the stationary phase of control and HHP-treated bacterial culture (1.5 mL) was extracted with the ExtractMe DNA Bacteria Kit (Blirt S.A.–DNA, Gdansk, Poland) as mentioned by Bucka-Kolendo et al. (2020). The quality and concentration of the obtained DNA were measured at the absorbance of 260 nm and the 260/280 nm ratio with a UV-Vis NanoDrop spectrophotometer (ThermoFisher Scientific, Waltham, MA, USA). The isolated DNA was stored at −20°C.

PCR amplification

PCR amplification of dnaK, ctsR, and hrcA sequences was performed using specific primers (Table 2). For dnaK, ctsR, and hrcA, a total PCR volume of 60 µL contained: 30 µL of Dream Taq PCR Master Mix (ThermoFisher Scientific, Waltham, MA, USA), 1 µl of each primer in the final concentration of 0.4 µM, and 10–20 ng of DNA. A peqSTAR 2X thermocycler (PeqLab, Germany) was used for the amplification run. PCR amplification was performed as described by Riccardi et al. (2012), with modifications. Reactions conditions were: initiating denaturation at 94°C for 5 min, followed by 35 cycles of denaturation at 94°C for 60 s, annealing at 58°C for 60 s (for dnaK and hrcA), annealing at 57°C for 60 s (for ctsR), and elongation at 72°C for 90 s, with final elongation at 72°C for 2 min.

Table 2. Primers designed to amplify the dnaK, ctsR, and hrcA genes’ sequences.

Primers Sequence Gene size (bp)
dnaK F: 5’- CGGTAGCGGTTCTTGAAGGT -3’
R: 5’- GCCTTTTCAACCGTGTCACC -3
295 bp
ctsR F: 5’- CGGACTCGGAGCATGTTGAA –‘3
R: 5’- GTATGAGGGCGTCCAACACA -‘3
204 bp
hrcA F: 5’- TCCGAGCGCTTCTATGTTGG -‘3
R: 5’- ACCCATCAGCCCAATCATCC -‘3
297 bp

Electrophoresis in 1.5% (w/v) gel agarose was performed to analyze the PCR product size, and BioImaging Systems 06-2d.1-G: BOX (Syngen, UK) was used to visualize the product. Sequencing of the PCR products was performed in the 96-capillary 3730xl DNA Analyzer (Applied Biosystems-Life Technologies), and the obtained sequences were analyzed in NCBI BLAST databases. The obtained sequences were used to evaluate the DNA polymorphism, and the results are shown in Tables 3–5, respectively.

Table 3. Polymorphism analysis on partial sequences of dnaK in LAB strains. Sites excluding gaps or missing data are shown in parentheses.

Strain Sites π S k H G+C% SS NSS Ω F
KKP 3565 239 0.05439 13 10 2 46.2 55.42 175.58 2.678 ∧ 1.8 fold
KKP 3566 271 (271) 0.03321 9 9 2 46.5 54.75 215.25 0.505 ∧ 1 fold
KKP 3568 291 (270) 0.02222 6 6 2 47 66.25 203.72 0 ∧ 1.3 fold
KKP 3569 290 (270) 0.02222 6 5 2 47.1 62.83 204.17 0 ∨ 1.8 fold
KKP 3570 290 (270) 0.02593 7 9 2 47.8 63.33 206.67 0.727 ∧ 0.8 fold
Total control strains 267 (260) 0.06077 32 15.8 5 46.7 49.87 199.13 0.956
Total HHP strains 294 (268) 0.01642 9 4.4 5 46.4 63.83 197.17 0.972

LAB, Lactic acid bacteria; HHP, high hydrostatic pressure; π, nucleotide diversity; S, number of polymorphism sites; k, the average number of nucleotide differences; H, number of Haplotypes; SS, synonymous sites; NSS, nonsynonymous sites; ω, Ka/Ks ratio of nonsynonymous substitutions to the synonymous substitutions; F, fitness effect (where ∧ overexpression, ∨ is underexpression).

Table 4. Polymorphism analysis on partial sequences of ctsR in LAB strains. Sites excluding gaps or missing data are shown in parentheses.

Strain Sites π S k H G+C% SS NSS ω F
KKP 3565 177 (167) 0 1 1 1 42.5 32.67 132.33 0 ∧ 1.27 fold
KKP 3566 179 (179) 0.01676 3 2 2 43.3 39.83 137.17 0.576 ∧ 0.8 fold
KKP 3568 176 0.02273 4 4 2 42.9 40.08 133.92 0.894 ∧ 0.73 fold
KKP 3569 175 0.01714 3 2 2 43.7 34.33 139.67 0 ∨ 1.16 fold
KKP 3570 172 0.01744 3 3 2 43.0 41.50 129.50 0.636 ∧ 0.74 fold
Total control strains 178 (171) 0.02515 8 4.3 5 42.4 32.33 132.67 0.971
Total HHP strains 178 (171) 0.00936 3 1.6 5 43.2 31.27 130.73 0

LAB, Lactic acid bacteria; HHP, high hydrostatic pressure; π, nucleotide diversity; S, number of polymorphism sites; k, the average number of nucleotide differences; H, number of Haplotypes; SS, synonymous sites; NSS, nonsynonymous sites; ω, Ka/Ks ratio of nonsynonymous substitutions to the synonymous substitutions; F, fitness effect (where ∧ overexpression, ∨ is underexpression).

Table 5. Polymorphism analysis on partial sequences of hrcA in LAB strains. Sites excluding gaps or missing data are shown in parentheses.

Strain Sites π S k H G+C% SS NSS ω F
KKP 3565 616 (281) 0.334 94 105 2 45.6 57.75 194.25 1.168 ∧ 1.65 fold
KKP 3566 606 (260) 0.33846 88 87 2 48.8 56.67 171.33 0.731 ∧ 0.77 fold
KKP 3568 299 (261) 0 0 0 1 49.0 62.17 192.83 0 ∨ 0.21 fold
KKP 3569 294 (270) 0.02963 8 9 2 48.5 65.33 204.67 2.587 ∨ 3.22 fold
KKP 3570 289 (268) 0.01119 3 3 2 48.7 65.33 201.67 0 ∧ 0.53 fold
Total control strains 607 (252) 0.2143 92 54 4 46 48.97 161.03 1.173
Total HHP strains 293 (265) 0.00868 6 2.3 5 49.1 63.67 197.33 0.1588

LAB, Lactic acid bacteria; HHP, high hydrostatic pressure; π, nucleotide diversity; S, number of polymorphism sites; k, the average number of nucleotide differences; H, number of Haplotypes; SS, synonymous sites; NSS, nonsynonymous sites; ω, Ka/Ks ratio of nonsynonymous substitutions to the synonymous substitutions; F, fitness effect (where ∧ overexpression, ∨ is underexpression).

Bacterial mRNA extraction and RT-qPCR analysis

The total mRNA extraction and relative gene expression analysis were performed according to Bucka-Kolendo et al. (2021). The custom TaqMan gene expression assays (ThermoFisher Scientific, TFS) were used for dnaK, ctsR, and hrcA. The relative expression levels of analyzed genes were standardized to an endogenous control 16S rRNA gene. Endogenous control’s stability was evaluated for control and pressurized LAB using the ΔCT algorithm. For relative expression ratios in target genes, the 2−ΔΔCT method was used. The analysis was performed as the mean of the two independent experiments.

Statistical data analysis

The sequences of each gene (dnaK, hrcA, ctsR) of each strain (control and pressurized) were trimmed, aligned, and analyzed. Multisequences’ alignment and phylogenetic analyses were performed using MEGA v. X (Kumar et al., 2016). Phylogenetical trees were created with the neighbor-joining evolution method based on sequences of 3 stress-related genes in control and pressurized strains. Evaluation of the number of polymorphic sites (S), nucleotide diversity (π) of the strains of variant obtained for control and HHP conditions, and the ratio of the nonsynonymous (Ks) to synonymous (Ka) mutations (ω) were calculated with the Dna SP. 5.1 (Rozas et al., 2017).

To evaluate known and predicted interactions between dnaK, ctsR, and hrcA proteins in Lactobacillus, the PPI network was created with the STRING database (string-db.org) (Szklarczyk et al., 2019). The PPI provided a new possibility for revealing molecular mechanisms.

Results

The DNA polymorphism analysis was performed to obtain preliminary insight into the effect of the HHP on the LAB and identify single mutations in stress-related dnaK, ctsR, and hrcA, (Tables 3–5, respectively). The gene sequences had a G+C% content of 46.5–47.8% (dnaK), 42.5–43.7% (ctsR), 45.6–49% (hrcA). The rate of Ka to the rate of Ks was calculated to determine the evolutionary pressure on protein-coding sequences. The Ka/Ks ratios (ω) for total control populations and total HHP treated vary between 0 and 1.12, with the number of mutations lower for pressurized strains (from 6 (hrcA) to 11 (dnaK)) than controlled (from 3 (ctsR) to 92 (hrcA)). ω ratios on the strains level were 0 to 2.678, where most were 0 or close to 1, suggesting little differences between synonymous and nonsynonymous mutations. DNA polymorphism analysis of bacteria under HHP suggests that strains are under neutral purifying selection pressure and acts on the genes in most cases. In all pressurized strains, the expression of genes was observed (Figures 2–4 respectively, for dnaK, ctsR, and hrcA), whereas, in control strains, the gene expression was detected in three of five strains. In most LAB, the SNPs do not significantly change gene expression or gene product function. Considering that all strains survived the pressurization (Bucka-Kolendo et al., 2021), it was assumed that although a low SNP level occurred in the genes, bacteria gained adaptive ability. However, defining the correlation between gene expression changes and the SS mutations’ fitness mutations was impossible. The expression levels for dnaK, ctsR, and hrcA under the 300 MPa/5’ did not correlate with the SNPs in those genes.

Figure 2. RT-qPCR analysis of dnaK in control and pressurized (300 Mpa/5’) LAB strains. Data are the mean of the two independent experiments. The standard deviations are indicated with vertical bars. * Statistically significant differences were estimated using student’s t-test (P < 0.05). LAB, Lactic acid bacteria.

Figure 3. RT-qPCR analysis of ctsR in control and pressurized (300 Mpa/5’) LAB strains. Data are the mean of the two independent experiments. The standard deviations are indicated with vertical bars. * Statistically significant differences were estimated using student’s t-test (P < 0.05). LAB, Lactic acid bacteria.

Figure 4. RT-qPCR analysis of hcrA in control and pressurized (300 Mpa/5’) LAB strains. Data are the mean of the two independent experiments. The standard deviations are indicated with vertical bars. * Statistically significant differences were estimated using student’s t-test (P < 0.05). LAB, Lactic acid bacteria.

PCA analysis

Our experiment showed that the pressurization in 300 MPa’5 significantly affected the cell’s survival (P < 0.05). For strain KKP 3570 Lacticaseibacillus rhamnosus, the decrease level was 2 log (CFU/mL). For other strains, the reduction was around 1 log (CFU/mL) (Figure 5).

Figure 5. Survival of LAB strains under HHP. Data are shown as the mean of the two independent experiments. Vertical bars indicate the standard deviations (a) control strain and (b) pressurized strain. Different letters over the bars are significantly different (P < 0.05). LAB, Lactic acid bacteria; HHP, high hydrostatic pressure.

Phylogenetic analysis of dnaK, ctsR, and hrcA

Maximum likelihood clustering of aligned sequences of each gene was performed to assess the changes in the relationship of genes across LAB for controlled and HHP treated strains. Alignment was performed using CLUSTAL W, and results for unrooted neighbor- joining phylogenetic trees are shown in Figures 6–8 for dnaK, ctsR, and hrcA, respectively. In addition, control and pressured strains were compared for each gene to express whether mutations were reflected in the phylogeny.

Figure 6. Phylogenetic trees of dnaK sequences in studied Lactobacillus. The evolutionary distances were computed with the neighbor-joining method. The trees are drawn to scale, where the branch lengths are evaluated in the number of substitutions per site. (A) control strains and (B) pressurized strains. The trees were obtained using the MEGAX software (Kumar et al., 2016).

Figure 7. Phylogenetic trees of ctsR sequences in studied Lactobacillus. The evolutionary distances were computed with the neighbor-joining method. The trees are drawn to scale, where the branch lengths are evaluated in the number of substitutions per site. (A) control strains and (B) pressurized strains. The trees were obtained using the MEGAX software (Kumar et al., 2016).

Figure 8. Phylogenetic trees of hrcA sequences in studied Lactobacillus. The evolutionary distances were computed with the neighbor-joining method. The trees are drawn to scale, where the branch lengths are evaluated in the number of substitutions per site. (A) control strains and (B) pressurized strains. The trees were obtained using the MEGAX software (Kumar et al., 2016).

All three genes were examined to display the discriminatory power and reproducibility in studied lactobacilli strains. For ctsR, the similarity was at 97% that mirrored in the trees’ topology of controlled and HHP-treated strains. For dnaK sequences obtained, tress had a high level of resemblance at 98% for controlled and 91% for pressurized strains. The strains grouped based on hrcA sequences showed the highest differentiation, wherein control strains L. backii KKP 3565 and KKP 3566 had 97% and were more distantly related (42%) to other studied strains. The hrcA sequences of strains after pressurization were highly related and had a 100% level of similarity. There was no correlation between the expression pattern and clustering of the strains based on the SNP sequences analysis, resulting from the small number of bacteria used or the complexity of the stress-response mechanisms in LAB. Phylogenetic trees of all three genes (dnaK, ctsR, hrcA) revealed a high degree of relatedness between LAB strains.

Protein-protein interactions and genes co-expression analysis

To forecast the network of PPI, the Search Tool for the Retrieval of Instances of Neighbouring Genes—STRING database was used (Szklarczyk et al., 2019). The PPI network is valuable for describing the molecular processes, and atypical PPI can be related to the many stress response changes. Therefore, genes with a high degree of connectedness were clustered into PPI based on the highest confidence of 0.9 and a maximum number of interactions ≤ 5. The study used a k-means clustering method with an average local grouping coefficient of 0.696. As a result, the network found three clusters containing eight nodes with an average node degree of 5.

Also, the co-expression analysis of the studied genes was performed with STRING (Figure 9). Results indicate the strength of data shown in experiments, describing the correlation of expression between two coding protein genes based on the patterns of RNA expression and protein coregulation. The PPI network analysis can give a network of molecular interactions formed between the protein products of the studied genes. As presented in Figure 9A, among studied genes, the highest degree of centrality had gene dnaK, indicating the importance of those proteins in the resulting PPI network. Figure 9B shows the functional relationship retrieved for genes belonging to the network, where the color indicates the similarity that correlates with the presence/absence of the interactions. The analysis suggests that the proteins can be biologically linked and interact more than expected (20 edges compared to 5 estimated).

Figure 9. Protein-protein interactions (PPI) networks analysis. (A) PPI network of dnaK, ctsR, hrcA protein cluster in Lactobacillus generated with STRING and (B) the similarities which target families occur patterns across genomes. The color intensity correlates with confidence in the expected functional interaction between genes and organisms.

Discussion

The present study investigated the polymorphism in three genes (dnaK, ctsR, hrcA) related to the HHP stress response in five LAB strains. The number of SNPs in our studies was low and did not display one specific outline of the tested conditions. Although SNPs and mutation rates of all three genes were comparable in all strains, it was observed that HHP trigged a strong response in LAB strains and induced stress-related genes expression (Figures 2–4), which can determine the HHP as a potent stress factor of those genes. Furthermore, our results showed the highest nucleotide diversity in the hrcA sequences (Table 5). However, we did not find significant diversity between control and pressurized LAB strains in dnaK, ctsR, and hrcA sequences (Tables 3–5, respectively).

Overall, the LAB strains adapt to the given pressure, as their decrease level was 2 log (CFU/mL) (Figure 5). Comparing the control with the pressurized strains, the impact of the HHP on the LAB strains based on the Ka/Ks ratio displays a low value. This suggests that NSS mutations have a small negative effect on fitness and can stay in the population long before being removed by natural selection. The diversity of the dnaK, ctsR, and hrcA sequences indicates small changes between studied LAB strains, which shows that strains were protected from functional mutations, or mutations that have occurred were almost neutral. These findings are in line with the reports of other authors (Bailey et al., 2014, 2021; Rocha, 2018). An observed high rate of nonsynonymous mutations suggests that adaptive selection favors different protein sequences depending on environmental changes. They can alert the function of the genes. Our studies indicate that the expression of the genes can be the defense mechanism of strains against HHP. However, among genes, dnaK, ctsR, and hrcA expression does not seem to play a significant role in stand bias and demonstrates genomic plasticity. As Bailey et al. (2014) noticed, the mutations indicate increased gene expression. However, in the less-performed codons, the molecular mechanisms are responsible for transcript regulation and are very important in LAB’s adaptation and evolutionary dynamics. The fitness effect of a SS mutation is not always due to the impact on the gene in which it occurs but rather through the changes in gene expression which has been studied in our previous work (Bucka-Kolendo et al., 2021).

Little is known regarding the fitness effect of these mutations and how they rely on stress-related factors like HHP. The changes in the environment, such as pressurization, may trigger the activation of adaptive mechanisms and thus maintain changes in the genome, proteome, transcriptome, or metabolome, leading to cell damage and death. Studies (Jeon et al., 2021) on SNPs related to cell wall synthesis in L. acidophilus under high temperatures suggested that changes in SNP can make the cell wall more rigid. Mutations that appear under natural selection increase the possibility of fixation and may often be adaptive. The cell morphology and division can be affected by HHP and impact the cytoskeleton and autolytic proteins (Molina-Hoppner et al., 2003). As Yang et al. (2021a) noticed, after the HHP treatment of 300–400 MPa, L. plantarum and L. curvatus acquire higher cell integrity, smoother cell surface, and uniformed protoplasm. These findings suggest that bacteria can increase their HHP resistance by modifying their structures.

The stress factor, like HHP, may impact the discriminatory power of the studied genes’ sequences. Thus, comparing the DNA polymorphism changes in stress-related genes with the phylogenetic trees assembled is required and can give a comprehensive understanding of the relationship between functional genes and phenotypic characteristics. The phylogenetic trees are highly valued tools that demonstrate the diversity of bacterial strains to develop more effective methods for their identification, prediction of gene functions, and underpin genetic research. The phylogenetic analysis of the unrooted trees demonstrates the relative relationship between strains and the impact of the HHP. Based on the maximum likelihood of aligned sequences (dnaK [Figure 6], ctsR [Figure 7], and hrcA [Figure 8]), the clustering did not reflect the phylogenetic relationships. It did not allow the separation into different clusters for studied LAB species. All strains had a high similarity. The PPI network (Figure 9) showed an insight into the molecular processes, and the functional relationship retrieved for the genes belonging to the network suggests that the proteins interactions can be biologically linked more than expected, and the whole process is more complicated. Although the studies showed significantly higher diversity of the hrcA sequences, it did not correspond with strain classification and their phylogenetical clustering. The highest diversity was obtained for the hrcA gene in control strains. Our studies confirmed the results of Guidone et al. (2015), on the alignment of the ctsR sequences, marking the gene’s taxonomic value for LAB classification. However, our results disagree with their results and that classification of the LAB species based on the hrcA gene sequence is a good taxonomic marker, as the hrcA sequences after pressurization had a 100% level of similarity. Contrary to the studies presented by Huang and Lee (2011), the dnaK gene is much more polymorphic than the housekeeping gene 16S rRNA in the LAB group and therefore has a discriminatory value in closely related species. Where Sharma et al. (2019), indicated in their research that, among other genes, dnaK showed no SNP while tracking in yogurt and probiotic powder, which confirmed the results of our work, where dnaK was characterized by low diversity.

Our current work, together with previous studies on proteome changes under HHP (Bucka-Kolendo et al., 2020), and changes in gene expression under HHP in LAB (Bucka-Kolendo et al., 2021), gives a more in-depth insight into the selected LAB strains stress response to HHP. Our findings are significant to the technological implications of LAB resistance to the food industry, host survival (in the case of probiotics), and bacteria stress behavior, where a complex regulatory network of genes regulates the bacteria’s response to HHP during food processing. Further studies on transcriptome and proteome are needed to confirm this hypothesis.

Additional studies on the larger group of Lactobacillus strains are needed to examine the impact of the HHP on the stress response mechanisms of those bacteria and how the HHP can contribute to the mutations in stress-related genes. As the mutations, particularly in the promoter region of the genes related to the stress, can have important implications for the fitness effect and the adaptive ability of the strains. The performance of the whole genome sequences (WGS) to determine the evolution of LAB genomes under HHP treatment would help to obtain knowledge of the bacteria with higher flexibility under stress conditions. With WGS methods, more insights can be gained, as it combines the determination of the strain’s similarities based on SNP and gene-by-gene approach. Thus, considering how the mutation impacts the whole strain, additional studies that affect other neighbor genes, especially those grouped in the operons and transcribed together (like dnaK and hrcA), can help understand SNP mechanisms adaptative evolution better.

REFERENCES

Bailey, S.F., Angela, L., Morales, A. and Kassen, R., 2021. Effects of synonymous mutations beyond codon bias: the evidence for adaptive synonymous substitutions from microbial evolution experiments. Genome Biology and Evolution 13(9): evab141. 10.1093/gbe/evab141

Bailey, S.F., Hinz, A. and Kassen, R., 2014. Adaptive synonymous mutations in an experimentally evolved Pseudomonas fluorescens population. Nature Communications 10(5): 4076. 10.1038/ncomms5076

Bangar, S.P., Sharma, N., Bhardwaj, A. and Phimolsiripol, Y., 2022. Lactic acid bacteria: a bio-green preservative against mycotoxins for food safety and shelf-life extension. Quality Assurance and Safety of Crops & Foods 14(2): 13–31. 10.15586/qas.v14i2.1014

Bron, P.A., Molenaar, D., de Vos, W.M. and Kleerebezem, M., 2006. DNA micro-array-based identification of bile-responsive genes in Lactobacillus plantarum. Journal of Applied Microbiology 100: 728–738. 10.4161/bbug.2.2.13910

Bucka-Kolendo, J., Juszczuk-Kubiak, E. and Sokołowska, B., 2021. Effect of high hydrostatic pressure on stress-related dnaK, hrcA, and ctsR expression patterns in selected Lactobacilli Strains. Genes 12: 1720. 10.3390/genes12111720

Bucka-Kolendo, J. and Sokołowska, B., 2017. Lactic acid bacteria stress response to preservation processes in the beverage and juice industry. Acta Biochimica Polonica 64: 459–464. 10.18388/abp.2017_1496

Bucka-Kolendo, J., Sokołowska, B. and Winiarczyk, S., 2020. Influence of high hydrostatic pressure on the identification of Lactobacillus by MALDI-TOF MS-preliminary study. Microorganisms 8: 813. 10.3390/microorganisms8060813

Chen, Y.-S., Wu, H.-C., Yu, C.-R., Chen, Z.-Y., Lu, Y.-C. and Yanagida, F., 2016. Isolation and characterization of lactic acid bacteria from Xi-Gua-Mian (fermented watermelon), a traditional fermented food in Taiwan. Italian Journal of Food Science 28(1): 9–14. 10.14674/1120-1770/ijfs.v451

De Angelis, M., Di Cagno, R., Huet, C., Crecchio, C., Fox, P.F. and Gobbetti, M., 2004. Heat shock response in Lactobacillus plan-tarum. Applied and Environmental Microbiology 70: 1336–1346. 10.1128/AEM.70.3.1336-1346.2004

Douillard, F.P., Ribbera, A., Xiao, K., Ritari, J., Rasinkangas, P., Paulin, L., Palva, A., Hao, Y., de Vos, W.M., 2016. Polymorphisms, chromosomal rearrangements, and mutator phenotype development during experimental evolution of Lactobacillus rhamnosus GG. American Society for Microbiology. Applied and Environmental Microbiology 82(13): 3783–3792. 10.1128/AEM.00255-16

Guidone, A., Parente, E., Zotta, T., Guinane, C.M., Rea, M.C., Stanton, C., Ross, R.P., Ricciadi, A., 2015. Polymorphisms in stress response genes in Lactobacillus plantarum: implications for classification and heat stress response. Annals of Microbiology 65: 297–305. 10.1007/s13213-014-0862-7

Han, J., Chen, D., Li, S., Li, X., Zhou, W.-W., Zhang, B. and Jia, Y., 2015. Antibiotic susceptibility of potentially probiotic Lactobacillus strains. Italian Journal of Food Science 27: 282–289. 10.14674/1120-1770/ijfs.v270

Huang, C.H. and Lee, F.L., 2011. The dnaK gene as a molecular marker for the classification and discrimination of the Lactobacillus casei group. Antonie Van Leeuwenhoek 99(2): 319–327. 10.1007/s10482-010-9493-6

Hunt, R., Sauna, Z.E., Ambudkar, S.V., Gottesman, M.M. and Kimchi-Sarfaty, C., 2009. Silent (synonymous) SNPs: should we care about them? Methods in Molecular Biology 578: 23–39. 10.1007/978-1-60327-411-1_2

Jeon, S., Kim, H., Choi, Y., Cho, S., Seo, M. and Kim, H., 2021. Complete genome sequence of the newly developed Lactobacillus acidophilus strain with improved thermal adaptability. Frontiers in Microbiology 12: 697351. 10.3389/fmicb.2021.697351

Kumar, S., Stecher, G. and Tamura, K., 2016. MEGA7: molecular evolutionary genetics analysis version 7.0. Molecular Biology and Evoluation 33: 1870–1874. 10.1093/molbev/msw054

Lebeuf-Taylor, E., McCloskey, N., Bailey, S.F., Hinz, A. and Kassen, R., 2019. The distribution of fitness effects among synonymous mutations in a gene under directional selection. eLife 8: e45952. 10.7554/eLife.45952

Liu, W., Su, X., Duo, N., Yu, J., Song, Y., Sun, T., Zha, M., Menghe, B., Zhang, H. and Sun, Z., 2019. A survey of the relationship between functional genes and acetaldehyde production characteristics in Streptococcus thermophilus by multilocus sequence typing. Journal of Dairy Science 102(11): 9651–9662. 10.3168/jds.2018-16203

López-González, M.J., Escobedo, S., Rodríguez, A., Neves Rute, A., Janzen, T. and Martínez, B., 2018. Adaptive evolution of industrial Lactococcus lactis under cell envelope stress provides phenotypic diversity. Frontiers in Microbiology 9: 2654. 10.3389/fmicb.2018.02654

Mahmmodi, P., Khoshkoo, Z., Basti, A.A., Shotorbani, P.M. and Khanjari, A., 2021. Effect of Bunium persicum essential oil, NaCl, bile salt, and their combinations on the viability of Lactobacillus acidophilus in probiotic yogurt. Quality Assurance and Safety of Crops & Foods 13(1): 37–48. 10.15586/qas.v13i1.858

Molina-Hoppner, A., Sato, T., Kato, C., Ganzle, M.G. and Vogel, R.F., 2003. Effects of pressure on cell morphology and cell division of lactic acid bacteria. Extremophiles 7: 511–516. 10.1007/s00792-003-0349-0

Papadimitriou, K., Alegría, Á., Bron, P.A., De Angelis, M., Gobbetti, M., Kleerebezem, M., Lemos, J.A., Linares, D.M., Ross, P., Stanton, C., Turroni, F., Van Sinderen, D., Varmanen, P., Ventura, M., Zuniga, M., Tsakalidou, E., Kok, J., 2016. Stress physiology of lactic acid bacteria. American Society for Microbiology. Microbiology and Molecular Biology Reviews 80(3): 837–890. 10.1128/MMBR.00076-15

Parente, E., Ciocia, F., Ricciardi, A., Zotta, T., Felis, G.E. and Torriani, S., 2010. Diversity of stress tolerance in Lactobacillus plantarum, Lactobacillus pentosus and Lactobacillus paraplantarum: a multivariate screening study. International Journal of Food Microbiology 144: 270–279. 10.1016/j.ijfoodmicro.2010.10.005

Price, M.N., Arkin, A.P. and Alm, E.J., 2006. The life-cycle of operons. PLoS Genetics 2(6): e96. 10.1371/journal.pgen.0020096

Riccardi, A., Parente, E., Guidone, A., Ianniello, R.G., Zotta, T., Abu Sayem, S.M. and Varcamonti, M., 2012. Genotypic diversity of stress response in Lactobacillus plantarum, Lactobacillus paraplantarum and Lactobacillus pentosus. International Journal of Food Microbiology 157: 278–285. 10.1016/j.ijfoodmicro.2012.05.018

Rocha, E.P.C., 2018. Neutral theory, microbial practice: challenges in bacterial population genetics. Molecular Biology and Evolution 35(6): 1338–1347. 10.1093/molbev/msy078

Rozas, J., Ferrer-Mata, A., Sánchez-Del Barrio, J.C., Guirao-Rico, S., Librado, P., Ramos-Onsins, S.E. and Sánchez-Gracia, A., 2017. DnaSP 6: DNA sequence polymorphism analysis of large data-sets. Molecular Biology and Evolution 34: 3299–3302 Editor Sudhir Kumar, Oxford University Press. 10.1093/molbev/msx248

Salvador-Castell, M., Oger, P. and Peters, J., 2020. Chapter 8 – high-pressure adaptation of extremophiles and biotechnological applications. In: Physiological and biotechnological aspects of extremophiles. Academic Press, pp. 105–122. 10.1016/B978-0-12-818322-9.00008-3

Serrano, L.M., Molenaar, D., Wels, M., Teusink, B., Bron, P.A., de Vos, W.M. and Smid, E.J., 2007. Thioredoxin reductase is a key factor in the oxidative stress response of Lactobacillus plantarum WCFS1. Microbial Cell Factories 6: 29. 10.1186/1475-2859-6-29

Sharma, A., Kaur, J., Lee, S. and Park, Y.S., 2019. Tracking of intentionally inoculated lactic acid bacteria strains in yogurt and pro-biotic powder. Microorganisms 8(1): 5. 10.3390/microorganisms8010005

Szklarczyk, D., Gable, A.L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., Simonovic, M., Doncheva, N.T., Morris, J.H., Bork, P., Jensen, L.J., von Mering, C., 2019. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Research 47: D607–D613. 10.1093/nar/gky1131

Tsuda, H., Okuda, S., Haraguchi, T. and Kodama, K., 2019. Influence of exopolysaccharide on the growth of lactic acid bacteria. Italian Journal of Food Science 31: 233–242. 10.14674/IJFS-1317

Van de Guchte, M., Serror, P., Chervaux, C., Smokvina, T., Ehrlich, S.D. and Maguin, E., 2002. Stress responses in lactic acid bacteria. Antonie Van Leeuwenhoek 82: 187–216. 10.1007/978-94-017-2029-8_12

Yaman, H., Sarica, E. and Coskun, H., 2020. A comparative study on the effect of high hydrostatic pressure on the ripening of Turkish white cheese from different milk species. Italian Journal of Food Science 32: 589–595. 10.14674/IJFS-1712

Yang, D., Zhang, Y., Zhao, L., Wang, Y., Rao, L. and Liao, X., 2021a. Pressure-resistant acclimation of lactic acid bacteria from a natural fermentation product using high pressure. LWT–Food Science and Technology 69: 102660. 10.1016/j.ifset.2021.102660

Yang, H., He, M. and Wu, C., 2021b. Cross protection of lactic acid bacteria during environmental stresses: stress responses and underlying mechanisms. LWT–Food Science and Technology 144: 111203. 10.1016/j.lwt.2021.111203

Zapaśnik, A., Sokołowska, B. and Bryła, M., 2022. Role of lactic acid bacteria in food preservation and safety. Foods 11(9): 1283. 10.3390/foods11091283