Effects of supplemental mixed bile acids on growth performance, body composition, digestive enzyme activities, skin color, and flesh quality of juvenile large yellow croaker (Larimichthys crocea) in soybean oil based diet

Now the replacement of fish oil (FO) with vegetable oils (VOs) has been broadly applied in aquatic feed, but studies reported that there were negative effects on growth, skin color, flavor and muscle texture. A 10-week feeding trial was conducted to investigate the effects of dietary mixed bile acids (BA) on growth performance, body composition, digestive enzyme activities, skin color and flesh quality of juvenile large yellow croaker (initial weight, 13.10 ± 0.18 g). Four iso-nitrogenous and iso-lipidic experimental diets were formulated and designated as soybean oil (SO), SO supplemented with 300 (BA300), 600 (BA600), and 1200 (BA1200) mg/kg bile acids. Two hundred and forty fish were randomly allocated into 12 floating net cages (1 m × 1 m × 1.5 m) that were located at marine fishing rafts. Each diet was randomly distributed in triplicate (three replicates per treatment, 20 fish per replicate). Fish were hand-fed to apparent satiation twice daily (05:30 and 17:30) for 10 weeks. Results showed that specific growth rate presented quadratic pattern with supplemental bile acids level, peaking at BA600 group. Meanwhile, feed conversion rate of cultured fish was significantly improved in BA600 group (P < 0.05). The redness (a*) of dorsal and lateral line skin and the yellowness (b*) of abdominal skin showed significantly quadratic pattern with the increase of supplemental BA level (P < 0.05), peaking at BA600 group. In terms of dorsal muscle texture, springiness showed a decreasing trend in significantly linear pattern with the increase of supplemental BA level, bottoming at BA600 group. However, cohesiveness and gumminess were significantly linear increased with the increase of supplemental BA level (P < 0.05). No significant differences were observed in lipase, amylase, and trypsin activities (P > 0.05). Total proportion of muscle n-3LCPUFA showed a linearly increasing trend with the increase of supplemental BA level. Muscle TG content was linearly increased with the increase of supplemental BA level (P < 0.05), and significant differences occurred as compared to fish fed diets with SO when supplemental BA level were reached 600 mg/kg and 1200 mg/kg (P < 0.05). The gene expression of acyl-CoA oxidase (aco) was significantly decreased in fish fed diets supplemented with 600 mg/kg bile acids compared with 300 mg/kg other than those fed diets with SO. In conclusion, these results suggested that mixed bile acids supplemented to SO-based diets could improve feed conversion ratio, skin color and flesh quality to some extent.

Now the replacement of fish oil (FO) with vegetable oils (VOs) has been broadly applied in aquatic feed, but studies reported that there were negative effects on growth, skin color, flavor and muscle texture. A 10-week feeding trial was conducted to investigate the effects of dietary mixed bile acids (BA) on growth performance, body composition, digestive enzyme activities, skin color and flesh quality of juvenile large yellow croaker (initial weight, 13.10 ± 0.18 g). Four isonitrogenous and iso-lipidic experimental diets were formulated and designated as soybean oil (SO), SO supplemented with 300 (BA300), 600 (BA600), and 1200 (BA1200) mg/kg bile acids. Two hundred and forty fish were randomly allocated into 12 floating net cages (1 m × 1 m × 1.5 m) that were located at marine fishing rafts. Each diet was randomly distributed in triplicate (three replicates per treatment, 20 fish per replicate). Fish were hand-fed to apparent satiation twice daily (05:30 and 17:30) for 10 weeks. Results showed that specific growth rate presented quadratic pattern with supplemental bile acids level, peaking at BA600 group. Meanwhile, feed conversion rate of cultured fish was significantly improved in BA600 group (P < 0.05). The redness (a*) of dorsal and lateral line skin and the yellowness (b*) of abdominal skin showed significantly quadratic pattern with the increase of supplemental BA level (P < 0.05), peaking at BA600 group. In terms of dorsal muscle texture, springiness showed a decreasing trend in significantly linear pattern with the increase of supplemental BA level, bottoming at BA600 group. However, cohesiveness and gumminess were significantly linear increased with the increase of supplemental BA level (P < 0.05). No significant differences were observed in

Introduction
In recent years, capture fisheries are gradually declining, which induces a steady increase in aquaculture production for meeting global demands for seafood (Watanabe, 2002;Turchini et al., 2010;Fawole et al., 2021). Thus, the demand for commercial aquatic feeds continue to increase. As an important lipid source in aquatic feeds, fish oil (FO) is abundant in n-3 polyunsaturated fatty acids, such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) (Yu et al., 2017), which are indispensable for growth and development of farmed marine animals. However, with the rapid expansion of aquaculture and the sharp increase of aquatic feed, the contradiction of FO between supply and demand becomes more and more prominent (Nasopoulou and Zabetakis, 2012). Compared with FO, vegetable oils (VOs) have many advantages including easy availability, relatively stable output and acceptable price (Bell et al., 2005;Tan et al., 2016), therefore they are considered as the best alternative lipid sources in aquafeed (Montero et al., 2015). Some studies reported that partial replacement of FO with VOs had no effects on growth and feed utilization in large yellow croaker (Duan et al., 2014), yellowtail (Aoki et al., 2000), Atlantic salmon (Bell et al., 2003), gilthead seabream (Menoyo et al., 2004;Fountoulaki et al., 2009), turbot (Regost et al., 2003b), cobia (Trushenski et al., 2011). However, high level VOs replacing FO had significantly negative effects on growth performance, body composition, fatty acid profiles, skin color, flavor and muscle texture in turbot (Regost et al., 2003a), gilthead seabream (Fountoulaki et al., 2009), red spotted grouper (Wang et al., 2016), large yellow croaker (Mu et al., 2018). Therefore, it is necessary to alleviate or eliminate the negative effects of FO substituted with VOs by nutrition regulation.
Bile acids are the main component of bile and synthesized in liver from cholesterol (Houten et al., 2006;Hofmann and Hagey, 2008). They have been well recognized to be able to emulsify lipids (Wang and Hartsuck, 1993;Alrefai and Gill, 2007;Romano et al., 2020), activate intestinal lipase (Adhami et al., 2017;Romano et al., 2020;Li et al., 2021b) and combine to long chain fatty acids to form a water-soluble complex (Einarsson et al., 1991) so as to improve lipids efficiency. Moreover, bile acids play vital role in emulsification, digestion and absorption of lipids. Meanwhile, they are considered as signaling molecules to regulate metabolism of triglyceride, cholesterol and maintain energy and glucose homeostasis (Insull, 2006;Nguyen and Bouscarel, 2008;Pols et al., 2011;Swann et al., 2011;Burrin et al., 2013;Di Ciaula et al., 2017;Molinaro et al., 2018). Previous studies in mammals and aquatic animals have shown that bile acids could effectively promote growth, improve lipids metabolism and reduce fat accumulation such as broiler (Lai et al., 2018;Ge et al., 2019), grass carp (Zhou et al., 2018), tilapia (Jiang et al., 2018), large yellow croaker (Du et al., 2017;Ding et al., 2020), largemouth bass (Yu et al., 2019;Yin et al., 2021) and improve body quality such as broiler (Lai et al., 2018), Wagyu heifer (Irie et al., 2011), grass carp (Zeng et al., 2017) and large yellow croaker (Yi et al., 2013). Given the functional properties of bile acids, dietary supplementary bile acids are gaining increasing concern in fish and providing the feasibility for improving the negative effects described above of FO replaced by VOs.
Large yellow croaker (Larimichthys crocea)is an important marine fish species with high economic value which has been widely farmed in China (Duan et al., 2014;Du et al., 2017;Ding et al., 2020;He et al., 2022). Previous studies indicated that VOs replacing FO had negative effects on growth, skin color and flesh quality etc. (Yi et al., 2013;Mu et al., 2018), which hindered the development of aquaculture industry and consumer market of large yellow croaker to some extent. In the previous study of our laboratory, supplemental chenodeoxycholic acid (CDCA) in diet replacing fish oil with soybean oil could significantly improve growth performance (Du et al., 2017). However, on the one hand, the price of CDCA is higher than mixed bile acids and it hasn't been extensively applied in aquaculture. On the other hand, no study was conducted on flesh quality of large yellow croaker in abovementioned study. Considering this, the purpose of the study was to comprehensively investigate the effects of the exogenous supplementation of mixed bile acids on growth performance, skin color and flesh quality of large yellow croaker.

Feed ingredients and diets formulation
Four iso-nitrogenous (42% crude protein) and iso-lipidic (12% crude lipid) experimental diets were formulated, including SO (formulation diet with soybean oil as the lipid source), BA300 (formulation diet with SO supplemented with 300 mg/kg mixed bile acids), BA600 (formulation diet with SO supplemented with 600 mg/kg mixed bile acids), and BA1200 (formulation diet with SO supplemented with 1200 mg/kg mixed bile acids). White fish meal, krill meal, and soybean meal were used as the main protein sources. Soybean oil was used as the main oil source. Ingredients were grinded to fine powder and mixed thoroughly to make pellets using pellet mill (EL-220, Haiyang City Huatong Machinery Co., Ltd., Shandong, China). After drying, diets were sieved into a proper size (diameter: 2 mm, length: 8 mm). Analysis of dry matter (105 C, 24 h), crude protein (Kjeldahl nitrogen × 6.25) and crude lipid (ether extraction by Soxhlet method) in experimental diets were performed according to the standard methods of Association of Official Analytical Chemists (AOAC, 2005). Feed formulation and proximate composition are presented in Table 1. All formulated diets were packed in separate bags and stored at -20°until further use.

Fish and experimental procedures
Hatchery produced juvenile large yellow croaker were purchased from Xiangshan Harbor Aquatic Seeds Company (Ningbo, China). Before the feeding trail, fish were reared for 2 weeks to acclimate the experimental diet and rearing condition. After acclimation, two hundred and forty fish (initial weight, 13.10 ± 0.18 g) were randomly allocated into 12 floating net cages (1 m × 1 m × 1.5 m) that were located at marine fishing rafts of Ningbo Marine Fishery Science and Technology Innovation Base. Each diet was randomly distributed in triplicate (three replicates per treatment, 20 fish per replicate). Fish were hand-fed to apparent satiation twice daily (05:30 and 17:30) for 10 weeks. During the feeding period, the water salinity, temperature, and dissolved oxygen were maintained at 25-28‰, 22-29°, 6.2-7.5 mg/L, respectively.

Sample collection
At the end of the experiment, all fish were fasted for 24 h and anesthetized with diluted MS-222 (Sigma-Aldrich, USA) before sampling. All fish were counted and weighed for measuring the survival rate、specific growth rate, respectively. Three fish from each cage were randomly sampled for the analysis of the body composition. The liver and viscera of six fish from each cage were weighed for measuring the hepatosomatic index and viscerosomatic index. Blood samples of three fish from each cage were obtained from the caudal vein using heparinized syringes and centrifuged at 3,000×g at 4°C for 10 min. The plasma was collected and stored at -80°C for subsequent analysis. Dorsal muscle and mid-intestine samples were flash-frozen in liquid nitrogen for determination of muscle quality and enzyme activities, respectively.

Body composition analysis
Proximate compositions of diets and whole fish were both analyzed by the standard methods of Association of Official Analytical Chemists (AOAC, 2005). Moisture was measured by drying at 105°C until constant weight in the oven. Crude protein (nitrogen × 6.25) was determined by the micro-Kjeldahl method using an Auto Kjeldahl System (FOSS KT260, Switzerland). Crude lipid was measured by solvent extraction using a Soxtec System (Soxtec System HT6; Tecator, Höganäs, Sweden). Ash was analyzed by combustion at 550°C for 4 h.

Skin color measurement
Skin color measurements were performed with a portable Minolta Chroma Meter CR-400 (Minolta, Osaka, Japan). According to CIE (CIE, 1976), the color parameters were L* for lightness, a* for redness or greenness, and b* for yellowness or blueness, respectively.

Dorsal muscle texture analysis
Texture analysis (hardness, adhesiveness, springiness, cohesiveness, gumminess, chewiness, resilience) was performed by a texture analyzer (TMS-PRO, FTC, USA). A texture profile analysis (TPA) was conducted using an 8-mm cylinder probe. The test parameters were set as follows: 0.1 N with the initial force, 30 mm/min of compression speed, and 60% of deformation of the original length .

Digestive enzyme activities assay
The intestine of fish was weighed (0.2 g) and then homogenized in 2 mL ice-cold PBS followed by centrifugation at 3,000×g for 10 min. The supernatant was subject to digestive enzyme activities assay. Lipase, amylase, and trypsin activities were measured by colorimetric enzymatic methods using commercial kits (Nanjing Jiancheng Bioengineering Institute, China). The total protein was measured by using commercial kits (Nanjing Jiancheng Bioengineering Institute, China).

Real-time quantitative polymerase chain reaction
In the present study, the sequences of primers were designed based on nucleotide sequences from GeneBank using software Primer Premier 5, and the method of total RNA extraction and complementary DNA (cDNA) synthesis were as described in the previous publication (Li et al., 2021a). qPCR measurements were performed in a total volume of 20 mL, containing 1 mL of each primer, 2 mL of cDNA product, 10 mL of SYBR-Premix ExTaqII (Takara, Japan) and 6 mL of RNA-free water and was carried out in a CFX Connect Real-Time System (Bio-Rad). The amplification procedure was employed as follows: 95°C for 2 min, followed by 39 cycles of 95°C for 10 s, 58°C for 10 s, and 72°C for 20 s. Primer of each gene for RT-qPCR was designed according to mRNA sequences of large yellow croaker, amplification efficiencies of all genes ranged from 0.95 to 1.05, and b-actin was chosen as the internal control. The expression level of genes was calculated and normalized using the 2 −DDCt method (Livak and Schmittgen, 2001). The primers were referred in previous study and listed in Table 2.

Calculations and statistical analysis
Where N t was the number of fish in each cage at the end of this experiment and N 0 was the number of fish in each cage at the beginning of this experiment.
Where W t was the final wet body weight (g), W 0 was the initial wet body weight and d was the experimental period in days.
One-way ANOVA was used for data statistics (SPSS 25.0, IBM, United States), and Tukey test was used for significance comparison between groups. Before Tukey's multiple range test, data were examined for normality and homogeneity of variances. Polynomial contrasts analysis was also performed in the pattern of linear and quadratic. P < 0.05 indicated significant difference. All results were expressed as the mean ± SEM.

Survival, growth performance and body index
No significant differences were observed in survival rate, final weight, specific growth rate, feed intake, hepatosomatic index, viscerasomatic index, and condition factor among different dietary treatments (P > 0.05), but specific growth rate showed quadratic pattern with supplemental mixed bile acids (BA) level, peaking at fish fed diets supplemented with 600 mg/kg (Table 3). Compared with fish fed diets with SO, FCR of those fed diets supplemented with 600 mg/kg BA was significantly improved (P < 0.05).

Whole body proximate composition and muscle TG content
Whole body proximate composition and muscle TG content of juvenile large yellow croaker are shown in Table 4. In terms of whole body proximate composition, no significant differences were observed in moisture, crude protein, crude fat, and crude ash among different dietary treatments (P > 0.05), but moisture exerted a decreasing trend in linear pattern with the increase of supplemental BA level. Muscle TG content was linearly increased with the increase of dietary BA level (P < 0.05), and significant differences occurred as compared with fish fed diets with SO when supplemental BA level were reached 600 mg/kg and 1200 mg/kg (P < 0.05), respectively.

Skin color
Skin color of juvenile large yellow croaker is shown in Table 5. No significant differences were observed in the lightness (L*) of dorsal, lateral line, and abdominal skin among different dietary treatments (P > 0.05). The b* of dorsal skin showed no significant difference among different dietary treatments. The redness (a*) of dorsal and lateral line skin and the b* of abdominal skin showed significantly quadratic pattern with the increase of supplemental BA level (P < 0.05), peaking at fish fed diets supplemented with 600 mg/kg. The b* of lateral line skin was significantly linearly decreased with the increase of supplemental BA level (P < 0.05), and significant differences occurred in fish fed diets supplemented with 600 mg/kg and 1200 mg/kg BA compared with those fed diets with SO (P < 0.05). The a* of abdominal skin showed both significantly linear and quadratic decreasing pattern with the increase of supplemental BA level (P < 0.05), and it was significantly higher in fish fed diets supplemented with 1200 mg/kg BA compared with those supplemented with 300 mg/kg or 600 mg/kg (P < 0.05).

Dorsal muscle texture
Muscle texture of juvenile large yellow croaker fed with different diets is shown in Table 6. No significant differences were observed in hardness, adhesiveness, chewiness, and resilience among different dietary treatments (P > 0.05). Springiness showed a significantly decreasing trend in significantly linear and quadratic pattern with the increase of supplemental BA level, bottoming at fish fed diets supplemented with 600 mg/kg. However, cohesiveness and gumminess were significantly linear increased with the increase of supplemental BA level (P < 0.05), and significant differences were found in the fish fed diets supplemented with 1200 mg/kg BA compared with those fed diets with SO and 300 mg/kg BA (P < 0.05).

Plasma metabolites
Effect of dietary BA on plasma metabolites of juvenile large yellow croaker is shown in Figure 1

Muscle fatty acids profiles
Fatty acids profile in the muscle of large yellow croaker is measured using a HP6890 gas chromatograph (Table 7). Muscle fatty acids profile showed no statistical difference among all treatments (P > 0.05). However, Total proportion of n-3LCPUFA showed a linearly increasing trend with the increase of supplemental BA level.

Real-time quantitative polymerase chain reaction
Gene expression of cd36, fas, scd, and cpt1 showed no statistical difference among all treatments (Figures 3A-D) (P > 0.05). However, gene expression of aco was significantly decreased in fish fed diets supplemented with 600 mg/kg BA compared with 300 mg/kg other than those fed diets with SO ( Figure 3E) (P < 0.05).

Discussion
In the present study, feed conversion rate (FCR) was significantly improved with the supplementation of mixed bile acids (BA) compared with fish fed diets with soybean oil (SO). The optimum dosage was 600 mg/kg and not in agreement with a previous finding in juvenile large yellow croaker (optimum dosage: 300 mg/kg) (Du et al., 2017;Ding et al., 2020). This discrepancy might be attributed to the differences of feed formations or ingredients of bile acids. It was reported that dietary bile acids could effectively increase feed utilization (Yamamoto et al., 2007;Jiang et al., 2018) and improve the growth performance of fish including largemouth bass (Micropterus salmoides) (Yu et al., 2019), rainbow trout (Oncorhynchus mykiss) (Yamamoto et al., 2007), and large yellow croaker (Larimichthys crocea) (Du et al., 2017). The reason why bile acids are important in improving FCR is that it can effectively improve the emulsification and transportation of lipids (Romanski, 2007).
Lipase activity can be recognized as an indicator of lipid utilization for animals. Previous studies showed dietary bile acids supplementation improved lipase activity as described in European eels (Zhai et al., 2020), rainbow trout (Adhami et al., 2017), large yellow croker (Ding et al., 2020), tilapia (Jiang et al., 2018), tongue sole (Li et al., 2021b), prenant's schizothoracin (Zeng et al., 2016;Xiang et al., 2019). Generally, bile acids are composed of a convex hydrophobic face and a concave hydrophilic face, this amphipathic structure is facilitated to lipid emulsification and breaks down large lipid molecules into small globules (Jiang et al., 2018), which provides lipase a higher surface area to digest lipids by catalyzing the hydrolysis of ester bonds (Romano et al., 2020). However, in this study, although lipase activity showed an upward trend, the difference was not significant. So we could attribute the improving FCR to the promoting emulsification of exogenous bile acids.
From the perspective of mixed bile acids, the main compositions are HDCA and CDCA. However, these components have seldom been studied in aquatic animals. Liao et al. (2020) reported that dietary bile acid (HDCA=69.9%) supplementation did not affect the growth performance of juvenile tiger puffer fed normal or high-lipid diets. Yin et al. (2021) found that dietary CDCA attenuated the adverse effects induced by high-fat diet on growth performance of largemouth bass, they attributed it to improving intestinal health, which conduced to the absorption and utilization of lipids. Du et al. (2017) reported that growth was significantly increased by dietary CDCA supplementation in large yellow croaker fed a soybean oil diet, they speculated that dietary CDCA might improve the growth performance by reducing inflammatory and improving health. Further studies are needed to explore the underlying mechanisms.
In addition to the improving feed conversion ratio, bile acids could also improve skin color and flesh quality of large yellow croaker in this study. Skin color is adopted as an index for assessing the quality of large yellow croaker, because it can directly affect the purchase of customers (Jiang et al., 2021). In this study, compared with fish fed diets with SO, supplementation of 600 mg/kg BA can significantly improve the redness (a*) of dorsal skin and the a* of lateral line skin. As a kind of fish characterized by yellow skin, the increase of skin color, especially yellowness (b*) of abdominal skin of large yellow croaker would improve commercial value. Moreover, muscle texture is also an important index for assessing the flesh quality of large yellow croaker. In the present study, the springness was decreased in fish  Li et al. 10.3389/fmars.2023.1149887 Frontiers in Marine Science frontiersin.org fed diets supplemented with BA compared with those fed diets with SO. The reason should be attributed to the increase of triglyceride in muscle, this was in consistent with previous study in largemouth bass (Geng et al., 2018). Meanwhile, the increase of cohesiveness and gumminess could contribute to chewiness, suggesting that the taste of fish muscle was more tender, palatable and not sticky. As to the reason for the improvement of skin color and muscle texture, it was speculated that bile acids promoted the absorption of lipids and lipid-soluble substances such as pigments (Insull, 2006;Islam et al., 2011;Yang et al., 2022). This process was achieved through solubilizing dietary lipids and lipid-soluble substances as micelle by bile acids to facilitate their easily diffuse across the intestinal mucosa (de Aguiar Vallim et al., 2013). Synthesis and catabolism of triglyceride are very pivotal for lipid accumulation in a specific tissue (Li et al., 2013). The present study found that triglyceride (TG) content in muscle was markedly increased when supplemented 600 mg/kg and 1200 mg/kg bile acids. Similarly, previous studies have verified that dietary bile acids play an important role in regulating triglyceride metabolism (Ding et al., 2020). Consequently, we examined the gene expression related to fat metabolism in muscle. ACO is a key enzyme of fatty Values are represented as mean ± SEM (n = 3). Li et al. 10.3389/fmars.2023.1149887 acid oxidation, the gene expression level of aco in muscle of fish fed diets 600 mg/kg BA was down-regulated in the present study. Therefore, we inferred that supplemented dietary bile acids enhanced TG deposit in muscle by inhibiting fatty acids oxidation, and thereby improve the flesh quality of large yellow croaker. Inclusion of bile acids could activate FXR (Du et al., 2017). FXR plays a central role in the regulation of bile acid synthesis, excretion and transport (Chiang, 2009), and fxr expression has been reported to reduce bile acids synthesis by inhibiting the expression of cyp7a1 in the liver (Zhou et al., 2018). Moreover, the activation of FXR could improve ppara expression and promote the b-oxidation of fatty acids (Xu et al., 2014). Therefore, in this study, the reason for lack of significant change in total bile acids (TBA) levels could be due to the downregulation of cyp7a1 expression and up-regulation of ppara expression, resulting in a reduction of BA synthesis and increase of fatty acids oxidation. The fatty acids profile of diets could affect the fatty acids composition of fish muscle, which was well documented for gilthead sea bream (Regost et al., 2003a;Grigorakis, 2007) and large yellow croaker . The substitution of fish oil by vegetable oils affected the fatty acids composition of fish, leading to increase of n-6PUFA and decrease of n-3PUFA in muscle and affecting flesh quality (Mu et al., 2018), In present study, compared with fish fed diets with SO, the level of total n-3LCPUFA showed a linearly increasing trend with the increase of supplemental bile acids level, which was consistent with study in large yellow croaker (Ding et al., 2020). n-3LCPUFA, especially, DHA and EPA play very important role in improving immunity, decreasing the risk of cardiovascular disease . Therefore, promoting the increase of n-3LCPUFA by dietary bile acids inclusion will be beneficial for health in human.

Conclusion
In summary, the results of the present study showed that dietary 600 mg/kg mixed bile acids supplementation could improve feed conversion ratio, skin color and muscle texture of juvenile large yellow croaker. These findings would contribute to the long-term sustainable development of large yellow croaker farming. We will 3 n-6 PUFA, n-6 poly-unsaturated fatty acids. 4 n-3 PUFA, n-3 poly-unsaturated fatty acids. 5 n-3PUFA/n-6PUFA, the ratio of n-3 poly-unsaturated fatty acids to n-6 poly-unsaturated fatty acids. 6 n-3 LC-PUFA, n-3 long chain poly-unsaturated fatty acids. 7 PUFA, poly-unsaturated fatty acids.
intensively study the effects of bile acids on large yellow croaker through more comprehensive work such as adding fish oil group as control, increasing the quantity of sampled fish or selecting other different bile acids etc. in the future.

Data availability statement
The original contributions presented in the study are included in the article/supplementary materials, further inquiries can be directed to the corresponding author.

Ethics statement
The animal study was reviewed and approved by Committee on the Ethics of Animal Experiments of Ocean University of China. Written informed consent was obtained from the owners for the participation of their animals in this study.  Gene expression of lipid metabolism (A-E) in the muscle of juvenile large yellow croaker in different treatments. Values are represented as mean ± SEM (n = 3). Bars with different lowercase letters represent significant difference (P < 0.05). cd36, cluster of differentiation 36; fas, fatty acid synthase; scd, stearoyl-coenzyme A desaturase; cpt1, carnitine palmitoyl transferase 1; aco, acyl-CoA oxidase.

Conflict of interest
Author JW was employed by Shandong Longchang Animal Health Product Co., Ltd.
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.