Introduction
Mutton is gradually becoming a superior substitute for the domestic market, as it is a good source of essential nutrients, high in protein and low in fat (Ramos et al., 2021). In the last few years, the consumption of high quality mutton products has risen dramatically. Chaka sheep as Qinghai Plateau wool and meat with semi-fine wool sheep is a breed of Tibetan sheep, which grows in the area of Chaka Salt Lake in Wulan County, is a cultivar of Chaka sheep (Chen et al., 2019a). In 1971, Wulan County introduced the British and New Zealand Romney-Marsh long-wool type of breeding sheep from Hudong Original Breeding Field of Qinghai Province and Daliuliu Original Breeding Field of Anhui Province to be cross-bred (Chen et al., 2019b). As the breed was improved and bred in the Chaka area, it was named “Chaka sheep”, which has been granted the protection of the Registration of Geographical Indications for agricultural products (Guo et al., 2022). Studies have shown that there are some differences in the nutritional composition, sensory quality, and processing quality of different parts of the muscle (Mohanty et al., 2010). Different parts of muscle have unique collagen content, structure, and properties (Zuo et al., 2016). Kohn et al. (2023) found significant differences in muscle fiber types and metabolic profiles in different parts of African black ostrich muscles.
Lipids, an important nutritional component in lamb, are critical biological small molecules that act as components of cell membranes, signal transduction, energy sources, and signaling pathway intermediates (Lugo Charriez et al., 2021). The lipid profile of meat is influenced by many factors, one of which is the type of muscle (Boselli et al., 2008). The lipid composition of tissues varies in different parts of the body (Yin et al., 2023). Li et al. (2022a) identified seventy-nine of 842 lipids in different muscles of Hu lambs. The content and type of lipids vary significantly in different meat products, and potential markers of lipids have gradually become one of the methods to analyze the quality and type of meat products. However, studies on the lipid profile and characteristics of different parts of the muscles of Chaka sheep are still incomplete. Thus, the lipidomic profile and its components of Chaka sheep need to be elucidated.
Therefore, this study determined and analyzed the nutrient composition, lipid content, and differential lipid molecules in the longissimus dorsi (LD), psoas major (PM), and biceps femoris (BF) of the Chaka sheep (Fig. 1). The aims of this research are to define the basic data about the composition and nutritional quality of three different muscles of the Chaka sheep and to provide insight into the utilization of lipidomics in the animal products industry.
Fig. 1. Experimental and analytical processes followed in this study.
Materials and Methods
Sample collection
The samples were obtained from Jintai Ranch, Chaka Town, Wulan County, Haixi Mongolian and Tibetan Autonomous Prefecture, Qinghai Province, China. We selected nine experimental sheep, all 1 year old, weighing about 45 kg, female. They were identified by experts as Chaka sheep. The Chaka sheep were manually slaughtered by puncturing the occipital atlantoaxial joint of the sheep with a sharp knife, severing the spinal nerves, then cutting the carotid artery to bleed out the blood, and then skinning the sheep, and after removing the head and hooves, 9 replicates each of LD (muscles of the back of the sheep between the 9th and 11th ribs), BF, and PM were collected. The LD, PM and BF of Chaka sheep were cut and packed into lyophilized vessels and then immediately stored in liquid nitrogen for lipidomic study, and the remaining meat samples were wrapped in cling film and placed in a thermal box and returned to the laboratory for nutrient analysis. Three replicates were set up for nutritional quality determination and six replicates for lipidomics.
Determination of the chemical components and characterization of the muscle fibers
An OPTO-LAB meat color analyzer (Konica Minolta, Tokyo, Japan) was employed to measure the color of the muscle. A PH-STAR carcass muscle pH direct tester (MATTHAUS, Portmeuse, Germany) was utilized for pH measurement of all samples. Crude protein (CP), total ash (Ash), ether extract (EE), and moisture content (MSTR) were measured according to Cohen’s methods (Cohen, 1971).
Conventional hematoxylin and eosin (H&E) staining was used to measure muscle fiber characteristics. Samples were scanned using a Leica RM-235 slide scanner (Leica, Nussloch, Germany) and the images were acquired by using CaseViewer 6.0 software (3D-Histech, Budapest, Hungary; Xu et al., 2023).
Determination of amino acids and fatty acids
A Sykam S-433D amino acid analyzer (Sykam, Munich, Germany) was used for amino acid determination. After weighing 50 mg of dry meat powder from a single sample and hydrolysis with hydrochloric acid solution (36% hydrochloric acid: water=1:1, v/v) for 24 h, 2 mL of the diluted specimens were transferred to the vial condenser. 10 mL of sample diluents were poured in and stirred. The compounded sample was withdrawn with a hypodermic syringe, filtered using a needle filter, and injected into the vial to assay (Ma et al., 2021).
The fatty acids were determined by gas chromatography-mass spectrometry (GC-MS). 5 g of the sample was taken, drawn and filtered into the injection vial with a stirred syringe. After adding 8 mL of a sodium hydroxide-methanol solution at 20 g/L relative density, a condensor for fat degradation and fatty acid methyl esterification was attached, n-heptane was poured in, shaken for 2.5 min, and a sodium chloride aqueous buffer solution was poured in. Supernatant was collected, a volume of sodium sulfate hydroxide was added, stirred thoroughly, and supernatant was transferred to vials for automated analysis (Schiavon et al., 2016).
Determination of volatile flavor compounds
The combined technology of solid phase microextraction (SPME) and GC-MS was applied to determine the volatile aroma compounds. The muscle was minced and ground, 4 g meat was added to a 20 mL headspace bottle, 20% NaCl solution and 5 µL of 2-methyl-3-heptanone solution was mixed, stirred thoroughly, and water bath (85℃) for 55 min, then chilled to 55℃. The SPME fiber probe was inserted into the GC-MS injector connection, and the adsorber head was in place for extraction for 40 min. The SPME fiber probe was inserted into the GC-MS injector connection, and the adsorber head was in place during extraction for 40 min, and the fiber was withdrawn for desorbing into the GC-MS insertion port for 3 min by GC-MS. The chromatographic parameters were as follows: the carrier gas was helium, the flow rate was 1.0 mL/min, the injection port temperature was 250℃, the column temperature was started at 40℃, maintained for 5 min, and then increased to 230℃ (5℃/min) for 8 min. The mass spectrometry conditions were as follows: the ion source temperature was 200℃, the ionization mode was EI, the electron energy was 70 eV, and the scanning mass range was 35–500 amu (He et al., 2020).
Lipid extraction
The samples were spiked with the appropriate amount of internal lipid standard, and then homogenized with 200 µL water and 240 µL methanol. Then 800 µL of methyl t-butyl ether was added to the mixture. The mixture was sonicated for 20 min at 55℃ and then allowed to stand for 30 min at room temperature. The upper layer of organic solvent was removed by centrifugation at 14,000×g for 15 min. The samples were dried under nitrogen.
Liquid chromatography with tandem mass spectrometry for the analysis of the lipidome in muscle
A CSH C18 column (1.7 µm, 2.1 mm×100 mm, Waters, Milford, MA, USA) was applied for LC separation. Lipid extracts were dissolved in 200 µL of 90% isopropanol-acetonitrile and incubated at 14,000 g for 15 min. Finally, 3 µL of the specimen was injected. The first mobile phase was 40% acetonitrile-isopropanol (1:9, v/v) at a flow rate of 300 µL per min. This was maintained for 3.5 min, followed by a linear increase to 75% acetonitrile-isopropanol (1:9, v/v) in 9 min. A linear increase to 99% acetonitrile-isopropanol for 6 min and then equilibration at 40% acetonitrile-isopropanol for 5 min.
The mass spectra were collected by the Q-Exactive Plus in the positive mode and in the negative mode, respectively. The ESI parameters were adjusted and optimized for all experiments. The experimental parameters were as follows: the source temperature was 300℃, the capillary temperature was 350℃, the ion spray voltage was 3,000 V, and the instrument scan range was 200–1,800 m/z (Tang et al., 2023).
Statistical analysis
Results of the study were presented as mean±SD in this assay with three replicates for each sample. SPSS 22.0 (IBM, Armonk, NY, USA) was used to analyze the regular nutritional composition data statistically. Orthogonal partial least squares discriminant analysis (OPLS-DA), variable importance in projection (VIP)>1, p-value<0.05 was defined as the criterion for significant differences.
Results and Discussion
Analysis of chemical components and muscle fiber parameters
There was no significant difference in the results of routine nutrient content measurements (Table 1) at pH45min, moisture, ash, and EE, but the CIE a*, CIE b*, and CIE L* of PM were significantly higher than LD and BF (p<0.05). The CIE a* is associated with meat hemoglobin content, CIE b* is negatively correlated with muscle freshness, and CIE L* is related to muscle whiteness (Liang et al., 2023). Therefore, we hypothesized that PM has the highest hemoglobin content and LD has the best freshness, which may be more popular with consumers in terms of meat color. The rank order of CP levels in three parts of Chaka sheep was: LD>BF>PM.
Table 1. The routine nutrient composition of different parts of muscles
The MSTR, Ash, CP and EE are on a dry matter basis.
a–c Means within a row with different subscripts differ when p-value<0.05.
LD, longissimus dorsi; PM, psoas major; BF, biceps femoris; MSTR, moisture; Ash, ash; EE, ether extract; CP, crude protein.
Previous studies showed that muscle fiber properties, such as type, size and total number of fibers in different parts were different, the higher muscle fiber density, the higher economic value of meat (Qin et al., 2020). The results of myofiber measurements are shown in Fig. 2, the total area of myofibers did not differ significantly between LD, PM, and BF, the total number of myofibers was ranked from most to least: PM>LD>BF, the density of myofibers was ranked: PM>LD>BF, the diameter of myofibers was the largest in BF. Our study demonstrated significant differences in total myofiber number, myofiber density, and myofiber diameter in LD, PM, and BF (p<0.05). The myofiber density in LD (1,557.93 n/mm2) and PM (1,588.74 n/mm2) were significantly higher than in BF (781.11 n/mm2), so it may suggested that LD and PM have a higher commercial value than BF.
Fig. 2. Characteristics of muscle fibers in different parts of the muscles. (A) Hematoxylin and eosin (H&E) staining, bar=50 μm, (B) total number, (C) total area, (D) the density, and the diameter (E) of muscle fibers in selected filed of view. a,b Means within a row with different subscripts differ when p-value<0.05. LD, longissimus dorsi; PM, psoas major; BF, biceps femoris.
Determination of amino acids and fatty acids
The results of amino acid measurement are presented in Table 2, the total amino acid and essential amino acid levels in different parts of Chaka sheep were not significant, but our previous study has indicated that there were differences in CP content. This inconsistency may be caused by calculating CP based on nitrogen content. However, in addition to protein, meat also includes other nitrogenous compounds such as phospholipids (PLs; Gao et al., 2022), so the level of CP was not exactly equal to TAA. There were significant differences in the levels of aspartic acid, glycine, histidine, and proline in different parts of the muscle. The contents of aspartic acid and histidine were the highest in LD, followed by BF, and the lowest in PM. The contents of glycine and proline were in the following order high to low: BF>LD>PM. Studies have indicated that the free precursor amino acids affect the flavor of the meat, which is essential in evaluating the quality of the meat (Tian et al., 2021), glycine often imparts a sweet flavor to mutton, aspartate gives the meat and umami taste, and favorited by consumers (Ribeiro de Araújo Cordeiro et al., 2020). Histidine is essential for producing histamine, a chemical messenger essential for immune reactions, digestion, and sleep/wake cycles (Zhang et al., 2022a). LD had the highest levels of aspartic acid and histidine, and BF had the highest levels of total glycine and proline, from which we hypothesized that the flavors of LD and BF might be more likely to be preferred by consumers.
Table 2. The amino acid content in different parts of Chaka sheep
The content of amino acids (%) in this study is expressed as the percentage of amino acids per unit mass of the sample.
EAA=Threonine+Valine+Methionine+Isoleucine+Leucine+Phenylalanine+Histidine+Lysine+Arginine.
NEAA=TAA–EAA.
a,b Means within a row with different subscripts differ when p-value<0.05.
LD, longissimus dorsi; PM, psoas major; BF, biceps femoris; EAA, essential amino acids; NEAA, nonessential amino acids; TAA, total amino acids.
Fatty acids are closely related to the eating quality of meat by determining the melting point and quality of meat fats (Dai et al., 2024). The contents of fatty acids in different parts of the muscle are shown in Table 3. A study has shown that intramuscular fat content was positively correlated with SFA (Zhang et al., 2019), which was consistent with the composition of crude fat in PM, LD, and BF in our findings. The content of linoleic acid (C18:2n6c) was, in decreasing order, ranked as follows: PM>BF>LD. Haag’s study found that eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) are metabolic byproducts of linoleic acid and linolenic acid, which are vital to the human body in the maintenance of cardiovascular health, brain function, and the reduction of pain in rheumatoid arthritis (Haag, 2003; Oh, 2005). Therefore, we hypothesized that the higher linoleic acid in PM would favor the production and accumulation of EPA and DHA. The total of ω-6 polyunsaturated fatty acid (PUFA) was ranked as: PM>BF>LD. Studies have found that ω-6 PUFA cannot be produced by the body, but must be consumed in the diet, and meat is an essential source of ω-6 PUFA. It is recommended that people should get these nutrients through a balanced diet to maintain the body’s ω-6 PUFA level (Simopoulos, 2008). Overall, our results indicated that there were differences in the fatty acid content of muscles from Chaka sheep, the PM exhibited the highest nutritional value in the fatty acids profile.
Table 3. The fatty acid content in different parts of Chaka sheep
∑SFA=C4:0+C6:0+C8:0+C10:0+C11:0+C12:0+C13:0+C14:0+C15:0+C16:0+C17:0+C18:0+C20:0+C21:0+C22:0+C23:0+C24:0.
∑MUFA=C14:1+C15:1+C16:1+C17:1+C18:1n9t+C18:1n9c+C22:1n9+C20:1+C24:1.
∑PUFA=C18:2n6t+C18:2n6c+C18:3n6+C18:3n3+C20:2+C20:3n3+C20:3n6+C20:4n6+C20:5n3+C22:2+C22:6n3.
∑UFA=∑MUFA+∑PUFA.
ω-3 PUFA=C18:3n3+C20:3n3+C20:5n3+C22:6n3.
ω-6 PUFA=C18:2n6t+C18:2n6c+C18:3n6+C20:3n6+C20:4n6.
a–c Means within a row with different subscripts differ when p-value<0.05.
LD, longissimus dorsi; PM, psoas major; BF, biceps femoris; SFA, saturated fatty acids; MUFA, monounsaturated fatty acid; PUFA, polyunsaturated fatty acid; UFA, unsaturated fatty acids; ND, not detected.
Determination of volatile flavor compounds
As presented in Table 4, aldehyde, alcoholic, ketonic, carboxylic, esteric, and other volatile aroma compounds have been identified in Chaka sheep muscle. There were marked differences in the alcohols, ketones, and esters of the Chaka sheep, with PM having the highest levels of total alcohols, ketones, and esters. Alcohols are derived from the breakdown of lipids and have a high olfactory threshold and a generally clear odor similar to that found in fruits and vegetables. Ketones, a category of volatile organic compounds with floral, fruity, and creamy aromas, have an impact on the overall flavor of meat products (Ge et al., 2023). Ester compounds, most of which have an aromatic flavor, are formed by the interesterification of acids and alcohols (Yang et al., 2024). We found that the PM exhibited the highest nutritional value in the fatty acids profile, so we speculated that the higher fatty acid oxidation may lead to the accumulation of flavor substances in PM.
Table 4. The content of volatile flavor substances in three parts of muscle
a–c Means within a row with different subscripts differ when p-value<0.05.
LD, longissimus dorsi; PM, psoas major; BF, biceps femoris; ND, not detected.
Lipid profiles analysis
The lipid composition analysis in different muscle parts
Lipids play an essential role in the production of the odor and the flavor of meat via autoxidation or thermal oxidation (Li et al., 2021). A detailed classification of the categories is crucial for the analysis of the metabolic pathways that underlie the transformation of lipids (Enser et al., 1999). A total of 2,639 lipid molecule species were detected in the Chaka sheep. The lipids were composed of 42 categories (Fig. 3), and the percentages of content, in descending order, were: triacylglycerol (TG: 18.91%), phosphatidylcholine (PC: 18.34%), phosphatidylethanolamine (PE: 14.25%), cardiolipid (CL: 6.14%), diacylglycerol (DG: 5.99%), monosaccharide ceramide (Hex1Cer: 5.53%), phosphatidylserine (PS: 4.96%), sphingomyelin (SM: 3.68%), ceramide (Cer: 3.64%) and others. Intramuscular fat is composed of PLs, TGs, and cholesterol, and the balance between triglyceride synthesis, breakdown, and absorption is reflected in IMF content (Roy et al., 2024). The effect of the lipid content on aroma retention is greater than the influence of the lipid type (Ayed et al., 2018). TGs and PLs have been shown to be the predominant lipids for predominant lipids for binding and generating flavor compounds (Fu et al., 2022). Liu’s study suggested that the highest types and concentrations of TG may contribute primarily to aroma retention in roasted mutton (Liu et al., 2022). Therefore, our hypothesis was that the high TG levels in Chaka sheep might be correlated with the taste of the meat. PC, PE, and CL being the major components of PLs. Zhou’s study found that the different PL forms (PC and PE) not only can regulate membrane fluidity and permeability but also are important components of PLs (Zhou et al., 2014). PLs, which accounted for almost half of the total number of molecules of lipid classes, were predominant in these three samples in our study. A study showed that PLs are abundant in PUFAs and are an essential resource of PUFAs (Zhou et al., 2013), and the flavor of meat is associated with the composition of PLs (Wu and Wang, 2019). Wang et al. (2009) found that the PLs were the main factor in forming the flavor of duck meat products, which also was one of the major determinants of the taste of duck meat products. Therefore, we hypothesized that the flavor of Chaka sheep would be associated with the levels of PLs.
Fig. 3. The molecular distribution of total lipids in Chaka sheep. TG, triacylglycerol; SM, sphingomyelin; PS, phosphatidylserine; PE, phosphatidylethanolamine; PC, phosphatidylcholine; DG, diacylglycerol; CL, cardiolipid.
The total number of lipid molecules in the different parts of the Chaka sheep is summarized in Fig. 4. The comparative results were as follows: PM>BF>LD. The difference between the groups was not significant. The composition of lipid subclasses in different parts of the Chaka sheep is shown in Fig. 5. PC, PE, phosphatidylinositol (PI), PS, sphingosine (SPH), TG, and WE were major subclasses lipids in Chaka, and all of them were present at levels >1,000 μg/g. TG levels varied significantly among the three comparison groups, with the highest levels in the BF group, followed by the PM group, and the lowest levels in the LD group (p<0.05). PC, PE, PI, PS, and SPH were belong to PLs. PC and PE are the major PL classes in porcine intramuscular fat, and PEs are the predominant glycerophospholipid in ruminant cell membranes (Fernandez et al., 1999). In conclusion, PLs were the highest total percentage of Chaka sheep, and PLs are also essential constituents of cell membranes. Hazel’s study indicated that long-chain PUFA is more highly mobilized at low temperatures, and the higher the PL content, the higher the long-chain PUFA content, which may correlate with better plateau adaptation in yak (Hazel and Williams, 1990). Lin’s study showed that fish muscle PLs were sensitive to decreasing temperature (Lin et al., 2019). Therefore, we hypothesized that the higher levels of PLs molecules in the muscles of Chaka sheep would facilitate their adaptation to the extreme environment of the Tibetan Plateau, which is characterized by high levels of cold and hypoxia.
Fig. 4. The total content of lipid molecules in different parts of Chaka sheep.
Fig. 5. The composition of lipid subclasses in different parts of the Chaka sheep. a−c Indicate significant differences between the groups compared in the figure (p<0.05). LD, longissimus dorsi; PM, psoas major; BF, biceps femoris; LPC, lysophosphatidylcholine; PA, phosphatidic acid; PC, phosphatidylcholine; PE, phosphatidylethanolamine; PI, phosphatidylinositol; PS, phosphatidylserine; SM, sphingomyelin; SPH, sphingosine; TG, triacylglycerol; WE, wax esters.
The differential lipid molecular analysis
To study the significant differences in lipids between different parts of Chaka sheep, OPLS-DA model was used to evaluate differential lipids. The OPLS-DA model maximizes the between-group differences in the horizontal coordinate, while the within-group variation is reflected in the orthogonal principal components in the vertical coordinate, and it is clear in Fig. 6 that the OPLS-DA model can distinguish the two groups of samples, and the model evaluation parameter Q2<0.5, indicating that the model has a great deal of stability, and the test results are real and reliable (Zhang et al., 2022b).
Fig. 6. The OPLS-DA displacement test in positive and negative ion mode. BF, biceps femoris; LD, longissimus dorsi; PM, psoas major; OPLS-DA, orthogonal partial least squares discriminant analysis.
The significantly different lipid molecules (VIP>1, p<0.05) screened in this experiment were visualized in the form of bubble plots in Fig. 7. Bubbles in graph represent significantly different lipid species, and the vertical coordinates show each subclass of lipid, marked by different colors, with the size of the bubbles representing the significance of the differences. In BF, TG, PS, PI, PG, PE, PC, lysophosphatidylethanolamine (LPE), lysophosphatidylcholine (LPC), and Hex1Cer were the highly significantly different lipid molecules compared with LD. In BF vs. PM, TG, PS, PE, PC, LPE, LPC, and Hex1Cer were the main significantly different lipid molecules. In PM vs. LD, TG, SM, PS, PI, PC, LPC, Hex1Cer, and DG were the main lipids with highly significant differences. PS, PE, PC, LPC, LPE in PLs, TG in triglycerides, and Hex1Cer in SMs were the three major co-differentiating lipid molecules. In Sun’s study, PLs were found to be the main factor distinguishing the muscle lipidome between cattle-yak and yak (Sun et al., 2018; Vahmani et al., 2020), and our findings were consistent with these results. Jia’s study mainly detected AcCa, Cer, DG, LPC, PC, PE, SM, and TG, less ChE, LPE, MG, PI, and SPH, and no Co, Hex1Cer, Hex2Cer, LSM, PG in Tan sheep meat (Jia et al., 2021). However, in our study, we detected lipid molecules such as Co, Hex1Cer, Hex2Cer, LSM, and PG. These differences may be related to diet, feeding habits, and other factors(Dannenberger et al., 2017). The top 20 differential lipids in the different comparison groups are shown in the Supplementary Table S1.
Fig. 7. The molecular bubble diagrams of differential lipids in different parts of Chaka sheep. BF, biceps femoris; LD, longissimus dorsi; PM, psoas major; TG, triacylglycerol; PS, phosphatidylserine; PI, phosphatidylinositol; PG, phosphatidylglycerol; PE, phosphatidylethanolamine; PC, phosphatidylcholine; LPE, lysophosphatidylethanolamine; LPC, lysophosphatidylcholine.
We used venn diagrams (Fig. 8) to show specific differential lipid molecules, and the results of the study showed that the numbers of significantly different lipid molecules between the different comparison groups were: 30 (BF vs. LD), 38 (BF vs. PM), and 88 (PM vs. LD). Among them, PE (18:2/18:2) and PC (25:0/11:3) were LD group-associated differential lipids. PS (20:3e/20:4), LPC (18:3), PE (8:1e/12:3), TG (18:0e/16:0/18:1), TG (18:0/18:0/18:0), TG (18:0e/18:0/18:1), and TG (18:0e/18:1/18:1) were PM group-associated differential lipids. LPC (16:0), LPC (18:1), LPE (18:1), PC (15:0/22:6), PE (18:1/18:1), Hex1Cer (d24:1/18:1), and PC (10:0e/6:0) were BF group-associated differential lipids. Dashdorj et al. (2015) found that PE (16:0/20:4), PE (18:0/19:1), and PE (18:1/22:5) were the top three lipid species, suggesting that they are the major contributors to shank/flank muscle liposome differentiation. In our study, PE (18:2/18:2) was mainly found in LD, PE (8:1e/12:3) was predominantly detected in PM, and PE (18:1/18:1) was mostly identified in BF. It can be concluded that PE is an important biomarker for the differentiation of the different parts of the muscle. Li et al. (2022b) showed that PLs containing PC and PE may be key lipids in the formation of aromatic substances in roasted lamb. Zhou et al. (2014) showed that PC is more easily oxidized to aldehydes due to its high content of unsaturated fatty acids. Guo’s study showed that PC and PE are prone to the formation of the aroma compounds (E,E)-2,4-decadienal, 1-octen-3-one and (E)-2-decene (Guo et al., 2021). Our study found significant differences in the content of the flavor compounds alcohols, ketones and esters among the three comparison groups. Therefore, we hypothesized that the cause of the differences in flavor compounds is related to the lipids PE (18:2/18:2), PE (8:1e/12:3), PE (18:1/18:1), PC (25:0/11:3), and PC (15:0/22:6). In addition, TG had a distinct influence on the volatile components and on the oral release of the most important hydrophilic volatile components. The fatty acid species contained in TG were found to affect the oxidation of volatile compounds (Frank et al., 2017). We found that volatile flavor alcohols, aldehydes, and esters were the most abundant in PM, while TG (18:0e/16:0/18:1), TG (18:0/18:0/18:0), TG (18:0e/18:0/18:1), and TG (18:0e/18:1/18:1) were the characteristic lipid molecules in PM. Alcohols have a higher taste threshold and less impact on the taste of meat compared to aldehydes (Huang et al., 2022). Studies have shown that aldehydes, usually derived from C18 PUFA such as linoleic and linolenic acids and C20 arachidonic acid, are the major products of lipid oxidation. Aldehydes have a low olfactory threshold and are important contributors to flavor (Varlet et al., 2007). Aldehydes were the highest percentage of flavor substances in our study and were significantly higher in the PM group than in the LD and BF groups. Therefore, we hypothesized that TG-like lipids are associated with the formation of aldehydes.
Fig. 8. The Venn diagram of differential lipid molecules in different parts of Chaka sheep. PE, phosphatidylethanolamine; PC, phosphatidylcholine; PM, psoas major; LD, longissimus dorsi; BF, biceps femoris; PS, phosphatidylserine; LPC, lysophosphatidylcholine; PE, phosphatidylethanolamine; TG, triacylglycerol; LPE, lysophosphatidylethanolamine.
The cluster analysis and correlation of differential lipid molecules
To evaluate the plausibility of differential lipids and to visualize more clearly the association between specimens and the differential expression patterns of lipids in different samples, we conducted hierarchical clustering (Fig. 9) using the expression levels of highly differential lipids (VIP>1, p-value<0.05) for each group of specimens. Differential lipids were analysed by clustering and heatmaps demonstrated significant differences in lipid profiles between the BF, LD and PM, indicating significant differences in lipid molecules in different parts of the muscle of Chaka sheep.
Fig. 9. The heat map for clustering analysis of differential lipid molecules. BF, biceps femoris; LD, longissimus dorsi; PM, psoas major.
To show the relationship of lipid coregulation more intuitively, the chord diagram (Fig. 10) was utilized to show the correlation between lipid sub-classes, where the beginning point of the link in the inner circle represented the significantly different lipid molecules, and the bows on the outer circles indicated the lipid sub-classes. In BF vs. LD, our study reflected a strong correlation within lipid subclasses PC, PE, and Hex1Cer. By comparing the BF group with the PM group, we obtained some intramolecular correlations of lipid subclasses such as PC, PE, PS, TG, and CL. By comparing PM and LD, our results suggested that there were strong intramolecular correlations between the glycerol ester lipid subclass molecules TG and DG. McMaster’s study showed that the majority of PE is synthesized from ethanolamines, while PE is a choline source. After the synthesis of PE, it can be converted to PC via sequential methylation steps (McMaster, 2018), which further validates the correlation between PC and PE in this study. Zhou’s study found that there was a strong correlation between PS and PI (Zhou et al., 2021). Hou et al. (2020) found a correlation between the expression of PSs and PIs in subcutaneous adipose tissue. However in our study, we did not find a strong correlation between PI and other lipid molecules, the relevance of other lipid molecules needs to be further verified.
Fig. 10. The plot of the differential lipid molecule correlation analysis. BF, biceps femoris; LD, longissimus dorsi; PM, psoas major; PE, phosphatidylethanolamine; PS, phosphatidylserine; TG, triacylglycerol; LPC, lysophosphatidylcholine; LPE, lysophosphatidylethanolamine; CL, cardiolipid; PI, phosphatidylinositol; DG, diacylglycerol; PC, phosphatidylcholine; PE, phosphatidylethanolamine; SM, sphingomyelin.
Conclusion
In our study, we found differences in the composition of nutrients and lipid profile in three parts of the muscles of the Chaka sheep. Specifically, the LD had the highest CP content and the highest levels of aspartic acid and histidine, the BF had the highest contents of glycine and proline, and the PM had the highest abundance of ω-6 PUFA, linoleic acid, and volatile flavor compounds. PLs were the predominant lipids, which accounted for almost half of the total lipids. The lipid profiles were different in the LD, PM and BF of the Chaka sheep. There were obvious intermolecular relationships between PC, PE, Hex1Cer, PS, TG, DG, and CL. This work is aimed at characterizing the nutritional composition and lipid profile of the Chaka sheep. Our findings may open new horizons in the utilization of the Chaka sheep for nutritional or biomedical products, thus enhancing the commercial value of these products.
Supplementary Materials
Supplementary materials are only available online from: https://doi.org/10.5851/kosfa.2024.e47.
Conflicts of Interest
The authors declare no potential conflicts of interest.
Acknowledgements
This research was supported by the Chief Scientist Program of Qinghai Province (2024-SF-102), the Second Comprehensive Scientific Expedition to the Qinghai-Tibet Plateau (2019QZKK040104), and the Joint Special Project of Sanjiangyuan National Park (LHZX-2022-02).
Author Contributions
Conceptualization: Xu X. Data curation: Guo T, Zhang Q. Formal analysis: Wang X, Li N. Methodology: Wang Y, Wei L. Software: Xu X, Liu H. Validation: Hu L. Investigation: Xu S. Writing - original draft: Xu X. Writing - review & editing: Xu X, Guo T, Zhang Q, Liu H, Wang X, Li N, Wang Y, Wei L, Hu L, Xu S.
Ethics Approval
The animal procedures used in this study were performed according to the Guidelines for the Care and Use of Laboratory Animals in China. This study was approved by the Experimental Animal Use Ethics Committee of the Northwest Institute of Plateau Biology, CAS (Permit No. NWIPB20160302).
참고문헌
- Ayed C, Martins SIFS, Williamson AM, Guichard E. 2018. Understanding fat, proteins and saliva impact on aroma release from flavoured ice creams. Food Chem 267:132-139.
- Boselli E, Pacetti D, Curzi F, Frega NG. 2008. Determination of phospholipid molecular species in pork meat by high performance liquid chromatography-tandem mass spectrometry and evaporative light scattering detection. Meat Sci 78:305-313.
- Chen H, Cheng J, Cao X, Hu L, Huang Y, Lan X. 2019a. Detecting copy number variation (CNV) marker related to growth trait of Chaka sheep used in molecular marker assisted selection breeding of Chaka sheep includes amplifying lysine methyltransferase 2D (KMT2D) gene and identifying. Northwest A&F University, Xianyang, China.
- Chen H, Yang Z, Hu L, Lan X, Huang Y, Hu S. 2019b. Detecting CNV marker of Chaka sheep BAG4 gene useful in molecular marker assisted selection breeding of Chaka sheep comprises e.g. using real-time fluorescent quantitative PCR, using genomic DNA of Chaka sheep. Northwest A&F University, Xianyang, China.
- Cohen EH. 1971. Comparison of the official AOAC method with rapid methods for the analysis of moisture in meats. J Assoc Off Anal Chem 54:1432-1435.
- Dai Z, Feng M, Feng C, Zhu H, Chen Z, Guo B, Yan L. 2024. Effects of sex on meat quality traits, amino acid and fatty acid compositions, and plasma metabolome profiles in White King squabs. Poult Sci 103:103524.
- Dannenberger D, Nuernberg G, Nuernberg K, Will K, Schauer N, Schmicke M. 2017. Effects of diets supplemented with n-3 or n-6 PUFA on pig muscle lipid metabolites measured by non-targeted LC-MS lipidomic profiling. J Food Compos Anal 56:47-54.
- Dashdorj D, Amna T, Hwang I. 2015. Influence of specific taste-active components on meat flavor as affected by intrinsic and extrinsic factors: an overview. Eur Food Res Technol 241:157-171.
- Enser M, Scollan ND, Choi NJ, Kurt E, Hallett K, Wood JD. 1999. Effect of dietary lipid on the content of conjugated linoleic acid (CLA) in beef muscle. Anim Sci 69:143-146.
- Fernandez X, Monin G, Talmant A, Mourot J, Lebret B. 1999. Influence of intramuscular fat content on the quality of pig meat: 1. Composition of the lipid fraction and sensory characteristics of m. longissimus lumborum. Meat Sci 53:59-65.
- Frank D, Kaczmarska K, Paterson J, Piyasiri U, Warner R. 2017. Effect of marbling on volatile generation, oral breakdown and in mouth flavor release of grilled beef. Meat Sci 133:61-68.
- Fu Y, Cao S, Yang L, Li Z. 2022. Flavor formation based on lipid in meat and meat products: A review. J Food Biochem 46:e14439.
- Gao H, Zhang Y, Liu K, Fan R, Li Q, Zhou Z. 2022. Dietary sodium butyrate and/or vitamin D3 supplementation alters growth performance, meat quality, chemical composition, and oxidative stability in broilers. Food Chem 390:133138.
- Ge Y, Gai K, Li Z, Chen Y, Wang L, Qi X, Sheng X, Xing K, Wang X, Xiao L, Ni H, Guo Y, Chen L, Sheng X. 2023. HPLC-QTRAP-MS-based metabolomics approach investigates the formation mechanisms of meat quality and flavor of Beijing You chicken. Food Chem X 17:100550.
- Guo D, Wan P, Liu J, Chen DW. 2021. Use of egg yolk phospholipids to boost the generation of the key odorants as well as maintain a lower level of acrylamide for vacuum fried French fries. Food Control 121:107592.
- Guo T, Wang X, Zhang Q, Wei L, Liu H, Zhao N, Hu L, Xu S. 2022. Comparative analysis of the composition of fatty acids and metabolites between Black Tibetan and Chaka sheep on the Qinghai-Tibet Plateau. Animals 12:2745.
- Haag M. 2003. Essential fatty acids and the brain. Can J Psychiat 48:195-203.
- Hazel JR, Williams EE. 1990. The role of alterations in membrane lipid composition in enabling physiological adaptation of organisms to their physical environment. Prog Lipid Res 29:167-227.
- He Y, Zhou M, Xia C, Xia Q, He J, Cao J, Pan D, Sun Y. 2020. Volatile flavor changes responding to heat stress-induced lipid oxidation in duck meat. Anim Sci J 91:e13461.
- Hou B, Zhao Y, He P, Xu C, Ma P, Lam SM, Li B, Gil V, Shui G, Qiang G, Wee Liew C, Du G. 2020. Targeted lipidomics and transcriptomics profiling reveal the heterogeneity of visceral and subcutaneous white adipose tissue. Life Sci 245:117352.
- Huang Q, Dong K, Wang Q, Huang X, Wang G, An F, Luo Z, Luo P. 2022. Changes in volatile flavor of yak meat during oxidation based on multi-omics. Food Chem 371:131103.
- Jia W, Li R, Wu X, Liu S, Shi L. 2021. UHPLC-Q-Orbitrap HRMS-based quantitative lipidomics reveals the chemical changes of phospholipids during thermal processing methods of tan sheep meat. Food Chem 360:130153.
- Kohn TA, Anley MJ, Magwaza SN, Adamson L, Hoffman LC, Brand TS. 2023. Muscle fiber type and metabolic profiles of four muscles from the African black ostrich. Meat Sci 200:109156.
- Li J, Yang Y, Tang C, Yue S, Zhao Q, Li F, Zhang J. 2022a. Changes in lipids and aroma compounds in intramuscular fat from Hu sheep. Food Chem 383:132611.
- Li J, Yang Y, Zhan T, Zhao Q, Zhang J, Ao X, He J, Zhou Q, Tang C. 2021. Effect of slaughter weight on carcass characteristics, meat quality, and lipidomics profiling in longissimus thoracis of finishing pigs. LWT-Food Sci Technol 140:110705.
- Li M, Ren W, Chai W, Zhu M, Man L, Zhan Y, Qin H, Sun M, Liu J, Zhang D, Wang Y, Wang T, Shi X, Wang C. 2022b. Comparing the profiles of raw and cooked donkey meat by metabonomics and lipidomics assessment. Front Nutr 9:851761.
- Liang Y, Jiao D, Du X, Zhou J, Degen AA, Ran F, Sun G, Ji K, Wu X, Cheng X, Ma X, Qian C, Yang G. 2023. Effect of dietary Agriophyllum squarrosum on average daily gain, meat quality and muscle fatty acids in growing tan lambs. Meat Sci 201:109195.
- Lin Y, Wang H, Rao W, Cui Y, Dai Z, Shen Q. 2019. Structural characteristics of dietary fiber (Vigna radiata L. hull) and its inhibitory effect on phospholipid digestion as an additive in fish floss. Food Control 98:74-81.
- Liu Z, Zhao M, Wang X, Li C, Liu Z, Shen X, Zhou D. 2022. Investigation of oyster Crassostrea gigas lipid profile from three sea areas of China based on non-targeted lipidomics for their geographic region traceability. Food Chem 386:132748.
- Lugo Charriez K, Soledade Lemos L, Carrazana Y, Rodriguez-Casariego JA, Eirin-Lopez JM, Hauser-Davis RA, Gardinali P, Quinete N. 2021. Application of an improved chloroform-free lipid extraction method to staghorn coral (Acropora cervicornis) lipidomics assessments. Bull Environ Contam Toxicol 107:92-99.
- Ma C, Li X, Zheng C, Zhou B, Xu C, Xia T. 2021. Comparison of characteristic components in tea-leaves fermented by Aspergillus pallidofulvus PT-3, Aspergillus sesamicola PT-4 and Penicillium manginii PT-5 using LC-MS metabolomics and HPLC analysis. Food Chem 350:129228.
- McMaster CR. 2018. From yeast to humans: Roles of the Kennedy pathway for phosphatidylcholine synthesis. FEBS Lett 592:1256-1272.
- Mohanty TR, Park KM, Pramod AB, Kim JH, Choe HS, Hwang IH. 2010. Molecular and biological factors affecting skeletal muscle cells after slaughtering and their impact on meat quality: A mini-review. J Muscle Foods 21:51-78.
- Oh R. 2005. Practical applications of fish oil (Ω-3 fatty acids) in primary care. J Am Board Fam Pract 18:28-36.
- Qin X, Zhang T, Cao Y, Deng B, Zhang J, Zhao J. 2020. Effects of dietary sea buckthorn pomace supplementation on skeletal muscle mass and meat quality in lambs. Meat Sci 166:108141.
- Ramos LMG, Bezerra LR, de Oliveira JPF, de Souza MP, da Silva AL, Pereira ES, Mazzetto SE, Pereira Filho JM, Oliveira RL. 2021. Effects of feeding growing-finishing lambs with cashew nut shell liquid on the growth performance, physicochemical attributes, lipid peroxidation and sensorial parameters of burger. Small Rumin Res 202:106468.
- Ribeiro de Araujo Cordeiro AR, de Medeiros LL, Alencar Bezerra TK, Bertoldo Pacheco MT, Galvao MS, Madruga MS. 2020. Effects of thermal processing on the flavor molecules of goat by-product hydrolysates. Food Res Int 138:109758.
- Roy BC, Coleman P, Markowsky M, Wang K, She Y, Richard C, Proctor SD, Bruce HL. 2024. Muscle fiber, connective tissue and meat quality characteristics of pork from low birth weight pigs as affected by diet-induced increased fat absorption and preferential muscle marbling. Food Sci Anim Resour 44:51-73.
- Schiavon S, Pellattiero E, Cecchinato A, Tagliapietra F, Dannenberger D, Nuernberg K, Bittante G. 2016. The influence of different sample preparation procedures on the determination of fatty acid profiles of beef subcutaneous fat, liver and muscle by gas chromatography. J Food Compos Anal 50:10-18.
- Simopoulos AP. 2008. The importance of the omega-6/omega-3 fatty acid ratio in cardiovascular disease and other chronic diseases. Exp Biol Med 233:674-688.
- Sun N, Chen J, Wang D, Lin S. 2018. Advance in food-derived phospholipids: Sources, molecular species and structure as well as their biological activities. Trends Food Sci Technol 80:199-211.
- Tang J, Zhang B, Liu D, Gao K, Dai Y, Liang S, Cai W, Li Z, Guo Z, Hu J, Zhou Z, Xie M, Hou S. 2023. Dietary riboflavin supplementation improves meat quality, antioxidant capacity, fatty acid composition, lipidomic, volatilomic, and proteomic profiles of breast muscle in Pekin ducks. Food Chem X 19:100799.
- Tian Y, Xu T, Li Y, Liu Y, Liu J. 2021. An untargeted LC-MS metabolomics approach to the metabolic profiles of bottom cultured scallops (Mizuhopecten yessoensis) subjected to mechanical shock in early post-harvest handling. Aquaculture 533:736061.
- Vahmani P, Ponnampalam EN, Kraft J, Mapiye C, Bermingham EN, Watkins PJ, Proctor SD, Dugan MER. 2020. Bioactivity and health effects of ruminant meat lipids. Invited review. Meat Sci 165:108114.
- Varlet V, Prost C, Serot T. 2007. Volatile aldehydes in smoked fish: Analysis methods, occurrence and mechanisms of formation. Food Chem 105:1536-1556.
- Wang DY, Zhu YZ, Xu WM. 2009. Comparative study of intramuscular phospholipid molecular species in traditional Chinese duck meat products. Asian-Australas J Anim Sci 22:1441-1446.
- Xu X, Liu H, Wang X, Zhang Q, Guo T, Hu L, Xu S. 2023. Evaluation of the Longissimus thoracis et Lumborum muscle quality of Chaka and Tibetan sheep and the analysis of possible mechanisms regulating meat quality. Animals 13:2494.
- Yang Z, Hou Y, Zhang M, Hou P, Liu C, Dou L, Chen X, Zhao L, Su L, Jin Y. 2024. Unraveling proteome changes of Sunit lamb meat in different feeding regimes and its relationship to flavor analyzed by TMT-labeled quantitative proteomic. Food Chem 437:137657.
- Yin M, Chen M, Matsuoka R, Song X, Xi Y, Zhang L, Wang X. 2023. UHPLC-Q-Exactive Orbitrap MS/MS based untargeted lipidomics reveals fatty acids and lipids profiles in different parts of capelin (Mallotus villosus). J Food Compos Anal 116:105096.
- Zhang X, Han L, Hou S, Raza SHA, Wang Z, Yang B, Sun S, Ding B, Gui L, Simal-Gandara J, Shukry M, Sayed SM, Al Hazani TMI. 2022a. Effects of different feeding regimes on muscle metabolism and its association with meat quality of Tibetan sheep. Food Chem 374:131611.
- Zhang Y, Rui X, Vaugeois R, Simpson BK. 2022b. Seal meat enzymatic hydrolysates and its digests: A comparison on protein and minerals profiles. LWT-Food Sci Technol 157:113072.
- Zhang Y, Zhang J, Gong H, Cui L, Zhang W, Ma J, Chen C, Ai H, Xiao S, Huang L, Yang B. 2019. Genetic correlation of fatty acid composition with growth, carcass, fat deposition and meat quality traits based on GWAS data in six pig populations. Meat Sci 150:47-55.
- Zhou J, Zhang Y, Wu J, Qiao M, Xu Z, Peng X, Mei S. 2021. Proteomic and lipidomic analyses reveal saturated fatty acids, phosphatidylinositol, phosphatidylserine, and associated proteins contributing to intramuscular fat deposition. J Proteomics 241:104235.
- Zhou L, Khalil A, Bindler F, Zhao M, Marcic C, Ennahar S, Marchioni E. 2013. Effect of heat treatment on the content of individual phospholipids in coffee beans. Food Chem 141:3846-3850.
- Zhou L, Zhao M, Bindler F, Marchioni E. 2014. Comparison of the volatiles formed by oxidation of phosphatidylcholine to triglyceride in model systems. J Agric Food Chem 62:8295-8301.
- Zuo H, Han L, Yu Q, Niu K, Zhao S, Shi H. 2016. Proteome changes on water-holding capacity of yak longissimus lumborum during postmortem aging. Meat Sci 121:409-419.