ARTICLE

Level Optimization of Beet Powder and Caramel Color for Beef Color Simulation in Meat Analogs before and after Cooking

Seonmin Lee1https://orcid.org/0000-0002-5713-1795, Kyung Jo1https://orcid.org/0000-0002-3006-5396, Seul-Ki-Chan Jeong1https://orcid.org/0000-0002-2163-8340, Yun-Sang Choi2https://orcid.org/0000-0001-8060-6237, Samooel Jung1,*https://orcid.org/0000-0002-8116-188X
Author Information & Copyright
1Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea
2Research Group of Food Processing, Korea Food Research Institute, Wanju 55365, Korea
*Corresponding author : Samooel Jung, Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134, Korea, Tel: +82-42-821-5774, Fax: +82-42-825-9754, E-mail: samooel@cnu.ac.kr

© Korean Society for Food Science of Animal Resources. This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Apr 20, 2023 ; Revised: Jul 24, 2023 ; Accepted: Aug 10, 2023

Published Online: Sep 01, 2023

Abstract

In this study, concentration levels of beet powder (BP) and caramel color (CC) were optimized to simulate beef color in meat analogs before and after cooking. The central composite design of response surface methodology (RSM) was used to set the levels of BP and CC, and the CIE L*, CIE a*, and CIE b* were selected as the responses for RSM. After optimization, myoglobin-free beef patties were prepared with three optimized levels of BP and CC. When raw, all the patties had the same color as natural beef; however, CIE L*, CIE a*, and CIE b* were statistically different from those of beef after cooking (p<0.05). Moreover, the use of BP and CC induced “browning” after the cooking process, with no excessive yellow color. Therefore, based on the overall desirability in the color optimization using RSM, the combination of BP (1.32%) and CC (1.08%) with the highest overall desirability can be used to simulate the color change of beef in meat analogs.

Keywords: meat analog; meat color; beet powder; caramel color

Introduction

The rapid increase in the world population and the limited production of food sources has led to the livestock industry becoming an unsustainable practice (Kumar et al., 2021). As more consumers become aware of environmental and animal welfare issues associated with the food they use, global demand for healthy alternative options has increased (Hwang et al., 2021). This trend has triggered changes in the food culture to develop sustainable meat analogs that mimic the sensorial and nutritional attributes of animal meats. The two major groups of meat analogs are cultured and plant-based. Cultured meats are generally described as “artificial,” “laboratory-grown,” or “synthetic” substitutes for meat, while plant-based analogs are composed of plant-origin proteins (Sun et al., 2021).

Among the various properties of meat and meat products, color is often regarded as the most important attribute for the consumers’ intent to purchase. Fresh meats are represented by diverse colors depending on the species, distribution of muscle fibers, and exposure to oxygen. Meat color mainly depends on the types and amount of a sarcoplasmic protein called myoglobin (King and Whyte, 2006). Upon exposure to high levels of oxygen, myoglobin converts into oxymyoglobin, which is bright cherry red in color. However, after slaughter, myoglobin exists as deoxymyoglobin exhibiting a purple-red color (Tomasevic et al., 2021). During the cooking process, the color of meat changes to brown due to the denaturation of proteins, primarily myoglobin (Sakai et al., 2022). It is desired that meat analog products undergo a similar color change before and after cooking. Cultured and plant-based meat analogs generally have a pale and beige or yellow color, respectively, due to the absence of myoglobin. Efforts have been made to mimic the red color of meat by the addition of natural pigments such as beet juice and soy leghemoglobin during the culturing of the tissues (Fraeye et al., 2020). However, it is still necessary to search for the ingredients and conditions that can simulate the thermally induced color changes of meat in its analogs are still required.

In a recent study, Sakai et al. (2022) showed that the use of beet pigments—often used to simulate red color in meat analogs—is associated with the appearance of a yellow-like brown color after thermal treatment when used in plant-based analogs. Herein, we simulated the color of raw and cooked beef patty samples by combining beet powder (BP) and caramel color (CC). We hypothesized that CC could partially simultaneously offset the high brightness in raw samples and the excessive CIE b* in cooked samples that was caused by BP. We selected beef as reference material and utilized myoglobin-removed beef to produce patties incorporating BP and CC to investigate whether BP and CC can potentially replicate beef color even in the meat analogs based on animal tissues. This approach enabled us to assess the impact of BP and CC on the color characteristics of the patties in the absence of myoglobin. The novelty of this work is a multivariate statistical approach using response surface methodology (RSM) to optimize the levels of BP and CC to simulate the color of beef before and after cooking by considering the interactions between the two pigments.

Materials and Methods

Sample preparation

Fresh beef (semitendinosus) and food grade BP and CC were purchased from a local market (Daejeon, Korea). Total of five carcasses (five batches) were obtained in five different days for the preparation of the source of beef with removed myoglobin (BRM). The excessive fats and connective tissues were removed, and the muscle was pulverized with a food mixer (C4 vv, Sirman SpA, Padova, Italy). To remove myoglobin (Lee et al., 2021), the ground beef was then mixed with an ice-cold buffer containing 0.1 M KCl, 2 mM MgCl2, 1 mM ethylene glycol tetraacetic acid, and 10 mM potassium phosphate (pH 7.0) in a 1:3 (w/w) ratio. The mixture was homogenized at 13,000 rpm for 30 s and centrifuged at 5,000×g for 10 min to separate the pellet. This procedure was repeated three times. Then, 0.1 M NaCl was added to the pellet in a ratio of 1:3 (w/w), and the homogenization and centrifugation were conducted three more times using the same conditions as described above. The pellet was washed with distilled water and used as BRM. The color attributes of beef before and after the removal of myoglobin are represented in Table 1.

Table 1. The color attributes (CIE L*, CIE a*, and CIE b*) of raw and cooked beef before and after removing myoglobin
Color attribute Treatments Before cooking After cooking SEM
CIE L* (+) Myoglobin 34.00Bb 45.05Ba 1.127
(–) Myoglobin 69.42Ab 71.81Aa 0.217
SEM 1.039 0.488
CIE a* (+) Myoglobin 21.76Aa 6.69Ab 0.721
(–) Myoglobin 2.89Ba 0.15Bb 0.079
SEM 0.705 0.170
CIE b* (+) Myoglobin 18.33a 14.62Bb 0.900
(–) Myoglobin 17.20a 16.70Ab 0.104
SEM 0.753 0.503

(+) Myoglobin and (–) myoglobin indicate the beef before and after removing myoglobin, respectively.

A,B Different uppercase letters indicate significant differences between means (p<0.05).

a,b Different lowercase letters indicate significant differences between means (p<0.05).

Download Excel Table
Design of experiments

The central composite design of RSM was used to set the optimum levels and the interactions between the independent factors influencing responses. First, we designated two independent factors (BP and CC) to a response (CIE L*, CIE a*, and CIE b* before and after cooking). Then, a five-level two-factor central composite design was employed, including five replicates of the central points. The coded levels of the independent variables are represented in Supplementary Table S1.

Preparation of the beef patties

The BRM was used to produce beef patties with appropriate levels of BP and CC to replicate the color of beef (Supplementary Table S1). The BRM and the colors were thoroughly mixed using a food processor (FPM250, De’Longhi-Kenwood Appliances, Treviso, Italy) for 1 min. Then, the mixture was stuffed in a stainless can (5 cm in diameter) and vacuum-packed, followed by heating at 80°C to reach the core temperature of 75°C. Then, the patty was cooled at 25°C, and the drip was removed to analyze the color. Five patties were prepared in independent days for each level to provide replicates for the color measurements and statistical analysis. Figs. 1A and B show the appearance of the beef patties with five levels of BP and CC.

kosfa-43-5-889-g1
Fig. 1. Appearance of the patties before and after cooking. (A) and (B) represent the appearance of the patties before and after cooking, representatively; (C) and (D) represent the appearance of the patties that were made by the optimized composition according to the overall desirability before and after cooking, respectively; (E) shows the appearance of the oven-cooked beef patties that were made by the optimized composition.
Download Original Figure
Color measurements

The color attributes (CIE L*, CIE a*, and CIE b*) were measured using a spectrophotometer (CM-5, Konica Minolta, Osaka, Japan) at two different parts of the cut inner surface (illumination area of 30 mm, illuminant D65, and standard observer of 10°). Blooming proceeded for 30 min before the color measurement. The average value of the two observation values was used as results for the analysis using Spectra Magic Software (SpectramagicTM NX, Konica Minolta).

Data analysis

Optimization of the concentration levels of the BP and CC was performed by RSM using Minitab 19 Statistical Software (Minitab, State College, PA, USA). Polynomial regression models were used to predict the responses, as presented in Table 2. In the model, the regression coefficients of the intercept, linear, quadratic, the involved interaction, and the probability level of p<0.05 were employed using one-way analysis of variance (ANOVA) to statistically test their effects. The response surface and contour plots were obtained by the Minitab software.

Table 2. The fitted polynomial equations for the response surfaces of the color (CIE L*, CIE a*, and CIE b*) of the patties (before and after cooking) with the optimized composition of BP and CC
Equation r2 Adjusted r2 Lack-of-fit AAD (%)
Before cooking
L*value=56.6689.358X116.622X2+0.8057X12+3.74X22+1.347X1X2 0.989 0.982 0.060 5.11
a*value=19.195+10.156X113.368X21.3606X12+5.075X221.579X1X2 0.984 0.972 0.059 10.88
b*value=16.984.735X1+10.513X2+0.5755X122.802X221.07X1X2 0.973 0.953 0.089 10.46
After cooking
L*value=66.972.62X116.98X2+0.449X12+1.827X22+1.787X1X2 0.992 0.986 0.767 2.43
a*value=4.841+3.087X13.714X20.0935X12+0.092X221.567X1X2 0.990 0.983 0.645 10.47
b*value=18.001.487X1+6.32X2+0.151X121.884X220.638X1X2 0.819 0.690 0.092 3.96

BP, beet powder; CC, caramel color; AAD, absolute average deviation.

Download Excel Table

The models were validated using the coefficient of determination (r2), residual analysis, lack-of-fit test, and absolute average deviation (AAD). The first three analyses were conducted in Minitab software, and AAD was calculated using the following Eq. (1):

AAD = { [ i = 1 p ( | y i exp y i cal | / y i exp ) ] / p } × 100
(1)

where p, yiexp, and yical indicate the experiment number, experimental response, and calculated response, respectively (Yolmeh and Jafari, 2017).

The levels (%) of BP and CC were optimized by the overall desirability features by setting nominal the best characteristics as the color of beef (CIE L*, CIE a*, and CIE b* before and after cooking). T-tests were carried out to determine the significance of the differences between the color of beef and the optimized color of the patties.

Results and Discussion

Color of beef before and after cooking

Table 3 shows the changes in the beef color before and after cooking based on the measured CIE L*, CIE a*, and CIE b*. Moreover, the appearance of the beef before and after cooking is represented in Figs. 1C and D. After cooking, the CIE L* increased from 34.00 (before cooking) to 45.05 (p<0.05) while the CIE a* (21.76 before cooking, and 6.69 after cooking) and CIE b* (18.33 before cooking, and 14.62 after cooking) decreased (p<0.05). These results indicate that the cooking process increased the CIE L* and decreased the CIE a* and CIE b*.

Table 3. Color attributes (CIE L*, CIE a*, and CIE b*) of beef semitendinosus and the patty (before and after cooking) that were made by the optimized composition according to the overall desirability
Color attribute Treatments Before cooking After cooking SEM
CIE L* Beef 34.00b 45.05Ba 1.127
Solution 1 36.46b 48.97Aa 0.409
SEM 1.039 0.598
Beef 34.00b 45.05Ba 1.127
Solution 2 35.75b 50.01Aa 0.288
SEM 1.037 1.127
Beef 34.00b 45.05a 1.127
Solution 3 33.72b 46.04a 0.309
SEM 1.038 0.536
CIE a* Beef 21.76a 6.69Bb 0.721
Solution 1 20.84a 11.04Ab 0.330
SEM 0.706 0.361
Beef 21.76a 6.69Bb 0.721
Solution 2 23.68a 11.38Ab 0.350
SEM 0.705 0.381
Beef 21.76a 6.69Bb 0.721
Solution 3 19.62a 11.60Ab 0.194
SEM 0.702 0.254
CIE b* Beef 18.33a 14.62Bb 0.900
Solution 1 19.39b 25.18Aa 0.440
SEM 0.758 0.654
Beef 18.33a 14.62Bb 0.900
Solution 2 18.25b 25.62Aa 0.351
SEM 0.751 0.607
Beef 18.33a 14.62Bb 0.900
Solution 3 20.51b 26.64Aa 0.357
SEM 0.751 0.611

A,B Different uppercase letters indicate significant differences between means (p<0.05).

a,b Different lowercase letters indicate significant differences between means (p<0.05).

Download Excel Table

Myoglobin, a pigment responsible for the meat color, is an oxygen-binding protein (Tomasevic et al., 2021). When meat is exposed to oxygen after being cut, oxymyoglobin is formed under the presence of oxygen at normal atmospheric pressure, favoring the bright red color of fresh meat (King and Whyte, 2006). The surface of the inner cut of the beef patties was exposed to air for 30 min before the color measurement, resulting in the formation of oxymyoglobin as the predominant form of myoglobin in the analyzed parts of the raw beef. Conversely, the temperature of the beef patty’s core reached 75°C during the cooking process. At this temperature, most oxy- and deoxymyoglobin contents denatured and transitioned into ferrohemochrome and then ferrihemochrome with a characteristic brown color (Suman and Joseph, 2013), resulting in a decrease in the CIE a*. Notably, although we consider the thermally induced color changes as the main cause of “browning,” the denaturation of myoglobin and other meat proteins can also form a gray color, which can increase the CIE L* and decrease the CIE b*. Although myoglobin accounts for only 0.5% of the total weight of red meat (King and Whyte, 2006), this protein pigment apparently changes meat color before and after cooking. The corresponding changes in the color should also be simulated in the meat analog.

Central composite design results and validation of the fitted model

Herein, the concentration levels of BP and CC were optimized using central composite design. We set the highest levels of BP and CC to 4.5% and 1.7%, respectively, based on our pilot study. Color measurements were carried out by setting a goal for the levels of BP and CC to reach the levels corresponding to raw and cooked natural beef. Table 2 shows the fitted polynomial equations for the responses (CIE L*, CIE a*, and CIE b* before and after cooking). The fitted models were all quadratic by considering both the second-order and interaction equations of the dependent variables. All the models had p values less than 0.05, indicating that the models are significant. All terms in the models were also significant (p<0.05, data not shown).

In this study, r2, p value of the lack-of-fit test, residual analysis, and AAD were used to validate the fitted models (Table 2). High r2 and adjusted r2 values were observed in all models, indicating the high accuracy of the applied models (Yolmeh and Jafari, 2017). Moreover, the predicted r2 values were 0.96, 0.92, and 0.93 for CIE L*, CIE a*, and CIE b*, respectively, before cooking, and 0.98, 0.95, and 0.62 for CIE L*, CIE a*, and CIE b*, respectively, after cooking. The difference between the predicted and adjusted r2 was less than 0.2 in all cases, indicating the existence of a reasonable agreement. As Table 2 shows, all the p values in the lack-of-fit test were higher than 0.05, confirming that the polynomial equation can adequately describe the data.

In the residual analysis, normal probability plot and the residual errors versus run plot were obtained (Supplementary Figs. S1 and S2) to confirm whether the residuals follow the normal distribution or not. In the normal probability plot, if the residual points follow a linear trend, the normal distribution of data is confirmed. Moreover, in the residual versus run plot, the residuals are arranged based on the order of the experimental runs; therefore, all data points should be randomly scattered to follow the normal distribution (Deepanraj et al., 2021). Herein, all the normal probability and residual versus run plots for the equations indicated that our data is normally distributed (Supplementary Figs. S1 and S2).

Finally, we assessed AAD and statistical variability (Yolmeh and Jafari, 2017). All the calculated AAD values were less than 11%, implying that the models had high accuracy as Zirebwa et al. (2014) reported. Therefore, according to the tests for model validation, the fitted models in Table 2 could be successfully employed for the optimization process.

Figs. 2 and 3 show the response surface and contour plots for CIE L*, CIE a*, and CIE b* before and after cooking. The response surface plots showed ridge surfaces. Observations showed that the level of BP had a direct correlation with CIE a* before and after cooking (Figs. 2B, 2E, 3B, and 3E) and CIE b* after cooking (Figs. 2F and 3F). CC had a direct relationship with CIE a* before cooking (Figs. 2B and 3B) and CIE b* before and after cooking (Figs. 2C, 2F, 3C, and 3F). The CIE L* showed an inverse relationship with both the levels of BP and CC (Figs. 2A, 2D, 3A, and 3D).

kosfa-43-5-889-g2
Fig. 2. Surface plots of the responses (CIE L*, CIE a*, and CIE b* of the patties) to the independent variables (levels of beet powder and caramel color). (A), (B), and (C) represent the surface plots of the responses as CIE L*, CIE a*, and CIE b* of the patties before cooking, respectively; (D), (E), and (F) represent the surface plots of the responses as CIE L*, CIE a*, and CIE b* of the patty after cooking, respectively.
Download Original Figure
kosfa-43-5-889-g3
Fig. 3. Contour plots of the responses (CIE L*, CIE a*, and CIE b* of the patty) to the independent variables (levels of beet powder and caramel color). (A), (B), and (C) represent the contour plots of the responses as CIE L*, CIE a*, and CIE b* of the patties before cooking, respectively; (D), (E), and (F) represent the surface plots of the responses as CIE L*, CIE a*, and CIE b* of the patty after cooking, respectively.
Download Original Figure
Optimization of the levels of beet powder (BP) and caramel color (CC) in beef patties

In RSM, desirability analysis was used to perform the multi-objective optimization considering a total of six responses: before and after cooking CIE L*, CIE a*, and CIE b*. The independent variables BP and CC were employed to optimize the responses by setting the nominal best characteristics as the color of beef. We obtained three solutions with the highest overall desirability. The setting satisfying the highest desirability values close to 1 is considered to be the optimum condition. Table 4 shows three combinations of the BP and CC with the highest overall desirability. Although the overall desirability of 1 was set as the goal (Yolmeh and Jafari, 2017), the six responses that we optimized to simulate BP before and after cooking resulted in overall desirability values between 0.60 and 0.65 for the combined objectives.

Table 4. Optimum levels of the variables to achieve the desired color of response variables (CIE L*, CIE a*, and CIE b* of beef semitendinosus before and after cooking) with maximum overall desirability
Solutions BP (%) CC (%) Overall desirability
Solution 1 1.32 1.08 0.65
Solution 2 1.80 1.01 0.61
Solution 3 1.30 1.51 0.60

BP, beet powder; CC, caramel color.

Download Excel Table

Next, the validation of the optimum values was conducted to evaluate the performance properties of the input parameters. A set of experiments was conducted by setting the optimum levels of BP and CC using the same experimental steps to obtain the model, and the output color responses were measured (Table 3). Results showed that the AAD of the patties made by the optimized solutions was below 14%, compared to the predicted values (Supplementary Table S2), indicating that the fitted model successfully predicted the responses. There were no significant color differences between the patties made by the levels of three optimized solutions and real raw beef (p>0.05). This indicates that the optimized levels of BP and CC could well simulate the color of raw beef. However, after cooking, all the color parameters of the three patties had significant differences compared to real cooked beef (p<0.05) except for solution 3 where there was no significant difference between the CIE L* of the patty and that of the real beef (p>0.05). This result implies that the combination of BP and CC could not accurately simulate the color of cooked beef. This might be related to the relatively low desirability (Table 4) for which we had to use six responses during the optimization. Nevertheless, the optimized levels of the three solutions showed similar trends in the CIE L* (increasing) and the CIE a* (decreasing) after cooking. However, the CIE b* of the patties increased after cooking, while the CIE b* of the real beef decreased after cooking. Moreover, as Figs. 1C and D show, the patties that were made by the levels suggested by solutions 1, 2, and 3 exhibited a color appearance similar to real beef. The combination of BP and CC also successfully simulated the color change from red to brown upon cooking, judged by the reduction of the excessive yellow-like-brown color to some degree. A lower CIE a* in the patties compared to beef was probably due to the color offset caused by the increase in the CIE b*, simulating the thermally induced “browning” of beef. In Fig. 1E, the appearance of oven-cooked patties is represented. Solution 3 had the most color similarity to beef, and its CIE L* after cooking showed no significant difference from that of beef. Therefore, considering the results in both the instrumental color and the appearance of the patties, 1.30% of BP and 1.51% of CC can be the optimum condition to simulate beef color.

The use of BP or beet extract to induce the color change of meat analog is limited due to their excessive expression of CIE b*. In Fig. 1, the patties added with BP but no CC also exhibit a purple color before cooking and yellow- or pink-like brown after cooking in both before and after cooking. Hence, a combination of the BP and CC is used to overcome the limitation of beet pigments.

Conclusion

This study aimed at stimulating the thermally induced color changes that beef undergoes during cooking. Combinations of BP and CC were used to replicate the “browning” of beef. Five levels of BP and CC were obtained by the central composite design of RSM. The obtained models for the six responses (CIE L*, CIE a*, and CIE b* before and after cooking) had great accuracy in predicting the response when we monitored the coefficient of determination, p value in lack-of-fit test, residual analysis, and AAD. Using the desirability function, we obtained three solutions for optimized levels of BP and CC: 1.32% of BP and 1.08% of CC, 1.80% of BP and 1.01% CC, and 1.30% of BP and 1.51% of CC. The patties of which color was designed by combining the optimized levels of BP and CC were able to successfully simulate the color of raw beef. However, the cooked patties had different colors when compared with the cooked beef. Nevertheless, the combination of the BP and CC could well simulate the thermally induced increase in the CIE L* and the decrease in the CIE a* caused by cooking. Moreover, the appearances of the patties were also comparable with the beef before and after cooking. Therefore, the combination of 1.32% of BP and 1.08% CC, with the highest overall desirability in this study, can be used to simulate the changes in the beef color from red to brown during cooking.

Meat color depends on several factors, including the animal species, parts or cuts, distribution of the intermuscular fats, muscle fiber types, and quality characteristics of meat (pH and resulting water-holding capacity). Therefore, future studies are required to simulate the color of meat in the meat analogs depending on the target characteristics of the products.

Supplemenrary Materials

Supplementary Materials

kosfa-43-5-889-suppl1.pdf

Conflicts of Interest

The authors declare no potential conflicts of interest.

Acknowledgements

This work was supported by Korean Institute of Planning and Evaluation for Technology in Food, Agriculture, and Forestry (IPET) through High Value-added Food Technology Development Program, funded by Ministry of Agriculture, Food and Rural Affairs (MAFRA) (grant number: 321025-05).

Author Contributions

Conceptualization: Lee S, Jung S. Data curation: Lee S. Formal analysis: Lee S, Jo K, Jeong SKC, Choi YS. Writing - original draft: Lee S. Writing - review & editing: Lee S, Jo K, Jeong SKC, Choi YS, Jung S.

Ethics Approval

This article does not require IRB/IACUC approval because there are no human and animal participants.

References

1.

Deepanraj B, Senthilkumar N, Ranjitha J, Jayaraj S, Ong HC. 2021; Biogas from food waste through anaerobic digestion: Optimization with response surface methodology. Biomass Convers Biorefin. 11:227-239

2.

Fraeye I, Kratka M, Vandenburgh H, Thorrez L. 2020; Sensorial and nutritional aspects of cultured meat in comparison to traditional meat: Much to be inferred. Front Nutr. 7:35

3.

Hwang J, Lee S, Jo M, Cho W, Moon J. 2021; The effect of sustainability-related information on the sensory evaluation and purchase behavior towards salami products. Food Sci Anim Resour. 41:95-109

4.

King NJ, Whyte R. 2006; Does it look cooked? A review of factors that influence cooked meat color. J Food Sci. 711:R31-R40

5.

Kumar P, Sharma N, Sharma S, Mehta N, Verma AK, Chemmalar S, Sazili AQ. 2021; In-vitro meat: A promising solution for sustainability of meat sector. J Anim Sci Technol. 63:693-724

6.

Lee S, Jo K, Yong HI, Choi YS, Jung S. 2021; Comparison of the in vitro protein digestibility of Protaetia brevitarsis larvae and beef loin before and after defatting. Food Chem. 338:128073

7.

Sakai K, Sato Y, Okada M, Yamaguchi S. 2022; Synergistic effects of laccase and pectin on the color changes and functional properties of meat analogs containing beet red pigment. Sci Rep. 12:1168

8.

Suman SP, Joseph P. 2013; Myoglobin chemistry and meat color. Annu Rev Food Sci Technol. 4:79-99

9.

Sun C, Ge J, He J, Gan R, Fang Y. 2021; Processing, quality, safety, and acceptance of meat analogue products. Engineering. 7:674-678

10.

Tomasevic I, Djekic I, Font-i-Furnols M, Terjung N, Lorenzo JM. 2021; Recent advances in meat color research. Curr Opin Food Sci. 41:81-87

11.

Yolmeh M, Jafari SM. 2017; Applications of response surface methodology in the food industry processes. Food Bioprocess Technol. 10:413-433

12.

Zirebwa FS, Sammie B, Kapenzi A, Mugandani R, Madanzi T. 2014; An evaluation of the performance and subsequent calibration of two reference evapotranspiration estimation models for Gweru, Zimbabwe. Acad J Agric Res. 2:167-174.