Skip to main content

ORIGINAL RESEARCH article

Front. Physiol., 23 August 2023
Sec. Aquatic Physiology
This article is part of the Research Topic Nutrients and Metabolism in Fish View all 9 articles

Integrated transcriptome and metabolome analysis reveals a possible mechanism for the regulation of lipid metabolism via vitamin A in rice field eel (Monopterus albus)

Huanhuan Huo,&#x;Huanhuan Huo1,2Chonghua Hu&#x;Chonghua Hu3Qiubai Zhou,
Qiubai Zhou1,2*Liufeng Xiong,Liufeng Xiong1,2Mo Peng,Mo Peng1,2
  • 1College of Animal Science and Technology of Jiangxi Agricultural University, Nanchang, China
  • 2Key Laboratory of Featured Hydrobios Nutrition Physiology and Healthy Breeding, Nanchang, China
  • 3Ganzhou Animal Husbandry and Fisheries Research Institute, Ganzhou, China

To understand the effects of vitamin A on lipid deposition in rice field eels, integrated liver transcriptome and metabolome were conducted and the changes in the genes and metabolites were assessed. Three groups of rice field eel were fed with 0, 200, and 16,000 IU/kg vitamin A supplementations in their diets for 70 days. The total lipid content in the whole body of the rice field eels was significantly increased with the vitamin A supplementations (p < 0.05). Comparative transcriptome analysis revealed 14 pathways and 46 differentially expressed genes involved in lipid metabolism. Sphingolipid metabolism, glycerolipid metabolism, primary bile acid biosynthesis and steroid hormone biosynthesis were significantly enriched pathways. In these pathways, three differential genes phospholipid phosphatase 1a (PLPP1a), phospholipid phosphatase 2b (PLPP2b), cytochrome P450 21a2 (CYP21a2) were consistent with the change trend of lipid content, and the other three differential genes aldo-keto reductase family 1 member D1 (AKR1D1), uridine diphosphate glucuronic acid transferase 1a1 (UGT1a1), cytochrome P450 1a (CYP1a) were opposite. Metabolomic analysis revealed that primary bile acid biosynthesis, sphingolipid metabolism, steroid hormone biosynthesis and biosynthesis of unsaturated fatty acids were all critical for rice field eel metabolic changes in response to vitamin A. Six important differential metabolites (eicosapentaenoic acid, sphinganine, 11-beta-hydroxyprogesterone, hydroxyeicosatetraenoic acid, cholic acid, and glycochenodeoxycholate) were identified and have provided new insights into how vitamin A regulates lipid deposition. Integrated transcriptome and metabolome analyses revealed that primary bile acid biosynthesis was the only remarkably enriched pathway in both the transcriptome and metabolome while that sphingosine was the main metabolite. Based on the above results, we have concluded that vitamin A promotes lipid deposition in the rice field eel through the primary bile acid synthesis pathway, and lipid deposits are widely stored in cell membranes, mainly in the form of sphingosine. These results will provide reference data to help improve our understanding of how vitamin A regulates lipid metabolism.

1 Introduction

Vitamin A (VA) is an essential nutrient for fish that helps to maintain growth and development as it is involved in various physiological reactions in the body including lipid metabolism (Roels, 1969). Studies have shown that the effects of VA on lipid metabolism have species-dependent effects. For example, the crude lipid contents in herring (Lian et al., 2017), largemouth bass (Chen et al., 2020) and sea bass (Zhang et al., 2015) were found to increase with increasing dietary VA levels, and then become stable. In contrast, the whole-body crude lipid content levels of Atlantic salmon (Ørnsrud et al., 2002), groupers (Mohamed et al., 2003) and juvenile flounder (Hernandez et al., 2005) were found to significantly decrease with increases in dietary VA. Furthermore (Thompson et al., 1995) found that VA intake had no significant effect on the body composition of rainbow trout and (Yu et al., 2019) found that VA had no significant effect on the crude lipid content of whole tilapia fish.

There are many mechanisms by which VA regulates lipid metabolism, including gene transcription factors that regulate lipid synthesis and lipid oxidation, signaling pathways related to lipid metabolism, number of adipocytes, secretion of adipocytokines, and epigenetic modifications in mammals (Wang and Yan, 2017). Peroxisome proliferator activated receptors (PPARs) are a class of lipid-activated transcription factors that regulate lipid metabolism. VA deficiency can decrease transcription levels of PPARβ, resulting in increased levels of polyunsaturated fatty acids, such as linoleic, linolenic, arachidonic, and docosahexaenoic (Yang et al., 2005). The activation of PPARα can be inhibited by the liver X receptor -sterol regulatory element-binding protein-1c pathway which downregulates the expression of fat synthesis genes, and the activation of liver X receptor inhibits PPARα-induced fatty acid oxidation (Tomohiro et al., 2003). Studies have suggested that VA promotes fat synthesis, which may be associated with epigenetic modifications. When VA is present, polycomb repression complexes can rapidly induce the separation of the zinc finger protein 423 (Zfp423) promoter, demethylation of the Zfp423 histone, and the expression of Zfp423 to promote fat synthesis (Wang et al., 2016). Kim et al. (2013) showed that the Wnt/β-catenin signaling pathway is involved in the differentiation of 3T3-L1 preadipocytes and that all-trans-retinoicacid inhibits the differentiation of 3T3-L1 preadipocytes through the transcriptional activation of β-catenin, which affects fat synthesis. The effects of VA on animal lipid metabolism are not consistent, as the mechanism is complex and there are organizational differences, and further research is thus required.

The Monopterus albus belongs to the genus eels of synbranchiidae and it is an economically important species in China. In 2022 years, its annual production was 334215 tons. The regulation of lipid levels by VA in rice field eels, however, has not yet been reported. Therefore, this study aims to explore the possible mechanism of VA-mediated regulation of lipid metabolism in rice field eels through comprehensive analysis of transcriptomics and metabolomics.

2 Materials and methods

2.1 Diet preparation

The fundamental feed formula is shown in Table 1. Fish meal was used as main protein sources. VA acetate (500,000 IU/g, Maclin Biochemical Co., Ltd., Shanghai, China) is the supplemental source of VA. Three isonitrogenous and isolipidic diets were formulated adding VA 0, 200 and 16000 IU/kg (3 201, 4 730 and 16890 IU/kg, actually) denoted as LA, MA and HA groups respectively. All raw materials were screened for 80 mesh and thoroughly mixed. After adding 16% moisture, the pellets were made by the pelletizer at 120°C. After air drying, the feed were stored at −20°C.

TABLE 1
www.frontiersin.org

TABLE 1. Composition and nutrient levels of the basal diet(DM basis).

2.2 Chemicals and reagents

Concentrated sulfuric acid, boric acid, hydrochloric acid, sodium hydroxide, petroleum ether, methanol, acetonitrile acetone, sodium sulfate and anhydrous ethanol purchased from Nanchang Jinsha chemical raw materials Co. LTD. Trizol reagent, electrophoretic buffer, Reverse transcription kit, loading buffer and SYBR Premix Ex Taq kite purchased from Jiangxi Guanyin Biotechnology Co. LTD.

2.3 Experimental animals and sample collection

The juvenile Monopterus albus were artificially bred by our team. All fish were reared temporarily for 14 days to acclimate to the conditions. After fasting for1 day, 480 fish (7.3 ± 0.02 g) were randomly assigned to 12 separate tanks of 3 groups (HA, MA, and LA, four repetitions per group) with 40 fish per tank. Each tank was filled with the right amount of fresh grass, Eichhornia crassipes. Fish were fed 3–5% of the body weight once daily at 18:00 p.m. A tuck net removed the feces and residual feeds. The breeding experiment lasted for 10 weeks and replaced with clean water weekly. The environmental conditions were: water depth 0.4 m; temperature, 28°C ± 1°C; dissolved oxygen >5.0 mg/L.

After the feeding trial, the fish were anaesthetized using MS222. Three fish from each tank were sampled randomly. The liver from twelve fish each group were dissected into 1.5 mL tubes and stored in −80°C for further RNA extraction. Three biological replicates were used for transcriptome sequencing and six biological replicates were used for metabolome detection. The remaining three fish each group were used to detect body composition.

2.4 Effects of dietary VA level on Monopterus albus

The crude protein, crude lipid, moisture and ash in whole fish body were measured following the methods of the Association of Official Analytical Chemists. The crude protein was determined by Kjeldahl method (FOSS, Kjeltec TM 8200). The crude lipid was estimated by using Soxhlet extraction method (FOSS, Soxtec 2050). The moisture was determined by drying samples in an oven at 105°C until constant weight. The ash was determined by combustion in a muffle furnace at 550°C for 4 h.

2.5 RNA extraction and transriptome sequencing

Liver RNA was extracted using Trizol reagent (Invitrogen, CA, United States) according to manufacturer’s instructions. Agarose gel electrophoresis was used to monitor RNA degradation. The Nano-Drop 2000 spectrophotometer was used to detect RNA concentration and purity. Further, the Agilent 2100 RNA 6000 Nano kitwas used to assess RNA integrity. Higher quality RNA was selected for high-throughput sequencing. Three RNA replicates of each sample were used to build RNA-seq library. After total RNA was extracted, the first strand cDNA was synthesized by random hexers and reverse transcriptase, and the second strand cDNA was synthesized by DNA polymerase I and RNase H. Then, a cDNA library was constructed for sequencing. The cDNA library was sequenced using Illumina HiSeqTM4000 by Gene Denovo Biotechnology Co. (Guangzhou, China). The transcriptome results were compared pair-to-pair.

2.6 Transriptome experimental validation using qPCR

Six genes (Table 2) were selected randomly for validation of RNA-Seq data by qPCR using a SYBR Premix Ex Taq kit (Invitrogen) according to the manufacturer’s instructions. The same RNA samples were used for both Illumina library synthesis and the qPCR verification assay. The first strand cDNA was obtained from 2 μg of total RNA using a PrimeScript first strand cDNA synthesis kit (Takara, Dalian, China). Melting curve analyses were performed following amplification. The specific primers used for qPCR are listed in Additional Table 2, and 18 s gene was used as an endogenous control. The thermal profile for SYBR Green qPCR was 95°C for 90 s, followed by 40 cycles of 95°C for 5 s, 60°C for 15 s, 72°C for 20 s.

TABLE 2
www.frontiersin.org

TABLE 2. The information of primers used to qPCR.

2.7 Metabolites extraction LC-MS/MS analysis

The Liver samples were thawed at 4°C, appropriate samples were added into the pre-cooled methanol/acetonitrile/water solution (2: 2, v/v), vortex mixing, ultrasonic for 30 min, standing at −20°C for 10 min. After 14,000 g at 4°C for 20 min, supernatant vacuum drying, adding 100 μL acetonitrile solution during mass spectrometry (acetonitrile: Water = 1:1, v/v) redissolved, vortex, 14,000 g centrifuged at 4°C for 15 min. Then, supernatant was taken for analysis. Analysis was performed using an UHPLC (1290 Infinity LC, Agilent Technologies) coupled to a quadrupole time-of-flight (AB Sciex TripleTOF 6600) in Shanghai Applied Protein Technology Co.,Ltd. The metabolome results were compared in pairs.

2.8 The integrated analysis of transcriptomic and metabolomics data

Firstly, differential genes and differential metabolites related to lipid metabolism were identified by KEGG pathway analysis. Secondly, O2PLS (two-way orthogonal projections to late structures) model was established to analyze the relationship between them. Finally, pearson correlation coefficients were used to assess the correlation between differential genes and differential metabolites.

2.9 Statistical analysis

The results were expressed as mean ± SD. (standard deviation of the mean) and were analyzed using one-way ANOVA in SPSS 19.0 software. Any significant differences were further investigated by comparing the group means using Turkey’s test. Statistical significance was considered if p < 0.05.

3 Results

3.1 Effects of dietary VA levels on Monopterus albus

The crude lipid content in the whole body was significantly different (p < 0.05) among the different treatments, and that in the whole body of the LA group was significantly lower than that of the MA and HA groups (Table 3). The crude lipid content of the whole body increased with increasing dietary VA levels. In contrast, the differences in crude protein, moisture, and ash in the whole body were not significant.

TABLE 3
www.frontiersin.org

TABLE 3. Effects of dietary VA levels on the nutritional composition of whole fish.

3.2 Transcriptomic analysis

To validate the transcriptome data, 6 DEGs were selected randomly for qPCR. Althought the test DEGs displayed different expression levels (Figure 1), in general, the qPCR results showed a positive correlation with transcriptome, indicating the reliability and accuracy of the transcriptome analysis.

FIGURE 1
www.frontiersin.org

FIGURE 1. Validation of transcriptome data by using qPCR.

After removing the low-quality reads, 50.52 million clean reads were obtained. A total of 32.0 Gb of data were generated by Illumina sequencing. The results of transcriptome have been submitted to NCBI (SRA accession no. SRP398570).Comparative transcriptome analysis identified 1323 differentially expressed genes (DEGs) among the three groups (Figure 2). Based on the KEGG annotation, lipid metabolism-related DEGs and pathways were further analyzed. A total of 46 DEGs were enriched in 14 lipid metabolism-related pathways (Table 4). Among them, LAvsMA contained 10 pathways and 16 DEGs; MAvsHA contained 14 pathways and 36 DEGs; LAvsHA contained 14 pathways and 46 DEGs. Significance analysis revealed that sphingolipid metabolism (ko00600), glycerol metabolism (ko00561), steroid hormone biosynthesis (ko00140), and primary bile acid biosynthesis (ko00120) were significantly enriched. There were 16 DEGs in the four significantly enriched pathways. Due to the crude lipid content was continuously rising, further trend analysis of the screened DEGs was necessary. Trend analysis revealed that PLPP1a, PLPP2b, and CYP21a2 were continuously upregulated, AKR1D1, UGT1a1, and CYP1a were continuously downregulated. The relationships between six genes and the four significantly enriched pathways are shown in Table 5.

FIGURE 2
www.frontiersin.org

FIGURE 2. The venn of DEGs distribution in three groups.

TABLE 4
www.frontiersin.org

TABLE 4. The fourteen lipid metabolism-related pathways.

TABLE 5
www.frontiersin.org

TABLE 5. Significantly enriched pathways and DEGs.

3.3 Metabolite extraction and LC-MS/MS analysis

According to analyses of the metabolites, 10878 metabolites and 8540 metabolites were obtained in positive and negative modes, respectively. Based on the variable important in projection (VIP) values of OPLS-DA, the significantly different metabolites were screened based on a VIP ≥1 and p ≤ 0.05. In negative ion modes, a total of 201 different metabolites were detected in the three groups (Figure 3). To explore the lipid metabolism pathways, we selected different lipid-relevant metabolites for KEGG enrichment analysis. The five lipid metabolism pathways were mainly related to arachidonic acid metabolism, primary bile acid biosynthesis, unsaturated fatty acid biosynthesis, sphingolipid metabolism, and steroid hormone biosynthesis. There were six metabolites identified in the five pathways: eicosapentaenoic acid, sphinganine, 11-beta-hydroxyprogesterone, hydroxyeicosatetraenoic acid, cholic acid, and glycochenodeoxycholate (Table 6). In positive ion modes, no metabolites were detected.

FIGURE 3
www.frontiersin.org

FIGURE 3. The venn of different metabolites distribution in three groups.

TABLE 6
www.frontiersin.org

TABLE 6. Significantly enriched pathways and differential metabolites.

3.4 Integrative transcriptome and metabolome analysis

To screen for DEGs and different metabolites associated with lipid metabolism, integrated transcriptome and metabolome analyses were performed using the O2PLS, correlation coefficient models and KEGG pathway. The common KEGG pathways were analyzed for the differential genes and metabolites. Integrative analysis revealed that primary bile acid and steroid hormone biosynthesis were enriched in both the transcriptome and metabolome.

Transcriptomic and metabolomic data were correlated analysis by O2PLS model. Components of the O2PLS model were calculated and reflected the contribution of the transcriptome and metabolome to the total variation. Therefore, the O2PLS model was very reliable for explaining the total variation in the comprehensive analysis of transcriptome and metabolome.The combined loading diagram of metabolites and transcripts was constructed according to the loading values shown in Figure 4. High absolute loading values indicated a strong correlation between the differential genes and differential metabolites. DEGs and differential metabolites farther from the origin were more closely related to each other. The correlation between metabolites and genes were evaluated by calculating pearson coefficients. The correlation matrix of the heat map revealed negative (blue) and positive (red) correlations between lipid-related differential metabolites and genes (Figure 5). Correlation network maps showed lipid-related differential genes or metabolites at important associated positions (Figure 6). Sphinganine was the main associated metabolite and it was correlated with five DEGs.

FIGURE 4
www.frontiersin.org

FIGURE 4. Loading plots for the differential metabolites and differential genes. Note: loading value represents the explanatory ability of a variable (metabolite/gene) in each component (that is, the contribution to the difference between groups). The positive and negative loading value represents the positive or negative correlation with another group. The greater the absolute value of the load, the stronger the correlation.

FIGURE 5
www.frontiersin.org

FIGURE 5. Heat plot showing the correlations between differential metabolites and differential genes. Note: red and blue indicate positive and negative correlations between the transcriptomics and metabolomics data, respectively.

FIGURE 6
www.frontiersin.org

FIGURE 6. Network showing the correlations between metabolites and genes. Note: circles represent metabolites and squares represent genes. The absolute value of pearson coefficient ≥0.5.

4 Discussion

Currently, the effects of VA on lipid metabolism in fish have been reported, however there are few studies on its regulatory mechanism. Therefore, this study aims to explore the possible mechanism of VA-mediated regulation of lipid metabolism through comprehensive analysis of transcriptomics and metabolomics.

4.1 VA and bile acid biosynthesis

The primary bile acid biosynthesis pathway was significantly enriched in both the transcriptome and metabolome, indicating that there is a complex interaction between VA and bile acids. Retinoic acid receptor (RAR), farnesoid X receptor (FXR), and retinoid X receptor (RXR) are the main regulators of VA and bile acid homeostasis (Zoubek et al., 2017). In the liver, VA and its metabolites activate RAR and RXR, generate RAR/RXR heterodimers, and regulate bile acid synthesis and metabolism (Li et al., 2021). In addition, retinoic acid can inhibit the genes involved 5in bile acid synthesis (Mamoon et al., 2014), in which RAR and RXR may be involved. Studies have shown that VA inhibits the expression of the bile salt synthesis rate-limiting enzyme CYP7A1, thereby reducing bile generation (Jahn et al., 2016). In this study, AKR1D1 was inhibited instead of CYP7A1. AKR1D1 plays an important role in steroid metabolism (Gathercole et al., 2021), but biologically, it is better known for its involvement in bile acid biosynthesis, where it catalyzes a key step that introduces the 5β-configuration into primary bile acids (Russell, 2009). This decrease in cholic acid levels may be due to the inhibition of AKR1D1. Another reason is that cholic acid is used to produce 12α-hydroxylated bile acid, which lead to hepatic lipid accumulation (Maegawa et al., 2021). In contrast, the glycochenodeoxycholate content was found to be increased in the metabolome. Glycochenodeoxycholate is a hydrophobic bile salt and that can cause lipid accumulation in the liver (Chen and Chen, 2014). Simultaneously, glycochenodeoxycholate can increase lipid solubility and promote intestinal absorption (Wu et al., 2016). Notably, high concentrations of glycochenodeoxycholate induced hepatocyte apoptosis (Wang et al., 2005). This may also explain why high levels of liver fat are disadvantageous in fish.

4.2 VA and fatty acid synthesis

Fatty acids are components of various molecules, including phospholipids, sphingolipids, and esters. The regulation of fatty acid metabolism by VA has been confirmed in other studies (Chen and Chen, 2014; Blaner, 2019). The genes for fatty acid extension, desaturation and synthesis pathways can be regulated by VA. This is achieved through retinoic acid (RA), a product of VA metabolism.

The first step in fatty acid synthesis is the conversion of acetyl-CoA to malonyl-CoA by acetyl-CoA carboxylase. It is then catalyzed by fatty acid synthase to form palmitic acid. Palmitic acid can be prolonged and desaturated to produce fatty acids with longer chains and double bonds (Kihara, 2012; Zhang et al., 2016). Elongases add two carbon atoms each time to create longer-chain fatty acids. Desaturases introduce double bonds into saturated and unsaturated fatty acids, resulting in the production of monounsaturated and polyunsaturated fatty acids, respectively. The mechanism of VA regulate fatty acid synthesis has been extensively studied. Studies have shown that when VA deficiency, the mRNA levels of acetyl-CoA carboxylase, fatty acid synthase and stearoyl-CoA desaturase 1 were reduced; when RA treatment, the mRNA levels were increased (Yang et al., 2021).

Ayuso et al. (2015) found that VA increases the content of polyunsaturated fatty acids in pork. Similar to the results of our study, the levels of the polyunsaturated fatty acids, eicosapentaenoic and hydroxyeicosatetraenoic acid, increased significantly. This is probably because of the increase in stearoyl-CoA desaturase 1. Stearoyl-CoA desaturase 1 was responsible for the formation of double bond and regulation of lipogenesis in the process of fatty acid synthesis (Koeberle et al., 2016; ALJohani et al., 2017). RA could induce stearoyl-CoA desaturase 1 mRNA expression (Miller et al., 1997) and desaturase index in adipose tissues (Gorocica-Buenfil et al., 2008).

4.3 VA and phospholipid synthesis

Oliveros et al. (2007) found that VA deficiency induced a hypolipidemic effect by decreasing the total phospholipid content in the liver, owing to low phosphatidylcholine synthesis and enhanced fatty acid oxidation in rats, and this was consistent with the results of our study. In the present study, PLPP1a and PLPP2b were significantly upregulated with increasing VA content. PLPP1a and PLPP2b are phospholipid phosphatases that belong to the lipid phosphate phosphatase family. Lipid phosphate phosphatase has wide substrate selectivity and can hydrolyze a variety of lipid phosphoric acids containing single lipid bonds, including phosphatidic acid, lysophosphatidic acid, diacylglycerol pyrophosphate, and sphingosine 1-phosphate, and it participates in the synthesis of glycerol (Toke et al., 1998; Li et al., 2017; Tang and Brindley, 2020). The significant increase in sphinganine levels in the metabolome may be related to the hydrolysis of sphingosine 1-phosphate.

4.4 VA and steroid hormone synthesis

VA is essential for the synthesis and secretion of steroid hormones. The expression of 11β-hydroxysteroid dehydrogenase type 1 was increased by VA deficiency in the hypothalamus (C62.5%) and hippocampus (C104.7%) (Marissal-Arvy et al., 2013). RA treatment could reverse the effects of VA deficiency on 11β-hydroxysteroid expression. The inhibitory effect of retinoids on 11β-hydroxysteroid dehydrogenase type 1 expression has been demonstrated peripherally in the liver and visceral fat of WNIN/Ob obese rats (Sakamuri et al., 2011), in the adipose tissue of mice (Sakuta et al., 2006), and in myotubes in vitro (Aubry and Odermatt, 2009). In contrast, in this study, the level of 11β-hydroxyprogesterone was increased by VA, but this discrepancy may be species related. Within the cell, 11β-hydroxysteroid dehydrogenase type 1 catalyzes the regeneration of active glucocorticoids, thereby amplifying their action (Marissal-Arvy et al., 2013). Similarly, in our study, we found that CYP21a2 was associated with glucocorticoid synthesis (Heidelberg, 2008) which was significantly upregulated.

In this study, comparative metabolomics and transcriptomics were applied to investigate lipid metabolite changes in rice field eel when fed different levels of VA. Comparative transcriptome analysis revealed that VA affected multiple lipid metabolic pathways, especially sphingolipid and glycerin metabolism, and steroid hormone and primary bile acid biosynthesis. Metabolite analysis revealed that the levels of lipids, such as eicosapentaenoic acid, sphinganine, 11-beta-hydroxyprogesterone, hydroxyeicosatetraenoic acid, and glycochenodeoxycholate, increased significantly. Integrative analyses of the transcriptome and metabolome revealed a significant correlation between differential metabolites and differential genes, suggesting that these are the main metabolic pathways. These findings provide reliable models at the transcriptional and metabolic levels which will facilitate further investigation into the effects of VA on lipid metabolism.

5 Conclusion

In conclusion, VA is found to promote lipid deposition in rice field eel. The most important pathway is primary bile acid synthesis and AKR1D1 plays a key role in this process. Integrated transcriptome and metabolome analyses have revealed that sphingosine is the main form of lipid deposits. Future work can be focused on the primary bile acid synthesis pathway, AKR1D1 gene and sphingosine.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

The animal study was approved by the Committee on Research Ethics of the Department of Laboratory Animal Science, Jingxi Agricultural University. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

HH: Conceptualization, Data curation, Writing–original draft. CH: Data curation, Project administration, Writing–original draft. QZ: Conceptualization, Supervision, Writing–review and editing. LX: Formal Analysis, Resources, Writing–review and editing. MP: Investigation, Methodology, Validation, Writing–review and editing.

Funding

The research was supported by the Youth Scientific Research Funds of Science and Technology Department of Jiangxi Province (20202BABL215026), the earmarked fund for CARS-46 and the Jiangxi Agricultural Researches System (JXARS-03).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Aljohani, A. M., Syed, D. N., and Ntambi, J. M. (2017). Insights into stearoyl-CoA desaturase-1 regulation of systemic metabolism. Trends Inendocrinology Metabolism 28, 831–842. doi:10.1016/j.tem.2017.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Aubry, E. M., and Odermatt, A. (2009). Retinoic acid reduces glucocorticoid sensitivity in C2C12 myotubes by decreasing 11beta-hydroxysteroid dehydrogenase type 1 and glucocorticoid receptor activities. Endocrinology 150, 2700–2708. doi:10.1210/en.2008-1618

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayuso, M., Fernández, A., Isabel, B., Rey, A., Benítez, R., Daza, A., et al. (2015). Long term vitamin A restriction improves meat quality parameters and modifies gene expression in Iberian pigs. J. Animal Sci. 93, 2730–2744. doi:10.2527/jas.2014-8573

PubMed Abstract | CrossRef Full Text | Google Scholar

Blaner, W. S. (2019). Vitamin A signaling and homeostasis in obesity, diabetes, and metabolic disorders. Pharmacol. Ther. 197, 153–178. doi:10.1016/j.pharmthera.2019.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S. J., Wu, C. L., and Ye, J. Y. (2020). Effect of dietary vitamin A on growth, serum biochemical index, digestive enzyme activities and glucose and lipid metabolism in juvenile Mylopharyngodon piceu. J. Fish. China 44 (1), 85–98. doi:10.11964/jfc.20190511785

CrossRef Full Text | Google Scholar

Chen, W., and Chen, G. X. (2014). The roles of vitamin A in the regulation of carbohydrate, lipid, and protein metabolism. J. Clin. Med. 3 (2), 453–479. doi:10.3390/jcm3020453

PubMed Abstract | CrossRef Full Text | Google Scholar

Gathercole, L. L., Nikolaou, N., Arvaniti, A., Harris, S. E., Poolman, T. M., Haziehurst, J. M., et al. (2021). Akr1d1-/- mice have a sexually dimorphic metabolic phenotype with reduced fat mass, increased insulin sensitivity and hypertriglyceridemia in males. Cold Spring Harb. Lab. doi:10.1101/2021.02.02.429227

CrossRef Full Text | Google Scholar

Gorocica-Buenfil, M. A., Fluharty, F. L., and Loerch, S. C. (2008). Effect of vitamin A restriction on carcass characteristics and immune status of beef steers. J. Animal Scince 86, 1609–1616. doi:10.2527/jas.2007-0241

PubMed Abstract | CrossRef Full Text | Google Scholar

Heidelberg, S. B. (2008). P450 mono-oxygenase system [M]. Springer Berlin Heidelberg.

Google Scholar

Hernandez, L. H. H., Teshima, S. I., Ishikawa, M., Alam, S., Koshio, S., and Tanaka, Y. (2005). Dietary vitamin A requirements of juvenile Japanese flounder Paralichthys olivaceus. Aquac. Nutr. 11 (1), 3–9. doi:10.1111/j.1365-2095.2004.00317.x

CrossRef Full Text | Google Scholar

Jahn, D., Sutor, D., Dorbath, D., Weiß, J., Götze, O., Schmitt, J., et al. (2016). Farnesoid X receptor-dependent and -independent pathways mediate the transcriptional control of human fibroblast growth factor 19 by vitamin A. Biochim. Biophys. Acta 1859 (2), 381–392. doi:10.1016/j.bbagrm.2015.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Kihara, A. (2012). Very long-chain fatty acids: elongation, physiology and related disorders. J. Biochem. 152, 387–395. doi:10.1093/jb/mvs105

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, D. M., Choi, H. R., Park, A., Shin, S. M., Bae, K. H., Lee, S. C., et al. (2013). Retinoic acid inhibits adipogenesis via activation of Wnt signaling pathway in 3T3-L1 preadipocytes. Biochem. Biophysical Res. Commun. 434 (3), 455–459. doi:10.1016/j.bbrc.2013.03.095

PubMed Abstract | CrossRef Full Text | Google Scholar

Koeberle, A., Löser, K., and Thürmer, M. (2016). Stearoyl-CoA desaturase-1 and adaptive stress signaling. Biochimica Biophysica Acta 1861, 1719–1726. doi:10.1016/j.bbalip.2016.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Cai, S. Y., and Boyer, J. L. (2021). The Role of the Retinoid receptor, RAR/RXR heterodimer, in liver physiology. Biochimica Biophysica Acta (BBA)-Molecular Basis Dis. 1867 (5), 166085. doi:10.1016/j.bbadis.2021.166085

CrossRef Full Text | Google Scholar

Li, Y. L., Zhang, Q. Q., and Hu, W. Q. (2017). Roles and regulation of phosphatidic acid phosphatase in lipid metabolism and signaling. Plant Physiology J. 53 (6), 897–904.

Google Scholar

Lian, X. Y., Chen, N. S., Wang, M. L., Yan, C., and Ding, G. (2017). Dietary vitamin A requirement of largemouth bass (Micropterus salmoides). Chin. J. Animal Nutr. 29 (10), 3819–3830. doi:10.3969/j.issn.1006-267X.2017.10.046

CrossRef Full Text | Google Scholar

Maegawa, K., Koyama, H., Fukiya, S., Yokota, A., Ueda, K., and Ishizuka, S. (2021). Dietary raffinose ameliorates hepatic lipid accumulation induced by cholic acid via modulation of enterohepatic bile acid circulation in rats. Br. J. Nutr. 127 (1), 1621–1630. doi:10.1017/S0007114521002610

PubMed Abstract | CrossRef Full Text | Google Scholar

Mamoon, A., Subaute, A., Subaute, M. C., and Subauste, J. (2014). Retinoic acid regulates several genes in bile acid and lipid metabolism via upregulation of small heterodimer partner in hepatocytes. Gene 550 (2), 165–170. doi:10.1016/j.gene.2014.07.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Marissal-Arvy, N., Hamiani, R., Richard, E., Moisan, M. P., and Pallet, V. (2013). Vitamin A regulates hypothalamic-pituitary-adrenal axis status in LOU/C rats. J. Endocrinol. 219 (1), 21–27. doi:10.1530/JOE-13-0062

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, C. W., Waters, K. M., and Ntambi, J. M. (1997). Regulation of hepatic stearoyl-CoA desaturase gene 1 by vitamin A. Biochem. Biophysical Res. Communiction 231, 206–210. doi:10.1006/bbrc.1997.6070

PubMed Abstract | CrossRef Full Text | Google Scholar

Mohamed, J. S., Sivaram, V., Roy, T. S. C., Peter Marian, M., Murugadass, S., and Raffiq Hussain, M. (2003). Dietary vitamin A requirement of juvenile greasy grouper (Epinephelus tauvina). Aquaculture 219 (1-4), 693–701. doi:10.1016/s0044-8486(02)00665-8

CrossRef Full Text | Google Scholar

Oliveros, L. B., Domeniconi, M. A., Vega, V. A., Gatica, L. V., Brigada, A. M., and Gimenez, M. S. (2007). Vitamin A deficiency modifies lipid metabolism in rat liver. Br. J. Nutr. 97, 263–272. doi:10.1017/S0007114507182659

PubMed Abstract | CrossRef Full Text | Google Scholar

Ørnsrud, R., Graff, I. E., Høie, S., Totland, G. K., and Hemre, G. I. (2002). Hypervitaminosis A in first-feeding fry of the Atlantic salmon (Salmo salar L). Aquac. Nutr. 8 (1), 7–13. doi:10.1046/j.1365-2095.2002.00185.x

CrossRef Full Text | Google Scholar

Roels, O. A. (1969). The fifth decade of vitamin A research. Am. J. Clin. Nutr. 22 (7), 903–907. doi:10.1093/ajcn/22.7.903

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, D. W. (2009). Fifty years of advances in bile acid synthesis and metabolism. J. Lipid Res. 50, S120–S125. doi:10.1194/jlr.R800026-JLR200

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakamuri, V. P. S. S., Ananthathmakula, P., Veettil, G. N., and Ayyalasomayajula, V. (2011). Vitamin A decreases pre-receptor amplification of glucocorticoids in obesity: study on the effect of vitamin A on 11beta-hydroxysteroid dehydrogenase type 1 activity in liver and visceral fat of WNIN/ob obese rats. Nutr. J. 10 (1), 70. doi:10.1186/1475-2891-10-70

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakuta, T., Uchiyama, T., and Kanayama, T. (2006). Topical ER36009, a RARgamma-selective retinoid, decreases abdominal white adipose tissue and elicits changes in expression of genes related to adiposity and thermogenesis. Endocrine 30 (1), 113–119. doi:10.1385/ENDO:30:1:113

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, X., and Brindley, D. N. (2020). Lipid phosphate phosphatases and cancer. Biomolecules 10 (9), 1263. doi:10.3390/biom10091263

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, I., Choubert, G., Houlihan, D. F., and Secombes, C. (1995). The effect of dietary vitamin A and astaxanthin on the immunocompetence of rainbow trout. Aquaculture 133 (2), 91–102. doi:10.1016/0044-8486(95)00024-v

CrossRef Full Text | Google Scholar

Toke, D. A., Bennett, W. L., Dillon, D. A., Wu, W. I., Chen, X., Ostrander, D. B., et al. (1998). Isolation and characterization of the saccharomyces cerevisiae DPP1 gene encoding diacylglycerol pyrophosphate phosphatase. J. Biol. Chem. 273 (6), 3278–3284. doi:10.1074/jbc.273.6.3278

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomohiro, I., Hitoshi, S., Tomohiro, Y., Yahagi, N., Amemiya-Kudo, M., Matsuzaka, T., et al. (2003). Cross-talk between peroxisome proliferator-activated receptor (PPAR) alpha and liver X receptor (LXR) in nutritional regulation of fatty acid metabolism. I. PPARs suppress sterol regulatory element binding protein-1c promoter through inhibition of LXR signaling. Mol. Endocrinol. 17 (7), 1240–1254. doi:10.1210/me.2002-0190

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, B., Yang, Q., Harris, C. L., Nelson, M. L., Busboom, J. R., Zhu, M. J., et al. (2016). Nutrigenomic regulation of adipose tissue development-role of retinoic acid: a review. Meat Sci. 120, 100–106. doi:10.1016/j.meatsci.2016.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Brems, J. J., Gamelli, R. L., and Ding, J. (2005). Reversibility of caspase activation and its role during glycochenodeoxycholate-induced hepatocyte apoptosis. J. Biol. Chem. 280 (25), 23490–23495. doi:10.1074/jbc.M411607200

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., and Yan, S. M. (2017). Regulatory effects of vitamin A on lipid metabolism of animals and its mechanism. Chin. J. Animal Nutr. 29 (5), 7. doi:10.3969/j.issn.1006?267x.2017.05.002

CrossRef Full Text | Google Scholar

Wu, F. Z., Yang, L. J., and Zhang, Q. (2016). Effect of bile salts and lipids in intestine on solubility of artemether as a poorly water-soluble drug. Chin. J. Exp. Traditional Med. Formulae 22 (20), 6–9. doi:10.13422/j.cnki.syfjx.2016200006

CrossRef Full Text | Google Scholar

Yang, F. C., Xu, F., Wang, T. N., and Chen, G. X. (2021). Roles of vitamin A in the regulation of fatty acid synthesis. World J. Clin. Cases 9 (18), 4506–4519. doi:10.12998/wjcc.v9.i18.4506

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Q., Graham, T. E., Mody, N., Preitner, F., Peroni, O. D., Zabolotny, J. M., et al. (2005). Serum retinol binding protein 4 contributes to insulin resistance in obesity and type 2 diabetes. Nature 436, 356–362. doi:10.1038/nature03711

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, L. J., Wen, H., Zhou, S. H., Wu, F., and Jiang, M. (2019). Effect of dietary vitamin A on growth performance, serum biochemical indices and meat quality of GIFT, Oreochiromis niloticus. J. Northwest A F Univ. Nat. Sci. Ed. 47 (2), 16–24. doi:10.13207/j.cnki.jnwafu.2019.02.003

CrossRef Full Text | Google Scholar

Zhang, J. Y., Kothapalli, K. S., and Brenna, J. T. (2016). Desaturase and elongase-limiting endogenous long-chain polyunsaturated fatty acid biosynthesis. Curr. Opin. Clin. Nutr. Metabolic Care 19 (2), 103–110. doi:10.1097/MCO.0000000000000254

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Li, J., Tan, F. F., Mai, K., Zhang, C., Li, H., et al. (2015). Effects of different dietary vitamin A levels on growth and serum biochemical parameters for Japanese seabass (Lateolabrax japonicus). J. Fish. China 39 (1), 88–96. doi:10.3724/SP.J.1231.2015.59406

CrossRef Full Text | Google Scholar

Zoubek, M. E., Trautwein, C., and Strnad, P. (2017). Reversal of liver fibrosis: from fiction to reality. Best Pract. Res. Clin. gastroenterology 31 (2), 129–141. doi:10.1016/j.bpg.2017.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Monopterus albus, vitamin A, transcriptome, metabolome, lipid metabolism

Citation: Huo H, Hu C, Zhou Q, Xiong L and Peng M (2023) Integrated transcriptome and metabolome analysis reveals a possible mechanism for the regulation of lipid metabolism via vitamin A in rice field eel (Monopterus albus). Front. Physiol. 14:1254992. doi: 10.3389/fphys.2023.1254992

Received: 08 July 2023; Accepted: 07 August 2023;
Published: 23 August 2023.

Edited by:

Dong Han, Chinese Academy of Sciences (CAS), China

Reviewed by:

Min Jin, Ningbo University, China
Ye Yuan, Shantou University, China

Copyright © 2023 Huo, Hu, Zhou, Xiong and Peng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Qiubai Zhou, zhouqiubai@163.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.