Next Article in Journal
The Use of Removed Mesquite Brush as a Fiber Replacement in Silage Production
Next Article in Special Issue
Effects of Bee Pollen on Growth Performance, Intestinal Microbiota and Histomorphometry in African Catfish
Previous Article in Journal
Measuring Biochemical Variables and Serum Amyloid A (SAA) in Working Mules in Central Italy
Previous Article in Special Issue
Effect of Tributyrin on Growth Performance and Pathway by which Tributyrin Regulates Oligopeptide Transporter 1 in Juvenile Grass Carp (Ctenopharyngodon idellus)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 12
Issue 20
10.3390/ani12202794
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Dietary Supplementation with Aurantiochytrium sp. on Zebrafish Growth as Determined by Transcriptomics

by
Hao Yang
1,2,
Yanlin Huang
1,
Zhiyuan Li
1,
Yuwen Guo
1,
Shuangfei Li
3,
Hai Huang
2,*,
Xuewei Yang
3,
Guangli Li
1 and
Huapu Chen
1,2,4,*
1
Guangdong Research Center on Reproductive Control and Breeding Technology of Indigenous Valuable Fish Species, Guangdong Provincial Key Laboratory of Pathogenic Biology and Epidemiology for Aquatic Economic Animals, Fisheries College, Guangdong Ocean University, Zhanjiang 524088, China
2
Key Laboratory of Utilization and Conservation for Tropical Marine Bioresources of Ministry of Education, Hainan Key Laboratory for Conservation and Utilization of Tropical Marine Fishery Resources, Yazhou Bay Innovation Institute, Hainan Tropical Ocean University, Sanya 572022, China
3
Shenzhen Key Laboratory of Marine Bioresource and Eco-Environmental Science, College of Life Sciences and Oceanography, Shenzhen University, Shenzhen 518060, China
4
School of Biomedical Sciences, The Chinese University of Hong Kong, Shatin, Hong Kong 999077, China
*
Authors to whom correspondence should be addressed.
Animals 2022, 12(20), 2794; https://doi.org/10.3390/ani12202794
Submission received: 4 September 2022 / Revised: 7 October 2022 / Accepted: 14 October 2022 / Published: 16 October 2022
(This article belongs to the Special Issue Early Development and Growth of Fishes)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

The high demand and price of fish oil and fishmeal, combined with the global strain on fish stocks, has attracted attention to environmentally friendly and sustainable microalgae. In this project, we fed zebrafish by feed supplementation for 56 days and found that Aurantiochytrium sp. extract did have a growth-promoting effect on zebrafish, with enhanced lipid metabolism and improved antioxidant capacity, reducing oxidative damage from lipid oxidation.

Abstract

The marine protist Aurantiochytrium produces several bioactive chemicals, including EPA (eicosapentaenoic acid), DHA (docosahexaenoic acid), and other critical fish fatty acids. It has the potential to improve growth and fatty acid profiles in aquatic taxa. This study evaluated zebrafish growth performance in response to diets containing 1% to 3% Aurantiochytrium sp. crude extract (TE) and single extract for 56 days. Growth performance was best in the 1% TE group, and therefore, this concentration was used for further analyses of the influence of Aurantiochytrium sp. Levels of hepatic lipase, glucose-6-phosphate dehydrogenase, acetyl-CoA oxidase, glutathione peroxidase, and superoxide dismutase increased significantly in response to 1% TE, while malic enzyme activity, carnitine lipid acylase, acetyl-CoA carboxylase, fatty acid synthase, and malondialdehyde levels decreased. These findings suggest that Aurantiochytrium sp. extract can modulate lipase activity, improve lipid synthesis, and decrease oxidative damage caused by lipid peroxidation. Transcriptome analysis revealed 310 genes that were differentially expressed between the 1% TE group and the control group, including 185 up-regulated genes and 125 down-regulated genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analyses of the differentially expressed genes revealed that Aurantiochytrium sp. extracts may influence liver metabolism, cell proliferation, motility, and signal transduction in zebrafish.

1. Introduction

Over the last two decades, global aquaculture production more than tripled. However, farmed fish oil production plummeted by a third, while prices more than doubled [1]. Fish oil is crucial for aquaculture because it contains critical fatty acids for fish growth, particularly n-3 highly unsaturated fatty acids (HUFA) [2]. Fish lack specific desaturases for the synthesis of n-3 and n-6 HUFA from short-chain fatty acids. They can only convert linoleic acid (LA) and linolenic acid (ALA) into 20-22C HUFA [3,4]. However, as a conventional source of unsaturated fatty acids, the extraction of n-3 HUFA from fish oil is costly and inefficient. As an alternative, marine microalgae are becoming well-known owing to their ability to synthesize long-chain unsaturated fatty acids [5].
Aurantiochytrium is a marine protist that produces carotenoids, squalene, sterols [6,7,8,9,10,11,12,13], and other terpenoids in addition to long-chain unsaturated fatty acids, such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) [14,15,16]. Aurantiochytrium has a wide range of applications in aquaculture. For example, nutritional supplementation with Aurantiochytrium sp. increased juvenile Nile tilapia (Oreochromis niloticus) growth at suboptimal low temperatures (22 °C), improved the ratio of unsaturated fatty acids, and was beneficial for human consumption [17]. Guimares et al. observed that substituting Aurantiochytrium sp. meal for fish oil improves shrimp development and feed conversion and increases shrimp muscle DHA levels. Some studies have indicated that Aurantiochytrium sp. bioactive compounds can increase the quality of aquatic goods, improve meat quality, reduce organic pollutant deposition, and enhance disease resistance [18,19,20]. Most research on Aurantiochytrium sp. as a dietary supplement has focused on its growth-promoting and other biological effects; little is known about the mechanism by which it stimulates fish growth [17,21,22,23,24].
Unsaturated fatty acids are the primary components of Aurantiochytrium sp., and the liver is crucial in lipid metabolism, oxidation, and HUFA production [25]. Transcriptome sequencing technology is widely used in studies of the molecular basis of fish growth [26,27]. It can be used for comprehensive analyses of the type, structure, and expression levels of all transcripts in a tissue or cell under different conditions, revealing the molecular regulatory mechanisms underlying specific biological processes.
Therefore, zebrafish (Danio rerio) was utilized as a research model to assess the influence of treatment with dietary Aurantiochytrium sp. for 56 days on zebrafish growth performance during the critical period. The mechanism by which Aurantiochytrium sp. affects zebrafish growth was investigated by enzyme-linked immunosorbent assay (ELISA) and transcriptomics to detect changes in liver biochemical indicators and genes.

2. Materials and Methods

2.1. Sample Source and Diet Preparation

Shenzhen University provided the Aurantiochytrium sp. Szu445 extract and the culture and extraction methods were described previously [28]. Table S1 lists the detailed composition of the Aurantiochytrium extract (determined by GC/MS). Various amounts of Aurantiochytrium sp. Szu445 extract as well as its single extract were incorporated into the diet, established based on prior exploratory experiments: 0% w/w (Control), 1% w/w TE (T1), 2% w/w TE (T2), and 3% w/w TE (T3); 1% w/w DHA, 2% w/w DHA, and 3% w/w DHA; 1% w/w EPA, 2% w/w EPA, and 3% w/w EPA. In the above amounts, extracts were added to commercial feed, which was then crushed and thoroughly mixed in a pelletizer to make feed pellets. The feed pellets were spread out and dried for 72 h at 30 °C. The feed was stored at −20 °C in a refrigerator. The commercial feed contained 50% crude protein, 10% crude fat, 88% dry matter, 16% ash, and 3% crude fiber.

2.2. Zebrafish Treatment and Sample Collection

The China Zebrafish Resource Center supplied AB-type wild-type zebrafish, cultivated in a zebrafish culture system (Shanghai Haisheng Biological Experiment Equipment Co., Ltd., Shanghai, China). The photoperiod was maintained at 14:10, the temperature was maintained at 28 ± 0.5 °C, the pH was between 7.0 and 8.0, and the dissolved oxygen was held between 6 and 7 mg/L by the aeration system. After parental breeding, one-month-old wild-type (AB) zebrafish were randomly chosen as experimental fish. Before the experiment, the body length and weight were measured (body weight, 0.007 0.001 g; body length, 0.90 0.029 cm).
The 600 zebrafish were randomly allocated into ten groups of three tanks (10 L) each, each tank comprising 20 zebrafish. Fish were fed 5% of their body weight daily (9 am and 4 pm). Throughout the eight-week feeding experiment, the fish were weighed every two weeks to adjust the feed weight. During the 56-day rearing experiment, nine zebrafish were chosen at random every two weeks to measure growth indices, such as mean body length and mean body weight. After 56 days, fish were selected at random, anesthetized with 20 mg/L eugenols, and then euthanized on ice. The samples were snap-frozen and kept at −80 °C until RNA extraction and biochemical analysis. Guidelines for the Care and Use of Laboratory Animals were meticulously adhered to in all animal handling procedures. This protocol was approved by the Animal Research and Ethics Committee of Guangdong Ocean University (NIH Pub. No. 85-23, revised 1996).

2.3. Detection of Liver Biochemical Indicators

Biochemical indicators tested in this study included hepatic amylase (AMS), glucose-6-phosphate dehydrogenase (G6PD), malic enzyme (ME), pyruvate kinase (PK), acetyl-CoA oxidase (ACO), hepatic lipase (LPS), fatty acid synthase (FAS), carnitine lipid acyltransferase (CACT), acetyl-CoA (ACC), superoxide dismutase (SOD), glutathione peroxidase (GSH-Px), catalase (CAT), and malondialdehyde (MDA). The liver samples were carefully weighed, diluted tenfold with phosphate-buffered saline (PBS) (Solarbio, Beijing, China), and homogenized with a tissue crusher in an ice-water bath. After homogenization, samples were centrifuged for 10 min at 4 °C with a 3000× g centrifuge, and the supernatant was kept. These parameters were assayed using a microplate reader according to the manufacturer’s instructions (Multiskan MK3, Thermo Fisher Scientific, Waltham, MA, USA). Product numbers can be found in Table S2.

2.4. RNA-Seq and Bioinformatics Analysis

For RNA-seq library preparation, zebrafish total RNA was extracted. Each group utilized three replicates, with each replicate including the livers of three fish. RNA was isolated using TRIzol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The NanoDrop 2000 was utilized to evaluate the concentration and purity of RNA (Thermo Fisher Scientific). RNA integrity was determined using the Agilent Bioanalyzer 2100 System and the RNA Nano 6000 Assay Kit (Agilent Technologies, Santa Clara, CA, USA). According to the manufacturer’s instructions, the NEBNext Ultra RNA Library Prep Kit for Il-lumina (NEB, Ipswich, MA, USA) was used to construct sequencing libraries, and index codes were added to assign sequences to each sample. The index-coded samples were clustered using a cBot Cluster Generation System equipped with a TruSeq PE Cluster Kit v4-cBot-HS (Illumina). The libraries were sequenced on the Illumina platform, and paired-end reads were generated after cluster creation. After cluster creation, libraries were sequenced on the Illumina platform. Perl scripts were used to remove adaptors and determine the Q20 and Q30 values, GC-content, and duplication rates in order to generate clean reads. The clean reads were then aligned with the genome sequence reference (GRCz11 release100). Using the reference genome as a guide, only reads containing a perfect match or one mismatch were annotated. HISAT2 was applied in order to map readings to the reference genome. Using RSEM v1.2.21, fragments per million reads per kilobase (FPKM) were utilized to assess gene expression levels [29]. Using the GOseq R package [30] and Wallenius’ non-central hypergeometric distribution, differentially expressed genes (DEGs) were analyzed for Gene Ontology (GO) enrichment. Using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and KOBAS [31], an enrichment analysis of pathways was performed.

2.5. Real-Time Quantitative PCR (RT-qPCR) Validation

TRIzol was used to extract total RNA from liver samples according to the manufacturer’s instructions (Invitrogen). Using a Nanodrop 2000 spectrophotometer, RNA sample concentration was measured (Thermo Scientific, Waltham, MA, USA). According to the manufacturer’s instructions for the RevertAid First-strand cDNA Synthesis Kit, the same quantity of RNA was used for the reverse transcription of each sample (Thermo Scientific). LightCycler 480 (Roche, Basel, Switzerland) and SYBR Premix Ex Taq II were utilized for RT-qPCR (TaKaRa Bio Inc., Shiga, Japan). The protocol for the thermal cycler was 95 °C for 30 s, followed by 40 cycles of 95 °C for 5 s, 56 °C for 20 s, and 72 °C for 20 s. Biological triplicates and three method repetitions were performed on each sample. The reference gene β-actin was employed as an internal control to standardize mRNA levels [32]. Table S3 offers a comprehensive list of primer pairs. Using the 2-ΔΔCt method, the relative gene expression levels were calculated [33].

2.6. Calculation of Growth Performance and Statistical Analysis

The following formulae were used to calculate growth performance:
S R   % = N t / N 0 × 100
R W G   % = W t − W 0 / W 0 × t × 100
S G R   % = ln W t − ln W 0 / t × 100
C F   % = W t / L t 3 × 100
where Nt: final number of fish, N0: initial number of fish, Wt: average fish weight after feeding (g), W0: average fish body weight at the start of feeding (g), Lt: average fish length after feeding (cm), t: number of days of feeding (d), SR: survival rate, RWG: relative growth rate for body weight, SGR: specific growth rate, and CF: condition factor.
All data were presented as mean ± SEM (standard error of mean). IBM SPSS Statistics 24.0 was utilized for the data’s one-way ANOVA (SPSS Inc., Chicago, IL, USA). At the 95% confidence level (p < 0.05), the least significant difference (LSD) test was employed to determine whether or not there were significant differences between group means.

3. Results

3.1. Effect of Dietary Aurantiochytrium sp. Crude Extract and Single Extract on Growth Performance in Zebrafish

Figure 1 depicts the growth of zebrafish in various groups. Before 14 days, there were no significant differences in the changes in body weight (Figure 1A) or body length (Figure 1B) among groups of zebrafish. After 28 days, zebrafish in the Aurantiochytrium sp. crude extract (TE) groups showed better growth rates, in terms of both body weight and body length, than those of the single extract group and control group.
The addition of Aurantiochytrium sp. crude extract and its single extract to the diet had no significant influence on the CF of zebrafish (p > 0.05), as shown in Table 1. However, the relative weight increase rates were 20.21% and 20.27% greater in fish fed diets supplemented with 1% and 2% TE than in the control group, respectively (p < 0.01). The 3% TE group did not show a significant difference with the control in growth performance, which may be due to lipid peroxidation due to the high lipid content. Similar results were obtained when using a high concentration in the preliminary experiment. There were no statistically significant differences between the single extract groups with 1% EPA, 2% EPA, and 2% DHA and the control group. It is worth mentioning that 1% DHA and 3% EPA had inhibitory effects on growth. The experimental results suggest that a low concentration of Aurantiochytrium sp. extract has a more significant impact on zebrafish growth than that of a single extract. As a result, 1% TE was chosen for further analyses.

3.2. Effects of 1% Aurantiochytrium sp. Extract on Biochemical Indicators in the Zebrafish Liver

The effects of Aurantiochytrium sp. extract on digestive enzyme activity, energy metabolism, lipid metabolism, and antioxidant levels in the liver were evaluated. With respect to digestive enzyme activity (Figure 2A), there was no significant difference in AMS activity, and LPS activity was significantly higher in the experimental group than in the control group (p < 0.01). In analyses of energy metabolism (Figure 2B), G6PD levels were significantly higher in the experimental group than in the control group (p < 0.01), while the opposite results were obtained for ME (p < 0.01), and PK showed no difference between groups. In addition, ACO activity was significantly higher in the group fed 1% TE, and the activities of CACT (p < 0.05), ACC (p < 0.01), and FAS (p < 0.01) were higher in the control group than in the experimental group. The effects of TE on antioxidant activity (Figure 2D) were also highly significant. MDA (p < 0.05) and CAT levels (p < 0.01) were significantly higher in the control group than in the 1% TE group, while GSH-Px and SOD levels were significantly higher in the 1% TE group than in the control group (p < 0.01).

3.3. RNA-Seq of the Liver

The Illumina HiSeq high-throughput sequencing platform was used to study the zebrafish liv transcriptomes of different groups after feeding for 56 days (C, control group; T, 1% Aurantiochytrium sp. extract group). Following quality control, 41.60 Gb of clean data was collected, with Q30 values (i.e., defined as the proportion of sequences with a sequencing error rate of less than 0.1%) for each sample of ≥93.81%. The GC content was approximately 48.19–48.57%. A detailed summary of sequencing statistics was provided in Table 2.
The zebrafish genome was utilized as a reference to generate a sequence alignment and in subsequent analyses. Table 3 summarizes the results of the HISAT2 comparative analysis. The mapping rates to the reference genome were between 90.49% and 92.15%. The data were high quality and can be utilized to investigate DEGs further. The reads from the StringTie alignment were used for assembly and quantification.

3.4. Differential Gene Expression Analysis

Gene expression levels in the zebrafish liver differed significantly between the control group and the 1% Aurantiochytrium sp. extract group. There were 185 up-regulated genes and 125 down-regulated genes in the experimental group. Each gene was clustered according to the difference in expression (q-value) and fold change in expression, and a volcano plot was generated for visualization (Figure 3).
Sixteen DEGs were selected for validation by RT-qPCR (Figure 4A). Expression levels obtained by RNA-Seq and RT-qPCR were significantly correlated (R2 = 0.96088), confirming the reliability and accuracy of gene expression levels quantified by the transcriptome analysis (Figure 4C). A hierarchical clustering analysis of the screened DEGs was performed, assigning genes with the same or similar expression patterns to clusters, as shown in Figure 4B. The repeatability within groups was good, and the overall difference between groups was significant.

3.5. GO and KEGG Pathway Enrichment Analyses

The DEGs were subjected to a GO functional enrichment analysis. DEGs were allocated to terms within the three major functional categories of cellular component (CC), molecular function (MF), and biological process (BP).
Centrosome (GO: 0005813), extracellular space (GO: 0005576), and methionine adenosyltransferase complex (GO: 0048269) have especially substantial enrichment coefficients in CC (Figure 5A). Rho-dependent protein serine/threonine kinase activity (GO: 0072518), heme binding (GO: 0020037), and GTP-Rho binding (GO: 0017049) had the highest enrichment coefficients in MF (Figure 5B). In addition, the most significant GO secondary terms enriched by BP (Figure 5C) were circadian regulation of gene expression (GO: 0032922), cortical actin cytoskeleton organization (GO: 0030866), regulation of cell junction assembly (GO: 1901888), negative regulation of circadian rhythm (GO: 0042754), mitotic cytokinesis (GO: 0042755), and regulation of cell junction assembly (GO:1 901888) (GO: 0000281). In addition, a number of GO terms associated with cell growth and division are enriched, including mitotic DNA replication initiation (GO: 1902975), pre-replicative complex assembly involved in nuclear cell cycle DNA replication (GO: 0006267), mitochondrial cell cycle phase transition (GO: 0044772), positive regulation of cell cycle G2/M phase transition (GO: 1902751), and regulation of cyclin-dependent protein serine/threonine kinase activity (GO: 0000079).
A KEGG enrichment analysis of liver DEGs revealed 103 enriched pathways. The top 20 pathways with the most reliable enrichment were: drug metabolism—other enzymes, cell cycle, drug metabolism—cytochrome P450, metabolism of xenobiotics by cytochrome P450, pentose and glucuronate interconversions, retinol metabolism, primary bile acid biosynthesis, porphyrin, chlorophyll metabolism, vascular smooth muscle contraction, riboflavin metabolism, ABC transporters, steroid hormone biosynthesis, p53 signaling pathway, and pyrimidine metabolism.
These pathways and the typical genes involved in the pathways can be classified into substance metabolism (carbohydrate metabolism, lipid metabolism, and metabolism of cofactors and vitamins), signal transduction, and cell growth and movement. DEGs were involved in the following metabolism-related KEGG pathways: ascorbate and aldarate metabolism, pentose and glucuronate interconversions, primary bile acid biosynthesis, steroid hormone biosynthesis, glutathione metabolism, retinol metabolism, and riboflavin metabolism (Figure 6A). These functional categories will guide to reveal differences in liver metabolism in response to Aurantiochytrium sp. in zebrafish. Figure 6A shows the top 20 up- and down-regulated pathways (smallest q-values) and the number of enriched genes; DEGs in each KEGG pathway were annotated according to the pathway type and summarized in Figure 6B. Pathways were enriched for cellular processes, environmental information processing, and metabolism, revealing the potential effects of Aurantiochytrium sp. extract on metabolism, cell growth and motility, and signal transduction in the zebrafish liver.

4. Discussion

Aurantiochytrium sp. produces several biologically active compounds, including n-3 unsaturated fatty acids and carotenoids. These nutrients have the potential to influence fish physiology. The most direct indicators of fish growth and health status are growth performance measures, such as RWG and SGR. Our results reveal that a diet supplemented with 1% TE Aurantiochytrium sp. considerably improved RWG and SGR as well as growth performance in zebrafish over those in the control group. Similar findings have been reported in Nile tilapia (Oreochromis niloticus) [17], catfish (Ictalurus punctatus) [34], and golden pomfret (Trachinotus ovatus) [35].
The fish liver has high enzymatic activity, which is essential for fish growth and development. The activity of digestive enzymes influences nutrition absorption and utilization in aquatic species, which contributes to growth performance [36]. We observed a significant increase in LPS activity in the experimental group, presumably because Aurantiochytrium sp. is rich in n-3 polyunsaturated fatty acids, which promote lipid emulsification in the digestive tract and increase the contact area with digestive enzymes to promote lipase digestion. This was consistent with the observed enrichment for the primary bile acid biosynthetic pathway (ko00120) in the transcriptome. Uncharacterized protein LOC393433 isoform X1, 24-hydroxycholesterol 7-alpha-hydroxylase precursor, 3 beta-hydroxysteroid dehydrogenase type 7 isoform X1, and cytochrome P450 were involved in this regulatory pathway.
Glucose 6-phosphate dehydrogenase is the first enzyme in the pentose phosphate metabolic pathway; its main physiological role is to create NADPH and ribose 5-phosphate, aid in fatty acid synthesis, and maintain glutathione in the reduced form [37]. In our study, G6PD activity increased in the group with 1% Aurantiochytrium sp. supplementation, whereas malic enzyme (NADP-ME) levels were suppressed. ME also catalyzes the generation of NADPH, connecting the glycolytic pathway with the citric acid cycle. However, the pentose phosphate pathway is essential for the maintenance of a normal NADPH/NADP ratio, and the loss of G6PD results in high NADP levels, resulting in enhanced compensatory ME flux [38]. Starch and sucrose metabolism (ko00500), fructose and mannose metabolism (ko00051), galactose metabolism (ko00052), and pentose and glucuronic acid interconversion (ko00040) were all up-regulated to various degrees based on the transcriptome analysis. These results indicate that Aurantiochytrium sp. supplementation increased lipid uptake and facilitated the metabolism of differently structured glycoconjugates in zebrafish and that these pathways and processes may be associated with improved growth performance in zebrafish.
The liver, as the core site of lipid metabolism, plays a vital role in the lipid cycle [25]. FAS, CACT, and ACC all contribute to the control of lipid synthesis [39,40,41]. In our study, FAS, CACT, and ACC activity were down-regulated by Aurantiochytrium sp. ACO, the first rate-limiting enzyme in the fatty acid oxidation pathway, has the maximum activity towards medium-chain fatty acyl-CoA, and its activity diminishes with increasing chain length [42]. ACO levels were elevated in the experimental group, and related lipid cycling pathways were also enriched in the transcriptome, including fatty acid breakdown (ko00071), fatty acid biosynthesis (ko00061), and fatty acid metabolism (ko01212). These findings imply that Aurantiochytrium sp. may influence zebrafish growth by modulating fatty acid production and β-oxidation pathways.
Fish health is related to their antioxidant capacity. MDA is a common lipid peroxidation indicator [43]. The MDA content was dramatically reduced in the 1% Aurantiochytrium sp. group, showing that Aurantiochytrium sp. reduces lipid peroxidation in zebrafish. GSH-Px and SOD both improved to various degrees. GSH-Px functions in the removal of lipid hydroperoxides. When combined with the transcriptome analysis, our results revealed that the glutathione metabolism enzymes gpx4a and mgst1.1 are involved in the control of GSH-Px. These findings are consistent with prior research showing that dietary supplementation with Aurantiochytrium sp. can increase the antioxidant capacity of juvenile mirror carp (Cyprinus carpio var. specularis) [44]. Aurantiochytrium carotenoids are a significant component of the natural antioxidant system. Improving the antioxidant capacity not only protects against the generation of lipid peroxidation, but also reduces the impairment of development induced by oxidative stress [45,46].
RNA-seq is a strong tool for examining the link between genotype and phenotype, providing insight into the underlying pathways and mechanisms that determine cell fate, development, and disease progression [47]. In the current investigation, zebrafish diets were supplemented with a 1% extract of Aurantiochytrium sp. for 56 days for RNA-seq analysis of the liver. Zebrafish DEGs were highly enriched in the GO terms catalytic activity, metabolic process, single-organism process, extracellular area, signaling, and membrane after the consumption of diets containing 1% Aurantiochytrium sp. extract. Aurantiochytrium sp. includes a variety of fatty acids, which regulate lipid metabolism and may play a role in signal transduction in zebrafish [48].
Analysis of KEGG pathways showed that DEGs were involved in environmental formation processing, genetic information processing, human disease, metabolism, and biological systems. In particular, metabolism-related pathways accounted for the majority of the top 20 enriched pathways. The pentose and glucuronide interconversion pathway, retinol metabolism pathway, primary bile acid biosynthesis pathway, porphyrin and chlorophyll metabolism pathway, ascorbate and aldose metabolism pathway, and steroid hormone biosynthesis pathway were significantly up-regulated. The differential genes of glutathione metabolism, riboflavin metabolism, and pyrimidine metabolism pathways were either up-regulated or down-regulated. They were involved in amino acid metabolism, carbohydrate metabolism, sugar biosynthesis, lipid metabolism, cofactor, and vitamin metabolism, covering most metabolic types of organic systems. In addition, some signaling and membrane transport pathways were similarly enriched, for example, p53 signaling pathway (down), ABC transporters (down), and FoxO signaling pathway (up and down). These results further suggested that the regulatory effects of Aurantiochytrium sp. on zebrafish growth were not mediated by a single lipid metabolic pathway.
We selected metabolic genes that may be involved in the control of zebrafish growth for qPCR validation. mgst1.1 could reduce lipid hydroperoxides directly in membranes and plays a crucial role in hematopoietic development [49,50]. gpx4a could react with hydrogen peroxide and a variety of lipid hydroperoxides, including those derived from cholesterol and cholesteryl esters [51]. odc1 is a metabolic enzyme closely related to polyamine biosynthesis associated with the initiation of apoptosis [52]. The control of glutathione metabolism suggested a role for Aurantiochytrium sp. extracts in zebrafish antioxidant enhancement. The liver was essential for maintaining glucose homeostasis in terms of carbohydrate metabolism. Extracts of Aurantiochytrium sp. modulated insulin resistance (ko04931), the insulin signaling pathway (ko04910), starch and sucrose metabolism (ko00500), and other associated pathways. In the insulin resistance pathway, protein O-GlcNAcase (oga), which was involved in glycoprotein metabolism and can suppress the development of type II diabetes is up-regulated [53]. In the insulin signaling pathway, protein phosphatase 1 regulatory subunit 3C-B (ppp1r3cb) and UDP-glucose pyrophosphorylase 2a (ugp2a) and ectonucleotide in starch and sucrose metabolism pyrophosphatase/phosphodiesterase family member 1 isoform X1 (enpp1) were also differently elevated. They contribute to glucose homeostasis and glycogen production or catabolism. Several mechanisms in lipid metabolism were also controlled. In the primary bile acid biosynthesis pathway, the cytochrome P450, family 46, subfamily A, polypeptide 1, and tandem duplicate 2 (cyp46a1.2) genes were increased. In the fatty acid breakdown pathway, acyl-CoA synthetase long-chain family member 1 isoform X1 (acl1b) was down-regulated. Within the glycerophospholipid metabolic pathway, lysophospholipid acyltransferase (lpcat4) was downregulated. acl1b catalyzes the conversion of long-chain fatty acids to acyl-coenzyme A; its downregulation reduces the oxidation of long-chain unsaturated fatty acids, hence increasing retention and deposition in muscle [54]. These data demonstrated altered expression of metabolically significant genes in the liver.
We discovered that dietary Aurantiochytrium sp. affects the FoxO signaling pathway, which regulates gene expression in cellular physiological events, such as apoptosis, cell cycle control, glucose metabolism, and oxidative stress resistance [49]. All these physiological events showed changes in the transcriptome. Aurantiochytrium sp. significantly affected genes downstream of the FoxO signaling pathway, including sgk2b (serine/threonine-protein kinase 2b), gadd45ba (growth arrest and DNA-damage-inducible, beta a), ccnb1 (Cyclin B1), ccnb2 (G2/mitotic-specific cyclin-B2,), plk1 (polo-like kinase), bnip3 (Wu:fj49c01), and cat (catalase). These genes played a role in the cell cycle, cell growth, and antioxidant defense. DHA (one of the main components of Aurantiochytrium sp.) negatively mediates FoxO to affect triacylglycerol metabolism in pigs [50]. However, similar results have rarely been reported in aquatic organisms. More research is needed to determine the relationship between the FoxO signaling pathway and Aurantiochytrium sp.

5. Conclusions

In summary, the liver transcriptome of zebrafish treated with dietary Aurantiochytrium sp. was successfully built and sequenced. Our results provide a detailed overview of the major metabolic pathways and genes that mediate the growth-promoting effects of Aurantiochytrium sp. Aurantiochytrium sp. extract has been shown to boost lipase digestion enzymes’ activity, improve lipid synthesis, and minimize oxidative damage caused by lipid peroxidation. In addition to the regulation of lipid metabolism, Aurantiochytrium sp. also mediates the metabolism and signal transduction of different structural sugar molecules in zebrafish. In particular, through the FoxO signaling pathway, Aurantiochytrium sp. may influence the metabolism of zebrafishIn. However, additional studies are needed to determine the underlying mechanism. The results of this study confirmed the growth-promoting effect of Aurantiochytrium sp. on zebrafish, provided insights into the mechanism by which Aurantiochytrium sp. regulates zebrafish growth, and provided a theoretical basis for the effective application of Aurantiochytrium sp. in aquaculture.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani12202794/s1, Table S1: Aurantiochytrium Extract Main Ingredients; Table S2: The Product Numbers of Commercial Assay kits; Table S3: Primers Used for Detection of DEGs.

Author Contributions

Data curation, H.Y. and Y.H.; formal analysis, H.Y., S.L., and Y.G.; funding acquisition, H.C.; investigation, H.Y., Y.H., and Z.L.; methodology, X.Y., H.Y., and H.H.; project administration, G.L. and H.C.; resources, H.C.; supervision, H.C.; writing—original draft, H.Y.; writing—review and editing, H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This article was supported by South China Aquaculture Breeding Funding (2022RHDKFKT01), the National Natural Science Foundation of China (32273131, 32160861 and 32172971), the Key Research and Development Program of Guangdong (2021B202020002), the talent team tender grant of Zhanjiang marine equipment and biology (2021E05035), the Guangdong Basic and Applied Basic Research Foundation (2019A1515010958), The Major Science and Technology plan of Hainan Province (ZDKJ2021017), Hong Kong Innovation and Technology Commission TCFS (GHP/087/18SZ) and Shenzhen Science and Technology Application Demonstration Project (KJYY20180201180253571).

Institutional Review Board Statement

Animal experiments were carried out in strict ac-cordance with the recommendations in the Guide for the Care and Use of Laboratory Animals. The protocol was approved by the Animal Research and Ethics Committee of Guangdong Ocean University (NIH Pub. No. 85–23, revised 1996).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the agreement with funding bodies.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Naylor, R.L.; Hardy, R.W.; Buschmann, A.H.; Bush, S.R.; Cao, L.; Klinger, D.H.; Little, D.C.; Lubchenco, J.; Shumway, S.E.; Troell, M. A 20-year retrospective review of global aquaculture. Nature 2021, 591, 551–563. [Google Scholar] [CrossRef] [PubMed]
  2. Tocher, D.R. Omega-3 long-chain polyunsaturated fatty acids and aquaculture in perspective. Aquaculture 2015, 449, 94–107. [Google Scholar] [CrossRef]
  3. Steffens, W.; Wirth, M. Freshwater fish-an important source of n-3 polyunsaturated fatty acids: A review. Fish. Aquat. Life 2005, 13, 5–16. [Google Scholar]
  4. Wang, S.; Monroig, O.; Tang, G.; Zhang, L.; You, C.; Tocher, D.R.; Li, Y. Investigating long-chain polyunsaturated fatty acid biosynthesis in teleost fish: Functional characterization of fatty acyl desaturase (Fads2) and Elovl5 elongase in the catadromous species, Japanese eel Anguilla japonica. Aquaculture 2014, 434, 57–65. [Google Scholar] [CrossRef]
  5. Abedi, E.; Sahari, M.A. Long-chain polyunsaturated fatty acid sources and evaluation of their nutritional and functional properties. Food Sci. Nutr. 2014, 2, 443–463. [Google Scholar] [CrossRef] [PubMed]
  6. Watanabe, K.; Arafiles, K.H.V.; Higashi, R.; Okamura, Y.; Tajima, T.; Matsumura, Y.; Nakashimada, Y.; Matsuyama, K.; Aki, T. Isolation of high carotenoid-producing Aurantiochytrium sp. mutants and improvement of astaxanthin productivity using metabolic information. J. Oleo Sci. 2018, 67, 571–578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Iwasaka, H.; Koyanagi, R.; Satoh, R.; Nagano, A.; Watanabe, K.; Hisata, K.; Satoh, N.; Aki, T. A possible trifunctional β-carotene synthase gene identified in the draft genome of Aurantiochytrium sp. strain KH105. Genes 2018, 9, 200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Kaya, K.; Nakazawa, A.; Matsuura, H.; Honda, D.; Inouye, I.; Watanabe, M.M. Thraustochytrid Aurantiochytrium sp. 18W-13a accummulates high amounts of squalene. Biosci. Biotechnol. Biochem. 2011, 75, 2246–2248. [Google Scholar] [CrossRef] [Green Version]
  9. Chen, G.; Fan, K.-W.; Lu, F.-P.; Li, Q.; Aki, T.; Chen, F.; Jiang, Y. Optimization of nitrogen source for enhanced production of squalene from thraustochytrid Aurantiochytrium sp. New Biotechnol. 2010, 27, 382–389. [Google Scholar] [CrossRef] [PubMed]
  10. Yoshida, M.; Ioki, M.; Matsuura, H.; Hashimoto, A.; Kaya, K.; Nakajima, N.; Watanabe, M.M. Diverse steroidogenic pathways in the marine alga Aurantiochytrium. J. Appl. Phycol. 2020, 32, 1631–1642. [Google Scholar] [CrossRef]
  11. Bagul, V.P.; Annapure, U.S. Isolation and characterization of docosahexaenoic acid-producing novel strain Aurantiochytrium sp. ICTFD5: A sterol with vitamin d-cholecalciferol, and cellulase and lipase producing thraustochytrid. Bioresour. Technol. Rep. 2021, 14, 100688. [Google Scholar] [CrossRef]
  12. Du, F.; Wang, Y.-Z.; Xu, Y.-S.; Shi, T.-Q.; Liu, W.-Z.; Sun, X.-M.; Huang, H. Biotechnological production of lipid and terpenoid from thraustochytrids. Biotechnol. Adv. 2021, 48, 107725. [Google Scholar] [CrossRef] [PubMed]
  13. Xie, Y.; Sen, B.; Wang, G. Mining terpenoids production and biosynthetic pathway in thraustochytrids. Bioresour. Technol. 2017, 244, 1269–1280. [Google Scholar] [CrossRef] [PubMed]
  14. Aasen, I.M.; Ertesvåg, H.; Heggeset, T.M.B.; Liu, B.; Brautaset, T.; Vadstein, O.; Ellingsen, T.E. Thraustochytrids as production organisms for docosahexaenoic acid (DHA), squalene, and carotenoids. Appl. Microbiol. Biotechnol. 2016, 100, 4309–4321. [Google Scholar] [CrossRef] [PubMed]
  15. Ji, X.-J.; Mo, K.-Q.; Ren, L.-J.; Li, G.-L.; Huang, J.-Z.; Huang, H. Genome sequence of Schizochytrium sp. CCTCC M209059, an effective producer of docosahexaenoic acid-rich lipids. Genome Announc. 2015, 3, e00819-15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Gupta, A.; Barrow, C.J.; Puri, M. Omega-3 biotechnology: Thraustochytrids as a novel source of omega-3 oils. Biotechnol. Adv. 2012, 30, 1733–1745. [Google Scholar] [CrossRef]
  17. Nobrega, R.O.; Batista, R.O.; Corrêa, C.F.; Mattioni, B.; Filer, K.; Pettigrew, J.E.; Fracalossi, D.M. Dietary supplementation of Aurantiochytrium sp. meal, a docosahexaenoic-acid source, promotes growth of Nile tilapia at a suboptimal low temperature. Aquaculture 2019, 507, 500–509. [Google Scholar] [CrossRef]
  18. Kousoulaki, K.; Berge, G.M.; Turid, M.; Aleksei, K.; Grete, B.; Trine, Y.; Mats, C.; John, S.; Bente, R. Microalgal Schizochytrium limacinum Biomass Improves Growth and Filet Quality When Used Long-Term as a Replacement for Fish Oil, in Modern Salmon Diets. Front. Mar. Sci. 2020, 7, 57. [Google Scholar] [CrossRef]
  19. Sprague, M.; Walton, J.; Campbell, P.J.; Strachan, F.; Dick, J.R.; Bell, J.G. Replacement of fish oil with a DHA-rich algal meal derived from Schizochytrium sp. on the fatty acid and persistent organic pollutant levels in diets and flesh of Atlantic salmon (Salmo salar, L.) post-smolts. Food Chem. 2015, 185, 413–421. [Google Scholar] [CrossRef]
  20. Shi, Y.; Cao, X.; Ye, Z.; Xu, Y.; Wang, Y.; Li, Z.; Hang, W.; He, N. Role of dietary Schizochytrium sp. in improving disease resistance of zebrafish through metabolic and microbial analysis. Aquaculture 2021, 539, 736631. [Google Scholar] [CrossRef]
  21. Fukada, H.; Kitagima, R.; Shinagawa, J.; Morino, H.; Masumoto, T. Effects of complete replacement of fish oil with plant oil mixtures and algal meal on growth performance and fatty acid composition in juvenile yellowtail Seriola quinqueradiata. Fisheries Science 2020, 86, 107–118. [Google Scholar] [CrossRef]
  22. Kumar, V.; Habte-Tsion, H.M.; Allen, K.M.; Bowman, B.A.; Thompson, K.R.; El-Haroun, E.; Filer, K.; Tidwell, J.H. Replacement of fish oil with Schizochytrium meal and its impacts on the growth and lipid metabolism of Pacific white shrimp (Litopenaeus vannamei). Aquac. Nutr. 2018, 24, 1769–1781. [Google Scholar] [CrossRef]
  23. Osmond, A.T.; Arts, M.T.; Hall, J.R.; Rise, M.L.; Bazinet, R.P.; Armenta, R.E.; Colombo, S.M. Schizochytrium sp. (T18) Oil as a Fish Oil Replacement in Diets for Juvenile Rainbow Trout (Oncorhynchus mykiss): Effects on Growth Performance, Tissue Fatty Acid Content, and Lipid-Related Transcript Expression. Animals 2021, 11, 1185. [Google Scholar] [CrossRef] [PubMed]
  24. Li, S.; Wang, B.; Liu, L.; Song, Y.; Lv, C.; Zhu, X.; Luo, Y.; Cheng, C.H.; Chen, H.; Yang, X. Enhanced Growth Performance Physiological and Biochemical Indexes of Trachinotus ovatus Fed With Marine Microalgae Aurantiochytrium sp. Rich in n-3 Polyunsaturated Fatty Acids. Front. Mar. Sci. 2021, 7, 609837. [Google Scholar] [CrossRef]
  25. Nguyen, P.; Leray, V.; Diez, M.; Serisier, S.; Bloc’h, J.L.; Siliart, B.; Dumon, H. Liver lipid metabolism. J. Anim. Physiol. Anim. Nutr. 2008, 92, 272–283. [Google Scholar] [CrossRef] [PubMed]
  26. Lu, X.; Chen, H.-M.; Qian, X.-Q.; Gui, J.-F. Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast-and slow-growing fish. Comp. Biochem. Physiol. Part D Genom. Proteom. 2020, 35, 100688. [Google Scholar] [CrossRef]
  27. Caballero-Solares, A.; Xue, X.; Parrish, C.C.; Foroutani, M.B.; Taylor, R.G.; Rise, M.L. Changes in the liver transcriptome of farmed Atlantic salmon (Salmo salar) fed experimental diets based on terrestrial alternatives to fish meal and fish oil. BMC Genom. 2018, 19, 796. [Google Scholar] [CrossRef] [Green Version]
  28. Liu, L.; Hu, Z.; Li, S.; Yang, H.; Li, S.; Lv, C.; Zaynab, M.; Cheng, C.H.K.; Chen, H.; Yang, X. Comparative Transcriptomic Analysis Uncovers Genes Responsible for the DHA Enhancement in the Mutant Aurantiochytrium sp. Microorganisms 2020, 8, 529. [Google Scholar] [CrossRef] [Green Version]
  29. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [Green Version]
  30. Young, M.D.; Wakefield, M.J.; Smyth, G.K.; Oshlack, A. Gene ontology analysis for RNA-seq: Accounting for selection bias. Genome Biol. 2010, 11, R14. [Google Scholar] [CrossRef] [Green Version]
  31. Mao, X.; Cai, T.; Olyarchuk, J.G.; Wei, L. Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics 2005, 21, 3787–3793. [Google Scholar] [CrossRef] [PubMed]
  32. Tang, R.; Dodd, A.; Lai, D.; McNabb, W.C.; Love, D.R. Validation of zebrafish (Danio rerio) reference genes for quantitative real-time RT-PCR normalization. Acta Biochim. Biophys. Sin. 2007, 39, 384–390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2 −ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  34. Li, M.H.; Robinson, E.H.; Tucker, C.S.; Manning, B.B.; Khoo, L. Effects of dried algae Schizochytrium sp., a rich source of docosahexaenoic acid, on growth, fatty acid composition, and sensory quality of channel catfish Ictalurus punctatus. Aquaculture 2009, 292, 232–236. [Google Scholar] [CrossRef]
  35. Xie, J.; Fang, H.; Liao, S.; Guo, T.; Yin, P.; Liu, Y.; Tian, L.; Niu, J. Study on Schizochytrium sp. improving the growth performance and non-specific immunity of golden pompano (Trachinotus ovatus) while not affecting the antioxidant capacity. Fish Shellfish Immunol. 2019, 95, 617–623. [Google Scholar] [CrossRef]
  36. Thongprajukaew, K.; Kovitvadhi, U.; Kovitvadhi, S.; Engkagul, A.; Rungruangsak-Torrissen, K. Evaluation of growth performance and nutritional quality of diets using digestive enzyme markers and in vitro digestibility in Siamese fighting fish (Betta splendens Regan, 1910). Afr. J. Biotechnol. 2013, 12, 6400–6410. [Google Scholar] [CrossRef] [Green Version]
  37. Ciftci, M.; Turkoglu, V.; Coban, T.A. Effects of some drugs on hepatic glucose 6-phosphate dehydrogenase activity in Lake Van Fish (Chalcalburnus Tarischii Pallas, 1811). J. Hazard. Mater. 2007, 143, 415–418. [Google Scholar] [CrossRef]
  38. Chen, L.; Zhang, Z.; Hoshino, A.; Zheng, H.D.; Morley, M.; Arany, Z.; Rabinowitz, J.D. NADPH production by the oxidative pentose-phosphate pathway supports folate metabolism. Nat. Metab. 2019, 1, 404–415. [Google Scholar] [CrossRef]
  39. Xiong, Q.; Feng, D.; Wang, Z.; Ying, Y.; Xu, C.; Wei, Q.; Zeng, S.; Yang, L. Fatty Acid Synthase Is the Key Regulator of Fatty Acid Metabolism and Is Related to Immunotherapy in Bladder Cancer. Front. Immunol. 2022, 13, 836939. [Google Scholar] [CrossRef]
  40. Lu, K.-L.; Xu, W.-N.; Wang, L.-N.; Zhang, D.-D.; Zhang, C.-N.; Liu, W.-B. Hepatic β-oxidation and regulation of carnitine palmitoyltransferase (CPT) I in blunt snout bream Megalobrama amblycephala fed a high fat diet. PLoS ONE 2014, 9, e93135. [Google Scholar] [CrossRef] [Green Version]
  41. Pang, Y.; Xu, X.; Xiang, X.; Li, Y.; Zhao, Z.; Li, J.; Gao, S.; Liu, Q.; Mai, K.; Ai, Q. High Fat Activates O-GlcNAcylation and Affects AMPK/ACC Pathway to Regulate Lipid Metabolism. Nutrients 2021, 13, 1740. [Google Scholar] [CrossRef] [PubMed]
  42. Morais, S.; Knoll-Gellida, A.; André, M.; Barthe, C.; Babin, P.J. Conserved expression of alternative splicing variants of peroxisomal acyl-CoA oxidase 1 in vertebrates and developmental and nutritional regulation in fish. Physiol. Genom. 2007, 28, 239–252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Moore, K.; Roberts, L.J. Measurement of lipid peroxidation. Free Radic. Res. 1998, 28, 659–671. [Google Scholar] [CrossRef] [PubMed]
  44. Xiao, F.; Xing, J.; Li, H.; Xu, X.; Hu, Z.; Ji, H. Effects of the defatted Schizochytrium sp. on growth performance, fatty acid composition, histomorphology and antioxidant status of juvenile mirror carp (Cyprinus carpio var. specularis). Aquac. Res. 2021, 52, 3062–3076. [Google Scholar] [CrossRef]
  45. Tsuduki, T.; Honma, T.; Nakagawa, K.; Ikeda, I.; Miyazawa, T. Long-term intake of fish oil increases oxidative stress and decreases lifespan in senescence-accelerated mice. Nutrition 2011, 27, 334–337. [Google Scholar] [CrossRef] [PubMed]
  46. Yang, Q.; Lu, T.; Yan, J.; Li, J.; Zhou, H.; Pan, X.; Lu, Y.; He, N.; Ling, X. Regulation of polyunsaturated fatty acids synthesis by enhancing carotenoid-mediated endogenous antioxidant capacity in Schizochytrium sp. Algal Res. 2021, 55, 102238. [Google Scholar] [CrossRef]
  47. Qian, X.; Ba, Y.; Zhuang, Q.; Zhong, G. RNA-Seq technology and its application in fish transcriptomics. Omics J. Integr. Biol. 2014, 18, 98–110. [Google Scholar] [CrossRef] [Green Version]
  48. Calviello, G.; Serini, S.; Palozza, P. n-3 polyunsaturated fatty acids as signal transduction modulators and therapeutical agents in cancer. Curr. Signal Transduct. Ther. 2006, 1, 255–271. [Google Scholar] [CrossRef]
  49. Glauser, D.A.; Schlegel, W. The emerging role of FOXO transcription factors in pancreatic β cells. J Endocrinol. 2007, 193, 195–207. [Google Scholar] [CrossRef] [Green Version]
  50. Chen, Y.-J.; Chen, C.-C.; Li, T.-K.; Wang, P.-H.; Liu, L.-R.; Chang, F.-Y.; Wang, Y.-C.; Yu, Y.-H.; Lin, S.-P.; Mersmann, H.J.; et al. Docosahexaenoic acid suppresses the expression of FoxO and its target genes. J Nutr Biochem. 2012, 23, 1609–1616. [Google Scholar] [CrossRef]
  51. Maiorino, M.; Thomas, J.P.; Girotti, A.W.; Ursini, F. Reactivity of phospholipid hydroperoxide glutathione peroxidase with membrane and lipoprotein lipid hydroperoxides. Free Radic Res Commun. 1991, 12, 131–135. [Google Scholar] [CrossRef] [PubMed]
  52. Pignatti, C.; Tantini, B.; Stefanelli, C.; Flamigni, F. Signal transduction pathways linking polyamines to apoptosis. Amino Acids 2004, 27, 359–365. [Google Scholar] [CrossRef] [PubMed]
  53. Webster, D.M.; Teo, C.F.; Sun, Y.; Wloga, D.; Gay, S.; Klonowski, K.D.; Wells, L.; Dougan, S.T. O-GlcNAc modifications regulate cell survival and epiboly during zebrafish development. BMC Dev Biol. 2009, 9, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Cheng, H.-L.; Chen, S.; Xu, J.-H.; Yi, L.-F.; Peng, Y.-X.; Pan, Q.; Shen, X.; Dong, Z.-G.; Zhang, X.-Q.; Wang, W.-X. Molecular cloning and nutrient regulation analysis of long chain acyl-CoA synthetase 1 gene in grass carp, Ctenopharyngodon idella L. Comp Biochem Physiol B Biochem Mol Biol. 2017, 204, 61–68. [Google Scholar] [CrossRef] [PubMed]
Animals 12 02794 g001 550
Figure 1. Effect of different diets on body weight and length of zebrafish. Growth data of dietary EPA group (A,B), dietary DHA group (C,D), and dietary TE group (E,F). Error bars indicate SEM (n = 9). EPA: eicosapentaenoic acid; DHA: docosahexaenoic acid; TE: Aurantiochytrium sp. extract.
Figure 1. Effect of different diets on body weight and length of zebrafish. Growth data of dietary EPA group (A,B), dietary DHA group (C,D), and dietary TE group (E,F). Error bars indicate SEM (n = 9). EPA: eicosapentaenoic acid; DHA: docosahexaenoic acid; TE: Aurantiochytrium sp. extract.
Animals 12 02794 g001
Animals 12 02794 g002 550
Figure 2. Comparative assessment of liver enzyme activities related to digestion (A), energy metabolism (B), lipid metabolism (C), and antioxidant defense (D) in the control group (C) and 1% TE group (T). Error bars indicate SEM (n = 9). Asterisks indicate significant differences between the two groups (* p < 0.05, ** p < 0.01), ns: not significant (p > 0.05) α-AMS: hepatic amylase; LPS: hepatic lipase; G-6-PD: dehydrogenase; NADP-ME: malic enzyme; PK: pyruvate kinase; FAS: fatty acid synthase; ACO: acetyl-CoA oxidase; CACT: carnitine lipid acyltransferase; ACC: acetyl-CoA; MDA: malondialdehyde; CAT: catalase; GSH–PX: glutathione peroxidase; SOD: superoxide dismutase.
Figure 2. Comparative assessment of liver enzyme activities related to digestion (A), energy metabolism (B), lipid metabolism (C), and antioxidant defense (D) in the control group (C) and 1% TE group (T). Error bars indicate SEM (n = 9). Asterisks indicate significant differences between the two groups (* p < 0.05, ** p < 0.01), ns: not significant (p > 0.05) α-AMS: hepatic amylase; LPS: hepatic lipase; G-6-PD: dehydrogenase; NADP-ME: malic enzyme; PK: pyruvate kinase; FAS: fatty acid synthase; ACO: acetyl-CoA oxidase; CACT: carnitine lipid acyltransferase; ACC: acetyl-CoA; MDA: malondialdehyde; CAT: catalase; GSH–PX: glutathione peroxidase; SOD: superoxide dismutase.
Animals 12 02794 g002
Animals 12 02794 g003 550
Figure 3. Volcano plot of differential gene expression in the zebrafish liver. Green dots represent down-regulated genes in the 1% Aurantiochytrium sp. extract group, red dots represent up-regulated genes, and black dots represent non-differentially expressed genes. Significance criteria were q < 0.05 and fold change ≥ 2.
Figure 3. Volcano plot of differential gene expression in the zebrafish liver. Green dots represent down-regulated genes in the 1% Aurantiochytrium sp. extract group, red dots represent up-regulated genes, and black dots represent non-differentially expressed genes. Significance criteria were q < 0.05 and fold change ≥ 2.
Animals 12 02794 g003
Animals 12 02794 g004a 550Animals 12 02794 g004b 550
Figure 4. Validation of DEGs and clustering heatmap. (A) Validation of differential expression by RT-qPCR. Error bars indicate SEM (n = 3). (B) Cluster map of gene expression patterns demonstrating general patterns in distinct liver samples. (C) Correlation analysis of log2-fold change values between RNA-Seq data (x-axis) and RT-qPCR data (y-axis).
Figure 4. Validation of DEGs and clustering heatmap. (A) Validation of differential expression by RT-qPCR. Error bars indicate SEM (n = 3). (B) Cluster map of gene expression patterns demonstrating general patterns in distinct liver samples. (C) Correlation analysis of log2-fold change values between RNA-Seq data (x-axis) and RT-qPCR data (y-axis).
Animals 12 02794 g004aAnimals 12 02794 g004b
Animals 12 02794 g005a 550Animals 12 02794 g005b 550
Figure 5. GO functional enrichment analysis of differentially expressed genes in the liver. Histograms of GO term enrichment in the (A) cellular component category, (B) molecular function category, and (C) biological process category.
Figure 5. GO functional enrichment analysis of differentially expressed genes in the liver. Histograms of GO term enrichment in the (A) cellular component category, (B) molecular function category, and (C) biological process category.
Animals 12 02794 g005aAnimals 12 02794 g005b
Animals 12 02794 g006 550
Figure 6. KEGG enrichment and classification. (A) Scatter plot of KEGG pathway enrichment results for differentially expressed genes: the ordinate represents the pathway name, the abscissa represents the enrichment factor, and the color represents the q-value (i.e., the p-value after correction for multiple hypothesis testing). The size of the circle indicates the number of genes involved in the pathway. (B) KEGG classification diagram: the left-hand vertical axis shows the KEGG pathway name, and the right-hand vertical axis shows the category corresponding to each pathway. A shared column color indicates the same category.
Figure 6. KEGG enrichment and classification. (A) Scatter plot of KEGG pathway enrichment results for differentially expressed genes: the ordinate represents the pathway name, the abscissa represents the enrichment factor, and the color represents the q-value (i.e., the p-value after correction for multiple hypothesis testing). The size of the circle indicates the number of genes involved in the pathway. (B) KEGG classification diagram: the left-hand vertical axis shows the KEGG pathway name, and the right-hand vertical axis shows the category corresponding to each pathway. A shared column color indicates the same category.
Animals 12 02794 g006
Table
Table 1. Effects of various diets on zebrafish development.
Table 1. Effects of various diets on zebrafish development.
Growth ParametersControl1% EPA2% EPA3% EPA
RWG (%/day)60.52 ± 4.6359.02 ± 3.2859.10 ± 3.3537.81 ± 2.32 **
SGR (%/day)6.33 ± 0.146.29 ± 0.106.30 ± 0.105.53 ± 0.11 **
CF (g/cm3)2.07 ± 0.102.03 ± 0.142.14 ± 0.082.11 ± 0.13
SR (%)93.33 ± 0.2790.00 ± 0.4790.00 ± 0.8291.67 ± 0.27
Growth ParametersControl1% DHA2% DHA3% DHA
RWG(%/day)60.52 ± 4.6347.44 ± 2.99 *61.52 ± 2.1176.26 ± 3.12 *
SGR (%/day)6.33 ± 0.145.92 ± 0.11 *6.37 ± 0.066.75 ± 0.07 *
CF (g/cm3)2.07 ± 0.102.02 ± 0.051.90 ± 0.142.19 ± 0.23
SR (%)93.33 ± 0.2788.33 ± 0.2790.00 ± 0.4793.33 ± 0.54
Growth ParametersControl1% TE2% TE3% TE
RWG (%/day)60.52 ± 4.6380.73 ± 2.33 **80.80 ± 3.74 **67.35 ± 4.80
SGR (%/day)6.33 ± 0.146.84 ± 0.05 *6.84 ± 0.08 *6.52 ± 0.12
CF (g/cm3)2.07 ± 0.102.23 ± 0.032.19 ± 0.092.00 ± 0.04
SR (%)93.33 ± 0.2791.67 ± 0.7295.00 ± 0.4793.33 ± 0.27
RWG: relative growth rate for body weight, SGR: specific growth rate, CF: condition factor, and SR: survival rate. (* p < 0.05, ** p < 0.01).
Table
Table 2. Sequencing data statistics.
Table 2. Sequencing data statistics.
SamplesClean ReadsClean BasesGC Content≥Q30%
C-120,237,3396,043,571,20648.47%93.81%
C-222,833,7296,803,261,10248.57%94.47%
C-329,402,4928,774,618,35048.43%94.45%
T-122,681,1846,757,938,49248.19%94.47%
T-222,571,8336,730,129,40648.49%94.35%
T-321,741,7086,489,332,01648.26%94.21%
Table
Table 3. Statistical summary of sequencing data and mapping to the zebrafish reference genome.
Table 3. Statistical summary of sequencing data and mapping to the zebrafish reference genome.
SampleTotal ReadsMapped ReadsUniq Mapped ReadsMultiple Map ReadsReads Map to ‘+’Reads Map to ‘−’
C-140,474,67836,863,629 (91.08%)34,099,049 (84.25%)2,764,580 (6.83%)18,396,027 (45.45%)18,428,782 (45.53%)
C-245,667,45842,081,877 (92.15%)39,375,737 (86.22%)2,706,140 (5.93%)20,994,119 (45.97%)21,028,413 (46.05%)
C-358,804,98453,960,399 (91.76%)50,048,173 (85.11%)3,912,226 (6.65%)26,914,470 (45.77%)26,978,224 (45.88%)
T-145,362,36841,432,683 (91.34%)37,606,267 (82.90%)3,826,416 (8.44%)20,648,798 (45.52%)20,715,844 (45.67%)
T-245,143,66640,850,423 (90.49%)37,875,132 (83.90%)2,975,291 (6.59%)20,380,508 (45.15%)20,420,021 (45.23%)
T-343,483,41639,832,673 (91.60%)36,965,709 (85.01%)2,866,964 (6.59%)19,868,815 (45.69%)19,905,921 (45.78%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, H.; Huang, Y.; Li, Z.; Guo, Y.; Li, S.; Huang, H.; Yang, X.; Li, G.; Chen, H. Effects of Dietary Supplementation with Aurantiochytrium sp. on Zebrafish Growth as Determined by Transcriptomics. Animals 2022, 12, 2794. https://doi.org/10.3390/ani12202794

AMA Style

Yang H, Huang Y, Li Z, Guo Y, Li S, Huang H, Yang X, Li G, Chen H. Effects of Dietary Supplementation with Aurantiochytrium sp. on Zebrafish Growth as Determined by Transcriptomics. Animals. 2022; 12(20):2794. https://doi.org/10.3390/ani12202794

Chicago/Turabian Style

Yang, Hao, Yanlin Huang, Zhiyuan Li, Yuwen Guo, Shuangfei Li, Hai Huang, Xuewei Yang, Guangli Li, and Huapu Chen. 2022. "Effects of Dietary Supplementation with Aurantiochytrium sp. on Zebrafish Growth as Determined by Transcriptomics" Animals 12, no. 20: 2794. https://doi.org/10.3390/ani12202794

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

Article Metrics

Back to TopTop