Liver Transcriptomic Analysis of the Effects of Dietary Fish Oil Revealed a Regulated Expression Pattern of Genes in Adult Female Spotted Scat (Scatophagus argus)

Despite the significance of n-3 long-chain polyunsaturated fatty acid (n-3 LC-PUFA) in fish oil in promoting the maturation of female broodstocks, the detailed mechanism of the effect on the expression of hepatic reproduction-related genes is still unclear. In this study, transcriptome sequencing was used to analyze the effect of the higher dietary n-3 LC-PUFA level on gene expression in the liver of adult females spotted scat. Two-year-old female spotted scat (average weight, 242.83 ± 50.90 g) were fed with diets containing 8% fish oil (FO) or 8% soybean oil (SO) for 40 days. The fatty acid profile in the serum, liver, and ovary was analyzed, and high proportions of n-3 LC-PUFA were observed in the FO group. The final average fish body weight and gonadosomatic index were similar between the FO and SO groups. The serum vitellogenin (Vtg) and hepatosomatic index (HSI) of the FO group were significantly higher and lower than that of the SO group, respectively. Comparatively, the liver transcriptome analysis showed 497 upregulated and 267 downregulated genes in the FO group. Among them, the expression levels of three estrogen-regulated genes (i.e., Vtga, Vtgb, and Zp4) were significantly higher in the FO than in the SO group. This expression pattern could be related to the upregulation of Hsd17b7 (the key gene for the synthesis of liver steroid hormone) and the downregulation of the Hsp90 (the estrogen receptor chaperone). The expression levels of Foxo1a and Lep, which are involved in the lipid metabolism, decreased significantly in the FO group, which may be related to the lower level of HSI in the FO group. The genes related to liver LC-PUFA absorption and transport, Fabp2 and Mfsd2ab, were significantly upregulated in the FO group, indicating that fish actively adapt to high-fish-oil diets. In brief, high-fish-oil diets can influence the expression of genes related to liver n-3 LC-PUFA metabolism and reproduction, inhibit the accumulation of liver fat, and promote the liver health and gonad development. This study will contribute to clarifying the mechanism of dietary n-3 LC-PUFA on promoting reproductive development in teleost fish.


INTRODUCTION
The effects of essential fatty acids in feed on the reproduction of broodstocks are among the main research areas of artificial fish breeding in recent years. Long-chain polyunsaturated fatty acids (LC-PUFA) supplied in fish oil, including docosahexaenoic acid (DHA,, eicosapentaenoic acid (EPA, 20:5n-3), and arachidonic acid (ARA,, are essential nutrients for the normal growth and development of fish and can play key roles in the lipid metabolism, synthesis of sex steroid hormones, and reproductive regulation (Zhou et al., 2011;Glencross et al., 2015;Xu et al., 2017). For example, in yellowfin sea bream (Acanthopagrus latus) and flame angelfish (Centropyge loriculus), the relative fecundity, hatching rate, and larval survival rate were increased by the supplementing higher level of n-3 LC-PUFA in diet (Zakeri et al., 2011;Callan et al., 2012). When Japanese flounder (Paralichthys olivaceus) were fed with diets containing three different levels of n-3 LC-PUFA (i.e., 2.1, 4.8, and 6.2%), egg production increased in 6.2% n-3 LC-PUFA group. At the same time, the hatching rate and percentage of normal larvae were higher in 2.1% n-3 LC-PUFA group, indicating that excessive n-3 LC-PUFA in the diet would reduce the egg quality (Furuita et al., 2002). Consistently, in Plectorhynchus cinctus broodstock fed with four different levels of n-3 LC-PUFA, the moderate level of n-3 LC-PUFA improved the reproductive performance. In comparison, too high or too low n-3 LC-PUFA levels caused a reduction in reproductive performance (Li et al., 2004).
The liver is a key metabolic organ that is widely involved in the bioconversion of fatty acid and the regulation of lipid homeostasis in fish (Lazzarotto et al., 2018;Esmaeili et al., 2021). The liver is also the most important active tissue for gonadal development in oviparous vertebrates. In teleost fish, most of the yolk is synthesized and secreted by liver cells in the form of a protein precursor: the vitellogenin (Vtg) (Bemanian et al., 2004). In addition, some zona pellucida (Zp) genes encoding Zp glycoproteins, which constitute the protective coat surrounding oocytes, were also expressed in the teleost fish liver (Wu et al., 2018). Hence, the liver is a vital protein synthesis "factory" for the normal ovary development in fish. Previous studies have shown that the addition of fish oil to increase the dietary n-3 LC-PUFA level can promote the expression of Vtg in the liver of female fish, including the spotted scat (Scatophagus argus) (Zhang et al., 2013;Peng et al., 2015). However, the regulatory mechanism of the reproductive-related gene caused by the dietary n-3 LC-PUFA level in the liver of broodstock is still unclear, even though very pertinent in breeding programs.
In artificial breeding, the primary step is to get broodstocks with healthy oocytes and sperm in females and males, respectively. Also, the usual way to obtain mature fish with a better reproductive ability is by nutrient supplementation in the fish diet. Nutrigenomics utilizes high-throughput genomics approaches to unravel how dietary nutrients regulate gene expression at the mRNA/protein level and ultimately modulate cellular and organism metabolism. The most widely applied technique of nutrigenomics is the transcriptome analysis via next-generation sequencing (gene expression patterns at the mRNA level), also known as RNA-seq (Afman and Müller, 2012).
Studies of the transcriptional effects of LC-PUFA on liver gene expression in fish have been carried out in postsmolts Atlantic salmon (Salmo salar) (Glencross et al., 2015;Xue et al., 2020) and 2-month-old golden pompano (Trachinotus ovatus) (Lei et al., 2019). Nonetheless, evaluating the effects of n-3 LC-PUFA on livers of female broodstocks using the transcriptional method is lacking. This study focused on the hepatic transcriptome analysis to comprehensively understand the detailed molecular mechanisms of the effects of fish oil on reproductive-related genes in the liver of adult female spotted scat.
Spotted scat is an important food fish mainly produced in the South and South-East Asian countries. Due to its good nutrient quality, high protein content, and delicious taste, the fish has gained high demand and become an important aquaculture species. It is also a popular aquaria fish due to its colorful appearance, strong adaptability, and easiness to be cultured indoors (Gupta, 2016). For these reasons, studies focusing on its reproduction biology are on the rise, including reproduction endocrinology Zhang et al., 2018;Mandal et al., 2020), sex determination and differentiation (Mustapha et al., 2018He et al., 2019), and artificial breeding (Cai et al., 2010;Gandhi et al., 2014;Mandal et al., 2021). Its artificial breeding is successful in several countries, including China (Cai et al., 2010) and India (Mandal et al., 2021). To obtain matured spotted scat broodstock, appropriate nutritional factors should be maintained to enhance gonad development and induce reproduction (Cai et al., 2010). Although fish oil provides adequate nutrition for broodstock, the regulatory mechanism is unclear. Also, there is a delay in female gonad maturation compared with that of male (Gandhi et al., 2014). Therefore, it is important to advance studies on the effects of fish oil on the expressions of liver genes to gain insight into gene regulatory mechanisms in female spotted scat. The findings will provide essential knowledge for further functional studies on fish oil in promoting teleost fish reproduction.

Ethical Considerations
This experiment was undertaken in accordance with the guidance, care, and use of laboratory animals in China with approval from the Committee on the Ethics of Animal Experiments of Guangdong Ocean University.

Fish Rearing
Two-year-old female spotted scats were bought from the Zhuhai HengDa Cultivation Base (Guangdong, China). The sex ratio of spotted scats is approximately 1:1 . The female spotted scats were empirically selected from the mixed-sex fish by the sexual morphological difference in the head, anal fin, body color, body weight, and genital pore (Barry and Fast, 1992). The fish were transferred to the Donghai Island experimental base, Guangdong province, China. The fish were reared in a concrete pond (12 m × 5 m × 2 m) and adapted to the rearing condition for 2 weeks. Aeration was supplied to the tank through air stones, and the average water temperature was 31.1 ± 0.8 • C.
The salinity was 8 , while pH and dissolved oxygen ranged from 7.5-8.5 and 6.0-7.0 mg/L, respectively. Fish were fed 2.5% of their body weight with commercial fish pellets consisting of 43% crude protein and 12% crude lipid. Fish were fed two times a day at 06:00 and 18:00 h throughout the adaptation period. The pond water was not recirculated. Approximately, 20% of the water in the pond was drained off and replaced with fresh water once in 3 days.

Experimental Design
After the adaptation period, 100 females (average initial weights of 242.83 ± 50.90 g, and length of 19.48 ± 1.13 cm) were randomly distributed into two treatment groups (50 fish in each group) containing different levels of fish oil (0% and 8%; ingredients are listed in Table 1). The diet containing 8% fish oil was marked as the FO group, while the diet of the control group containing 8% soybean oil marked the SO group. Each treatment group was duplicated and contained 25 fish in each net cage (5 m × 3.5 m × 1.8 m). Two independent concrete ponds (12 m × 5 m × 2 m) were used; one concrete pond had two net cages for both FO and SO treatments. The fish were fed with the two diets for 40 days before sampling. The culturing condition was consistent with the adaptation stage.

Growth Measurements and Sampling
At the termination of the experiment, the spotted scats were sampled after 24 h starvation. Five fish per cage were sacrificed for sampling. The final number of fish, body weight, body length, gut weight, liver weight, and gonad weight were taken individually. The survival rate (SR), condition factor (CF), viscerosomatic index (VSI), hepatosomatic index (HSI), and gonadosomatic index (GSI) were calculated using the following formulas: SR (%) = final number of fish/initial number of fish × 100; CF = final body weight (g)/final body length (cm) 3 × 100; VSI (%) = final viscera weight (g)/final body weight (g) × 100; HIS (%) = final liver weight (g)/final body weight (g) × 100; GSI (%) = final gonad weight (g)/final body weight (g) × 100. The liver and ovary samples of the fish were excised and quickly frozen in liquid nitrogen. Blood samples were collected from the caudal vein of the fish. The blood samples were stored at 4 • C overnight and centrifuged (at 3,000 rpm for 10 min at 4 • C), and the upper layer serums were collected. All samples were subsequently stored at −80 • C for further analysis.
Fatty Acid Composition of the Experimental Diets, Fish Serum, Liver, and Ovary The total lipid from the ovary and liver (n = 6 per group) and feed were extracted using diethyl ether after acid hydrolysis, whereas lipid from the serum (n = 6 per group) was extracted with chloroform according to the procedure of the National Standards of the People's Republic of China (5009.168-2016). Extracted lipid was saponified with NaOH-methanol. Their constitutive fatty acids were methylated through the action of boron trifluoride-methanol, separated using Agilent 7890B gas chromatograph (Agilent Technologies, Santa Clara, CA, United States) on a capillary column (SPTM-2560; 100 m × 0.25 mm × 0.25 mm), and quantified using glyceryl triundecanoate as an internal standard.

RNA Extraction, Library Construction, and Sequencing
To evaluate the effects of fish oil on the genes expressed in the liver, the transcriptome analysis was carried out. First, total RNA (n = 6 per group) from the two groups was extracted from fish liver using RNeasy R Mini Kit reagents according to the instructions of the manufacturer (Qiagen GMBH, Hilden, Germany). The RNA quality was determined and quantified using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, United States) and checked using the RNase-free agarose gel electrophoresis. The mRNA was enriched by oligo (dT) beads; the mRNA was interrupted. The first strand of cDNA was synthesized in M-MuLV reverse-transcriptase system, and the second strand was synthesized from dNTPs in DNA polymerase I system. Then, the cDNA fragments were ligated to Illumina sequencing adapters. The ligation products were size selected using agarose gel electrophoresis, PCR amplified, and sequenced using novaseq 6000 by Guangzhou Genedenovo Biotechnology Co. (Guangzhou, China). All clean libraries of sequencing data were submitted to the NCBI Sequence Read Archive (SRA) database (Accession No: PRJNA765892).

Data Filtering, Reads Mapping, and Differential Gene Expression Analysis
To ensure that high-quality data are achieved, the original data were filtered before analyzing. The raw reads were filtered using fastp (version 0.18.0). The adapters read with more than 10% unknown nucleotides and over 50% low-quality nucleotides (Qvalue ≤ 20) were removed. The ribosome RNA (rRNA) reads were removed with the alignment tool Bowtie2 (version 2.2.8). Then, these clean reads were mapped to the spotted scat reference genome  with HISAT2.2.4. Genes with unknown transcripts or in intergenic regions were defined as novel genes. The fragments per kilobase of transcript per million mapped reads (FPKM) method was used to estimate the gene expression levels. Differentially expressed genes (DEGs) in two groups were identified using the EdgeR software package. In comparison, to avoid some important genes being filtered out, genes with | log2 (fold change)| > 1 and P < 0.05 were identified as significant DEGs. DEGs were annotated with the enrichment analysis of Gene Ontology (GO) functions and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The KEGG pathway with P < 0.05 was regarded as a significant enrichment pathway.

Validation by Real-Time PCR
To validate the gene expression profiles from the RNAseq analysis, 24 differentially expressed genes related to fat metabolism and ovarian development were selected for the realtime PCR analysis. Real-time PCR was performed using the PerfectStart TM Green qPCR SuperMix (TRAN, China). The PCR reaction was carried out in a total volume of 20 µl with 10 µl of 2 × PerfectStart TM Green qPCR SuperMix, 0.4 µl of each primer, and 2 µl cDNA. The thermal cycling program included an initial denaturation at 94 • C for 30 s, followed by 40 cycles of 94 • C for 5 s, 60 • C for 15 s, and 72 • C for 10 s. The reference gene, β-actin, was used to normalize the expression values. The 2 − Ct method was used to calculate the relative expression of the target genes. Each group contains 10 independent samples, and each sample was detected in triplicate. Primer sequences used for real-time PCR are listed in Supplementary Table 1.

Vitellogenin Protein Levels of Serum
The content of serum Vtg was detected using Grouper Vitellogenin ELISA Kit according to the instructions of the manufacturer (Cusabio, Wuhan, China). The serum samples from 10 fish in each group were analyzed.

Statistical Analysis
All data are expressed as means ± standard error (SE), and an independent sample t-test was used with the significance level at P < 0.05. The statistical analysis of the data was performed using the SPSS 17.0 software (SPSS Inc., Chicago, IL, United States).

Fatty Acid Composition of the Feed
The DHA, EPA, and ARA contents in FO diet are higher than that of the SO diet. On the contrary, the oleic acid (C18:1n-9) and linoleic acid (C18:2n-6) contents in FO diet are lower than that of the SO diet ( Table 2). The differences are mainly attributed to their different oil sources. 3 MUFA, mono-unsaturated fatty acids. 4 n-6 PUFA, n-6 polyunsaturated fatty acids. 5 n-3 PUFA, n-3 polyunsaturated fatty acids. 6 n-3 LC-PUFA, n-3 long-chain polyunsaturated fatty acid. 7 n-3 PUFA/n-6 PUFA, the ratio of n-3 polyunsaturated fatty acids to n-6 polyunsaturated fatty acids.

Survival Rate, Growth Performance, and Body Indices
At the end of the experiment, no fish died in both FO and SO groups. The final body weight, body length, CF, GSI, and VSI of Data are presented as means ± SE (n = 6). Values in the same row with different superscripts are significantly different (P < 0.05). 1 ND, non-detectable. 2 SFA, saturated fatty acids. 3 MUFA, mono-unsaturated fatty acids. 4 n-6 PUFA, n-6 polyunsaturated fatty acids. 5 n-3 PUFA, n-3 polyunsaturated fatty acids. 6 n-3 LC-PUFA, n-3 long-chain poly-unsaturated fatty acid. 7 n-3 PUFA/n-6 PUFA, the ratio of n-3 polyunsaturated fatty acids to n-6 polyunsaturated fatty acids.
the SO and FO fish are similar (P > 0.05, Table 4). The HSI of fish in FO group was significantly lower than that of the SO group (P < 0.05).

Raw Sequencing Reads and Quality Statistics
Twelve cDNA libraries were constructed from livers from the SO and FO groups. In this experiment, a total of 534,586,540 (i.e., 289,243,632 in the SO group and 245,342,908 in the FO group) clean reads were obtained from the transcriptome sequencing. The percentages of Q30 bases were more than 93% for all the samples, indicating a high-quality sequence (

Differentially Expressed Genes Verification Using Real-Time PCR
Twenty-four genes related to ovarian development and fat metabolism from RNA-seq were verified using real-time PCR (Figure 4). The results showed that the expression levels of the selected genes were consistent with RNA-seq, suggesting that the RNA-seq results are reliable.

Change of Serum Vitellogenin Levels
The mRNA expression levels of Vtga and Vtgb in livers of the FO group are higher than that of the SO group. The Vtg protein level in the serum was detected. The results showed that the serum Vtg level in the fish from the FO group was significantly higher than that of the SO group (P < 0.05; Figure 5). The results at the mRNA and protein levels are consistent.

Dietary n-3 Long-Chain Polyunsaturated Fatty Acids Effected Growth Performance
The growth of the fish was not affected by the oil source during the 40-day experimental stage. In sharpsnout seabream (Diplodus puntazzo), replacing fish oil with vegetable oil did not affect growth after 3 months of feeding (Piedecausa et al., 2007). However, in cobia (Rachycentron canadum), the proportional replacement of fish oil with soybean oil gradually increased (i.e., 0% SO, 33% SO, 67% SO, and 100% SO) and resulted in a gradual decrease in weight gain (Trushenski et al., 2011). The discrepancy could be related to the differences in fatty acid requirement among species (Xu et al., 2015). The current study suggests that spotted scat can be fed with a diet containing lower fish oil for normal growth. That is probably because spotted scat is omnivorous fish species. Of course, we could not conclude that the spotted scat does not need LC-PUPA for normal growth, as residual fish oil in fish meal could supply the LC-PUPA in the SO group. Decreased HSI in the FO group reveals that a high level of dietary n-3 PUFAs significantly reduced fat accumulation in the liver. It is well known that adult zebrafish suffer from hepatic adipose infiltration, which might influence reproduction and offspring (Virote et al., 2020). The liver with less fat tissue should be healthier, and the healthy liver might have a stronger ability to synthesize the proteins and lipids, such as the Vtg and Zp proteins required for gonadal development. In female spotted scat, the higher GSI values were observed in the reproduction season, indicating the maturing of the ovary (Gandhi et al., 2014). Consistently, the GSI in the FO group is relatively higher than that of SO group, although there was no statistical difference (P = 0.306). The possible reason for the large standard error between groups is the genetic background of the fish affected the regulation of lipid metabolism (Morais et al., 2011). These results suggest that a high level of dietary n-3 LC-PUFA might be beneficial to a healthy liver and promote ovarian development in female spotted scat.

Dietary n-3 Long-Chain Polyunsaturated Fatty Acids Affected Reproduction-Related Genes Expressions in Livers
The liver transcriptome analysis was carried out to confirm whether the high level of LC-PUPA would benefit the reproduction via influencing the liver genes. Vitellogenesis-and zonagenesis-related genes are induced by the estrogen in the liver to synthesize Vtg and zona radiata proteins (Zrp), respectively. Vtg and Zrp are secreted and transported from blood to the ovary and are taken up by the maturing oocytes (Arukwe and Goksøyr, 2003;Levi et al., 2009). Vtga, Vtgb, and Vtgc are all expressed increasingly during the ovary maturation from phase II to phase IV in the female liver of spotted scat (Cui et al., 2017). In this study, both transcriptome and real-time PCR data showed that Vtga and Vtgb are highly expressed in the FO group. Consistently, the serum Vtg protein is also higher in the FO group than that of the SO group. Zp4 is one of the genes controlling the making of Zp (in mammals) or vitelline envelope (in fish, amphibians, and birds), a glycoprotein layer surrounding the oocytes (Laine et al., 2019). In addition, Zp4 expression is significantly upregulated in the livers of FO fish. Taken together, the fish oil might promote ovary development by upregulating genes critical for vitellogenesis (i.e., Vtga and Vtgb) and zonagenesis (i.e., Zp4) in spotted scat.
The estrogen signaling pathway is mediated by the nuclear receptors. Three estrogen receptors, erα, erβ1, and erβ2, were characterized in many fish species, including spotted scat (Cui et al., 2017). Although, the three estrogen receptors are expressed in the female livers, their expressions were not affected by the high level of LC-PUPA in the presents study. In contrast, the liver estrogen level might be increased by the upregulation of the Hsd17b7. Hsd17b7 could transform estrone to estradiol in vitro (Ohnesorg et al., 2006). Similarly, fish oils upregulate 17β-hsd gene expression in the testicular tissues of rats (Mohammad et al., 2015). The regulatory mechanism of fish oil on Hsd17b7 expression remains to be elucidated.
Hsp90 is a ubiquitously expressed molecular chaperone of estrogen receptors and is involved in the signal transduction process of estrogen receptors (Chang et al., 2014). Herein, the Hsp90 expression was downregulated in the FO group. Whether less Hsp90 would either promote or inhibit the estrogen signaling pathway is still unclear. Considering the upregulation of Vtgs and Zp4 expressions in the FO group, the results suggested that less Hsp90 might enhance the estrogen signaling pathway. The possible mechanism is that decreasing Hsp90 would increase the number of free estrogen receptors. Those free receptors would combine with estrogen to be translocated into the nuclear to take action. The aryl hydrocarbon nuclear translocator (Arnt) has been found to be a coactivator of erα-and erβ-dependent transcription (Brunnberg et al., 2003). The Arntl2 gene has been considered a candidate gene in the molecular regulation of reproduction seasonality (Zhao W. et al., 2020). However, its expression was decreased in the FO group. So, the n-3 LC-PUPA might have a diverse influence on the expression of the genes involved in the network of the estrogen signaling pathway. The total output of the estrogen signaling pathway in the female liver should be enhanced in the group that contains higher level n-3 LC-PUPA.

Dietary n-3 Long-Chain Polyunsaturated Fatty Acids Affected Lipid Metabolism-Related Genes Expressions in Livers
The forkhead transcription factor Foxo1 is a critical regulator of hepatic lipid and glucose metabolism, and it functions in hepatic lipid homeostasis (Shin et al., 2012). The DHA supplemented diets reduced the expression of Foxo1 in the liver of pigs after 30-day feeding experiment (Chen et al., 2012). In contrast, EPA was more effective in reducing Foxo1 expression than DHA in rats after an 8-week feeding experiment (Hong et al., 2019). DHA or EPA alone or together could reduce the expression of Foxo1 in 3T3-L1 adipocytes in vitro (Martins et al., 2020). In mice, gain and loss of function proved that Foxo1 could promote the insulin sensitivity and lipid synthesis (Matsumoto et al., 2006). Foxo1a was highly expressed in insulin-responsive tissues (e.g., adipose tissue, liver, and muscle) in grass carp (Ctenopharyngodon idella) (Sun et al., 2017). In contrast, the addition of bovine insulin in vitro stimulates Vtg expression in hepatopancreatic explants from the mud crab (Scylla paramamosain) (Huang et al., 2017). The decreasing of the HIS in the FO group could be explained by the DHA/EPA in the fish oil decreasing the Foxo1a expression and reducing the lipid synthesis. Whether the increase of Vtg in FO group is affected by Foxo1-mediated insulin needs further verification. Leptin (Lep) is a protein hormone that regulates body weight and energy homeostasis (Paolucci et al., 2020). In fish, the liver seems to be the primary organ for Lep synthesis (Frøiland et al., 2010). In large yellow crookeder (Larimichthys crocea), Lep levels in blood were significantly lower in the FO group than in palm oil (Wang et al., 2019). Similarly, the liver Lep mRNA in the fish of FO group was significantly lower than that of SO group in the present study. The high level of LC-PUPA might decrease the Lep expression level in the liver in fish. However, there was a contrary relationship between liver lipid content and liver Lep mRNA expression in Arctic charr (Salvelinus alpinus) (Frøiland et al., 2010). The upregulation of the Lep mRNA expression in the SO group containing high HIS might be a negative feedback mechanism to regulate the liver lipid. The E2-induced Vtg mRNA transcript was influenced negatively by recombinant Lep in the rainbow trout (Oncorhynchus mykiss) in vitro (Paolucci et al., 2020). Thereafter, the induction of the Vtgs expression by higher n-3 LC-PUPA level might also depend on the decreasing of Lep.
The Ppar-related gene expression could be regulated by dietary fat that affects hepatic lipid storage and maintenance of metabolic health (Soni et al., 2015). The fatty acid-binding protein (Fabp) is the DEG enriched by the KEGG analysis in the PPAR signaling pathway (ko03320). The long-chain fatty acids are assisted by some fatty acid transporters to move through the cell membrane and are then transported by Fabps to the corresponding metabolic target sites (Lei et al., 2019). A Fabp2 −/− zebrafish decreased intestinal LC-PUPA and increased hepatic lipid accumulation compared with wild-type zebrafish. At the same time, transgenic overexpression Fabp2 significantly increased intestinal LC-PUPA . These results suggested that Fabp2 also promoted LC-PUPA absorption in fish. Consistently, Fabp2 expression was upregulated in the FO group, which contains high LC-PUFA than the SO group in the present study. In golden pompano (Trachinotus ovatus), the liver expressions of the Fabp -1, -4, and -6 were higher in the 2.1% LC-PUFA than that of 1.0% LC-PUFA feed group, suggesting that Fabps may be involved in LC-PUFA metabolism in this species (Lei et al., 2019). Taken together, Fabp gene family members may be essential for the LC-PUFA uptake and deposition in fish, while different subtypes of Fabp may be used for gene duplication and functional diversification in different species. In addition, the major facilitator superfamily domain-containing 2a (Mfsd2a) is confirmed as a transporter of DHA as lysophosphatidylcholine-DHA and is highly expressed in the liver (Pauter et al., 2017). Consistently, Mfsd2ab is also upregulated in the FO group in spotted scat. Therefore, the n-3 LC-PUFA supplementation affected the genes involved in LC-PUFA metabolism in spotted scat; hence, the regulator mechanism of these genes should be studied in the future.

CONCLUSION
The dietary supplementation of fish oil will significantly increase the deposition of n-3 LC-PUFA in spotted scat. The fish growth rates were similar between the FO and SO groups, while the HIS of the FO group was significantly lower than that of the SO group. The dietary fish oil supplementation increased the serum Vtg level, which is required for the ovary development. The liver transcription analysis showed that 764 genes were differentially expressed between the FO and SO groups. These data indicate that, the physiological response in the liver is associated with dietary oil. The expression of liver genes analyzed by RNA-seq gives us some new implications on how the dietary n-3 LC-PUFA supplied from the fish oil could promote the reproductive development of marine teleost.

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 below: https://www.ncbi.nlm. nih.gov/, PRJNA765892.

ETHICS STATEMENT
The animal study was reviewed and approved by Committee on the Ethics of Animal Experiments of Guangdong Ocean University.

AUTHOR CONTRIBUTIONS
TW and DJ designed and took part in the whole process of the experiment and wrote the draft of this manuscript. HS, SD, HC, and CZ co-conceived the experiment. ZL participated in the experiments. UM and WL revised the draft. GL supervised the experiments. All authors read and approved the final manuscript.