Next Article in Journal
Functional Anatomy of the Thoracic Limb of the Komodo Dragon (Varanus komodoensis)
Next Article in Special Issue
Improving the Aerobic Capacity in Fingerlings of European Sea Bass (Dicentrarchus labrax) through Moderate and Sustained Exercise: A Metabolic Approach
Previous Article in Journal
Effects of Concentrate Feeding Sequence on Growth Performance, Nutrient Digestibility, VFA Production, and Fecal Microbiota of Weaned Donkeys
Previous Article in Special Issue
Metabolic Rates of Japanese Anchovy (Engraulis japonicus) during Early Development Using a Novel Modified Respirometry Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Metabolome and Transcriptome Analyses Reveal the Efficacy of Steroidal Saponins for Glucose and Lipid Metabolism in Hybrid Grouper (♀Epinephelus fuscoguttatus × ♂Epinephelus lanceolatu) Fed Higher-Lipid Diets

1
Laboratory of Aquatic Nutrition and Feed, College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China
2
Guangzhou Fishtech Biotechnology Co., Ltd., Guangzhou 510640, China
3
Guangdong Engineering Technology Research Center of Aquatic Animals Precision Nutrition and High Efficiency Feed, Zhanjiang 524088, China
*
Author to whom correspondence should be addressed.
Animals 2023, 13(18), 2894; https://doi.org/10.3390/ani13182894
Submission received: 24 July 2023 / Revised: 31 August 2023 / Accepted: 5 September 2023 / Published: 12 September 2023
(This article belongs to the Special Issue Metabolic Adaptation and Regulation in Aquatic Animals)

Abstract

:

Simple Summary

The ingestion of higher-lipid diets by fish could lead to excessive lipid deposition and metabolic disorders in vivo. This experiment analyzed the effects of steroidal saponins on hepatic glucose and lipid metabolism in hybrid grouper fed higher-lipid diets using a multi-omics approach. Data from grouper fed diets with 0% and 0.1% steroidal saponins show that 0.1% steroidal saponins inhibited glycogen synthesis, gluconeogenesis, and the arachidonic acid metabolism pathway and activated glycogenolysis, glycolysis, and the fatty acid β-oxidation pathway. This study provides a deeper insight into the glycolipid metabolic processes in the liver of grouper fed higher crude-lipid diets.

Abstract

An analysis of the extent of the effect of steroidal saponin addition on glucose and lipid metabolism in hybrid grouper liver was performed at the transcriptomic and metabolomic levels. Feeds (52% crude protein, 14% crude lipid) were prepared containing 0% (S0), 0.1% (S0.1), and 0.2% (S0.2) steroidal saponins. After eight weeks of feeding trial, compared to the S0 group, the activities of serum albumin, alanine aminotransferase, and aspartate transaminase were significantly lower and the activities of lysozyme, acid phosphatase, and alkaline phosphatase were significantly higher in the S0.1 group (p < 0.05). The superoxide dismutase, catalase, and glutathione peroxidase activities in the livers of the S0.1 group were significantly higher than those of the S0 group, while the malondialdehyde content was significantly lower than that of the S0 group (p < 0.05). There were forty-two differentially expressed genes and thirty-two differential metabolites associated with glucose and lipid metabolism enriched using KEGG and GO. In the S0 group, the expression of prostaglandin-endoperoxide synthase 1, prostaglandin E synthase 1, and thromboxane-2 synthase mRNA was significantly higher than in the S0.1 group (p < 0.05). The expression levels of genes in the S0 group were significantly higher than those in the S0.1 group (p < 0.05), including for glycogen synthase kinase, glucose-6-phosphatase catalytic subunit 2, fructose-1,6-bisphosphatase, phosphoenolpyruvate carboxykinase, glucose transporter 4, and malate dehydrogenase. The expression of mRNA such as fatty acid synthase, acetyl-CoA carboxylase, and sterol regulatory element-binding protein 1 was significantly lower in the S0.1 group than in the S0 group, while the expression of carnitine acyltransferase 1, acyl-CoA synthetase, and acyl-CoA dehydrogenase genes was significantly higher in the S0 group (p < 0.05). In summary, glycogen synthesis, gluconeogenesis, and the arachidonic acid metabolism pathway were inhibited by 0.1% steroidal saponins, and glycogenolysis, glycolysis, the tricarboxylic acid cycle, and the fatty acid β-oxidation pathway were activated. This study aims to provide a reference for the formulation of grouper feeds with a higher crude-lipid level.

1. Introduction

In the past five years, China’s cultivated grouper production grew from 131,500 tons to 204,100 tons, an increase of 55.21% [1,2]. The development of compound diet technology is a key factor in ensuring grouper production. Further, fishmeal, which is an optimal protein source for grouper diets, is rich and balanced in nutrients, especially in protein and fatty acids [3]. Between 2017 and 2021, Chinese fishmeal production rose from 639,200 tons to 659,000 tons. However, in fact, the annual demand for fishmeal in China was about 1.7 million tons, which meant that the Chinese fishmeal imports had increased from 1,571,600 tons to 1,823,400 tons with an external dependence of over 70% [1,2,4,5]. In view of this, fishmeal influenced the formula composition and price of feed and limited the efficient development of aquaculture.
Appropriately increasing lipid content in diets is an effective measure to improve protein efficiency and feed utilization [6,7,8,9], but long-term intake of high-lipid diets can result in abnormal accumulation of hepatic lipids and induce metabolic diseases such as fatty liver, which leads to metabolic disorders and immune suppression in the organism, affecting fish growth and reducing the economic efficiency of culturing [10,11].
The Panax notoginseng saponin can enhance the antioxidant capacity of the organism by increasing superoxide dismutase and catalase enzyme activities and decreasing malondialdehyde content in the liver of grouper [12]. Yucca schidigera saponins and soya saponins were found to reduce malondialdehyde, total cholesterol, and triglyceride levels in the serum and liver of white shrimp (Litopenaeus vannamei) and turbot (Scophthalmus maximus L.) [13,14]. The steroid saponins belong to a steroidal class of saponins, including spirostanes, furanosteroids, and furanospiralosteroids, which have the ability to scavenge free radicals, such as -OH and O2− [15], to lower serum and hepatic levels of cholesterol and triglycerides [16] as well as the ability to promote the breakdown of hepatic glycogen [17]. The bitter melon saponins (Momordica charantia, steroids) have been shown to enhance glucose metabolism by preventing the absorption of glucose in the small intestine of mice and promoting the breakdown of glycogen in the liver [18,19]. Saikosaponin d has been confirmed to alleviate high-fat-diet-induced hepatic steatosis in hybrid grouper by targeting the AMPK/PPARα pathway [20]. Our previous study found that steroidal saponins improve the utilization of lipids and glucose in grouper liver when the fish ingested a high-lipid diet [21]. What is the exact mechanism?
The transcriptome is the pivot that links the genetic information of genes to functional proteins [22]. For example, it can explain that taurine can promote hepatic endogenous taurine synthesis, bile secretion, insulin secretion, and fatty acid β-oxidation efficiency to reduce fat accumulation in grouper ingesting high-lipid diets (15% crude lipid) [23]. As a bridge between genotype and phenotype, the metabolome ensures the maintenance of aerobic metabolism and activates gluconeogenesis as key metabolic strategies of turbot kidney against heat stress [24]. Association analysis of transcriptomic and metabolomic data provides a more detailed analysis of the relationship between critical genes and metabolites to explore the potential mechanisms of the organism [25]. The down-regulated key genes (lysophosphatidylcholine acyltransferase and glycerol-3-phosphate dehydrogenase) and increased metabolites (choline and taurine) were found to contribute to the alleviation of lipid metabolism disorders caused by elevated temperature stress in turbot [26]. In this study, we combine transcriptome and metabolome association analysis to investigate how steroidal saponins alter hepatic glycolipid metabolism in grouper.

2. Materials and Methods

2.1. Experimental Design

Three isoproteic (52% crude protein) and isolipidic (14% crude lipid) diets were formed by adding 0% (S0), 0.1% (S0.1), and 0.2% (S0.2) steroidal saponins (Guangzhou Fishtech Biotechnology Co., Ltd., Guangzhou, China) to diets with fishmeal, concentrated cottonseed protein, poultry by-product meal, and soybean meal as the main protein source and fish oil and soybean oil as the main lipid sources. Please refer to Supplementary Material (Tables S1 and S2). Experimental fish were fed the experimental diets to apparent satiation at 8:00 and 16:00 daily for 8 weeks. The hybrid grouper were separated into nine buckets (1 m3 Fiberglass farming buckets) with twenty-five fish in each bucket. During the experiment, the fish were continuously oxygenated every day, dissolved oxygen 5–6 mg/L, ammonia nitrogen <0.2 mg/L, temperature 30.5 ± 0.8 °C, salinity 28–32. Please refer to our previous study for the detailed feeding process [21].

2.2. Sample Collection and Chemical Analysis

Blood samples were obtained from the tail vein of hybrid grouper. The samples were centrifuged at 4 °C and 3500 rpm for 15 min to extract the serum and then immediately stored at −80 °C.
Meanwhile, the soybean-sized liver (approximately 600 mg) was obtained from each replicate on ice and then washed with saline. At once, the liver samples were stored at −80 °C for analysis of the gene expression of glucose lipid metabolism enzyme.
The grouper were anesthetized with MS-222 (1:10,000), and then the fresh soybean-sized liver was cut from living fish fed diets S0 and S0.1 and immediately washed with saline to clean the surface of blood remaining on the tissue [21]. The liver samples were numbered, stored in 2 mL freeze tubes, quickly put in liquid nitrogen, and subsequently transferred to −80 °C for metabolome and transcriptome analyses. The detailed methodology was described in Yang and Li et al. [9,13].
Serum parameters were measured using commercial kits (Nanjing Jiancheng Institute of Biological Engineering Co., Ltd., Nanjing, China). The parameters included total protein (TP, microplate method, A045-4-1), albumin (ALB, microplate method, A028-2-1), alanine aminotransferase (ALT, microplate method, C009-2-1), aspartate transaminase (AST, microplate method, C010-2-1), lysozyme (LZM, turbidimetric method, A050-1-1), acid phosphatase (ACP, microplate method, A060-2-2), and alkaline phosphatase (AKP, microplate method, A059-2-2).
Antioxidative enzymes in liver including superoxide dismutase (SOD, hydroxylamine method, A001-2-2), catalase (CAT, ammonium molybdate method, A007-1-1), glutathione peroxidase (GSH-PX, colorimetric method, A005-1-2), and hepatic malondialdehyde (MDA, TBA method, A003-1-2) contents were analyzed with commercial kits (Nanjing Jiancheng Institute of Biological Engineering Co., Ltd., Nanjing, China).

2.3. Metabolome Analysis

The 100 mg liver samples (six replicates per group) were ground separately in liquid nitrogen, and 500 μL of precooled 80% methanol was resuspended using well vortex to homogenize. The liver samples were ice-bathed for 5 min, followed by centrifugation at 15,000× g for 20 min at 4 °C. The supernatants were diluted with mass-spectrometry-grade water to a methanol content of 53% and then centrifuged (15,000× g, 20 min, 4 °C). Thus, the supernatants were injected into the UPLC–MS/MS system for analysis. UPLC–MS/MS analyses were conducted at Gene Denovo Ltd. (Guangzhou, China), and a Vanquish UPLC system (ThermoFisher, Berlin, Germany) and an Orbitrap Q ExactiveTM HF-X mass spectrometer (Thermo Fisher, Berlin, Germany) were used.
The raw data files were processed using UPLC–MS/MS with Compound Discoverer 3.1 (CD3.1, Thermo Fisher) performing peak alignment, peak extraction, and quantification for each metabolite. Peak was matched to mzCloud (https://www.mzcloud.org/, accessed on 16 February 2022), mz Vault, and Mass List database to obtain accurate qualitative and relative quantitative results. We used the statistical software R (R-3.4.3 version), Python (Python version 2.7.6), and CentOS (CentOS version 6.6) for statistical analysis. A normal transformation was attempted using the area normalization method when the data were not normally distributed.

2.4. Transcriptome Analysis

The total RNA of liver (about 300 mg) was extracted with Trizol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s protocol. The RNA quality integrity and concentration were assessed with agarose gel electrophoresis and Nanodrop 2000 assay. The input concentration of RNA sample per replicate in S0 group was 765.8, 576.9, and 406.3 ng/uL; in S0.1 group, 764.7, 537.0, and 792.5 ng/uL; in S0.2 group, 808.7, 675.4, and 855.0 ng/uL. The enriched mRNA was then fragmented with fragmentation buffer and reversed to cDNA using the NEBNext Ultra RNA Library Prep Kit for Illumina (NEB #7530, New England Biolabs, Ipswich, MA, USA). The cDNA libraries were sequenced with Gene Denovo Biotechnology (Guangzhou, China) using Illumina Novaseq6000.
The raw reads were quality controlled using fastp (https://github.com/OpenGene/fastp, accessed on 14 February 2022) software with default parameters to obtain clean reads. All assembled sequences were matched to the protein sequence databases Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.jp/, accessed on 14 February 2022) and Gene and Genomes Database (GO, http://geneontology.org, accessed on 14 February 2022) with blastx (E-value cut-off of 8 × 105). The following differential expression analysis and functional enrichment were performed on the Omicsmart cloud platform (https://www.omicsmart.com, accessed on 14 February 2022). The expression of each unigene was calculated using the RSEM program (http://deweylab.github.io/RSEM/, accessed on 14 February 2022). Differentially expressed unigenes (DEGs) were analyzed using DESeq2 software with |log2FC| > 1 and p value ≤ 0.05 (http://www.bioconductor.org/packages/release/bioc/html/DESeq.html, accessed on 14 February 2022).

2.5. Quantitative RT-PCR Analysis of Gene

Quantitative reverse transcription PCR (qRT-PCR) was performed to verify the relative expression levels of the fourteen selected DEGs so that the gene expression data obtained with RNA-Seq could be validated. Liver total RNA was extracted with Trizol reagent (Invitrogen, Carlsbad, CA, USA). The cDNA was composed with the Prime Script RT kit (Takara, Osaka, Japan), and the quantitative thermal cycler (Light Cycler 480II, Roche Diagnostics, Basel, Switzerland) was used for qRT-PCR with the SYBR Premix Ex Taq kit (Takara, Osaka, Japan) within a reaction volume of 10 μL, including 3.2 μL of sterilized double-distilled water, 1 μL of cDNA, 0.4 μL of each primer, and 5 μL of SYBR Premix Ex Taq (Takara, Osaka-shi, Japan). Requirements were a cycle of 30 s at 95 °C followed by thirty-five cycles at 95 °C for 5 s, 60 °C for 25 s, and 72 °C for 30 s. The β-actin (AY510710.2) was used as internal control. The results were calculated using the 2−ΔΔCt method [12]. A list of primers is given in Table 1.

2.6. Statistical Analysis

Mean ± standard error (SE) was used to present the results of the experiment. The statistic was analyzed using one-way analysis of variance (ANOVA) followed by Tukey’s multiple comparison tests with the software SPSS 21.0 (SPSS Inc., Michigan Avenue, Chicago, IL, USA). Significant differences were identified as p < 0.05. The Graph Pad Prism 8.0 software (8.0.2.263) was used to image experimental results.

3. Results

3.1. Activities of ALT, AST, and Levels of TP, ALB in Serum

The results of the serum biochemical indexes are shown in Table 2. The enzyme activity of ALT between the S0.1 and S0.2 groups was not significantly different, but both groups’ ALT activity was significantly lower than that in the S0 group (p < 0.05). The AST activity in the S0.1 group was significantly lower than that in the other two groups (p < 0.05). There was no significant difference in the TP content of serum among the three groups (p > 0.05). The ALB level in the S0.2 group was significantly higher than that in the S0.1 group but significantly lower than that in the S0 group (p < 0.05).

3.2. Levels of LZM, ACP, and AKP in Serum

The serum nonspecial immune indexes of the hybrid grouper are shown in Table 3. The enzyme activities of LZM, ACP, and AKP in the S0 group were not significantly different from those in the S0.2 group, while they were significantly lower than those in the S0.1 group (p < 0.05).

3.3. Liver Antioxidative Indexes

As shown in Table 4, the content of MDA in the S0 group was significantly higher than that in the experiment groups (p < 0.05). The SOD and CAT enzyme activities in the S0.1 group were significantly higher than those in the other groups (p < 0.05). The enzyme activity of GSH-PX in the S0.2 group was significantly higher than that in the S0 group but significantly lower than that in the S0.1 group (p < 0.05).

3.4. Transcriptome Analysis in S0 and S0.1 Group Fish Liver

To explore the potential mechanisms of steroidal saponins on glucose lipid metabolism in hybrid grouper liver, first, the high-throughput RNA-sequencing technique was used in this experiment. The goodness of reproducibility between samples within groups and the degree of variability of samples between groups were analyzed with PCA (Figure 1A). Then, differentially expressed genes (DEGs) were separated into the S0 group and the best-addition S0.1 group of steroidal saponins. A total of 5667 DEGs were identified under the log2|Fold change| > 2.0 threshold and with a p value < 0.05 (Figure 1C), including 2687 up-regulated DEGs and 2980 down-regulated DEGs (Figure 1B).

3.5. Signaling Pathway Network and Key Genes of Glucose and Lipid Metabolism

The DEGs were entered into the GO and KEGG databases to search their involved signaling pathways (Figure 2). The present experiment involved 46 key different genes (Table 5), including 18 down-regulated and 28 up-regulated key genes. Of these, 15 genes were significantly down-regulated, and 24 genes were significantly up-regulated. Then, 30 DEGs were selected, as shown in Table 5, and clustered into a heatmap depending on their expression level and classification (Figure 2C).
The interesting DEGs were enriched in GO, focusing on the thirty items in this study. To be exact, among the 30 enrichment items analyzed using the rich factor value, 24 items were associated with glucose metabolism, including glycogen metabolism (positive regulation of glycogen mediated process, positive regulation of glycogen biosynthetic process, UDP-glucose-4-epimerase activity, and UDP-glucose: glycoprotein glycosyltransferase activity), glycolysis (glucose catabolic process, glucose catabolic process to pyruvate, glucose catabolic process to lactate, and glucose catabolic process to lactate via pyruvate), glycogenesis (glycolytic process through glucose-6-phosphate, glucose-6-phosphate metabolic process, and DTDP-glucose-4,6-dehydratase activity), tricarboxylic acid cycle (glucose-6-phosphate isomerase activity), glucose metabolic process, regulation of glucose-mediated signaling pathway, invasive growth in response to glucose limitation, negative regulation of glucose-mediated signaling pathway, glucose-mediated signaling pathway, glucose homeostasis, regulation of invasive growth in response to glucose limitation, response to glucose, response to oxygen glucose deprivation, cellular response to oxygen glucose deprivation, cellular glucose homeostasis, and cellular response to glucose stimulus (Figure 2A). Meanwhile, six items were related to lipid metabolism, including phospholipid-hydroperoxide glutathione peroxidase activity, regulation of lipid storage, calcium-dependent phospholipid binding, positive regulation of sphingolipid biosynthetic process, fatty acid binding, and protein–lipid complex assembly (Figure 2A).
In addition, the focused DEGs were enriched into 30 interesting KEGG signaling pathways. The 30 signaling pathways showed that glucose metabolism (glycolysis/gluconeogenesis, oxidative phosphorylation, carbon metabolism, pentose phosphate pathway, fructose and mannose metabolism, glyoxylate and dicarboxylate metabolism, and thermogenesis), amino acids metabolism (biosynthesis of amino acids, proteasome, protein processing in endoplasmic reticulum, glutathione metabolism, cysteine and methionine metabolism, and protein digestion and absorption), lipid metabolism (metabolism of xenobiotics by cytochrome P450 and drug metabolism by cytochrome P450), and other signaling pathways (antigen processing and presentation, cardiac muscle contraction, apoptosis, estrogen signaling pathway, HIF-1 signaling pathway, drug metabolism by other enzymes, nucleocytoplasmic transport, phagosome, regulation of actin cytoskeleton, necroptosis, ferroptosis, mineral absorption, lysosome, RNA degradation, and toll and imd signaling pathway) were mainly involved in this experiment (Figure 2B).

3.6. Metabolism Analysis in S0 and S0.1 Group Fish Liver

After constructing the model with Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) and Partial Least Squares Discriminate Analysis (PLS-DA), the sample clustering results were described using OPLS-DA and PLS-DA score plots (Figure 3A,B). It was shown that metabolites may be clustered depending on differences in the presence or absence of steroidal saponins in high-lipid diets. For example, the S0.1 group was located on the positive side of the T-score and Partial Least Squares-1 (PLS1). The metabolite products were identified in this research, including prostaglandin G2/H2, thromboxane-2, malate, phosphoenolpyruvic acid, and fructose-6-phosphate. Moreover, the metabolites were enriched into 15 KEGG signaling pathways (Figure 3C). The results showed that lipid metabolism (bile secretion, ppar signaling pathway, adipocytokine signaling pathway, regulation of lipolysis in adipocytes, fatty acid metabolism, fatty acid elongation, fatty acid degradation, fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, arachidonic acid metabolism, and steroid hormone biosynthesis) and other metabolism pathways were mainly involved in this experiment.
In addition, thirty-two metabolites related to lipid metabolism are listed in Table 6, of which eight metabolites were significantly down-regulated, nine metabolites were significantly up-regulated, nine metabolites were down-regulated, and six metabolites were up-regulated. At the same time, there are eighteen key metabolites related to glucose metabolism listed in Table 6, of which three metabolites were significantly down-regulated, nine metabolites were down-regulated, two metabolites were significantly up-regulated, and four metabolites were up-regulated.

3.7. Combined Analysis of Transcriptomic and Metabolomic Data

In an attempt to explore the changes in metabolites and mRNA transcripts during glycose lipid metabolism, we performed correlation analysis using transcriptomic and metabolomic data (Figure 4). It was shown that the addition of steroidal saponins to high-lipid diets contributed to activating glycogen breakdown, glycolysis, the TCA cycle, and the fatty acid β-oxidation pathway while inhibiting glycogen synthesis, glycogenesis, fatty acids biosynthesis, cyclooxygenase, and the arachidonic acid signaling pathway.
Next, 14 DEGs were randomly selected from the samples of the S0 and S0.1 groups. Then, they were examined using qRT-PCR to verify the reliability of the RNA-Seq data. The results of qRT-PCR and RNA-seq were consistent, as shown in Figure 4, indicating a reliable result of data analysis.

3.8. Gene Expression of ptgs-1, ptges-1, and tbxas-2

There was no significant difference in liver ptgs-1 and tbxas-2 mRNA expression between the S0.1 group and the S0.2 group, but both groups’ expression levels were significantly lower than those in the S0 group (p < 0.05) (Figure 5). In the S0.2 group, the expression of ptges-1 mRNA was significantly higher than in the S0.1 group while significantly lower than in the S0 group (p < 0.05).

3.9. Gene Expression of Glucose Metabolism in Liver

As indicated in Figure 6, the expression of gsk, f-1,6-bp, and pepck mRNA in the S0.1 group was not significantly different from that in the S0.2 group (p > 0.05), but the three genes’ expression levels were significantly lower in the S0.1 and S0.2 groups than in the S0 group (p < 0.05). The g6pca2 mRNA expression level did not show differences between the S0 and S0.2 groups (p > 0.05), but it was significantly higher in the S0.1 group (p < 0.05). The expression of glut-4 mRNA in the S0 group was significantly higher than that in the S0.2 group but significantly lower than that in the S0.1 group (p < 0.05). The md mRNA expression levels in the S0 group and S0.2 group were not found to be significantly different, while both groups’ levels were significantly lower than that in the S0 group (p < 0.05).

3.10. Gene Expression of Lipid Metabolism in Liver

As can be seen in Figure 7, the fas, acc, and srebp-1 mRNA expression levels between the S0.1 and S0.2 groups were not significantly different (p > 0.05), but they were significantly lower than those in the S0 group (p < 0.05). Although the expression of cpt-1 and acsl mRNA in the S0 group was significantly higher than that in the S0.2 group, it was significantly lower than that in the S0.1 group (p < 0.05). The acadm mRNA expression in the S0 group was significantly lower than that in the S0.1 and S0.2 groups (p < 0.05).

4. Discussion

After ingestion of high-lipid diets (17% and 15% crude lipid), the area of lipid droplets in hybrid grouper liver increased, and the activities of the LZM, AKP, ACP, SOD, and CAT enzymes in the liver and serum decreased, while MDA, ALB, AST, and ALT levels were higher [7,27]. When high-lipid diets (15% crude lipid) were supplemented with saponins, the lipid droplet area decreased in the liver of hybrid grouper, and SOD, CAT, and GSH-PX activities significantly increased to enhance the antioxidant capacity of the fish [20]. In carp [17], white shrimp [28], and swimming crab (Portunus trituberculatus) [29], studies have demonstrated that saponins can decrease serum ALB, AST, and ALT enzyme activities and increase antioxidant and nonspecific immune enzyme activities in the liver to enhance the immune capacity of the body. A previous study demonstrated that 0.1% steroidal saponin significantly increased the activities and gene expression of SOD, CAT, GSH-PX, and GR in hybrid grouper liver and serum and decreased the MDA content [21]. In the present experiment, compared to the control group, supplementation of 0.1% steroidal saponins to a 14% crude lipid diet significantly increased LZM, ACP, and AKP enzyme activities and decreased ALB, ALT, and AST enzyme activities in the serum of hybrid grouper.
Previous studies on Chinese mitten crab (Eriocheir sinensis) [30] and olive flounder (Paralichthys olivaceus) [31] found that arachidonic acid can enhance the antioxidant capacity of the organism by increasing SOD, CAT, AKP, ACP, and LZM activities in the liver and serum. Three metabolic pathways—cytochrome P450, cyclooxygenase, and lipoxygenase—are all engaged in the metabolism of arachidonic acid [32], in which inhibition of the expression and content of key genes and metabolites of the cyclooxygenase metabolic pathway could improve the anti-inflammatory effects of the body [33]. Compared with the control group, the sod, cat, and gsh-px genes’ expression in hybrid grouper liver was up-regulated with the addition of 0.1% steroidal saponins, and the levels of PGG2/PGH2, PGA2 and TXA2 metabolites, ptgs-1 and ptgs-2, ptges-1 and ptges-2, and tbxas-2 mRNA significantly decreased. Both KEGG and GO showed that these genes and metabolites were enriched in the cyclooxygenase metabolic pathway. The upstream ptgs-1 and ptgs-2 genes can activate the cyclooxygenase metabolic pathway to release the intermediate products PGG2 and PGH2, which stimulate the downstream ptges-1, ptges-2, and tbxas-2 genes to produce the downstream products PGA2 and TXA2, mediating the inflammatory response [34,35]. The present study also demonstrated that the supplementation of 0.1% and 0.2% steroidal saponins in the diets contributed to a decrease in the expression of ptgs-1 and ptgs-2, ptges-1 and ptges-2, and tbxas-2 mRNA in the liver of hybrid grouper, inhibited the expression of key genes in the cyclooxygenase metabolic pathway, and played an anti-inflammatory role in the body’s immune system. Research in mice demonstrated that saponins can down-regulate ptgs-1 and ptgs-2, ptges-1 and ptges-2, and tbxas-2 mRNA expression in the liver and decrease the levels of PGG2, PGH2, PGA2, and TXA2 metabolites [36,37]. Consequently, steroid saponins can protect the liver by effectively removing lipid peroxidation products from the body by maintaining oxidative homeostasis through non-specific immune factors, thus improving the immunity of the fish and enabling the orderly metabolism of nutrients in the fish. However, compared to the S0.2 group, 0.1% steroidal saponin significantly decreased serum and liver ALB and AST enzyme activities and ptges-1 gene expression and significantly increased LZM, ACP, and AKP enzyme activities. Research on Atlantic salmon (Salmo salar L.) [38] found that the dietary addition of 0.2% soya saponins increased serum SOD and CAT activities, but supplementation of 0.2% or 0.32% soy saponin decreased carnivorous field eels’ (Monopterus albus) [39] and Japanese flounder’s [40] antioxidant capacity. Research on omnivorous carp [41] found that antioxidant capacity was significantly inhibited when dietary Momordica charantia saponins were provided at levels above 0.64%. Thus, high doses of saponins are toxic to fish and decrease antioxidant capacity.
A previous study showed that supplementation of 0.1% steroidal saponins to high-lipid diets significantly increased the activity and gene expression of hexokinase, pyruvate kinase, and phosphofructokinase in the liver and significantly decreased glucose and glycogen content in the liver and serum of grouper [21]. Those phenomena indicated that steroid saponins had a potential role in stimulating carbohydrate utilization by the organism by elevating the activity and gene expression of the enzymes mentioned above. In this experiment, analysis of transcriptomic and metabolomic data showed that 0.1% steroidal saponins increased the expression of genes such as gp, pk, and md and metabolites such as G-6-P, F-6-P, and PEP in hybrid grouper liver, which were enriched in the glycolysis pathway, and decreased the expression of genes such as gsk, g6pca2, f-1,6-bp, and pepck and metabolite levels such as PYR and OAA in the gluconeogenesis pathway. It was suggested that steroidal saponins can inhibit the gluconeogenesis pathway by activating the glycolysis pathway and gradually reducing intermediate products, such as pyruvate, a product of glycolysis and a substrate of the gluconeogenesis pathway. At the same time, the ingestion of diets containing 0.1% and 0.2% steroidal saponins up-regulated glut-4 and md gene expression in hybrid grouper liver, while down-regulating gsk, g6pca2, f-1,6-bp, and pepck mRNA expression, consistent with the analysis of omics data. An investigation on grouper and carp found that hepatic expression of hk, gk, and pk genes was up-regulated by saponins and down-regulated by g-6-p, f-1,6-bp, and pepck genes in the liver [42,43]. As a result, steroidal saponins can effectively promote the ability of hybrid grouper to utilize carbohydrates in the diet by activating the expression of genes and metabolite content in the glycolysis pathway while inhibiting the gluconeogenesis pathway.
The ingestion of high-lipid diets by grass carp (Ctenopharyngodon idella) (11% and 16% crude lipid) [44,45] and Nile tilapia (Oreochromis niloticus) (12% crude lipid) [46] caused an increase in serum and liver triglyceride and cholesterol levels. Studies on field eel [39], catfish (heterobranchus longifilis) [47], and largemouth bass (Micropterus salmoides) [48] revealed that saponin reduced fas, acc, ppar-γ, and srebp-1 mRNA expression in the liver while elevating the gene expression of cpt-1, cpt-2, acsl, and acadm. In this experiment, supplementation with 0.1% steroidal saponins significantly reduced triglyceride and cholesterol levels in the liver and serum of hybrid grouper [21], decreased the expression of genes such as fas and ppar-γ and metabolite levels such as eicosapentaenoic acid and palmitic acid in the liver, thereby inhibiting fatty acid biosynthesis while elevating the fatty acid β-oxidation pathway in cpt-1, cpt-2, acsl, and acadm gene expression and metabolite levels such as carnitine and citric acid. Meanwhile, the mRNA expression levels of fas, acc, and srebp-1 in the liver significantly decreased, while the gene expression levels of cpt-1, acsl, and acadm increased in the hybrid grouper who were fed a diet supplemented with 0.1% and 0.2% steroidal saponins. Steroid saponins can activate the fatty acid β-oxidation signaling pathway, increase the expression of key genes and metabolite levels, promote fatty acid transport to mitochondria, and, thus, reduce lipid deposition in grouper liver [49]. However, previous studies have found that saikosaponin alleviates hepatic steatosis in hybrid grouper by activating the AMPK/PPARα signaling pathway after the fish consume a high-lipid diet (15% crude fat) [20]. This was inconsistent with the findings of the present study. Saikosaponin belongs to the triterpenoid saponins, while our experimental steroid saponins belong to the steroid group of saponins. It was possible that the saponin species had different effects on hybrid grouper.

5. Conclusions

Compared with the control group, the 0.1% of steroidal saponins effectively scavenged lipid peroxidation products and enhanced the immune defense mechanism of fish by inhibiting the levels of metabolites and genes in the cyclooxygenase metabolic pathway. Further, 0.1% of steroidal saponins promoted more efficient utilization of carbohydrates and lipids in the diets by activating the expression of genes and metabolites levels in the glycolysis pathway and fatty acid β-oxidation pathway and inhibiting the gluconeogenesis pathway, facilitating the conversion and deposition of feed proteins in grouper.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani13182894/s1, Table S1: Ingredients (g/100 g diet) and proximate composition (% dry matter), Table S2: The growth performance of hybrid groupers.

Author Contributions

Conceptualization, B.T., J.Z., Q.Y., S.X. and W.L.; formal analysis, H.D.; funding acquisition, B.T. and S.C.; investigation, H.D. and S.C.; methodology, B.T., G.C., J.Z., Q.Y., X.D., S.X. and W.L.; project administration, S.C.; resources, X.D. and G.C.; software, S.X. and G.C.; writing—original draft, H.D.; writing—review and editing, S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the China Agriculture Research System of MOF and MARA (CARS-47); the Guangdong Provincial Higher Education Key Field Special Project (2020ZDZX1034); the Zhanjiang Science and Technology Plan Project (2020A02001); the Key laboratory construction project of Zhanjiang (2020A05003); and the Postgraduate Education Innovation Project of Guangdong Ocean University (202331).

Institutional Review Board Statement

The hybrid grouper used in this experiment were provided by the Fish Species Technology Extension Station (Zhanjiang, China) and approved by the Animal Research and Ethics Committee of Guangdong Ocean University (approval ID: GDOUIACUC-2021-A0207; approval date: 8 May 2021).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the main article.

Conflicts of Interest

The authors declared that they had no known competing financial interests or personal relationships that could have influenced the work reported in this paper. The author that is affiliated with a company that supplies the product, introduces the effect of product application in other species, and provides comments and suggestions for the design and implementation of this experiment.

References

  1. Fisheries Administration of the Ministry of Agriculture and Rural Areas; National Fisheries Technology Extension Center; China Society of Fisheries. Fishery Statistical Yearbook 2018; Chinese Agricultural Press: Beijing, China, 2018. (In Chinese) [Google Scholar]
  2. Fisheries Administration of the Ministry of Agriculture and Rural Areas; National Fisheries Technology Extension Center; China Society of Fisheries. Fishery Statistical Yearbook 2022; Chinese Agricultural Press: Beijing, China, 2022. (In Chinese) [Google Scholar]
  3. Siddaiah, G.M.; Kumar, R.; Kumari, R.; Chandan, N.K.; Debbarma, J.; Damle, D.K.; Giri, S.S. Dietary fishmeal replacement with Hermetia illucens (Black soldier fly, BSF) larvae meal affected production performance, whole body composition, antioxidant status, and health of snakehead (Channa striata) juveniles. Anim. Feed. Sci. Technol. 2023, 297, 115597. [Google Scholar] [CrossRef]
  4. Lin, S. Review of fish meal market in 2017 and prospect of situation in 2018. Guangdong Feed 2017, 27, 17–18. [Google Scholar]
  5. Lin, S. Review of fish meal market in 2021 and prospect of situation in 2022. Guangdong Feed 2021, 27, 20–22. [Google Scholar]
  6. Li, X.; Jiang, Y.; Liu, W.; Ge, X. Protein-sparing effect of dietary lipid in practical diets for blunt snout bream (Megalobrama amblycephala) fingerlings: Effects on digestive and metabolic responses. Fish Physiol. Biochem. 2012, 38, 529–541. [Google Scholar] [CrossRef]
  7. Xie, R.T.; Amenyogbe, E.; Chen, G.; Huang, J.S. Effects of feed fat level on growth performance, body composition and serum biochemical indices of hybrid grouper (Epinephelus fuscoguttatus × Epinephelus polyphekadion). Aquaculture 2021, 530, 735813. [Google Scholar] [CrossRef]
  8. Meng, Y.; Tian, H.; Hu, X.; Han, B.; Li, X.; Cangzhong, L.; Ma, R. Effects of dietary lipid levels on the lipid deposition and metabolism of subadult triploid rainbow trout (Oncorhynchus mykiss). Aquac. Nutr. 2022, 2022, 6924835. [Google Scholar] [CrossRef]
  9. Li, T.; Yan, X.; Dong, X.; Pan, S.; Tan, B.; Zhang, S.; Zhang, H. Choline Alleviates Disorders of Lipid Metabolism in Hybrid Grouper (♀ Epinephelus fuscoguttatus × ♂ E. lanceolatus) Caused by High-Lipid Diet. Aquac. Nutr. 2022, 2022, 998849. [Google Scholar]
  10. Fan, Z.; Li, J.; Zhang, Y.; Wu, D.; Zheng, X.; Wang, C.A.; Wang, L. Excessive dietary lipid affecting growth performance, feed utilization, lipid deposition, and hepatopancreas lipometabolism of large-sized common carp (Cyprinus carpio). Front. Nutr. 2021, 8, 694426. [Google Scholar] [CrossRef]
  11. Fei, S.; Xia, Y.; Chen, Z.; Liu, C.; Liu, H.; Han, D.; Xie, S. A high-fat diet alters lipid accumulation and oxidative stress and reduces the disease resistance of overwintering hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × P. vachelli ♂). Aquac. Rep. 2022, 23, 101043. [Google Scholar] [CrossRef]
  12. Sun, Z.; Tan, X.; Ye, H.; Zou, C.; Ye, C.; Wang, A. Effects of dietary Panax notoginseng extract on growth performance, fish composition, immune responses, intestinal histology, and immune related genes expression of hybrid grouper (Epinephelus lanceolatus ♂ × Epinephelus fuscoguttatus ♀) fed high lipid diets. Fish Shellfish Immunol. 2018, 73, 234–244. [Google Scholar]
  13. Yang, Q.H.; Tan, B.P.; Dong, X.H.; Chi, S.Y.; Liu, H.Y. Effects of different levels of Yucca schidigera extract on the growth and nonspecific immunity of Pacific white shrimp (Litopenaeus vannamei) and on culture water quality. Aquaculture 2015, 439, 39–44. [Google Scholar] [CrossRef]
  14. Yun, B.; Ai, Q.; Qian, X.; Mai, K. Effects of soya saponins on feed intake, growth performance, and cholesterol metabolism in juvenile turbot (Scophthalmus maximus L.). Isr. J. Aquac. Bamidgeh 2015, 67, 20734. [Google Scholar] [CrossRef]
  15. Yu, C.H.; Xie, G.; He, R.R.; Zhai, Y.J.; Li, Y.F.; Tsoi, B.; Yang, D.P. Effects of a purified saponin mixture from alfalfa on plasma lipid metabolism in hyperlipidemic mice. J. Health Sci. 2011, 57, 401–405. [Google Scholar] [CrossRef]
  16. Upadhyay, S.; Jeena, G.S.; Shukla, R.K. Recent advances in steroidal saponins biosynthesis and in vitro production. Planta 2018, 248, 519–544. [Google Scholar] [CrossRef]
  17. Dawood, M.A.; Gewaily, M.S.; Monier, M.N.; Younis, E.M.; Van Doan, H.; Sewilam, H. The regulatory roles of yucca extract on the growth rate, hepato-renal function, histopathological alterations, and immune-related genes in common carp exposed with acute ammonia stress. Aquaculture 2021, 534, 736287. [Google Scholar] [CrossRef]
  18. Zheng, S.C.; Zhao, Q.T.; Sang, X.Y.; Jiang, J.J.; Xun, L.Y. Study progress on active ingredients of Chinese materia medica and their mechanism to improve insulin resistance. J. Tradit. Chin. Med. 2021, 40, 94–99. [Google Scholar]
  19. You, L.L.; Chen, Y.H.; Liu, J.F.; Jin, X.M.; Feng, C. Influence of bitter melon saponin on the glucose consumption of insulin resistance HepG2 cell in vitro. Sci. Technol. Food Ind. 2014, 35, 338–340+345. [Google Scholar]
  20. Zou, C.; Du, L.; Wu, J.; Gan, S.; Li, Q.; Babu, V.S.; Lin, L. Saikosaponin d alleviates high-fat-diet induced hepatic steatosis in hybrid grouper (Epinephelus lanceolatus ♂ × Epinephelus fuscoguttatus ♀) by targeting AMPK/PPARα pathway. Aquaculture 2022, 553, 738088. [Google Scholar] [CrossRef]
  21. Deng, H.; Zhang, J.; Yang, Q.; Dong, X.; Zhang, S.; Liang, W.; Chi, S. Effects of Dietary Steroid Saponins on Growth Performance, Serum and Liver Glucose, Lipid Metabolism, and Immune Molecules of Hybrid Groupers (♀ Epinephelus fuscoguttatus × ♂ Epinephelus lanceolatu) Fed High-Lipid Diets. Metabolites 2023, 13, 305. [Google Scholar] [CrossRef]
  22. Cui, W.X.; Ma, A.J. Transcriptome analysis provides insights into the effects of myoinositol on the turbot Scophthalmus maximus. Fish Shellfish Immunol. 2020, 106, 691–704. [Google Scholar] [CrossRef]
  23. Chen, M.; Bai, F.; Song, T.; Niu, X.; Wang, X.; Wang, K.; Ye, J. Hepatic transcriptome analysis provides new insight into the lipid-reducing effect of dietary taurine in high-fat fed groupers (Epinephelus coioides). Metabolites 2022, 12, 670. [Google Scholar] [CrossRef] [PubMed]
  24. Yang, S.; Zhao, T.; Ma, A.; Huang, Z.; Liu, Z.; Cui, W.; Yuan, C. Metabolic responses in Scophthalmus maximus kidney subjected to thermal stress. Fish Shellfish Immunol. 2020, 103, 37–46. [Google Scholar] [CrossRef]
  25. Xiang, Q.; Yan, H.; Luo, X.; Kang, Y.; Hu, J.; Chen, L. Integration of transcriptomics and metabolomics reveals damage and recovery mechanisms of fish gills in response to nanosilver exposure. Aquat. Toxicol. 2021, 237, 105895. [Google Scholar] [CrossRef] [PubMed]
  26. Zhao, T.; Ma, A.; Yang, S.; Huang, Z. Integrated metabolome and transcriptome analyses revealing the effects of thermal stress on lipid metabolism in juvenile turbot Scophthalmus maximus. J. Therm. Biol. 2021, 99, 102937. [Google Scholar] [CrossRef]
  27. Pan, S.; Yan, X.; Dong, X.; Li, T.; Suo, X.; Tan, B.; Zhang, H. The positive effects of dietary inositol on juvenile hybrid grouper (♀ Epinephelus fuscoguttatus × ♂ E. lanceolatu) fed high-lipid diets: Growth performance, antioxidant capacity and immunity. Fish Shellfish Immunol. 2022, 26, 84–95. [Google Scholar] [CrossRef]
  28. Su, B.; Chen, J.C. Effect of saponin immersion on enhancement of the immune response of white shrimp and its resistance against Vibrio alginolyticus. Fish Shellfish Immunol. 2008, 24, 74–81. [Google Scholar] [CrossRef]
  29. Ng’ambi, J.W.; Li, R.; Mu, C.; Song, W.; Liu, L.; Wang, C. Dietary administration of saponin stimulates growth of the swimming crab Portunus trituberculatus and enhances its resistance against Vibrio alginolyticus infection. Fish Shellfish Immunol. 2016, 59, 305–311. [Google Scholar] [CrossRef] [PubMed]
  30. Miao, S.; Wan, W.; Hu, J.; Chang, E.; Zhou, Z.; Zhou, Y.; Sun, L. Dietary arachidonic acid affects the innate immunity, antioxidant capacities, intestinal health, and microbiota in Chinese mitten crab (Eriocheir sinensis). Aquaculture 2022, 548, 737635. [Google Scholar] [CrossRef]
  31. Medagoda, N.; Lee, K.-J. Effects of dietary arachidonic acid supplementation in high plant protein diets on growth, feed utilization, and immunity of olive flounder, Paralichthys olivaceus. Aquaculture 2023, 571, 739431. [Google Scholar] [CrossRef]
  32. Wang, B.; Wu, L.; Chen, J.; Dong, L.; Chen, C.; Wen, Z.; Wang, D.W. Metabolism pathways of arachidonic acids: Mechanisms and potential therapeutic targets. Signal Transduct. Target. Ther. 2021, 6, 94. [Google Scholar] [CrossRef]
  33. Ruan, Y.; Fan, Y.; Xie, Y.; Ma, C.; Mo, B.; Lai, Y.; Kuang, W. Modified Xiaoqing long decoction alleviates lipopolysaccharide-induced acute lung injury in mice by regulating arachidonic acid metabolism and exerting anti-apoptotic and anti-inflammatory effects. Anat. Rec. 2022, 305, 1672–1681. [Google Scholar] [CrossRef] [PubMed]
  34. Yoshimura, R.; Sano, H.; Masuda, C.; Kawamura, M.; Tsubouchi, Y.; Chargui, J.; Wada, S. Expression of cyclooxygenase-2 in prostate carcinoma. Cancer 2000, 89, 589–596. [Google Scholar] [CrossRef] [PubMed]
  35. Shoieb, S.M.; El-Ghiaty, M.A.; Alqahtani, M.A.; El-Kadi, A.O. Cytochrome P450-derived eicosanoids and inflammation in liver diseases. Prostaglandins Other Lipid Mediat. 2020, 147, 106400. [Google Scholar] [CrossRef] [PubMed]
  36. Lee, J.E.; Song, B.K.; Kim, J.H.; Siddiqi, M.Z.; Im, W.T. Production of Prosaikogenin F, Prosaikogenin G, Saikogenin F and Saikogenin G by the Recombinant Enzymatic Hydrolysis of Saikosaponin and their Anti-Cancer Effect. Molecules 2022, 27, 3255. [Google Scholar] [CrossRef]
  37. Jung, H.J.; Kim, S.G.; Nam, J.H.; Park, K.K.; Chung, W.Y.; Kim, W.B.; Park, H.J. Isolation of Saponins with the Inhibitory Effect on Nitric Oxide, Prostaglandin E2 and Tumor Necrosis Factor-α Production from Pleurospermum kamtschaticum. Biol. Pharm. Bull. 2005, 28, 1668–1671. [Google Scholar] [CrossRef]
  38. Gu, M.; Gu, J.N.; Penn, M.; Bakke, A.M.; Lein, I.; Krogdahl, A. Effects of diet supplementation of soya-saponins, isoflavones and phytosterols on Atlantic salmon (Salmo salar L.) fry fed from start-feeding. Aquac. Nutr. 2015, 21, 604–613. [Google Scholar] [CrossRef]
  39. Hu, Y.; Zhang, J.; Xue, J.; Chu, W.; Hu, Y. Effects of dietary soy isoflavone and soy saponin on growth performance, intestinal structure, intestinal immunity, and gut microbiota community on rice field eel (Monopterus albus). Aquaculture 2021, 537, 736506. [Google Scholar] [CrossRef]
  40. Chen, W.; Ai, Q.H.; Mai, K.M.; Xu, W.; Liufu, Z.G.; Zhang, W.B.; Cai, Y.H. Effects of dietary soybean saponins on feed intake, growth performance, digestibility, and intestinal structure in juvenile Japanese flounder (Paralichthys olivaceus). Aquaculture 2011, 318, 95–100. [Google Scholar] [CrossRef]
  41. Fan, Z.; Li, J.; Wu, D.; Li, C.; Cao, D.; Miao, L.; Wang, L. Preliminarily curative effectiveness of long-term bitter melon Momordica charantia saponins administration for the glucose homeostasis of juvenile common carp (Cyprinus carpio) fed a high-starch diet. Aquac. Rep. 2022, 25, 101232. [Google Scholar] [CrossRef]
  42. Tan, X.; Sun, Z.; Liu, Q.; Ye, H.; Zou, C.; Ye, C.; Lin, H. Effects of dietary ginkgo biloba leaf extract on growth performance, plasma biochemical parameters, fish composition, immune responses, liver histology, and immune and apoptosis-related genes expression of hybrid grouper (Epinephelus lanceolatus ♂ × Epinephelus fuscoguttatus ♀) fed high lipid diets. Fish Shellfish Immunol. 2018, 72, 399–409. [Google Scholar]
  43. Serrano, A.E., Jr. Effects of Quillaja saponins on growth, feed efficiency, digestive enzyme activities and metabolism of common carp (Cyprinus carpio L.). Aquac. Nutr. 2013, 19, 468–474. [Google Scholar] [CrossRef]
  44. Li, A.; Yuan, X.; Liang, X.F.; Liu, L.; Li, J.; Li, B.; Tao, Y.X. Adaptations of lipid metabolism and food intake in response to low and high fat diets in juvenile grass carp (Ctenopharyngodon idellus). Aquaculture 2016, 457, 43–49. [Google Scholar] [CrossRef]
  45. Tang, T.; Hu, Y.; Peng, M.; Chu, W.; Hu, Y.; Zhong, L. Effects of high-fat diet on growth performance, lipid accumulation and lipid metabolism related MicroRNA/gene expression in the liver of grass carp (Ctenopharyngodon idella). Comp. Biochem. Physiol. Part B Biochem. Mol. Biol. 2019, 234, 34–40. [Google Scholar] [CrossRef] [PubMed]
  46. Lv, H.B.; Ma, Y.Y.; Hu, C.T.; Lin, Q.Y.; Chen, L.Q.; Zhang, M.L.; Qiao, F. The individual and combined effects of hypoxia and high-fat diet feeding on nutrient composition and flesh quality in Nile tilapia (Oreochromis niloticus). Food Chem. 2021, 343, 128479. [Google Scholar] [CrossRef] [PubMed]
  47. Ekanem, A.P.; Inyang Etoh, A.P.; Inyang Etoh, P.C. Antibacterial efficacies of some plant extracts against aeromonas and pseudomonas diseases of farmed catfish (Heterobranchus longifilis). Vet. Sci. Dev. 2011, 1, 11. [Google Scholar] [CrossRef]
  48. Kwasek, K.; Wojno, M.; Patula, S.; Oliaro, F.; Cabay, C.; Pinnell, L.J. The effect of first feeding exposure of larval largemouth bass to a formulated soybean meal-based or soy saponin-supplemented diet on fish growth performance and gut microbiome. N. Am. J. Aquac. 2021, 83, 312–326. [Google Scholar] [CrossRef]
  49. Wen, M.; Fu, X.; Han, X.; Hu, X.; Dong, P.; Xu, J.; Wang, Y. Sea cucumber saponin echinoside A (EA) stimulates hepatic fatty acid β-oxidation and suppresses fatty acid biosynthesis coupling in a diurnal pattern. J. Nutr. Sci. Vitaminol. 2016, 62, 170–177. [Google Scholar] [CrossRef]
Figure 1. Transcriptome analysis in hybrid grouper liver. Note: (A) principal component analysis of livers in S0 and S0.1 (PCA); (B) number of up-regulated and down-regulated differentially expressed genes (DEGs); (C) volcano plot of DEGs between S0 group and S0.1 group.
Figure 1. Transcriptome analysis in hybrid grouper liver. Note: (A) principal component analysis of livers in S0 and S0.1 (PCA); (B) number of up-regulated and down-regulated differentially expressed genes (DEGs); (C) volcano plot of DEGs between S0 group and S0.1 group.
Animals 13 02894 g001
Figure 2. GO and KEGG enrichment of DEGs in fish liver. Note: (A) GO enrichment of DEGs; (B) KEGG enrichment of DEGs; (C) heatmap of DEGs.
Figure 2. GO and KEGG enrichment of DEGs in fish liver. Note: (A) GO enrichment of DEGs; (B) KEGG enrichment of DEGs; (C) heatmap of DEGs.
Animals 13 02894 g002
Figure 3. Metabolome analysis in hybrid grouper liver. Note: (A) the result of discriminant analysis of Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA); (B) identification analysis was performed with Partial Least Squares Discriminate Analysis (PLS-DA); (C) KEGG enrichment of metabolites.
Figure 3. Metabolome analysis in hybrid grouper liver. Note: (A) the result of discriminant analysis of Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA); (B) identification analysis was performed with Partial Least Squares Discriminate Analysis (PLS-DA); (C) KEGG enrichment of metabolites.
Animals 13 02894 g003
Figure 4. Combining transcriptome and metabolome analysis and quantitative reverse transcription PCR verification of DEGs in hybrid grouper liver glucose lipid metabolism. The dynamics of metabolites and gene expression are shown in the diagram above. The detailed DEGs and metabolite expression are shown in Table 3 and Table 4. Note: glycogen phosphorylase (gp), glycogen synthase kinase-1 (gsk-1), glucose (Glu), hexokinase (hk), glucokinase (gk), glucose-6-phosphatase (g-6-p), glucose-6-phosphate (G-6-P), fructose-6-phosphate (F-6-P), phosphofructokinase (pfk), fructose-1,6-bisphosphatase (f-1,6-bp), fructose-1,6-bisphosphate (F-1,6-BP), phosphoenolpyruvic acid (PEP), pyruvate kinase (pk), pyruvic acid (PYR), pyruvate carboxylase (pc), phosphoenolpyruvate carboxykinase (pepck), oxaloacetic acid (OAA), citrate synthase (cs), citrate (CIT), cis-aconitate (CIS), isocitrate (ISO), oxalosuccinate (OXA), isocitrate dehydrogenase (id), 2-oxoglutarate (2-OG), succinyl-CoA (SUC-CoA), succinyl-CoA synthetase (scs), succinic acid (SUCC), succinate dehydrogenase (sd), fumarate (FUM), malate (MAL), malate dehydrogenase (md), peroxisome proliferator-activated receptor γ (ppar-γ), fatty acid synthase (fas), prostaglandin-endoperoxide synthase (ptgs), prostaglandin G2/H2 (PGG2/PGH2), prostaglandin E synthase (ptges), prostaglandin A2 (PGA2), thromboxane-2 synthase (tbxas-2), thromboxane A2 (TXA2), acyl-CoA synthetase (acsl), carnitine acyltransferase 1 (cpt-1), carnitine acyltransferase 2 (cpt-2), and acyl-CoA dehydrogenase (acadm). For the comparison of the S0 group with the S0.1 group, the Student’s t-test was used in the experiment. Asterisks (** and ***) stand for significant differences with p < 0.01 and p < 0.001, respectively.
Figure 4. Combining transcriptome and metabolome analysis and quantitative reverse transcription PCR verification of DEGs in hybrid grouper liver glucose lipid metabolism. The dynamics of metabolites and gene expression are shown in the diagram above. The detailed DEGs and metabolite expression are shown in Table 3 and Table 4. Note: glycogen phosphorylase (gp), glycogen synthase kinase-1 (gsk-1), glucose (Glu), hexokinase (hk), glucokinase (gk), glucose-6-phosphatase (g-6-p), glucose-6-phosphate (G-6-P), fructose-6-phosphate (F-6-P), phosphofructokinase (pfk), fructose-1,6-bisphosphatase (f-1,6-bp), fructose-1,6-bisphosphate (F-1,6-BP), phosphoenolpyruvic acid (PEP), pyruvate kinase (pk), pyruvic acid (PYR), pyruvate carboxylase (pc), phosphoenolpyruvate carboxykinase (pepck), oxaloacetic acid (OAA), citrate synthase (cs), citrate (CIT), cis-aconitate (CIS), isocitrate (ISO), oxalosuccinate (OXA), isocitrate dehydrogenase (id), 2-oxoglutarate (2-OG), succinyl-CoA (SUC-CoA), succinyl-CoA synthetase (scs), succinic acid (SUCC), succinate dehydrogenase (sd), fumarate (FUM), malate (MAL), malate dehydrogenase (md), peroxisome proliferator-activated receptor γ (ppar-γ), fatty acid synthase (fas), prostaglandin-endoperoxide synthase (ptgs), prostaglandin G2/H2 (PGG2/PGH2), prostaglandin E synthase (ptges), prostaglandin A2 (PGA2), thromboxane-2 synthase (tbxas-2), thromboxane A2 (TXA2), acyl-CoA synthetase (acsl), carnitine acyltransferase 1 (cpt-1), carnitine acyltransferase 2 (cpt-2), and acyl-CoA dehydrogenase (acadm). For the comparison of the S0 group with the S0.1 group, the Student’s t-test was used in the experiment. Asterisks (** and ***) stand for significant differences with p < 0.01 and p < 0.001, respectively.
Animals 13 02894 g004
Figure 5. Efficacy of dietary steroidal saponins on arachidonic acid signaling pathway of experimental fish livers. Note: (a–c): Values from smallest to largest.
Figure 5. Efficacy of dietary steroidal saponins on arachidonic acid signaling pathway of experimental fish livers. Note: (a–c): Values from smallest to largest.
Animals 13 02894 g005
Figure 6. Efficacy of dietary steroidal saponins on glucose metabolism in hybrid grouper liver. Note: (a–c): Values from smallest to largest.
Figure 6. Efficacy of dietary steroidal saponins on glucose metabolism in hybrid grouper liver. Note: (a–c): Values from smallest to largest.
Animals 13 02894 g006
Figure 7. Effect of dietary steroidal saponins on lipid metabolism in the liver of hybrid grouper. Note: (a–c): Values from smallest to largest.
Figure 7. Effect of dietary steroidal saponins on lipid metabolism in the liver of hybrid grouper. Note: (a–c): Values from smallest to largest.
Animals 13 02894 g007
Table 1. Primer pair sequences for qRT-PCR.
Table 1. Primer pair sequences for qRT-PCR.
Target GeneNucleotide Sequence (5′-3′)Accession No.
Glycogen phosphorylase
(gp)
F: TCTCCCGTGTTCTTTACCC
R: GCCATTGCTGGATGAGTG
MT548584
Glycogen synthase kinase
(gsk)
F: TACACTGCCTGACCAAGACC
R: TAATGTGGCTGGAGACGAAT
MT548583
Hexokinase-1
(hk)
F: ACTTTGGGTGCAATCCTGAC
R: AGACGACGCACTGTTTTGTG
BC067330.1
Glucokinase
(gk)
F: TGGGTTTTACCTTCTCCTT
R: AGTCCCCTCGTCTCTTGAT
MH213270
Phosphofructokinase
(pfk)
F: AAACGCCCATGCAAACTAC
R: CAACCTCTCTGACAGCCAC
MH213271
Pyruvate kinase
(pk)
F: TCCTGGAGCATCTGTGTCTG
R: GTCTGGCGATGTTCATTCCT
BC152219.1
Glucose-6-phosphatase
(g-6-p)
F: GCCTGTGGATGCTAATGGG
R: GGAGGTCAAGAAGAGAGTCGTG
MT548587
Glucose-6-phosphatase catalytic subunit 2
(g6pca2)
F: CACAGTCCGTCCTCACAT
R: GCAAAACAGCGTCCATAA
MH213269
Fructose-1,6-bisphosphatase
(f-1,6-bp)
F: ACAGTCTGAATGAAGGCTAC
R: CTCATACAACAGCCTCAGCT
XM019099298.1
Pyruvate carboxylase
(pc)
F: GTTGCCCACAACTTCAGCAAG
R: ATGCTGGGCTG TGAAGTCATC
NM_022172.3
Phosphoenolpyruvate carboxykinase (pepck)F: CGTGCTGGACTGGATGTTC
R: CCAAGCCCTGGAGGTTCA
MT548586
Glucose transporter 4
(glut-4)
F: ATTCCTCTACATCGTCCGC
R: CCATCTTCCTCTTCTCCTC
XM_010737169.3
Citrate synthase
(cs)
F: GAAGGAACTCGGTGGTGAGG
R: GACCATAGCCTGGCACAACT
BC165743.1
Succinate dehydrogenase
(sd)
F: ACCTTGGACCTGCTGTGTTGATG
R: GCGGTAGAGGGAGAACGGAT
SRP152627
Malate dehydrogenase
(md)
F: ATGCCATGAATGGAAAGGAG
R: TCCTCGCCCTTCTTGATAGA
XM_038711981.1
Prostaglandin-endoperoxide synthase 1 (ptgs-1)F: AAGGTCTGGGTCACGGAGTG
R: CAAGAGAATCTCCATAAATGTGTCCA
GENE-020530-1
Prostaglandin-endoperoxide synthase 2 (ptgs-2)F: GAGTACTGGAAGCCGAGCAC
R: GATATCACTGCCGCCTGAGT
AM296029
Prostaglandin E synthase 1
(ptges-1)
F: TAATAACGGGACAAGTGAGGC
R: GGCAAACAAGTAGGCAACG
XM_010750942
Prostaglandin E synthase 2
(ptges-2)
F: TGGCTGTGTTTGGCGTCCTCA
R: TCTTGGTGTTCTGCGGTGTCCT
XM_019272321
Thromboxane-2 synthase
(tbxas-2)
F: AGCTGCATGATGGGATCTGTCAATC
R: AGGTGGATGATGCGATGTGTGAATC
AY398422
Peroxisome proliferator-activated receptor γ
(ppar-γ)
F: ACCGCAGCACGAAGAACAAC
R: TGGACGCCATAGTGAAACCC
KM052849.1
Peroxisome proliferator-activated receptor β
(ppar-β)
F: GCTACAGAGCAGCACGACA
R: CTCCTCATCTTCGCTTTCC
DQ232867.1
Fatty acid synthase
(fas)
F: GGCGGCATTGTAGGCATTA
R: CAATCAAAGTGTAGCCTCGGTAG
FJ196231.1
Acetyl-CoA carboxylase
(acc)
F: GAGAAGGCACCAGAAGATCATAG
R: CACAGTACCTGCACTCACATAG
KX066238.1
Sterol regulatory element-binding protein 1
(srebp-1)
F: TGGAGTTTGGAGGACTGTTTG
R: AAGATGAGGGTGGAGTTGGA
KU179485.1
Acyl-CoA synthetase
(acsl)
F: CCATACACATCCACACCGAGT
R: ACAGCTCTTCTATCTGGGGTC
GLEAN_10004723
Carnitine acyltransferase 1
(cpt-1)
F: TGCTCCACGGAAAGTGCTAC
R: GCAGTGACCCTCCTCAGTGTAT
HM037343.1
Carnitine acyltransferase 2
(cpt-2)
F: CTATCTGGAGTGACATCATG
R: CTCACTCACAGGTAGAGATG
XM019122007.1
Acyl-CoA dehydrogenase
(acadm)
F: ACAAGGTTTTGAGGGCAGGT
R: TTGGCACTAGCTTGAGCACT
NM_213010.2
β-actinF: GGCTACTCCTTCACCACCACA
R: TCTGGGCAACGGAACCTCT
AY510710.2
Table 2. Effects of dietary steroidal saponins on biochemical indexes of the hybrid grouper serum.
Table 2. Effects of dietary steroidal saponins on biochemical indexes of the hybrid grouper serum.
IndexesS0S0.1S0.2
ALT/(U/L)66.92 ± 2.29 b45.32 ± 1.05 a45.34 ± 1.71 a
AST/(U/L)8.02 ± 0.28 b6.32 ± 0.22 a7.55 ± 0.18 b
TP/(mg/mL)3.24 ± 0.143.84 ± 0.183.26 ± 0.23
ALB/(g/L)52.05 ± 0.68 c29.72 ± 0.31 a35.89 ± 0.29 b
Note: Values range from minimum to maximum. There are significant differences indicated with different letters (p < 0.05).
Table 3. Effects of dietary steroidal saponins on nonspecial immune indexes of fish serum.
Table 3. Effects of dietary steroidal saponins on nonspecial immune indexes of fish serum.
IndexesS0S0.1S0.2
LZM/(U/mL)48.97 ± 4.48 a77.81 ± 4.38 b52.69 ± 1.25 a
ACP/(King’s unit/L)3.60 ± 0.05 a5.96 ± 0.10 b4.36 ± 0.18 a
AKP/(King’s unit/L)3.04 ± 0.11 a4.66 ± 0.11 b3.22 ± 0.12 a
Note: Values range from minimum to maximum. There are significant differences indicated with different letters (p < 0.05).
Table 4. Efficacy of dietary steroidal saponins on antioxidative indexes of hybrid grouper livers.
Table 4. Efficacy of dietary steroidal saponins on antioxidative indexes of hybrid grouper livers.
IndexesS0S0.1S0.2
MDA/(nmol/mgprot)3.95 ± 0.27 b3.24 ± 0.11 a4.24 ± 0.21 a
SOD/(U/mgprot)137.35 ± 2.65 a164.22 ± 1.73 b133.57 ± 1.98 a
CAT/(U/mgprot)1.68 ± 0.06 a3.43 ± 0.15 b2.11 ± 0.04 a
GSH-PX/(U)4.73 ± 0.09 a7.61 ± 0.04 c6.54 ± 0.18 b
Note: Values range from minimum to maximum. There are significant differences indicated with different letters (p < 0.05).
Table 5. Different expressions of DEGs related to glucose lipid metabolism in liver transcriptome of hybrid grouper.
Table 5. Different expressions of DEGs related to glucose lipid metabolism in liver transcriptome of hybrid grouper.
Gene SymbolRelated DescriptionGene IDLog2FCp Value
pepckphosphoenolpyruvate carboxykinaseUnigene0066979−11.530.00
tbxas-2thromboxane synthaseUnigene0016223−11.320.00
scssuccinyl-CoA synthetaseUnigene0027404−11.210.00
scssuccinyl-CoA synthetaseUnigene0027402−10.940.00
scssuccinyl-CoA synthetaseUnigene0027400−10.490.00
tbxas-2thromboxane synthaseUnigene0025629−10.160.01
ptgesprostaglandin E synthaseUnigene0066035−10.010.00
scssuccinyl-CoA synthetaseUnigene0052415−9.900.00
gsk-1glycogen synthase kinaseUnigene0014315−9.720.03
scssuccinyl-CoA synthetaseUnigene0027399−9.670.00
scssuccinyl-CoA synthetaseUnigene0027397−9.580.00
f-1,6-bpfructose-1,6-bisphosphataseUnigene0026457−9.220.04
gsk-1glycogen synthase kinaseUnigene0014320−9.040.02
fasfatty acid synthaseUnigene0053299−8.940.11
ppar-γperoxisome proliferator-activated receptor gammaUnigene0038259−1.200.00
f-1,6-bpfructose-1,6-bisphosphataseUnigene0086471−0.750.01
pcpyruvate carboxylaseUnigene0012223−0.550.13
g-6-pglucose-6-phosphataseUnigene0029240−0.380.28
cpt-1carnitine palmitoyltransferase 1Unigene00554480.770.08
gkglucokinaseUnigene00457171.520.00
cpt-2carnitine palmitoyltransferase 2Unigene00478902.060.21
cscitrate synthaseUnigene00181032.140.00
sdsuccinate dehydrogenaseUnigene00691702.160.03
acadmacyl-CoA dehydrogenaseUnigene00745612.160.03
sod-2superoxide dismutase 2Unigene00773592.340.00
sdsuccinate dehydrogenaseUnigene00440062.900.03
mdmalate dehydrogenaseUnigene00526662.900.00
cortcarnitineUnigene00650473.080.03
idisocitrate dehydrogenaseUnigene00800073.720.00
pfkphosphofructokinaseUnigene00044883.790.03
pkpyruvate kinaseUnigene00599443.880.00
sdsuccinate dehydrogenaseUnigene00029484.390.01
gpglycogen phosphorylaseUnigene00437454.460.00
hkhexokinaseUnigene00716374.880.00
idisocitrate dehydrogenaseUnigene00079295.210.00
gshglutathione peroxidaseUnigene00043259.520.37
mdmalate dehydrogenaseUnigene00604239.570.01
ptgs-2prostaglandin G/H synthaseUnigene00439989.600.03
mdmalate dehydrogenaseUnigene001791110.030.03
idisocitrate dehydrogenaseUnigene002215110.110.04
acadmacyl-CoA dehydrogenaseUnigene002857710.210.01
catcatalaseUnigene002858410.350.08
acslvery long-chain-specific acyl-CoA dehydrogenaseUnigene007129710.380.03
acadmacyl-CoA dehydrogenaseUnigene002731510.410.02
ptgs-1prostaglandin G/H synthaseUnigene004399910.480.00
sod-1superoxide dismutase 1Unigene005237710.670.00
Table 6. Diverse levels of metabolites related to glucose lipid metabolism in liver metabolomic of hybrid grouper.
Table 6. Diverse levels of metabolites related to glucose lipid metabolism in liver metabolomic of hybrid grouper.
Metabolite NameMetabolite IndexLog2FCp Value
Lipid metabolism
Methyl palmitateCom_5883_pos−2.180.00
Pentadecanoic acidCom_12401_pos−1.950.00
11(Z),14(Z),17(Zz)-Eicosatrienoic acidCom_7643_pos−1.870.00
11-Deoxy prostaglandin F1βCom_422_pos−0.930.01
Cytidine 5′-diphosphocholineCom_1078_pos−0.400.04
LipoamideCom_4536_pos−0.390.02
Prostaglandin H1Com_5241_neg−0.360.53
15(R)-15-Methyl prostaglandin A2Com_2092_pos−0.340.04
Pyridoxal 5′-phosphate hydrateCom_5949_pos−0.280.04
Docosapentaenoic acidCom_21_pos−0.230.18
Prostaglandin F2α-1-glyceryl esterCom_3467_pos−0.20.26
7Z,10Z,13Z,16Z,19Z-docosapentaenoic acidCom_11503_pos−0.20.73
Thromboxane A2Com_7080_pos−0.10.73
all-cis-4,7,10,13,16-Docosapentaenoic acidCom_778_pos−0.090.67
d-threo-Isocitric acidCom_12240_pos−0.090.81
Glycerol-3-phosphateCom_483_neg−0.070.83
Prostaglandin A1 ethyl esterCom_1353_neg−0.010.95
Arachidonic acidCom_22_neg0.000.99
Citric acidCom_168_neg0.060.94
Palmitoyl sphingomyelinCom_1814_pos0.170.54
LevalbuterolCom_11294_pos0.240.66
1-PalmitoylglycerolCom_10_pos0.250.40
l-Palmitoyl carnitineCom_13685_pos0.380.27
N-Palmitoyl taurineCom_3290_pos0.450.01
2-Arachidonyl Glycerol etherCom_1450_pos0.560.00
Myristic acidCom_8920_pos1.100.00
OleamideCom_1_pos1.250.03
12-Oxo phytodienoic acidCom_669_pos1.460.03
4-acetyl-4-(ethoxycarbonyl)heptanedioic acidCom_2516_pos2.310.00
N-TetradecanamideCom_3122_pos2.990.00
(R)-3-Hydroxy myristic acidCom_2141_neg3.190.00
HexadecanamideCom_148_pos3.660.00
Glucose metabolism
Phenyl pyruvic acidCom_735_neg−1.280.45
Guanosine5′-diphosphate (GDP)Com_2804_neg−0.910.02
Trehalose-6-phosphateCom_12250_pos−0.630.33
Hydroxyglutaric acidCom_28_neg−0.600.02
Phosphoenolpyruvic acidCom_1634_neg−0.60.35
Citraconic acidCom_2432_pos−0.520.14
d-GalactosamineCom_904_pos−0.330.03
Malic acidCom_82_neg−0.200.49
Uridine 5′-monophosphateCom_2039_pos−0.150.46
Glycerol-3-phosphateCom_177_neg−0.150.68
Fumaric acidCom_322_neg−0.150.51
d-Glucose-1-phosphateCom_6290_pos−0.020.96
d-Glucose-6-phosphateCom_237_neg0.100.74
N-Acetyl-d-glucosamine 6-phosphateCom_7268_neg0.110.64
Glycerol-1-hexadecanoateCom_1112_pos0.290.04
Benzyl 6-O-beta-d-glucopyranosyl-beta-d-glucopyranosideCom_15594_pos0.340.04
6-phospho-d-glucono-1,5-lactoneCom_4851_pos0.340.47
d-Glucosamine-6-phosphateCom_10302_neg0.950.52
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Deng, H.; Chen, G.; Zhang, J.; Yang, Q.; Dong, X.; Xie, S.; Liang, W.; Tan, B.; Chi, S. Integrated Metabolome and Transcriptome Analyses Reveal the Efficacy of Steroidal Saponins for Glucose and Lipid Metabolism in Hybrid Grouper (♀Epinephelus fuscoguttatus × ♂Epinephelus lanceolatu) Fed Higher-Lipid Diets. Animals 2023, 13, 2894. https://doi.org/10.3390/ani13182894

AMA Style

Deng H, Chen G, Zhang J, Yang Q, Dong X, Xie S, Liang W, Tan B, Chi S. Integrated Metabolome and Transcriptome Analyses Reveal the Efficacy of Steroidal Saponins for Glucose and Lipid Metabolism in Hybrid Grouper (♀Epinephelus fuscoguttatus × ♂Epinephelus lanceolatu) Fed Higher-Lipid Diets. Animals. 2023; 13(18):2894. https://doi.org/10.3390/ani13182894

Chicago/Turabian Style

Deng, Hongjin, Guiqiong Chen, Jiacheng Zhang, Qihui Yang, Xiaohui Dong, Shiwei Xie, Weixing Liang, Beiping Tan, and Shuyan Chi. 2023. "Integrated Metabolome and Transcriptome Analyses Reveal the Efficacy of Steroidal Saponins for Glucose and Lipid Metabolism in Hybrid Grouper (♀Epinephelus fuscoguttatus × ♂Epinephelus lanceolatu) Fed Higher-Lipid Diets" Animals 13, no. 18: 2894. https://doi.org/10.3390/ani13182894

APA Style

Deng, H., Chen, G., Zhang, J., Yang, Q., Dong, X., Xie, S., Liang, W., Tan, B., & Chi, S. (2023). Integrated Metabolome and Transcriptome Analyses Reveal the Efficacy of Steroidal Saponins for Glucose and Lipid Metabolism in Hybrid Grouper (♀Epinephelus fuscoguttatus × ♂Epinephelus lanceolatu) Fed Higher-Lipid Diets. Animals, 13(18), 2894. https://doi.org/10.3390/ani13182894

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop