Principal component analysis
From the principal component analysis in Fig 1, it can be seen that there is a certain grouping trend of subcutaneous fat under different feeding conditions and the difference between the G24_SF and F24_SF groups is large and the difference with the G18_SF group is small, in which the QC samples are clustered together, which indicates that the samples have good reproducibility, the quality of the QC is high and the data are reliable and can be subjected to further analysis of variance.
Screening for differential metabolites
Venn diagrams show the metabolites that are shared or unique between metabolic sets. To obtain the differential metabolites between metabolic sets, differential metabolite analysis was performed. As shown in Fig 3, there were 97 differential metabolites between the G24_SF and G18_SF sets and 110 differential metabolites between the F24_SF and G24_SF sets (Fig 2).
Bioinformatics analysis
Classification of KEGG compounds
The analyzed metabolites were compared to the KEGG Compound database to obtain an overview of metabolite classification, as shown in Fig 3a and it was found that the metabolites differing between the F24_SF and G24_SF groups were classified into a total of seven major types, including vitamins and cofactors, peptides, nucleic acids; the different metabolites produced were mainly peptides and lipid compounds. It was found that the differential metabolites of the G18_SF and G24_SF groups were categorized into five major types, including peptides, nucleic acids, lipids, hormones and transmitters and carbohydrates; the differential metabolites produced were mainly peptides and carbohydrates (Fig 3b).
Functional analysis of the KEGG database
The experiment searched the KEGG database for the differential metabolites of each sample group and obtained 31 and 18 signaling pathways belonging to 6 and 5 classes, respectively and the results showed (Fig 4) that the differential metabolites of each sample group were significantly annotated to the metabolism of amino acid metabolism, carbohydrate metabolism and lipid metabolism, among others; Digestive system and Nervous system in the system of organismal systems,
etc.; Cancer: overview in the human disease,
etc.; Membrane transport in the environmental information processing,
etc.
Hierarchical cluster analysis of differential metabolites
Cluster analysis can provide a better understanding of the relationship between the differential metabolites screened above. Fig 5 Horizontal coordinate is the sample name, vertical coordinate is the differential metabolite, red color indicates high expression of the substance content and blue color indicates low expression of the substance content. As shown in Fig 5a, the G24_SF group showed high expression of phosphatidyl ethanolamine (PE) and phosphatidylcholine (PC) related metabolites and N-acetylneuraminic acid and low expression of LysoPE, LysoPC, Beta-D-Glucosamine and Methycholine. Low expression of docosapentaenoic acid, ophthalmic acid, N-acetylneuraminic acid and PE and PC related metabolites and high expression of LysoPE, LysoPC, Beta-D-Glucosamine and Methycholine in the F24_SF group.As can be seen in Fig 5b, in the G24_SF group, the expression of PE and PC related metabolites and LysoPE-related metabolites was low and the expression of Palmitoyl-L-carnitine and 2-Pentyl-3-phenyl-2-propenal was high. The expression of PE and PC related metabolites, LysoPE-related metabolism and Beta-D-Glucosamine was high in the G18_SF group;low expression of Palmitoyl-L-carnitine and 2-pentyl-3-phenyl-2-propene.
Volcanic composition
Each dot in Fig 6 represents a specific metabolite. The dots on the left are metabolites that are differentially down-regulated in expression and the dots on the right are metabolites that are differentially up-regulated in expression. In Fig 6a, 83 metabolites were down-regulated and 27 metabolites were up-regulated in the G24_SF relative to the F24_SF group and in Fig 6b, 41 metabolites were down-regulated and 56 metabolites were up-regulated in the G24_SF relative to the G18_SF group. Table 5 and Table 6 show that products with annotations were detected in the G24_SFvsF24_SF and G24_SFvsG18_SF groups, respectively.
KEGG pathway enrichment
By comparing the F24_SF with the G24_SF group (Fig 7a and Table 7), KEGG pathway enrichment was found to involve metabolic pathways such as Alanine, aspartate and glutamate metabolism, Arachidonic acid metabolism, D-Glutamine and D-glutamate metabolism and others;By comparing the G18_SF with the G24_SF group (Fig 7b and Table 8), KEGG pathway enrichment was found to involve metabolic pathways such as Alanine, aspartate and glutamate metabolism, Caffeine metabolism, Starch and sucrose metabolism and others.
In this study, yaks were grazed and housed during the dry grass period and yak subcutaneous fat was used as the object of study to investigate the effects of different feeding methods and different ages on yak fat metabolism by analyzing the metabolomic changes of yak subcutaneous fat at the ages of 18 months and 24 months.
Differential metabolite analysis
2-Pentyl-3-phenyl-2-propenal, N-Acetylneuraminic acid and Xanthine were the most abundant in F24_SF, followed by G24_SF and the least abundant in G18_SF.N-acetylneuraminic acid (Neu5Ac), also known as sialic acid, is an acetylated derivative of the acidic sugar neuraminic acid, Neu5Ac enhances intestinal absorption of minerals and vitamins and promotes bone development; and Neu5Ac has also been associated with neurotransmission, promotion of brain development, cosmetic whitening, promotion of the growth of probiotics, anti-tumor, anti-viral or bacterial infections and anti-Alzheimer’s disease, among other functions
(Bavaresco et al., 2008).Studies have shown that Neu5Ac enhances immunity in mice. Neu5Ac prevents HFD-induced hyperlipidemia in rats by regulating lipid-associated and coagulation-associated genes, which in turn induce changes in metabolites and proteins
(Yida et al., 2015). Hypoxanthine, also known as 6-hydroxypurine, which can be obtained biologically by the removal of one amino group from adenine, is an important alkaloidal purine that is widely distributed in the body and is involved in the regulation of the body’s physiological functions. Hypoxanthine may be used to assess health status (
Mahanty and Xi, 2020). Arginine is an essential amino acid for young animals, which can regulate protein synthesis and metabolism, promote growth and development, improve intestinal development and the number of intestinal villi and increase the immunity of animals (
Nieves and Langkamp-Henken, 2002;
Peranzoni et al., 2007). It was shown that feeding 1.15% L-arginine to ZDF rats significantly reduced abdominal and perirenal fat in rats (
Clapés and Rosa Infante, 2009;
Fu et al., 2005). Therefore, it is hypothesized that the nutritional value of yaks housed at 24 months of age may be somewhat higher; and the above differential metabolites may affect fat deposition in yaks, which can be promoted later by regulating them.
Dimethylglycine, niacinamide and glutamine were lower than in the housed group. Dimethylglycine (DMG) is a chemical complex that appears in the pathway of choline to glycine metabolism and has the emulsifying effect of a surfactant, which can effectively emulsify oil-water mixtures, break down fat particles in the mixtures, promote the chemical action of digestive enzymes in the intestinal tract and increase the contact area between the villi of the small intestine and the nutrients to promote the digestive absorption of the nutrients. Nicotinamide (NAM) is the amide form of vitamin B3.NAM can regulate fat metabolism, enhance body immunity and promote fat deposition in lambs; it was found that peri rumen nicotinamide can promote the growth and development of fattening lake sheep and increase the average daily weight gain; and Li suggested that NAM, by inhibiting the activity of SIRT1 dehydrogenase to promote lipid synthesis in animals
(Li et al., 2015). Glutamine (GIn) is a non-essential amino acid, on the one hand, it can promote the proliferation of intestinal cells and intestinal development, enhance the activity of small intestinal enzymes and promote the absorption of nutrients; on the other hand, it can promote the proliferation of immune cells and improve the immunity of the body.It was found that the addition of 0.8% GIn significantly increased the mean daily feed intake and mean daily weight gain of 1- and 4-week-old chicks; the addition of GIn to the diet increased the mean daily weight gain and apparent digestibility of weaned piglets
(Hwang et al., 2015). The above DMG, hypoxanthine, NAM, GIn and GIcN all promoted fat deposition, which was lower in the G24_SF group, presumably because yaks are in a cold environment and the body must provide sufficient energy to resist the cold, resulting in excessively low fat deposition compared to housed yaks.DMG, NAM, arginine, GIn, GIcN
(Shmagel et al., 2019; Bowman and Wolfgang, 2019), PRO, CAR
(Salem et al., 2015; Brown et al., 2014; Bao et al., 2015), AA
(Demetz et al., 2014; Zheng and Hasegawa, 2016) and Neu5Ac all improved the immunity of yaks and were lower in the G24_SF group. Lower temperatures in the cold season on the Tibetan Plateau, a decrease in drinkable water resources and an uneven distribution of pasture grasses resulted in a lower nutrient intake by yaks during the cold season, which reduced the immunity of yaks. Therefore, grazing may affect the expression levels of the above differential metabolites, which in turn affects fat deposition and immunity in yaks; secondly, the nutritional value of housed yak meat is a little bit better compared to grazed yak meat.
KEGG enrichment pathway analysis
The Alanine, aspartate and glutamate metabolism pathways are thanked as significantly different metabolic pathways common to G24_SF, F24_SF and G24_SF. The metabolites of alanine (Ala) can provide the body with carbon shelf, nitrogen source and energy. Aspartate (Asp) not only regulates the immune system, but also plays an important role in the proliferation of immune cells. Glutamate (Glu) regulates normal physiological activities such as cell proliferation, differentiation, secretion and contraction through glutamate receptors (
Skerry and Genever 2001;
Newsholme et al., 2003). Glutamine is synthesized from glutamate, NH
3 and ATP in skeletal muscle by the action of glutamine synthetase. It has been found that glutamine promotes the cytokine TNF to enhance immunity (
Calder and Yaqoob 1999). Differential metabolites enriched to this pathway may be related to fat deposition, immunocompetence and growth and development in yaks and are important for resolving the mechanisms of yak adipocyte proliferation and differentiation as well as immunoregulation.
In summary, these metabolites regulate fat metabolism and may promote fat deposition by inhibiting adipocyte proliferation, promoting intestinal development and increasing nutrient absorption. Therefore, different feeding practices as well as different months of age will affect fat deposition in yaks and these significantly enriched pathways and differential metabolites will be the main metabolites and metabolic pathways that are highly affected during yak growth and development.