- 1Key Laboratory of Animal Nutrition and Feed Science of Yunnan Province, Yunnan Agricultural University, Kunming, China
- 2Faculty of Animal Science, Xichang University, Xichang, China
- 3Animal Husbandry Station of Yunnan Province, Kunming, China
- 4Animal Husbandry Technology Extension Station of Nujiang Lisu Autonomous Prefecture, Lushui, China
- 5Yunnan Tengsen Investment Group Limited, Tengchong, China
- 6Livestock Work Station of Tengchong City, Tengchong, China
Objective: The aim of this study was to evaluate the influence of different production systems, including low-altitude indoor feeding, high-altitude indoor feeding, and high-altitude grazing, on the slaughter performance and meat quality of gayals.
Methods: We slaughtered 15 male gayals (five from each of three feedlots) that were randomly selected, with similar body weights within each feedlot. The gayals were raised under different production systems: indoor feeding at 1,325 m (TC, Tengchong feedlot; n = 5), indoor feeding at 2,240 m (FHS, Fenghuangshan feedlot; n = 5), and grazing at an altitude of 2,600 m (JMD, Jiumudang alpine pasture; n = 5). Gayals were slaughtered to evaluate slaughter performance, and the left longissimus dorsi (LD) and biceps femoris (BF) muscles were collected for meat quality analyses (pH45min, , , , muscle fiber diameter, water loss rate, shear force), chemical composition (moisture, crude protein, ether extract, cholesterol), amino acid composition, and fatty acid profile.
Results: The results indicated that compared with the other two groups, the TC group exhibited significantly greater slaughter performance, including dressing percentage, lean meat percentage, and meat: bone ratio (p < 0.05). The TC group also showed higher levels of monounsaturated fatty acids (MUFA; especially C16:1 and C18:1 c9), essential amino acids (EAA), n-6:n-3 polyunsaturated fatty acids (PUFA) ratio, and thrombogenic index (TI), but lower cholesterol content (p < 0.05). In contrast, compared with the TC group, the JMD group had higher contents of PUFAs (including C18:3n3, C20:5n3, C22:5n3, and C22:6n3), a lower n-6:n-3 PUFA ratio, and a lower TI, with the TI decreased by 29.35% in LD muscle (p < 0.05). Additionally, the JMD group had higher contents of sweet amino acids in the LD muscle and umami amino acids in the BF muscle (p < 0.05). The slaughter performance and meat quality of the FHS group were intermediate between the TC and JMD groups, with closer alignment to the JMD group.
Conclusion: Low-altitude indoor-fed gayals demonstrated superior slaughter performance and higher beef yield, whereas high-altitude grazed gayals produced beef with a more favorable fatty acid profile and enhanced flavor-related amino acids, despite exhibiting lower productivity. High-altitude indoor-fed gayals exhibited slaughter performance and meat quality indicators that were intermediate between the two systems.
1 Introduction
Gayal (Bos frontalis) is a distinctive bovid species characterized by its unique semi-wild and semi-domesticated attributes (1). The gayal diverged from the gaur approximately 994,000 years ago, forming a separate subspecies or even species, rather than a domestic form of gaur, and is likely domesticated from an extinct subspecies or species (2). Therefore, the gayal holds a unique and critical branch position in the evolution of bovine species. The gayal mainly inhabits the hill forest areas of China, India, Bangladesh, Myanmar, and Bhutan, and plays a significant role in promoting forest conservation and maintaining ecological balance (3).
In China, gayals are traditionally grazed in the alpine valleys and subtropical rainforests of the Dulong and Nujiang River basins, located in the high-altitude regions of Yunnan Province (4). In 2014, the gayal was officially recognized and integrated into the Nujiang Lisu Autonomous Prefecture of China's Records of Specialty Livestock and Poultry Resources. Nutritional analysis of gayal meat reveals a high protein content and low fat percentage, complemented by the presence of beneficial amino acids and fatty acids (5). Moreover, gayal meat possesses favorable meat quality traits, including a smaller muscle fiber diameter and increased tenderness (6). Consumer demand for this high-quality beef is increasing steadily. However, the current population of gayals in China is < 4,000 individuals, resulting in gayal beef commanding a market price nearly twice that of other beef types. As a consequence, gayal beef production plays a significant role in the local agricultural economy.
Different production systems have direct effects on slaughter performance and meat quality (7). Previous studies have demonstrated that grain-fed bison (8) and grain-fed yaks (9) exhibited significantly superior slaughter performance, including higher carcass weight and dressing percentage, compared to their grass-fed counterparts, and, in terms of meat quality and nutritional value, grain-fed beef had lower moisture levels, higher intramuscular fat (IMF) content, and significantly improved tenderness, whereas grass-fed beef had a more favorable fatty acid profile. Specifically, grass-fed beef had lower levels of cholesterol-raising saturated fatty acids (C12:0–C16:0) and higher total n-3 polyunsaturated fatty acids (PUFA) content, as well as increased levels of long-chain n-3 PUFAs such as C20:5n3, C22:5n3, and C22:6n3, therefore, conferring greater health benefits and providing enhanced protection for the cardiovascular system (10). Furthermore, beef from different altitudes exhibited significant differences in meat quality after slaughter. As altitude increased, the a* value, pH, and shear force of beef significantly increased, whereas the b* value, L* value, and water-holding capacity significantly decreased (11). Additionally, previous reports indicated that yaks subjected to low-altitude indoor feeding exhibited higher dressing percentages, lean meat percentages, and overall edible meat quality compared with yaks grazing at high altitudes (12).
In China, due to the low efficiency of beef production under the high-altitude grazing system, some farmers have relocated gayals to other feedlots for indoor feeding, which are located at different altitudes. Nonetheless, to date, no systematic evaluations of their slaughter performance and meat quality under different production systems have been conducted. Therefore, the current study was conducted to evaluate the slaughter performance and meat quality of gayals under different production systems, including low-altitude indoor feeding, high-altitude indoor feeding, and high-altitude grazing conditions. The objective was to provide a comprehensive reference for the efficient and sustainable farming of gayals.
2 Materials and methods
2.1 Ethics statement
This research was approved by the Laboratory Animal Management and Use Committee of Yunnan Agricultural University (Approval No. YNAU 20180020).
2.2 Experimental animals
The study was conducted by randomly selecting five male gayals (~3 years old) for slaughter from each of the three different feedlots, totaling 15 male gayals. The three feedlots employed distinct production systems. In the low-altitude indoor feeding system (TC group), five healthy gayals were randomly selected based on similar body weight from the Tengchong feedlot in Qushi Town, Yunnan Province, China (altitude 1,325 m). These gayals had grazed at altitudes ranging from 2,500 m to 3,200 m until the age of 2 years before being transferred to the feedlot, with the fattening period lasting for 317 days (initial weight 261.1 ± 30.2 kg). For the high-altitude indoor feeding system (FHS group), five healthy gayals with matched body weights were randomly selected from the Fenghuangshan feedlot located in Lushui City, Yunnan Province, China (altitude 2,240 m). These gayals (initial weight 294.4 ± 66.4 kg) were subjected to a fattening period of only 47 days at the feedlot, having previously grazed at an altitude of 2,600 m before entering the fattening facility. Both the TC and FHS groups were provided with a corn silage-based diet. In the high-altitude grazing system (JMD group), five healthy gayals were selected from the Jiumudang alpine pasture in Gongshan County, Yunnan Province, China (altitude 2,600 m) to ensure similar initial body weights. These gayals were grazed year-round in alpine valleys where the plants primarily consist of alpine arrow bamboo, reeds, artemisia, and weeds. The ingredients and chemical compositions of the diets for the three experimental groups are shown in Table 1.
2.3 Slaughter performance measurements
All slaughter procedures followed the China National Operation Procedure of Cattle Slaughtering Standards (GB/T 19477-2018). Animals were fasted for 24 h and deprived of water for 8 h before slaughter. They were stunned after weighing and then exsanguinated, with the viscera, digestive system, head, skin, and hooves removed. The meat, bone, viscera and adipose tissue were also weighed in order to evaluate the slaughter performance and carcass characteristics. Dressing percentage was calculated as carcass weight divided by slaughter weight. Lean meat percentage was determined by dividing the lean meat weight by the carcass weight. Meat–bone ratio was calculated as a ratio of lean meat weight to bone weight. The backfat thickness of the left carcasses was measured over the deepest part of the loin-eye muscle (13), specifically, after the entire carcass was split into left and right halves, the backfat thickness was measured at the cross-sectional area between the 12th and 13th ribs on the left side of the carcass using a digital vernier caliper, with the measurement taken perpendicular to the muscle surface. Subsequently, the color of fat was evaluated using a color chart (1 = pure white, 8 = yellow), and the fat color was categorized into eight levels based on color depth (14), with a larger number indicating a yellower color. The cross-sectional area of the longissimus dorsi (LD) muscle between the 12th and 13th ribs was outlined on sulfuric acid drawing paper, and then the loin eye area was calculated using a digital planimeter (LK-300, Kezhe Biotech Co., Ltd., Shanghai, China) (15). This involved placing the sulfuric acid drawing paper on the LK-300 measuring platform, scanning the contour using an automatic sensor, and automatically calculating the enclosed area as the loin eye area.
After the animals were slaughtered, the carcasses were segmented according to the China National Beef Carcass and Cuts Standards (GB/T 27643-2011). Subsequently, samples of the LD muscle at the 12th and 13th ribs and the biceps femoris muscle (BF, muscle from the lateral hind legs located along the femoral margin of the semitendinosus muscle) were taken from the left half of each carcass and divided into two portions. One portion of the LD and BF muscles was used for on-site analysis of meat quality traits, including pH (45 min), L* (45 min), a* (45 min), b* (45 min), muscle fiber diameter, and water loss rate. The other portion was immediately sealed in sterile vacuum bags and stored at −20 °C for subsequent analysis of shear force, chemical composition (moisture, crude protein, ether extract, cholesterol), amino acid composition, and fatty acid profile.
2.4 Meat quality and chemical composition measurements
The pH value was measured by directly inserting an electrode 45 min after slaughter (1.5 cm deep into the muscle). The pH meter (BJ-260, Leici Co., Ltd., Shanghai, China) had been previously calibrated with pH 4.0 and pH 7.0 buffers. Muscle fiber diameter was determined following established methods from previous studies (16). The samples were placed in a constant temperature water bath and reaching a center temperature of 70 °C for 20 min, the meat was cooled to room temperature. Shear force was measured using a texture analyser (Model C-LM 3, Harbin Instruments Inc., China) on three cooked meat cores per sample (cylindrical, diameter 1.27 cm, length 3 cm) (17), with the mean value recorded after repeating the process three times. Meat colorimetric traits (L* value, a* value, and b* value) were measured at 45 min postmortem using a chromometer (CR-300, Minolta Co., Ltd., Osaka, Japan) at three randomly selected locations, and the average values were calculated (18). The chromameter was calibrated with a standardized white plate before each measurement. Water loss rate was measured using a texture analyser (Stable Micro Systems Inc., Godalming, UK). The middle part of the left LD and BF muscles were trimmed to 30 mm × 35 mm × 2 mm, weighed. The meat samples were placed between 20 pieces of filter paper and pressed at 35 kg for 300 s. Then, the samples were weighed again. The water loss rate was calculated in accordance with the methodology outlined by Ge et al. (19). The moisture, crude protein and ether extract of the muscle were determined using the AOAC method (20). According to the China National Standards for Determination of Cholesterol in Foods (GB 5009.128-2016), HPLC was used to determine cholesterol levels and the external standard method was used for quantification. The samples were homogenized and 1 g of samples was placed in a tube, mixed with 5 ml (60% KON: ethanol = 1:2) solution, and bathed in water at 70 °C until no precipitation occurred. After cooling, 10 ml (diethyl ether: petroleum ether = 1:1) solution was added, shaken, and centrifuged (9,000 rpm, 5 min). The waste liquid in the lower layer was discarded. The remaining solution was nitrogen-blown until dry. After drying, 5 ml anhydrous ethanol was added. The sample was filtered through a 0.45 μm membrane into a sample bottle for computerized testing.
2.5 Fatty acid measurements
The determination of fatty acid levels was optimized based on the method of the China National Determination of Fatty Acids in Foods Standards (GB 5009.168-2016). After homogenizing the muscle, 30 mg of sample was transferred to a 50 ml EP tube and mixed with 500 μl of extraction solution (isopropanol: hexane = 2:3). The mixture was vortexed for 30 s and ground for 4 min at 40 Hz using a grinding instrument. After that, it was ultrasound-treated in an ice water bath for 5 min, repeated thrice. The samples were centrifuged at 12,000 rpm and 4 °C for 15 min, and the supernatant was then transferred to another EP tube. The combined supernatant was concentrated by nitrogen blowing, dried completely, and then treated with a methanol: trimethylsilyl diazomethane (1:2) solution was maintained at room temperature for 30 min before being blow-dried again with nitrogen. Finally, hexane was added to redissolve the sample before centrifuging at high speed and filtering through a membrane (0.45 μm) into the loading bottle. Fatty acid analysis was conducted using a gas chromatography-mass spectrometry (GC-MS) instrument (Trace 1310 ISQ, Thermo, USA). The atherogenic index (AI), thrombogenic index (TI) and hypocholesterolemic: hypercholesterolemic ratio (h:H) were calculated according to the following equations (7, 21):
2.6 Amino acid measurements
Following the methods outlined in the China National Determination of Amino Acids in Foods Standards (GB 5009.124-2016), ion-exchange chromatography was used with post-column derivatization of ninhydrin to conduct the amino acid analysis. The samples were hydrolyzed with 10 ml of 6 mol/L HCl and three to four drops of phenol at 110 ± 1 °C for 22 h. The mixed amino acid standard solution and the sample solution were injected into the amino acid analyzer in equal volumes. The concentration of amino acid in the sample determination solution was calculated using the peak area of the external standard method.
2.7 Statistical analysis
All analyses were conducted using the General Linear Model (GLM) procedure in IBM SPSS Statistics (version 25.0). The effect of the three production systems on all variables was analyzed using a statistical linear model: Yij = μ + M? + eij, where Yij is the observed value for the j-th animal within the i-th production system (i = TC, FHS, JMD), μ is the overall mean, Mi is the fixed effect of the i-th production system, and eij is the residual error for the j-th animal within the i-th production system. The level of significance between the groups was determined according to Tukey's HSD post-hoc test using p < 0.05 as the limit to identify significant differences. In addition, principal component analysis was conducted using Origin (version 9.8.0.200) to identify the components that absorb greater variability and separate the three production systems. Only the qualitative parameters that exhibited significant differences among the three groups were included in the analysis; these parameters refer to indicators of slaughter performance and meat quality of the LD muscle, such as dressing percentage, pH value, and amino acid composition.
3 Results
3.1 Slaughter performance
Data on slaughter performance are presented in Table 2. The dressing percentage, lean meat percentage, meat: bone ratio, and backfat thickness in the TC group were significantly higher than those in the FHS and JMD groups (p < 0.05). The loin eye area was also significantly larger in the TC group than in the JMD group (p < 0.05). The color of fat was lower in the FHS group than in the JMD group (p < 0.05). With regards to viscera tissue, compared to the TC group, the FHS group showed a higher lung weight proportion and a lower spleen weight proportion relative to slaughter weight (p < 0.05). The proportion of kidney weight relative to slaughter weight in the JMD group was significantly higher than that in the TC group (p < 0.05).
Table 2. Comparison of the differences in slaughter performance of gayals under different production systems.
3.2 Meat quality and chemical composition
The meat quality and chemical composition data are shown in Table 3. At 45 min postmortem, the pH values in the FHS and JMD groups were significantly higher than those in the TC group (p < 0.05). The LD muscle fiber diameter was larger in the TC and JMD groups compared to the FHS group (p < 0.05), and the BF muscle fiber diameter increased in the TC group than in the FHS and JMD groups (p < 0.05). Regarding color parameters, the BF muscle from the FHS group exhibited lower b* values than those in the TC and JMD groups (p < 0.05). In terms of chemical composition, the crude protein content was higher in the LD muscle of the TC group compared to the FHS group (p < 0.05), whereas the highest crude protein content in the BF muscle was observed in the JMD group (p < 0.05). Additionally, it was observed that the ether extract content was highest in the TC group, but the cholesterol level was lowest among the three groups (p < 0.05).
Table 3. Comparison of the differences in meat quality and chemical composition of gayals under different production systems.
3.3 Fatty acid composition
The results regarding the differences in fatty acid composition among various production systems are presented as follows in Table 4. While the fatty acid profiles were generally similar between the two muscle types, some notable differences were observed. The LD muscle exhibited higher contents of C16:0 and C18:1 c9 across all three groups compared to the BF muscle. In contrast, the BF muscle showed a tendency toward higher levels of C18:2n6c, C20:4n6, total PUFA content, and the n-6:n-3 PUFA ratio and the h:H ratio.
Table 4. Comparison of the differences in fatty acid composition of gayals under different production systems (%).
Significant differences in fatty acid composition were evident among the three production systems. In terms of saturated fatty acids (SFA), the JMD group exhibited higher SFA content in both muscles, with significantly elevated levels of C15:0, C17:0, C20:0, C21:0, C23:0, and C24:0 compared to the TC and FHS groups (p < 0.05). Only the C14:0 content in the LD muscle was higher in the TC group than in the other groups (p < 0.05). Regarding monounsaturated fatty acids (MUFA), the TC group showed significantly higher contents of C16:1, C18:1 c9, C18:1 c11, and total MUFA in the LD muscle, as well as elevated C18:1 c9 and C18:1 c11 levels in the BF muscle (p < 0.05). However, the C18:1 t11 content in the BF muscle was lowest in the TC group (p < 0.05). For PUFA, the FHS and JMD groups demonstrated significantly higher contents of multiple PUFAs, including C18:3n3, C20:5n3, C22:5n3, and C22:6n3 in both muscles compared to the TC group (p < 0.05). Additionally, in the LD muscle, C18:2n6t, C20:2, C20:3n6, C20:4n6, and total PUFA content were higher in the FHS and JMD groups than in the TC group (p < 0.05). The n-6:n-3 PUFA ratio was significantly lower in both muscles of the FHS and JMD groups compared to the TC group (p < 0.05). Furthermore, the TI in LD muscle was lower in the JMD group than in the TC group (p < 0.05).
3.4 Amino acid composition
Table 5 demonstrates that the muscle amino acid composition was significantly influenced by different production systems.
Table 5. Comparison of the differences in amino acid composition of gayals under different production systems (%).
Regarding nonessential amino acids (NEAA), the total NEAA content in both the LD and BF muscles was higher in the FHS group than in the TC group (p < 0.05). In the LD muscle, the content of aspartic acid was significantly higher in the TC and FHS groups than in the JMD group (p < 0.05), whereas serine and glycine were lower in the TC group than in the FHS group (p < 0.05). Additionally, alanine and cysteine contents were lowest in the TC group (p < 0.05). In the BF muscle, proline content was reduced in the JMD group compared to the TC group (p < 0.05), and both aspartic acid and serine contents were notably lower in the TC group than in the other groups (p < 0.05).
For essential amino acids (EAA), the total content was higher in the TC group than in the FHS group in both muscles (p < 0.05). In the LD muscle, valine, methionine, and isoleucine contents were higher in the TC group than in the JMD group (p < 0.05), and the lysine content was highest in the JMD group (p < 0.05). In the BF muscle, valine, isoleucine, and histidine contents were notably elevated in the TC group compared to the FHS group (p < 0.05).
Additionally, the results show that the contents of sweet amino acids in the LD muscle and umami amino acids in the BF muscle were significantly higher in the JMD group than in the TC group (p < 0.05). Moreover, bitter amino acids contents were significantly higher in the TC group than in the FHS group (p < 0.05).
3.5 Principal component analysis
Through principal component analysis, Figure 1 shows the overall variation in slaughter performance, nutritional value, and flavor of the meat. PC1 and PC2 explained 67.0% of the variability. Based on the loading values, dressing percentage, EAA, PUFA, sweet amino acids and bitter amino acids were the primary factors for distinguishing the TC group from the FHS and JMD groups (Figure 1A). The TC group was clearly separated from the FHS and JMD groups, while there was an overlap between the FHS and JMD groups (Figure 1B).
Figure 1. Principal component analysis loading plot and score plot for the major variability variables on the first two PCs were absorbed. (A) Loading plot of the investigated variables; (B) PC score plot of the first two principal components. CP, crude protein; EE, ether extract; MUFA, monounsaturatedfatty acids; PUFA, polyunsaturated fatty acids; TI, thrombogenic index; EAA, essential amino acids; Bitter amino acids: valine, methionine, isoleucine, leucine, tyrosine, phenylalanine, histidine, argnine; Sweet amino acids: glycine, serine, threonine, alanine, proline, lysine.
4 Discussion
4.1 Slaughter performance
Cattle production performance is commonly assessed through carcass characteristics, which can be influenced by the feeding system (22). Previous studies have shown that the buffaloes raised on pasture exhibited a lower dressing percentage compared to those raised under intensive management (7). Yuan et al. (12) reported that the dressing percentage and lean meat percentage of indoor-fed yaks at low-altitude were higher than those of grazing yaks. Our findings in this study were consistent with the above results. Additionally, the TC group showed the highest dressing percentage (63.11%), which surpassed the 57.39 and 54.97% observed in the other two groups. Our study indicated that under low-altitude indoor feeding conditions, an improved feeding environment and a more comprehensive nutrient intake contributed to enhanced carcass characteristics in gayals.
Different production system also affected fat deposition. Fat deposition dependents on residual net energy expenditure, with higher dietary energy levels and total carcass fat contents increasing (22), thus indoor feeding led to increased backfat thickness and IMF in the TC group compared to the JMD group (p < 0.05), however, no significant difference was found between the FHS and JMD groups (p > 0.05). The observed differences may be attributed to several factors, the FHS group was provided with feed of lower nutrient levels relative to the TC group, additionally, due to the brief transition period from grazing to indoor feeding in the FHS group, they had not fully acclimated to the environment of feedlot, moreover, the colder ambient temperature at high altitude required more energy for thermoregulation.
In terms of viscera tissue, in this study, the percentage of lung to slaughter weight increased in FHS and JMD group compared to the TC group. This result may be due to the fact that the lungs are the first organ to respond to hypoxic environments, and in order to adapt to high altitude environments, animals tend to have larger lungs (23). In contrast, the TC group (low-altitude indoor feeding system) faced less hypoxic stress, resulting in a relatively smaller lung proportion.
4.2 Meat quality and chemical composition
A diet rich in feed can lead to accumulation of more glycogen in the muscle, the conversion of glycogen to lactate after slaughter contributes to decreasing pH value (24), explaining why TC muscle pH value was significantly lower than that of JMD group (p < 0.05). Furthermore, the higher activity of lactate dehydrogenase in the LD muscle of cattle at high altitude (25) resulted in a minimal difference in muscle pH value between JMD and FHS groups, but both were higher than TC group (p < 0.05). This study demonstrated that the pH value of the BF muscle was higher than that of the LD muscle, which is consistent with previous findings by Ge et al. (19).
Carvalho et al. (26) suggested that diets rich in lutein and beta-carotene increased the b* value of meat. In our study, the BF muscle b* value in the JMD group was significantly higher than that in the FHS group (p < 0.05), indicating a substantial presence of lutein and beta-carotene in the JMD group. The significant increase in BF muscle b* value in the TC group compared to the FHS group (p < 0.05) may be attributed to the addition of more corn meal to the TC diet, which, as a feed component may lead to an increase in b* value (27). The change pattern of BF muscle b* values across the three groups was consistent with the trend in fat color reported in Table 2.
Shear force serves as an intuitive indicator of muscle tenderness (28). The size of muscle fiber diameter is a crucial determinant of tenderness, with tenderness diminishing as the fiber diameter increases (29). The LD muscle fiber diameter was larger in the TC group than in the FHS group, and the BF muscle fiber diameter increased in the TC group compared to the FHS and JMD groups (p < 0.05). However, IMF also significantly influences meat tenderness, with an increase in IMF leading to a decrease in shear force. In this study, the IMF of the TC group was significantly higher than that of the other two groups (p < 0.05). Therefore, the above synthesis explained why no significant difference in shear force was observed among the three groups.
Among the nutrients provided by beef, protein is the most valuable component for consumers. Ge et al. (19) conducted a study to measure the crude protein (CP) content of muscle in nine different breeds of beef cattle in China. The results indicated that there was variability in the CP content within the same breed of cattle across different muscle parts, specifically, the CP content of the LD muscle was the highest in Luxi cattle (23.27%), while the BF muscle of Yunling cattle exhibited the highest CP content (23.66%). Our results showed that the LD muscle of the TC group had the highest CP content (23.57%), followed by the BF muscle of the JMD group (22.45%). Our findings fall within the range of previously reported values, and the variation in CP content may arise from differences in muscle fiber composition and basal metabolic activity between the LD and BF muscles. Additionally, the CP content of the BF muscle was higher in the JMD group than in the TC group (p < 0.05), whereas the CP content of the LD muscle did not differ significantly between the TC and JMD groups (p > 0.05), consistent with the findings of Yuan et al. (12). This can be attributed to the inverse relationship between moisture content and crude protein in the LD and BF muscles: as moisture content increases, crude protein tends to decrease, and vice versa. Furthermore, no significant difference in moisture content was observed among the three groups in the LD muscle.
Meat serves as a primary dietary source of cholesterol. The dietary guidelines issued by the US Department of Agriculture recommend a daily cholesterol intake of < 300 mg (30). In our study, if 100 g of gayal beef were consumed daily, the cholesterol intake would represent only 7.79%−11.41% of the recommended maximum intake, indicating the potential benefits of consuming gayal beef due to its relatively low cholesterol content. Interestingly, the cholesterol levels in the JMD and FHS groups were significantly higher than those in the TC group (p < 0.05), acclimatization to high altitudes showed a significant positive correlation with increasing serum cholesterol, it was mainly due to the fact that cholesterol promotes hypoxic adaptation at high altitude in vivo and in vitro by inducing the expression of hypoxia-inducible factor (HIF) 1-alpha and inducible nitric oxide synthase (iNOS), thereby protecting hepatocytes against hypoxia-mediated cell death (31).
4.3 Fatty acid
Currently, health-conscious consumers are increasingly concerned about the content and composition of fatty acid in beef. A lower AI and TI are utilized to assess the healthy fatty acid composition of meat (7). The h:H ratio is another metric for evaluating meat healthiness, the value above two signifies a healthier fatty acid composition, which is associated with a reduced risk of cardiovascular disease (21). In our study, no statistically significant differences were observed in the AI or h:H ratio among the three groups (p > 0.05). This lack of significance may be partly attributed to the relatively small sample size (n = 5 per group), which could limit the statistical power to detect subtle differences. Nevertheless, compared with the TC group, the JMD group exhibited a trend toward a lower AI and consistently maintained a higher h:H ratio. More importantly, the TI in the LD muscle was significantly lower in the JMD group than in the TC group (p < 0.05), indicating that the meat from gayals in the grazing system had greater health benefits.
C18:3n3 and C18:2n6 are classified as essential fatty acids (EFA). C18:3n3 can be prolongated and desaturated to synthesize C20:5n3, C22:5n3, and C22:6n3 (32). C18:3n3, C20:5n3 and C22:6n3 are EFA that regulate cellular metabolism and function, reduce plasma triglyceride levels, enhance myocardial relaxation and contraction, and mitigate the risk of cardiovascular diseases (33). C18:2n6 is metabolically converted to C20:4n6 through processes of desaturation and chain elongation (34). C20:4n6 affects the activity of ion channels and the fluidity of cell membranes, therefore, it plays a pivotal role in the brain, retina, and central nervous system (35). In our study, C18:2n6, C18:3n3, C20:4n6, C20:5n3, C22:6n3, and C22:5n3 in the LD muscle and C18:3n3, C20:5n3, C22:6n3, and C22:5n3 in the BF muscle of JMD group were significantly higher than those in the TC group (p < 0.05), indicating that grazing gayals could provide healthier fatty acid for humans. Interestingly, C20:5n3, C22:6n3, C22:5n3, and C20:4n6 also showed higher levels in the FHS group, which may be attributed to the similar altitude between the FHS and JMD groups and the relatively brief transition period from grazing to indoor feeding for the FHS group.
Grass-fed beef contains a higher concentration of C18:1 t11, which serves as an important precursor for the de novo synthesis of conjugated linoleic acid (CAL: C18:2 c9, t11), a potent anti-carcinogenic fatty acid formed within body tissues (36). A linear correlation was observed between dietary C18:1 t11 levels and the rate of CAL synthesis in human subjects (37). In our study, the content of C18:1 t11 in the BF muscle of the JMD group was significantly higher than that in the TC group (p < 0.05). Therefore, we can speculate that the BF muscle of the JMD group contains more CAL. Previous research has also confirmed that grass-fed cattle produce beef with a higher content of CLA (38).
The UK Department of Health recommended a PUFA: SFA ratio above 0.45 (39), Gibson et al. (40) reported that n-6:n-3 PUFA ratio close to 1 or 2 and < 5 is optimal for health. In the current study, the PUFA: SFA ratio in the LD muscle of the JMD group was 0.93, which was significantly higher than that of the TC group (p < 0.05), and the n-6:n-3 PUFA ratio of LD and BF muscles of the JMD group were 2.82 and 2.92, respectively, which were significantly lower than those in the TC group (p < 0.05), these findings further confirmed the potential of the JMD group in balancing fatty acid composition to provide healthier meat products.
4.4 Amino acid
Meat serves as a crucial source of high-quality protein, providing humans with a substantial supply of amino acid, some of which play significant physiological roles in human health (41). Methionine and branched-chain amino acids (BCAA) are classified as EAA. In the LD muscle, the levels of methionine, valine, and isoleucine in the TC group were found to be significantly higher than those in the JMD group, while in BF muscle, the levels of valine and isoleucine in the TC group were significantly higher than those in the FHS group (p < 0.05), indicating that beef from the TC group was more advantageous for human health in terms of amino acid nutritional value.
Free amino acids are water-soluble compounds that contribute to flavor formation (42), the composition and proportion of amino acid in muscle tissue significantly contribute to the overall flavor profile of the muscle (43). Bitter amino acids, such as valine, methionine, isoleucine, leucine, tyrosine, and phenylalanine, contrast with sweet amino acids, which include glycine, serine, threonine, alanine, proline, and lysine, and umami amino acids consist of aspartic acid and glutamic acid (44). In our study, both the LD and the BF muscles in the TC group exhibited the highest content of EAA, however, the TC group also had the highest bitter amino acids levels, as most EAA were bitter amino acids. Furthermore, the contents of sweet amino acids in the LD muscle and umami amino acids in the BF muscle of the TC group were significantly lower than those in the JMD group (p < 0.05), indicating that although the TC group provide high-quality meat enriched in EAA, its poor flavor make the JMD group a more preferred choice for consumers focused on taste. Grazing has been shown to increase flavor amino acid levels in animal muscle tissue, with this elevated content closely associated with natural forage quality and grazing intensity (45). Previous studies have also reported similar findings, showing that grazing increases the contents of sweet and umami amino acids in the LD muscle of Tan lambs compared with the indoor feeding group, an effect that may be attributed to changes in key rumen bacteria, such as Schwartzia and Moryella, which may affect amino acid deposition in meat (46).
4.5 Principal component analysis
Principal component analysis (PCA) effectively visualized the relationships among the three production systems. As shown in the score plot (Figure 1B), the TC group exhibited a distinct separation from the other two groups, whereas significant overlap was observed between FHS and JMD groups. Based on the loading values (Figure 1A), dressing percentage, EAA, PUFA, sweet amino acids and bitter amino acids were the primary factors for distinguishing the TC group from the FHS and JMD groups. Specifically, the TC group displayed a significantly higher dressing percentage and richer content of EAA and bitter amino acids. In contrast, FHS and JMD groups showed elevated levels of PUFA and sweet-tasting amino acids. This result may be attributed to the similarity in altitude between the FHS and JMD groups, as well as the short transition time from grazing to indoor feeding conditions for the FHS group, which resulted in the FHS group appearing closer to the JMD group in terms of slaughter performance, meat quality, and muscle nutritional value.
5 Conclusions
This study demonstrates that the production system significantly influences the slaughter performance and meat quality of gayals. Gayals raised under low-altitude indoor feeding conditions (TC group) exhibited superior slaughter performance, including higher dressing percentage, lean meat percentage, and meat: bone ratio, which contributed to an increased yield of beef products. The meat from the TC group was characterized by a lower cholesterol content, higher levels of EAA, and an enriched content of MUFA, particularly C16:1 and C18:1 c9. However, this system also resulted in a higher IMF content and a less favorable fatty acid profile from a cardiovascular health perspective, as evidenced by a higher n-6:n-3 PUFA ratio and a higher TI value. In contrast, gayals subjected to high-altitude grazing (JMD group) displayed reduced slaughter performance, yet produced beef with a healthier fatty acid composition, including higher levels of PUFAs, such as C18:3n3, C20:5n3, C22:5n3, and C22:6n3, a lower n-6:n-3 PUFA ratio, and a lower TI value. This group also exhibited higher contents of flavor amino acids, specifically sweet amino acids in the LD muscle and umami amino acids in the BF muscle, suggesting an enhanced sensory quality. The high-altitude indoor feeding system (FHS group) yielded results intermediate to the TC and JMD groups for most parameters. However, due to the similar altitude and a relatively brief transition period from grazing to indoor feeding, the slaughter performance and meat quality in the FHS group were found to be intermediate between the TC and JMD groups, with characteristics more closely resembling the JMD group.
A major limitation of this study was the restricted sample size (n = 5 per group), which was largely constrained by the small population of gayals in China (< 4,000). Future research should aim to incorporate larger sample sizes where feasible and to include longer adaptation periods when transitioning animals between different production systems, to enable a more robust and comprehensive evaluation of how different production systems affect the slaughter performance and meat quality of gayals. Further research could also explore the effects of each production system on animal health, meat sensory attributes, and consumer acceptability, thereby providing foundational data for sustainable production strategies.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Ethics statement
This research was approved by the Laboratory Animal Management and Use Committee of Yunnan Agricultural University (Approval No. YNAU 20180020). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
LZ: Data curation, Formal analysis, Software, Visualization, Writing – original draft, Writing – review & editing. QL: Conceptualization, Methodology, Supervision, Writing – review & editing. ZW: Validation, Writing – review & editing. MZ: Software, Writing – review & editing. ShiH: Investigation, Writing – review & editing. HZ: Software, Writing – review & editing. SD: Validation, Writing – review & editing. LZ: Investigation, Writing – review & editing. SZ: Validation, Writing – review & editing. FJ: Investigation, Writing – review & editing. LW: Project administration, Writing – review & editing. SheH: Project administration, Writing – review & editing. DK: Project administration, Writing – review & editing. CD: Project administration, Writing – review & editing. MY: Project administration, Writing – review & editing. HM: Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by the Key Industry (Dairy Cattle) Expert Group Work for Agricultural Modernization (K2420240105), and Creation of Jingtian Science and Technology Innovation Demonstration Village in Huize County - Direct (K2420240061).
Conflict of interest
DK was employed by Yunnan Tengsen Investment Group Limited.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Gao H, Yu Y, Lv Y, Wang D, Li H, Li Z, et al. Metagenomic sequencing reveals the taxonomic and functional characteristics of rumen micro-organisms in gayals. Microorganisms. (2023) 11:1098. doi: 10.3390/microorganisms11051098
2. Wu DD, Ding XD, Wang S, Wójcik JM, Zhang Y, Tokarska M, et al. Pervasive introgression facilitated domestication and adaptation in the Bos species complex. Nat Ecol Evolu. (2018) 2:1139–45. doi: 10.1038/s41559-018-0562-y
3. Dorji T, Wangdi J, Shaoliang Y, Chettri N, Wangchuk K. Mithun (Bos frontalis): the neglected cattle species and their significance to ethnic communities in the Eastern Himalaya - a review. Anim Biosci. (2021) 34:1727–38. doi: 10.5713/ab.21.0020
4. Mukherjee S, Mukherjee A, Jasrotia RS, Jaiswal S, Iquebal MA, Longkumer I, et al. Muscle transcriptome signature and gene regulatory network analysis in two divergent lines of a hilly bovine species Mithun (Bos frontalis). Genomics. (2020) 112:252–62. doi: 10.1016/j.ygeno.2019.02.004
5. Lalchamliani, Chauhan G, Ramesh V, Shivanagowda GP, Khate K, Karabasanavar NS, et al. Carcass characteristics, meat quality, and nutritional profiles of Mithun (Bos frontalis) meat reared under a semi-intensive system. Anim Sci J. (2024) 95:e13942. doi: 10.1111/asj.13942
6. Sun Y, Zhang X, Xi D, Li G, Wang L, Zheng H, et al. Isolation and cDNA characteristics of MHC-DRA genes from gayal (Bos frontalis) and gaytle (Bos frontalis × Bos taurus). Biotechnol Biotechnol Equip. (2015) 29:33–9. doi: 10.1080/13102818.2014.986128
7. Contò M, Cifuni GF, Iacurto M, Failla S. Effect of pasture and intensive feeding systems on the carcass and meat quality of buffalo. Anim Biosci. (2022) 35:105–14. doi: 10.5713/ab.21.0141
8. Janssen J, Cammack K, Legako J, Cox R, Grubbs JK, Underwood K, et al. Influence of grain- and grass-finishing systems on carcass characteristics, meat quality, nutritional composition, and consumer sensory attributes of bison. Foods. (2021) 10:1060. doi: 10.3390/foods10051060
9. Liu YX, Ma XM, Xiong L, Wu XY, Liang CN, Bao PJ, et al. Effects of intensive fattening with total mixed rations on carcass characteristics, meat quality, and meat chemical composition of yak and mechanism based on serum and transcriptomic profiles. Front Vet Sci. (2020) 7:599418. doi: 10.3389/fvets.2020.599418
10. Nogoy KMC, Sun B, Shin S, Lee Y, Zi Li X, Choi SH, et al. Fatty acid composition of grain- and grass-fed beef and their nutritional value and health implication. Food Sci Anim Resour. (2022) 42:18–33. doi: 10.5851/kosfa.2021.e73
11. Xin K, Tian K, Yu Q, Han L, Zang Z. Effects of altitude on meat quality difference and its relationship with HIF-1α during postmortem maturation of beef. J Food Biochem. (2022) 46:e14470. doi: 10.1111/jfbc.14470
12. Yuan P, Zhang X, Guan J, An T, Li H, Tan W, et al. Effects of low altitude house-feeding on growth performance, slaughter performance and muscle quality of yak and cattle-yak. Chin J Anim Nutr. (2021) 33:4791–800.
13. Zhang Z, Wang Y, Chen Y, Zhang L, Zhang Y, Liu Y, et al. The dietary supplemental effect of nitroethanol in comparison with monensin on methane emission, growth performance and carcass characteristics in female lambs. Animals. (2021) 11:327. doi: 10.3390/ani11020327
14. Qiu S, Li K, He X, Gu M, Jiang X, Lu J, et al. The effects of composite alkali-stored spent Hypsizygus marmoreus substrate on carcass quality, rumen fermentation, and rumen microbial diversity in goats. Animals. (2024) 14:166. doi: 10.3390/ani14010166
15. Tian XZ, Li JX, Luo QY, Wang X, Xiao MM, Zhou D, et al. Effect of supplementation with selenium-yeast on muscle antioxidant activity, meat quality, fatty acids and amino acids in goats. Front Vet Sci. (2021) 8:813672. doi: 10.3389/fvets.2021.813672
16. Wu J, Liu Y, Xu N. Histological characteristics of musculus longissimus dorsi and their correlation with restriction fragment length polymorphism of the myogenin gene in Jinghua × Pietrain F(2) crossbred pigs. Meat Sci. (2009) 81:108–15. doi: 10.1016/j.meatsci.2008.07.010
17. Han L, Yu Y, Fu R, Fu B, Gao H, Li Z, et al. Impact of various ration energy levels on the slaughtering performance, carcass characteristics, and meat qualities of Honghe yellow cattle. Foods. (2024) 13:1316. doi: 10.3390/foods13091316
18. Wang Y, Wang Q, Dai C, Li J, Huang P, Li Y, et al. Effects of dietary energy on growth performance, carcass characteristics, serum biochemical index, and meat quality of female Hu lambs. Anim Nutr. (2020) 6:499–506. doi: 10.1016/j.aninu.2020.05.008
19. Ge F, Li J, Gao H, Wang X, Zhang X, Gao H, et al. Comparative analysis of carcass traits and meat quality in indigenous Chinese cattle breeds. J Food Compos Anal. (2023) 124:105645. doi: 10.1016/j.jfca.2023.105645
20. Latimer GW; AOAC International. Official Methods of Analysis of AOAC International. 19th ed. AOAC International (2012).
21. Silva J, Rodrigues LS, Lourenço-Júnior JB, Alfaia CM, Costa MM, Castro VCG, et al. Total lipids, fatty acid composition, total cholesterol and lipid-soluble antioxidant vitamins in the longissimus lumborum muscle of water buffalo (Bubalus bubalis) from different production systems of the Brazilian eastern amazon. Animals. (2022) 12:595. doi: 10.3390/ani12050595
22. Muir PD, Deaker JM, Bown MD. Effects of forage- and grain-based feeding systems on beef quality: a review. N Z J Agric Res. (1998) 41:623–35. doi: 10.1080/00288233.1998.9513346
23. Low AT, Medford AR, Millar AB, Tulloh RM. Lung function in pulmonary hypertension. Respir Med. (2015) 109:1244–49. doi: 10.1016/j.rmed.2015.05.022
24. Priolo A, Micol D, Agabriel J. Effects of grass feeding systems on ruminant meat colour and flavour. A review Anim Res. (2001) 50:185–200. doi: 10.1051/animres:2001125
25. Lin YQ, Wang GS, Feng J, Huang JQ, Xu YO, Jin SY, et al. Comparison of enzyme activities and gene expression profiling between yak and bovine skeletal muscles. Livest Sci. (2011) 135:93–7. doi: 10.1016/j.livsci.2010.06.134
26. Carvalho ME, Gasparin G, Poleti MD, Rosa AF, Balieiro JCC, Labate CA, et al. Heat shock and structural proteins associated with meat tenderness in Nellore beef cattle, a Bos indicus breed. Meat Sci. (2014) 96:1318–24. doi: 10.1016/j.meatsci.2013.11.014
27. Stefanello FS, Fruet APB, Trombetta F, da Fonseca PAF, dos Santos da Silva M, Stefanello S, et al. Stability of vacuum-packed meat from finishing steers fed different inclusion levels of brewer's spent grain. Meat Sci. (2019) 147:155–61. doi: 10.1016/j.meatsci.2018.09.004
28. Ding Y, Jiang X, Yao X, Zhang H, Song Z, He X, et al. Effects of feeding fermented mulberry leaf powder on growth performance, slaughter performance, and meat quality in chicken broilers. Animals. (2021) 11:3294. doi: 10.3390/ani11113294
29. Renand G, Picard B, Touraille C, Berge P, Lepetit J. Relationships between muscle characteristics and meat quality traits of young Charolais bulls. Meat Sci. (2001) 59:49–60. doi: 10.1016/S0309-1740(01)00051-1
30. Carson JAS, Lichtenstein AH, Anderson CAM, Appel LJ, Kris-Etherton PM, Meyer KA, et al. Dietary cholesterol and cardiovascular risk: a science advisory from the American heart association. Circulation. (2020) 141:e39–53. doi: 10.1161/CIR.0000000000000743
31. Miron N, Tirosh O. Cholesterol prevents hypoxia-induced hypoglycemia by regulation of a metabolic ketogenic shift. Oxid Med Cell Longevity. (2019) 2019:5829357. doi: 10.1155/2019/5829357
32. Scollan N, Hocquette JF, Nuernberg K, Dannenberger D, Richardson I, Moloney A. Innovations in beef production systems that enhance the nutritional and health value of beef lipids and their relationship with meat quality. Meat Sci. (2006) 74:17–33. doi: 10.1016/j.meatsci.2006.05.002
33. Mozaffarian D, Wu JHY. Omega-3 fatty acids and cardiovascular disease: effects on risk factors, molecular pathways, and clinical events. J Am Coll Cardiol. (2011) 58:2047–67. doi: 10.1016/j.jacc.2011.06.063
34. Burr GO, Burr MM. On the nature and rôle of the fatty acids essential in nutrition. J Biol Chem. (1930) 86:587–621. doi: 10.1016/S0021-9258(20)78929-5
35. Nejatian M, Darabzadeh N, Bodbodak S, Saberian H, Rafiee Z, Kharazmi MS, et al. Practical application of nanoencapsulated nutraceuticals in real food products; a systematic review. Adv Colloid Interface Sci. (2022) 305:102690. doi: 10.1016/j.cis.2022.102690
36. Bauman DE, Lock AL. Conjugated Linoleic Acid: Biosynthesis and Nutritional Significance. Springer US (2006).
37. Turpeinen AM, Mutanen M, Aro A, Salminen I, Basu S, Palmquist DL, et al. Bioconversion of vaccenic acid to conjugated linoleic acid in humans. Am J Clin Nutr. (2002) 76:504–10. doi: 10.1093/ajcn/76.3.504
38. Leheska JM, Thompson LD, Howe JC, Hentges E, Boyce J, Brooks JC, et al. Effects of conventional and grass-feeding systems on the nutrient composition of beef. J Anim Sci. (2008) 86:3575–85. doi: 10.2527/jas.2007-0565
39. Committee on Medical Aspects of Food Policy CRG. Nutritional Aspects of Cardiovascular Disease. 1th ed. HMSO Publications Centre (1994).
40. Gibson RA, Neumann MA, Lien EL, Boyd KA, Tu WC. Docosahexaenoic acid synthesis from alpha-linolenic acid is inhibited by diets high in polyunsaturated fatty acids. Prostaglandins, Leukotrienes Essent Fatty Acids. (2013) 88:139–46. doi: 10.1016/j.plefa.2012.04.003
41. Wu S, Luo X, Yang X, Hopkins DL, Mao Y, Zhang Y. Understanding the development of color and color stability of dark cutting beef based on mitochondrial proteomics. Meat Sci. (2020) 163:108046. doi: 10.1016/j.meatsci.2020.108046
42. Dinh TTN, Legako JF, Miller MF, Brooks JC. Effects of USDA quality grade and cooking on water-soluble precursors of beef flavor. Meat Sci. (2018) 146:122–30. doi: 10.1016/j.meatsci.2018.08.008
43. Yin J, Li Y, Zhu X, Han H, Ren W, Chen S, et al. Effects of long-term protein restriction on meat quality, muscle amino acids, and amino acid transporters in pigs. J Agric Food Chem. (2017) 65:9297–304. doi: 10.1021/acs.jafc.7b02746
44. Gan M, Shen L, Fan Y, Guo Z, Liu B, Chen L, et al. High altitude adaptability and meat quality in Tibetan pigs: a reference for local pork processing and genetic improvement. Animals. (2019) 9:1080. doi: 10.3390/ani9121080
45. Zhang Q, Que M, Li W, Gao S, Tan X, Bu D. Gangba sheep in the Tibetan plateau: validating their unique meat quality and grazing factor analysis. J Environ Sci. (2021) 101:117–22. doi: 10.1016/j.jes.2020.06.024
46. Wang B, Wang Y, Zuo S, Peng S, Wang Z, Zhang Y, et al. Untargeted and targeted metabolomics profiling of muscle reveals enhanced meat quality in artificial pasture grazing tan lambs via rescheduling the rumen bacterial community. J Agric Food Chem. (2021) 69:846–58. doi: 10.1021/acs.jafc.0c06427
Keywords: altitude, different production systems, gayal, meat quality, slaughter performance
Citation: Zhi L, Li Q, Wang Z, Zhang M, He S, Zhang H, Dai S, Zhang L, Zhao S, Jiang F, Wang L, He S, Kang D, Duan C, Yang M and Mao H (2025) Impact of different production systems on the slaughter performance and meat quality of gayal (Bos frontalis). Front. Vet. Sci. 12:1538622. doi: 10.3389/fvets.2025.1538622
Received: 03 December 2024; Accepted: 05 September 2025;
Published: 08 October 2025.
Edited by:
Yafeng Huang, Anhui Agricultural University, ChinaReviewed by:
Sandra Sofia Quinteiro Rodrigues, Centro de Investigação de Montanha (CIMO), PortugalLetian Kang, Inner Mongolia Agricultural University, China
Sibel Soycan Önenç, Namik Kemal University, Türkiye
Copyright © 2025 Zhi, Li, Wang, Zhang, He, Zhang, Dai, Zhang, Zhao, Jiang, Wang, He, Kang, Duan, Yang and Mao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huaming Mao, bWh1YW1pbmdAc2luYS5jb20=
Qing Li1