Optimization of the process for the preparation of peanut meal enzyme digest by enzymatic hydrolysis method
One-way test of hydrolysis conditions by composite alkaline protease
As shown in Figure 1, the most suitable one-factor reaction conditions of peanut protein hydrolysis in the one-factor test were screened to be: reaction time 4 h, initial reaction pH 9.0, reaction temperature 55 ℃ and enzyme addition of 800 U·g-1, considering the time cost and economic cost, respectively.
Response surface results and analysis of variance (ANOVA)
The corresponding test results under different conditions are shown in Table 6.
Table 6 Analysis of Box-Behnken response surface test results
serial
number
|
A
|
B
|
C
|
D
|
DH (%)
|
1
|
-1
|
-1
|
0
|
0
|
18.17
|
2
|
1
|
-1
|
0
|
0
|
23.25
|
3
|
-1
|
1
|
0
|
0
|
13.34
|
4
|
1
|
1
|
0
|
0
|
14.86
|
5
|
0
|
0
|
-1
|
-1
|
17.83
|
6
|
0
|
0
|
1
|
-1
|
17.22
|
7
|
0
|
0
|
-1
|
1
|
22.41
|
8
|
0
|
0
|
1
|
1
|
22.31
|
continuation sheet 6
serial
number
|
A
|
B
|
C
|
D
|
DH (%)
|
|
9
|
-1
|
0
|
0
|
-1
|
14.73
|
|
10
|
1
|
0
|
0
|
-1
|
18.24
|
|
11
|
-1
|
0
|
0
|
1
|
22.90
|
|
14
|
0
|
1
|
-1
|
0
|
13.41
|
|
15
|
0
|
-1
|
1
|
0
|
20.62
|
|
16
|
0
|
1
|
1
|
0
|
12.88
|
|
17
|
-1
|
0
|
-1
|
0
|
18.46
|
|
18
|
1
|
0
|
-1
|
0
|
22.68
|
|
19
|
-1
|
0
|
1
|
0
|
17.95
|
|
20
|
-1
|
0
|
1
|
0
|
22.58
|
|
21
|
0
|
-1
|
0
|
-1
|
16.95
|
|
22
|
0
|
1
|
0
|
-1
|
11.46
|
|
23
|
0
|
-1
|
0
|
1
|
21.20
|
|
24
|
0
|
1
|
0
|
1
|
16.28
|
|
25
|
0
|
0
|
0
|
0
|
23.12
|
|
26
|
0
|
0
|
0
|
0
|
24.95
|
|
27
|
0
|
0
|
0
|
0
|
24.01
|
|
28
|
0
|
0
|
0
|
0
|
23.54
|
|
29
|
0
|
0
|
0
|
0
|
23.80
|
A response regression model was established based on the experimental results to obtain the quadratic multiple regression model equation for the hydrolysis degree of peanut protein:
DH = - 669.0911 + 23.9143A + 89.6385B + 7.5216C + 0.1051D - 0.8900AB + 0.0205AC - 3.2000AD + 6.5000BC + 7.1250BD + 1.2750CD - 1.6016A2 - 5.0153B2 - 0.0710C2 - 5.7633D2
The analysis of variance (ANOVA) was performed on the quadratic multiple regression equation and the results are shown in Table 7, the P-value of the total model < 0.0001 indicates that the quadratic multiple regression model reached significance and the P-value of the dislocation term = 0.2611 > 0.05 indicates that at the level of dislocation term α = 0.05 the effect of non-testing factors on the results of the test is not significant. The coefficient of determination of the model, R2 = 0.9748, and the modification coefficient, Rajd2 = 0.9496, indicate that the variations in the experimental data can be explained by the quadratic multiple regression model and that the significance of the model is good. The above results indicate that the regression model established in this experiment fits well, the regression is significant and reliable, the model is valid, and the model can be used to predict the experimental results of the peanut protein hydrolysis process by alkaline protease.
The F-value in the table is proportional to the effect of response value, the larger the F-value, the greater the effect on the response value, the F-value of A is 41.19, the F-value of B is 159.99, the F-value of C is 0.65, the F-value of D is 109.18, i.e., B > D > A > C, which indicates that the strength of the effect of each factor on the degree of hydrolysis of peanut proteins is as follow: the initial reaction pH > amount of enzyme addition > reaction time > reaction temperature.
Table 7 Box-Behnken test regression model ANOVA
source of variance
|
squares sum
|
degrees of freedom (physics)
|
mean square
|
F-value
|
P-value
|
significance
|
general model
|
433.98
|
14
|
31
|
38.65
|
< 0.0001
|
**
|
A-reaction time
|
33.03
|
1
|
33.03
|
41.19
|
< 0.0001
|
**
|
continuation sheet 7
source of variance
|
squares sum
|
degrees of freedom (physics)
|
mean square
|
F-value
|
P-value
|
significance
|
|
B-initial pH
|
128.31
|
1
|
128.31
|
159.99
|
< 0.0001
|
**
|
|
C-reaction temperature
|
0.53
|
1
|
0.53
|
0.65
|
0.432
|
—
|
|
D-enzyme addition
|
88.13
|
1
|
88.13
|
109.88
|
< 0.0001
|
**
|
|
residual
|
11.23
|
14
|
0.8
|
|
|
|
|
lost proposal
|
9.37
|
10
|
0.94
|
2.01
|
0.2611
|
—
|
|
pure error
|
1.86
|
4
|
0.47
|
|
|
|
|
aggregate
|
445.21
|
28
|
|
|
|
|
Note:* denotes significant difference(P < 0.05); ** denotes highly significant difference(P < 0.01); —denotes insignificant difference(P > 0.05).
Response surface factor interactions
Response surface Design Expert 2013 software was used to fit the results of 29 groups of experiments to obtain the response surface three-dimensional plot of the two factors affecting the degree of hydrolysis of peanut protein, in order to explore the effect of the interaction between the reaction time, the initial reaction pH, the reaction temperature and the enzyme addition amount, the results are shown in Figure 2. Combined with the comparison of the significance analysis in Table 7, it shows that the initial reaction pH has the greatest effect on the degree of hydrolysis of peanut protein, while the reaction temperature has the least effect.
Validation of the regression model
The resulting regression equation model was subjected to DH Wangda analysis, and the enzymatic hydrolysis process was optimized by response surface analysis, which showed that the theoretical hydrolysis degree of peanut protein was 25.596% under the conditions of reaction time of 4.519 h, initial reaction pH of 8.635, reaction temperature of 54.855 °C, and enzyme addition of 900.773 U·g-1. The reaction conditions were adjusted according to the actual situation, and the modified optimal enzymatic hydrolysis process conditions were as follows: reaction time of 4.5 h, initial reaction pH of 8.6, reaction temperature of 55.0°C, and enzyme addition of 900.0 U·g-1. Three replicated validation tests were carried out with the optimized process, and the hydrolysis degrees were obtained as 24.96%, 25.13%, and 24.97%, respectively. The average hydrolysis degree was 25.02%, and the standard deviation (SD) was calculated to be 0.00095, indicating that the regression model can better predict the actual hydrolysis situation and has better practical value.
Determination of Enzymatic Solution Product Composition
The peanut meal enzyme digest was prepared using a complex alkaline protease under the optimal reaction conditions obtained by optimizing the hydrolysis process in a regression model, and the resulting product was subjected to compositional determination, which was converted to the results shown in Table 8.
Table 8 Composition of Enzymatic Solution (g·kg-1)
N
|
P2O5
|
K2O
|
Ca
|
Mg
|
Organic matter
|
5.10
|
5.30
|
9.70
|
0.92
|
0.67
|
50.40
|
Effects of peanut meal enzymatic solution and fermentation solution replacing part of nitrogen fertilizer with on potato growth and soil environment
Potato yield and yield components
As shown in Figure 3, compared with MCK treatment, the total potato yield of group M treatments were significantly increased, of which M10 treatment had the best yield increase effect, which was significantly increased by 26.63% compared with MCK treatment; the treatments that played a role in yield increase effect in group F treatments were F5 and F10 treatments, which had a significant increase in yield by 11.38% and 4.58%, respectively. However, the total potato yield of M10 treatment was significantly higher than that of F5 treatment.
As shown in Figure 4 (a) and (b), under the test conditions, with the increase of the proportion of peanut meal enzyme solution replacing inorganic nitrogen fertilizer, the effect on the increase of big potato yield showed a trend of increasing and then decreasing, of which the M10 treatment had the best effect on the increase of big potato yield, which was significantly increased by 30.90% compared with the MCK treatment; the M20 treatment had the highest rate of big potato, which was significantly increased by 5.97% compared with the control.
The accumulation of macro-elements and dry matter in potato tubers
The results are shown in Figure 5 (a), (b), (c) and (d). The effects of the four treatments in group M on the accumulation of nitrogen, phosphorus, potassium, and dry matter in potato tubers showed a trend of increasing and then decreasing, which were better than the MCK treatments as a whole, with the best effect of the M10 treatment, which was significantly more than the MCK treatments by 19.04%, 22.47%, 29.32%, and 31.86%, respectively. Compared with MCK treatment, F10 treatment had the best effect on nitrogen and phosphorus, which increased by 9.8% and 17.51%, respectively, and F5 treatment had the best effect on dry matter accumulation, which effectively increased by 15.62%. Both of them significantly increased the accumulation of potassium with small differences.
Potato quality
As shown in Figure 6, compared with MCK treatment, the protein and vitamin C contents of potato tubers from M10, M15 and M20 treatments in group M treatments were significantly increased, and the vitamin C content of M20 treatment was 253.95 mg·kg-1, which was significantly increased by 20.54%; the soluble sugar content of tubers was higher than that of MCK treatment, and reached the highest in M15 treatment, which was significantly higher than MCK treatment by 24.10%. Among the treatments in group F, the protein content of F15 treatment was significantly increased by 26.40% compared with that of MCK treatment, while in terms of vitamin C content of potato, only F15 treatment was significantly higher than that of MCK treatment, and there was no significant difference among the other treatments.
Soil physical and chemical properties
As shown in Figure 7(a), (b), (c), (d) and (e), soil pH, EC, effective phosphorus, quick-acting potassium and organic matter content reached their highest values at M20 and F20 treatments, which were significantly higher than MCK treatments by 5.70% and 4.43%, 40.06% and 23.22%, 31.65% and 22.70%, 14.28% and 10.53%, 13.48% and 17.83%, respectively; as shown in Figure 7(f), soil alkaline dissolved nitrogen content reached its highest value at M15 and F20 treatments, which were 106.26 mg·kg-1 and 97.23 mg·kg-1, respectively.
Sequencing information of soil bacteria from different treatments
The microbial sequencing information of soil samples under different treatments is shown in Table 9, the sequence number (Seq_num) of each treatment was in the range of 62601 ~ 84437, the number of sequence bases (Base_num) was in the range of 25823612 ~ 34861824, and the average optimized sequence length (Mean_length) was in the range of 412.113672 bp ~ 413.942411 bp, the shortest sequence length (Min_length) is in the range of 200 bp ~ 201 bp, and the longest average sequence (Max_length) is in the range of 474 bp ~ 515 bp.
Table 9 Sequencing information of soil microorganisms in different treatments
treatment
|
Seq_ num
|
Base_ num
|
Mean_ length
|
Min_ length
|
Max_ length
|
MCK_1
|
62601
|
25823612
|
412.511174
|
200
|
475
|
continuation sheet 9
treatment
|
Seq_ num
|
Base_ num
|
Mean_ length
|
Min_ length
|
Max_ length
|
|
MCK_2
|
66732
|
27533267
|
412.594662
|
201
|
475
|
|
MCK_3
|
73979
|
30513706
|
412.464429
|
201
|
486
|
|
M5_1
|
68257
|
28166848
|
412.658746
|
201
|
486
|
|
M5_2
|
63703
|
26319866
|
413.165251
|
201
|
475
|
|
M5_3
|
66614
|
27492707
|
412.716651
|
201
|
475
|
|
M10_1
|
68667
|
28397532
|
413.554284
|
201
|
475
|
|
M10_2
|
68520
|
28302419
|
413.053400
|
201
|
474
|
|
M10_3
|
73183
|
30222052
|
412.965470
|
201
|
475
|
|
M15_1
|
63710
|
26292872
|
412.696154
|
201
|
475
|
|
M15_2
|
64670
|
26693092
|
412.758497
|
201
|
487
|
|
M15_3
|
67496
|
27939457
|
413.942411
|
201
|
481
|
|
M20_1
|
67487
|
27906646
|
413.511432
|
201
|
486
|
|
M20_2
|
71325
|
29489591
|
413.453782
|
201
|
475
|
|
M20_3
|
68223
|
28125699
|
412.261246
|
201
|
491
|
|
F5_1
|
77530
|
31951173
|
412.113672
|
200
|
475
|
|
F5_2
|
78863
|
32560251
|
412.871068
|
201
|
475
|
|
F5_3
|
71458
|
29480399
|
412.555613
|
201
|
482
|
|
F10_1
|
82144
|
33938081
|
413.153499
|
201
|
475
|
|
F10_2
|
81387
|
33562282
|
412.378906
|
200
|
502
|
|
F10_3
|
83191
|
34286274
|
412.139222
|
201
|
477
|
|
F15_1
|
79725
|
32927029
|
413.007576
|
201
|
515
|
|
F15_2
|
78632
|
32418199
|
412.277432
|
201
|
492
|
|
F15_3
|
73793
|
30487328
|
413.146613
|
201
|
475
|
|
F20_1
|
78788
|
32517041
|
412.715655
|
201
|
475
|
|
F20_2
|
81659
|
33682505
|
412.477559
|
201
|
475
|
|
F20_3
|
84437
|
34861824
|
412.873788
|
201
|
475
|
Venn diagrams of soil bacterial communities at the OTU level for different treatments
As shown in Figure 8, the nine treatments had a total of 2,272 bacterial OTUs, or 15.32% of the total, with the F10 treatment having the largest number of unique bacterial OTUs at 910, or 6.13% of the total, and the M10 treatment having 855 unique bacterial OTUs, or 5.76% of the total.
Alpha diversity index of soil bacterial communities in different treatments
As shown in Table 10, the Coverage index of all treatments exceeded 96%, indicating that the sequencing results of soil samples from all treatments were reliable and could truly reflect the characteristics of soil bacterial communities. In group M and F treatments, the Ace index, Chao index and Sobs index of soil bacterial community of M10 and F10 treatments were significantly higher than the remaining three treatments in the group, indicating that the species of soil bacteria were more abundant and the status of dominant species was more prominent under M10 and F10 treatments, which were not significantly different.
Table 10 Alpha diversity indices of soil bacterial communities
treatment
|
Coverage Index
|
Ace Index
|
Chao Index
|
Sobs Index
|
MCK
|
0.969±0.0 abc
|
4709.5±78.5 cd
|
4529.9±47.5 cde
|
3562.3±71.0 bc
|
M5
|
0.970±0.0 ab
|
4651.7±88.0 d
|
4491.7±71.0 de
|
3527.3±68.8 bc
|
M10
|
0.967±0.0 bcd
|
5083.9±120.8 ab
|
4872.4±126.1 ab
|
3898.0±56.6 a
|
M15
|
0.969± 0.0 abc
|
4737.7±126.9 cd
|
4578.1±144.2 bcde
|
3540.3±76.1bc
|
M20
|
0.971±0.0 a
|
4442.3±10.4 d
|
4290.5±30.3 e
|
3370.0±16.1 c
|
F5
|
0.969±0.0 abc
|
4770.1±92.3 bcd
|
4614.2±76.1 bcd
|
3570.0±83.2bc
|
F10
|
0.966±0.0 d
|
5175.1±83.2 a
|
4962.3±63.5 a
|
3813.0±138.6 ab
|
F15
|
0.967±0.0 cd
|
5002.0±160.6 abc
|
4833.9±154.2 abc
|
3708.7±125.7 ab
|
F20
|
0.969±0.0abc
|
4751.2±98.6 cd
|
4554.6±62.5 cde
|
3545.3±72.2 bc
|
Note: The data in the table are the mean ± standard error of 3 repetitions, and those with the same letter at the end of the data in the same column indicate that the difference is not significant (P > 0.05).
Structural composition of soil bacterial communities at different treatment genus levels and analysis of significance tests for differences between groups
As shown in Figure 9, there were 10 genera that ranked in the top 20 in terms of abundance and could be accurately categorized, and the abundance of these 10 genera varied to different degrees under the treatments of Groups M and F. Among them, four genera with abundance greater than 1% and could be accurately categorized with highly significant or significant differences among the treatments of MCK, Groups M and F were Sinomonas (P ≤ 0.01), Terrabacter (P ≤ 0.05), Sphingomonas (P ≤ 0.05) and Arthrobacter (P ≤ 0.05), respectively, Terrabacter (P ≤ 0.05), Sphingomonas (P ≤ 0.05), and Arthrobacter (P ≤ 0.05). As shown in Figure 10, the abundance of Sinomonas increased by 149.8% in the F20 treatment compared to the M20 treatment, and the abundance of Terrabacter, Sphingomonas, and Arthrobacter increased by 36.1%, 91.8%, and 173.9%, respectively, in the M5 treatment compared to that at F5.
As shown in Figure 11, based on the results of the sample clustering tree analysis, the MCK treatment was more similar to the M10 treatment, the M20 treatment to the M15 treatment, and the M20 treatment to the F20 treatment in terms of soil dominant genus composition, respectively.
Heat map analysis of correlation between soil bacterial communities and environmental factors at different treatment genus levels
Five genera could be accurately categorized in the figure as shown in Figure 12. Arthrobacter and Terrabacter showed significant or highly significant negative correlation (P<0.05 and P<0.01) with soil effective phosphorus; fast potassium was significantly negatively correlated with Arthrobacter (P<0.05) and significantly positively correlated with Bacillus; Arthrobacter, Terrabacter and Sphingomonas were significantly negatively correlated with organic matter, pH and EC, while organic matter was significantly positively correlated with Sinomonas.