S. Adam Fuller1*†
Jason W. Abernathy2*†Nithin Muliya Sankappa2,3†Benjamin H. Beck2Steven D. Rawles1Bartholomew W. Green1
Kurt A. Rosentrater4Matthew E. McEntire1George Huskey1Carl D. Webster1*- 1USDA-ARS Harry K. Dupree Stuttgart National Aquaculture Research Center (HKDSNARC), Stuttgart, AR, United States
- 2USDA-ARS Aquatic Animal Health Research Unit (AAHRU), Auburn, AL, United States
- 3Oak Ridge Institute for Science and Education (ORISE), ARS Research Participation Program, Oak Ridge, TN, United States
- 4Iowa State University, Agricultural and Biosystems Engineering, Ames, IA, United States
White bass (Morone chrysops) are a popular sportfish throughout the southern United States, and one parent of the commercially-successful hybrid striped bass (M. chrysops ♂ x M. saxatilis ♀). Currently, white bass are cultured using diets formulated for other carnivorous fish, such as largemouth bass (Micropterus salmoides) or hybrid striped bass and contain a significant percentage of marine fish meal. Since there are no studies regarding the utilization of alternative proteins in this species, we evaluated the global gene expression of white bass fed diets in which fish meal was partially or totally replaced by various combinations of soybean meal, poultry by-product meal, canola meal, soy protein concentrate, wheat gluten, or a commercial protein blend (Pro-Cision™). Six isonitrogenous (40% protein), isolipidic (11%), and isocaloric (17.1 kJ/g) diets were formulated to meet the known nutrient and energy requirements of largemouth bass and hybrid striped bass using nutrient availability data for most of the dietary ingredients. One of the test diets consisted exclusively of plant protein sources. Juvenile white bass (40.2 g initial weight) were stocked into a flow-through aquaculture system (three tanks/diet; 10 fish/tank) and fed the test diets twice daily to satiation for 60 days. RNA sequencing and bioinformatic analyses revealed significant differentially expressed genes between all test diets when compared to fish meal control. A total of 1,260 differentially expressed genes were identified, with major ontology relating to cell cycle and metabolic processes as well as immune gene functions. This data will be useful as a resource for future refinements to moronid diet formulation, as marine fish meal becomes limiting and plant ingredients are increasingly added as a reliable protein source.
Introduction
Aquaculture has traditionally used marine fish meal (FM) for the majority of protein in prepared diets due to its excellent nutrient composition, especially for carnivorous fish. However, it is the most expensive macro-ingredient in the diet and is of finite supply. Thus, within the past 4 decades, much research has been conducted to partially or totally replace FM with alternative, more renewable protein ingredients. Both animal-source protein ingredients (such as poultry by-product meal) and plant-source protein ingredients (such as soybean meal and canola meal) have been added to aquaculture diets, partially or totally replacing FM depending upon the fish species. However, with this dietary shift away from FM, there has been evidence that variations in gene expression have been reported. Panserat et al. (2008) reported that a decrease in protein biosynthesis capacity and nitrogen metabolism occurred in rainbow trout (Oncorhynchus mykiss) in liver genes, while Panserat et al. (2009) reported that 2% of the hepatic transcriptome were differentially expressed between fish fed a diet containing FM and fish fed a plant-based diet. However, no major dysfunction of hepatic metabolism and no over-expression of hepatic genes involved in stress and health were reported.
Webster et al. (1999) stated that blends of animal and plant protein ingredients in aquafeeds often result in better growth performance in fish because they provide complimentary nutrient profiles, especially vitamins and essential amino acids (EAAs), that meet the animal’s nutrient requirements. Further, they may reduce potential negative palatability (diet acceptability) associated with using predominantly one protein source (Webster et al., 1992). Among plant-protein ingredients, soybean meal (SBM) is the most widely used in aquaculture diets because of its wide availability, good nutrient composition, and high digestibility (Akiyama et al., 1989; Gatlin et al., 2007; Lemos et al., 2009). However, SBM also has negative attributes that may limit its use at high percentages, or as the predominant protein in commercial aquaculture diets for some fish species. For example, SBM has lower overall EAA concentrations than FM (deficiencies in methionine—MET, lysine—LYS, and threonine -THR), and lack the essential n-3 fatty acids eicosapentaenoic acid (EPA; 20:5n-3), and docosahexaenoic acid (DHA: 22:6n-3; Fox et al., 2004; Gatlin et al., 2007; NRC, 2011). Soybean meal also contains anti-nutritional factors (ANFs), such as trypsin inhibitors, lectins, phytic acid, saponins, antivitamins, oligosaccharides, and high levels of non-starch polysaccharides that may affect nutrient digestibility and/or availability to fish (Francis et al., 2001; Gilani et al., 2005; Gatlin et al., 2007).
The second-most widely used plant protein ingredient is rapeseed (canola) meal (CM). Canola is the name for varieties of rapeseed (Brassica napus and B. campestris) that have low levels of erucic acid and glucosinolates. Canola meal possesses favorable EAA composition, bioactive peptides (Akbari and Wu, 2015), and phenolic compounds that have antioxidant activity (Mazumder et al., 2016). However, use of CM in aquaculture diets is limited by the presence of ANFs and high fiber levels (Higgs et al., 1982; Lim et al., 1998; Webster et al., 2000; Glencross et al., 2004; Ngo et al., 2016; Dossou et al., 2018).
Animal-protein ingredients, such as poultry by-product meal (PBM), have variable success as partial or total alternatives to FM in diets for various fish. For several carnivorous fish, PBM is an excellent replacement for moderate-to-high percentages of FM (Nengas et al., 1999; Tidwell et al., 2005; Wang et al., 2015; Wu et al., 2018). However, animal by-product meals can vary greatly in quality depending upon what is used to make the final product. Further, there are commercial protein blends available for use as FM substitutes in aquaculture diets, such as Aqua-Pak Elite™ series and Pro-Cision™ (H.J. Baker Bros., Inc., Tuscola, TX, United States) blends that are comprised of proprietary combinations of animal proteins that include PBM, meat and bone meal, feather meal, blood meal, and plant proteins that include corn gluten, wheat gluten (WG), distiller’s dried grains with solubles, soy protein concentrate (SPC), and soy protein isolate (SPI) (McCann et al., 2021).
Nutrigenomics, the integration of nutrition and genomic analysis, allows for the understanding of the effects, and affects, of diet on gene expression. This allows for a more thorough and complete assessment of ingredients used in diet formulations for fish. Replacement of FM by other animal and/or plant protein sources is a vital area of research for aquaculture as use of plant-protein ingredients has been shown to affect growth, body composition, and biochemical/histological composition of fish at the molecular level (Gómez-Requieni et al., 2004). Growth rates in fish have been reported to be influenced by growth hormone (GH) and insulin-like growth factor (IGF) axis (Company et al., 2001) and protein sources may affect the expression of GH and IGF-1 encoding genes (Kumar et al., 2011).
White bass (WB; Morone chrysops) are an under-utilized fish in the U.S., but has the potential to become a valuable food fish industry. Because it is carnivorous, WB typically are fed high-protein (>35% crude protein) diets formulated with high FM content similar to other commercially important carnivorous fish. Nutrient requirements or diet formulations for WB may be similar to hybrid striped bass (HSB; M. chrysops ♂ x M. saxatilis ♀ and M. saxatilis ♂ x M. chrysops ♀) (NRC, 2011) and largemouth bass (LMB; Micropterus salmoides) (Tidwell et al., 1996; Coyle et al., 2000; Portz et al., 2001; Tidwell et al., 2005), but to date, there are no published reports of nutrient requirements or diet formulations for WB and thus, no data on gene expression when WB are fed diets in which FM has been replaced. The objectives of this feeding trial were to evaluate global gene expression of white bass fed diets in which FM was partially or totally replaced by a combination of SBM and PBM, and other plant proteins that included CM, SPC, WG, or the commercial protein blend Pro-Cision™.
Materials and methods
Fish and culture systems
A total of 180 white bass (40.1 ± 1.83 g; mean ± SE) were randomly stocked into 200-L tanks (n = 18) at a density of 10 fish/tank (three replicate tanks per dietary treatment). The experimental protocol was approved by the HKD-SNARC Institutional Animal Care and Use Committee and conformed to USDA Agricultural Research Service Policies and Procedures 130.4 and 635.1. Water from a freshwater well was supplied at 2.0 L/min in a flow-through design. Fish were carefully hand-fed their respective pelleted test diets twice daily (0700 and 1530 h) all they would consume in 20 min for the 60-day feeding trial. No more diet was fed when the last pellet was consumed. Mortalities were replaced during the first week of the feeding trial. To mimic natural environment, light was provided by overhead fluorescent ceiling lights set on a photoperiod of 15 h:9 h light:dark cycle. All tanks were siphoned daily to remove any fecal waste. If present, mortalities were recorded and removed daily. All tanks were covered with polyethylene mesh to prevent fish loss due to escapement or jumping to adjacent tanks.
Dissolved oxygen (DO) and water temperature (°C) were measured daily (Model 556MPS, YSI, Inc., Yellow Springs, OH, United States), while alkalinity, hardness, chlorides, total ammonia nitrogen (TAN), nitrite nitrogen, nitrate, and pH were recorded weekly. TAN and nitrite nitrogen were analyzed by multiparameter colorimeter (Model DR 900, Hach, Co., Loveland, CO, United States) using the salicylate method (TAN) and azo-dye method (nitrite nitrogen). Average water quality parameters (±SE) during the feeding trial were 22.5°C ± 0.1°C water temperature, 7.99 ± 0.03 mg DO/L, 0.03 ± 0.01 mg TAN/L, 0.01 ± 0.00 mg nitrite/L, 191.5 ± 1.9 mg total alkalinity/L, 110.8 ± 1.4 mg hardness/L, 202.1 ± 3.1 mg chloride/L, and a pH of 8.14 ± 0.12. All parameters were within acceptable limits for fish growth and health (Boyd, 1979).
Diets, feeding, and sampling
Six isonitrogenous (40% protein), isolipidic (11%), and isocaloric (17.1 kJ/g) diets were formulated (Table 1) on an as-fed basis to meet the known nutrient and energy requirements of LMB (Portz and Cyrino, 2003; Portz and Cyrino, 2004). Diets were formulated on a digestible protein basis or an available amino acid basis when available for FM, SBM, PBM, and wheat as determined in LMB (Portz and Cyrino, 2004), WG (Joana et al., 1998) and soy protein concentrate (Kaushik, 2002) as determined in European seabass, and CM as determined in HSB (Gaylord et al., 2004). Amino acid availabilities for WG, SPC, and CM in white bass are not published and were assumed to be 90%. Nutrient requirements for LMB were used as it has been previously demonstrated that LMB can be fed a diet formulated to meet nutrient requirements of HSB (Tidwell et al., 2005) and LMB require some EAA at higher levels than HSB (NRC, 2011): THR (Rahman et al., 2021), LYS (Dairiki et al., 2007), and ARG (Zhou et al., 2012). EAA composition of each formulated diet is provided in Table 2. Lipid contribution from dry ingredients was balanced with canola oil and menhaden fish oil to maintain the diets isolipidic and isoenergetic, while the contribution of fish oil to total lipid was maintained constant in all diets by including 1.5%–2% menhaden fish oil in the fish meal free diets. Diet 1 (30% FM-control) was formulated as a high-quality, commercial carnivorous fish diet containing 30% FM and 37% SBM. Diets 2–5 were formulated to replace FM partially or totally with other commercial protein ingredients. Diet 6 (Pro-Cision™ + SBM) was formulated to replace FM completely with a commercial protein blend (AquaPak Pro-Cision™, H.J. Bakers Bros., Inc., Tuscola, TX, United States) at a 40% inclusion level to have similar formulated nutrient composition as the control (Diet 1).
TABLE 1. Composition (g/kg; as-fed basis) of six practical test diets fed to juvenile white bass. Values are means of three determinations per diet.
TABLE 2. Essential amino acid (EAA) composition (g/kg; as-fed basis) of six test diets fed to juvenile white bass. Values are means of two determinations per diet.
All ingredients were ground in a pilot-scale hammermill (Viking, CPM Roskamp Champion, Waterloo, IA, United States), and then blended in a rotating mixer (Model 043206, Kobalt/Monarch Industries, Winnipeg, Canada). Extrusion was performed using a single-screw extruder (Model 500, Insta-Pro, Inc., Des Moines, Iowa, United States) with a 45 mm diameter screw and a 20:1 screw length to diameter (L/D) ratio. Feed blends were manually fed into the extruder at a constant rate. The extruder was connected to a 7.5 HP motor and screw speed was set at 600 RPM. A circular die plate (with multiple 3 mm holes) was attached to the extruder. The mass flow rate was determined by collecting extrudate samples at 30 s intervals during extrusion processing and weighing the samples on an electronic balance. Observed mass flow rates ranged from 0.089 to 0.095 kg/s. The temperatures of the die and of the resulting pellets were recorded every 2 minutes using an infrared thermometer and determined to be 53°C ± 5°C for the die plate and 63°C–70°C for the extrudates. After extrusion processing, the extrudates were air dried for 24 h, then bagged and frozen. Diets were stored at −20°C in plastic containers until fed.
Data collection and sample analysis
Fish were not fed for 18-h prior to the conclusion of the study. Subsequently, fish in each tank were enumerated to determine survival (%), individually weighed to the nearest 0.01 g on an electronic scale (Model AB54-S; Mettler Toledo, Columbus, OH, United States), and measured for total length to the nearest 0.1 mm. Three fish from each of the three replicate tanks each diet (n = 9) were arbitrarily selected, euthanized according to American Veterinary Medical Association guidelines (Underwood and Anthony, 2020), and liver samples dissected and placed in RNAlater Stabilization Solution (ThermoFisher Scientific, Waltham, MA) and stored at −80°C for transcriptomic analysis.
RNA isolation from tissues and creation of sequencing libraries
Total RNA was extracted with the Qiagen RNeasy Mini Kit (Cat. No. 74106; Germantown, MD) per the manufacturer’s instruction. Total RNA was treated with Amplification Grade DNase I (Cat. No. AMPD1; Sigma-Aldrich, St. Louis, MO) per the package insert. Afterword, each sample was ethanol-precipitated and resuspended in nuclease-free water. Sample concentrations were then measured via spectrophotometry (Bio-Tek Synergy H1, Winooski, VT) and each sample standardized to 100 ng/μL by dilution with nuclease-free water. RNA quality was assessed using the Agilent BioAnalyzer with the RNA 6000 Nano Kit (Cat. No. 5067-1511; Santa Clara, CA), where RNA Integrity Numbers (RIN) > 9 were considered high-quality for processing the RNA further to be sequenced. Sequencing libraries were then prepared using the NEBNext Ultra II Directional RNA Library Prep Kit for Illumina (Cat. No. E7760L) with the NEBNext Poly(A) mRNA Magnetic Isolation Module (Cat. No. E7490L; New England Biolabs, Ipswich, MA). Barcodes used were the NEBNext Multiplex Oligos for Illumina (Index Primers Sets 1 and 2; Cat. Nos. E7335L and E7500L, respectively). Library insert size and quality were assessed on the Agilent BioAnalyzer using the High Sensitivity DNA Kit (Cat. No. 5067-4626). The NEBNext Library Quant Kit for Illumina (Cat. No. E7630L) was used to quantify each sample. Equimolar amounts (5 nM) of indexed samples were combined, and these pooled libraries were sent for 150 bp paired-end sequencing on an Illumina HiSeq X Ten (San Diego, CA) via a service provider (Novogene, Sacramento, CA).
Curation of the protein-coding white bass transcriptome
A total of 185,531 published white bass transcriptome contigs (Li et al., 2014) were re-analyzed and a curated set of protein-coding transcripts selected to detect differentially expressed genes (DEGs), herein referred to as white bass transcriptome. To create the WB transcriptome, all transcripts were analyzed by BLASTx to the non-redundant (nr) database and only those with blast score E-value < e-20 were retained. Then, several measures were taken to reduce redundancy. First, contigs with >90% similarity were clustered using CD-HIT-EST (Fu et al., 2012). Sequence similarity was further detected using BLAT (Kent, 2002), where duplicates at a >80% Similarity threshold were removed. Then, custom shell scripts were used to manually curate these clusters, where those with top-hit BLASTx returns of hypothetical/unknown protein identities were removed, and those with duplicate top-hit descriptions were reduced to include only the longest sequence. This produced a protein-coding WB transcriptome containing 20,914 sequences. BUSCO analysis (Simão et al., 2015) against the Actinopterygii_odb9 lineage dataset indicated a substantially completed transcriptome, with only 8.4% reported as missing.
Identification of differentially expressed genes
Raw RNAseq reads were processed for quality control using the TrimGalore! Software (Krueger et al., 2021). This software was used to trim low-quality ends from reads and for Illumina adapter removal. Bowtie2 software (Langmead and Salzberg, 2012) was used for alignment of reads to the curated WB transcriptome using the alignment parameters as recommended by the eXpress software (Roberts and Pachter, 2013). Alignment files were piped directly into the eXpress software using the default parameters to generate read-counts. Effective read counts output from the eXpress software were used to calculate pairwise fold-change using the DESeq2 package (Love et al., 2014) of R-bioconductor. A transcript that met the criteria of greater than a 1.5-fold-change difference with a p-value of p < 0.05 after adjusting for multiple comparisons testing was considered a significant differentially expressed gene (DEG).
Gene ontology and enrichment analysis
For Gene Ontology (GO) enrichment analyses, WB transcriptome sequences were uploaded into the OmicsBox program (OmicsBox, 2019). Transcripts were functionally re-annotated by BLASTx (E-value < 1E-3) searches to the nr protein sequence database at the NCBI and mapped to GO terms where available. This list provided a “reference set” for enrichment analyses. Significant DEG lists between each comparison were used as a “test set” for enrichment analysis to the “reference set” by two-sided Fisher’s Exact Test. Enrichment of all the three major GO categories (Biological Process (BP), Molecular Function (MF) and Cellular Component (CC)) was determined. Data was considered statistically significant using a False Discover Rate (FDR) of FDR <5%. Gene Set Enrichment Analysis (GSEA) was also performed as we have previously detailed (Lange et al., 2018; Lange et al., 2019), where significance was set at FDR <5% after 1,000 permutation tests.
Quantitative PCR validation of RNAseq
Reverse transcription quantitative PCR (RT-qPCR) was used to validate RNAseq analyses. Each RNA was quantified by spectrophotometry (Bio-Tek). High-quality RNA samples (RIN >9) from liver tissues extracted previously were used for reverse transcription for all the experimental diet samples. The NEB LunaScript RT SuperMix Kit (Ipswich, MA) was used for the synthesis of cDNA from RNA samples. Total reaction volume was 20 μL, including 4.0 μL of LunaScript RT SuperMix (5X), 200 ng template RNA, and the volume adjusted to 20 μL with nuclease-free water. Reaction conditions were as follows: 25°C for 2 min, followed by 55°C for 10 min and heat inactivation at 95°C for 1 min. After reverse transcription, cDNA was stored at
The expression levels of liver for different diets were evaluated by RT-qPCR using a 384-well Roche LightCycler 480 (Indianapolis, IN). A total of 14 genes were selected to evaluate the DEGs. Primer sets were either designed in-house using the Primer3 software (Untergasser et al., 2012) or were taken from the literature (Cooper et al., 2006; Childress et al., 2016), as indicated (Supplementary Table S1). RT-qPCR was performed in 10 μL reaction in triplicate for each biological replicate. All nine biological replicates were used in the evaluation. Each 10 µL RT-qPCR reaction consisted of 5 μL of the Luna universal qPCR Master Mix (2X), 0.5 μL of each primer (1 μM), 3 μL of nuclease-free water, and 1 μL of cDNA that had been diluted 1:10 with nuclease-free water. Thermal cycling conditions used were as suggested following the Luna NEB protocol: initial denaturation at 95°C for 15 s, followed by 45 cycles of 95°C for 15 s and 60°C for 30 s. Melting curve analyses was performed on all RT-qPCR runs. No RT enzyme controls, and no template controls were included. Two genes (18S rRNA and β-actin) were evaluated for use as housekeeping genes for data normalization. The 18S rRNA gene was selected, as the cycle threshold (Ct) values for all replicates was <1 Ct variation. For relative gene expression analyses, target samples were processed concurrently with the reference gene and the relative expression was given by using the 2−ΔΔCT method (Pfaffl, 2001). Since the standard curve was between 0.9 and 1.0, it was assumed that the accuracy of primer efficiency for computation was 100% for all targets and the reference gene.
Results
Growth and diet utilization
The highest final weight, weight gain, SGR, and lowest FCR were observed for fish fed Diet 4 (PBM + SBM). The lowest final weight, weight gain, SGR, and highest FCR were observed for fish fed Diet 6 (Pro-Cision™ + SBM). Although final individual weight, percentage weight gain, and SGR of fish fed Diet 6 was significantly lower (p ≤ 0.05) than those of fish fed Diet 4, fish fed the other diets (Diets 1–3; Diet 5) did not significantly differ (p > 0.05) from each other (Table 3). Fish fed Diet 5 (all-plant) had significantly lower final weight and SGR compared to fish fed Diet 4 (PBM + SBM), but were not different from fish fed the other test diets. FCR was not significantly different among dietary treatments and averaged 1.74 (Table 3).
TABLE 3. Growth and feed conversion of juvenile white bass (initial average weight, 40.1 ± 1.83 g) fed six diets containing various percentages of fish meal. Means (±SE) within a row with different letters are different (p < 0.05). Data adapted from Rawles et al. (2022).
Differential gene expression
We evaluated differential gene expression (DEG) of all dietary treatments, with all treatments being compared to Diet 1 (FM control), with all comparisons as follows: Diet1:Diet 2, Diet1:Diet 3, Diet1:Diet 4, Diet1:Diet 5, Diet1:Diet 6. We found a total of 1,260 DEGs that were differentially expressed greater than 1.5-fold magnitude in our dietary comparisons (Figures 1A, B; Table 4). For the Diet1:Diet 2 comparison, 293 DEGs were identified, including 148 down-regulated and 145 up-regulated genes (Supplementary Table S2). For the Diet1:Diet 3 comparison, 11 DEGs were identified, including 3 down-regulated and eight up-regulated genes (Supplementary Table S3). For the Diet1:Diet 4 comparison, 52 DEGs were identified, including 37 down-regulated and 15 up-regulated genes (Supplementary Table S4). For the Diet1:Diet 5 comparison, 651 DEGs were identified, including 311 down-regulated and 340 up-regulated genes (Supplementary Table S5). For the Diet1:Diet 6 comparison, 253 DEGs were identified, including 134 down-regulated and 119 up-regulated genes (Supplementary Table S6).
FIGURE 1. The number of differentially expressed genes (DEGs) in the liver of white bass fed a diet containing soybean meal (SBM-Diet 2), poultry by-product meal (PBM- Diet 3), soy protein concentrate (SPC- Diet 4), canola meal (CM- Diet 5), and commercial protein blend (Pro-Cision™- Diet 6) relative to white bass fed a diet composed predominantly of prime FM (Diet 1). (A) The number of differentially expressed genes between experimental diets, the blue filled are down regulated DEGs green filled are upregulated DEGs and the dark blue filled are total DEGs (FDR< 0.05; Fold-change >1.5), and (B) Venn diagram showing the number of unique and shared gene cluster in five experimental diets.
TABLE 4. Differently expressed genes between different dietary treatment combinations (Diets 1–6). All dietary treatments were compared to Diet 1 - fish meal (FM) control. Values indicate contigs/genes passing cutoff values of fold change ≥1.5 (p < 0.05) and read number per contig ≥5. Diet treatments include: Diet: 1 (30%FM-control); 2 ((15%FM + SBM)); 3 (15%FM + PBM + SBM); 4 (PBM + SBM); 5 (All-plant); 6 (Pro-Cision™ + SBM).
GO analysis was performed on the genes underlying each of these diet treatment comparisons previously listed. Genes were mapped to GO terms and then two-sided Fisher’s exact test (FET) was used to assess significance between comparison genes and the WB transcriptome. While significant (FDR <0.05) representation of ontology could not be found using FET for the Diet1:Diet2, Diet1:Diet3, and Diet1:Diet4 comparisons, upon further examination with GSEA, several genes from the Diet1:Diet2 comparison mapped to BP (Table 5) GO terms relating to cell cycle; Diet1:Diet4 for BP GO terms relating to metabolic and cell process and MF GO terms relating to catalytic, cofactor, anion, and small molecule binding (Table 5); Diet1:Diet6 for CC GO terms relating to mitochondrial inner membrane, organelle inner membrane and mitochondrion were downregulated along with MF GO terms of proton transmembrane transporter activity and BP GO terms of iron and transition metal ion transport. Cellular component (Table 5). In Diet1:Diet2, 14 genes were upregulated, and these genes were involved in the GO terms relating to cell cycle BP for the progression in vascular smooth muscle cell (Supplementary Table S7). BP and MF GO terms were upregulated in the Diet1:Diet4 for DNA replication (three genes), Catalytic activity (two genes), and organophosphate metabolic process (three genes) as provided in Supplementary Table S7. In Diet1:Diet6, 40 genes were downregulated with 30 genes corresponds to CC, six genes to BP and four genes to MF (Supplementary Table S7). Significant (FDR <0.05) perturbation of immune pathways (Table 5) based on FET was observed in the Diet1:Diet5 comparison; out of 28 genes underlying the significant GO, 18 genes were upregulated and primarily represented chemokines, erythropoietin, MHC II, toll-like receptors, and insulin-like growth factor, while ten genes were downregulated. GO analyses revealed genes underlying the Diet1:Diet6 comparison involved in the glutathione biosynthetic process and non-ribosomal peptide biosynthetic process (Table 5); out of twenty genes, six genes were upregulated and 14 genes were downregulated (Supplementary Table S7). These pathways are involved in cellular metabolism, and more specifically in amino acid biosynthesis and protein translation.
TABLE 5. Significant descriptions for GO categories within each of the dietary comparisons found to be significantly different than the reference transcriptome following GO mapping, Gene Set Enrichment Analysis (GSEA), and a two-sided Fisher’s exact test (FET) are listed below. A complete listing of all descriptions is provided in Supplementary Tables S2–S6. Diet: 1 (30%FM-control); 2 ((15%FM + SBM)); 3 (15%FM + PBM + SBM); 4 (PBM + SBM); 5 (All-plant); 6 (Pro-Cision™ + SBM). As per GSEA, Diet1:Diet3 and Diet1:Diet5 were not significant whereas in case of FET, Diet1: Diet2; Diet1:Diet3 and Diet1:Diet4 were not significant.
To validate our differential gene expression analyses using RNA sequencing, seven growth-related genes and seven immune related genes (Figures 2, 3) were analyzed by RT-qPCR for all liver samples. For growth-related genes, myogenin (myog) expression in all fish was upregulated in the Diet 1 (30% FM-control) group relative to all other dietary treatments, with similar results observed for 40S ribosomal protein S19 (RPS19) (Figure 2). Serine/threonine-protein kinase PLK2 isoform 1 (PLK2) was significantly upregulated in fish fed Diet 5 (all-plant) and downregulated in fish fed Diet 2 (15%FM + SBM) and Diet 6 (Pro-Cision™ + SBM). Dipeptidyl peptidase 4-like (DPP4), glutathione S-transferase omega-1 (GSTF6), WD repeat and HMG-box DNA-binding protein1 (WDHD1), and olfactory marker protein (OMP) were significantly upregulated when compared to fish fed Diet 1 (30% FM-control). Immune related genes such as hepcidin (HAMP), heat shock protein 70 (HSP70), interleukin-6 receptor subunit beta-like (IL6), immunoglobulin delta heavy chain (IGHD), MHC class II antigen beta chain (MHC2-Ab1), and C-X-C chemokine receptor type 2-like (CXC) were significantly upregulated in fish fed Diet 5 (all-plant protein) except insulin-like growth factor II (IGF2) (Figure 3). In fish fed Diet 6 (Pro-Cision™ + SBM), IGF2 and CXC were significantly downregulated whereas others were upregulated with a significant upregulation in HAMP. Overall, all RNA sequencing data and RT-qPCR results were correlated (Table 6).
FIGURE 2. Relative mRNA expression of growth-related genes such as myogenin (MYOG), dipeptidyl peptidase 4-like (DPP4), glutathione S-transferase omega-1 (GSTF6), WD repeat and HMG-box DNA-binding protein1 (WDHD1), olfactory marker protein (OMP), serine/threonine-protein kinase PLK2 isoform 1 (PLK2) and 40S ribosomal protein S19 (RPS19) in liver in white bass fed with different experimental diets as measured by quantitative real-time PCR. Data are presented as mean ± SE (n = 9). Different alphabetical letters indicate significant difference with p < 0.05.
FIGURE 3. Relative mRNA expression of immune related genes such as hepcidin (HAMP), insulin-like growth factor II (IGF2), heat shock protein 70 (HSP70), interleukin-6 receptor subunit beta-like (IL6), immunoglobulin delta heavy chain (IGHD), MHC class II antigen beta chain (MHC2-Ab1) and C-X-C chemokine receptor type 2-like (CXC) in liver in white bass fed with different experimental diets as measured by quantitative real-time PCR. Data are presented as mean ± SE (n = 9). Different alphabetical letters indicate the significant difference with p < 0.05.
TABLE 6. Comparison of RNA-seq and RT-PCR fold change in the liver for the diets (1–6).
Discussion
RNA sequencing offers insights into an organism’s gene expression for a specific time and can be a useful analytical tool for identifying gene pathways affected by dietary changes in aquaculture species (Chávez-Calvillo et al., 2010; Shiel et al., 2017; Schroeter et al., 2018; Meiler and Kumar, 2021; Lorenz et al., 2022). This is the first published study to evaluate partial to total replacement of FM in diets for WB and their corresponding effects on gene expression. Rawles et al. (2022) reported that growth performance of WB fed an all-plant diet and a diet containing a commercially-blended protein product were reduced compared to fish fed diets containing FM (control) and a FM-free diet containing SBM and poultry by-product meal (PBM). The aim of this study was to utilize RNA sequencing to evaluate possible differences in gene expression in WB to changes in diet. Mechanisms and diet-related gene pathways that may be responsible for any observed changes are discussed.
Differentially expressed genes compared to diet 1 (30%FM-control)
Changes in the liver transcriptome with replacement of FM by plant protein sources have been reported in previous studies on several fish species (Geay et al., 2011; Tacchi et al., 2012; Juárez et al., 2021). This is the first published report in which differential liver gene expression of WB fed diets containing different protein sources was analyzed. The study was focused on the liver because it is the main organ involved in the use of nutrients and is the main center of intermediary metabolism in animals (Panserat et al., 2009). In the present study, there were more differentially expressed genes (DEGs) in comparisons to Diet 1 (30%FM-control) when the test diet contained no FM and/or high levels of SBM (Diets 2, 5, and 6; 293, 651, 253 DEGs, respectively). The exception to this is Diet 4 (PBM + SBM), which replaced FM with PBM, where there were only 52 DEGs (Table 4).
To further analyze each of the dietary comparisons, we looked for DEGs potentially involved in metabolic pathways and oxidative stress, immune function, and apoptosis. While there were certainly many more functions impacted, given the nature of the study and the composition of test diets, it seemed most prudent to focus on these areas. For a more detailed look at all DEGs by test diet comparisons, refer to the Supplementary Tables S2–S6.
Metabolism
The capability of the liver to maintain the balance between the intake and release of glucose plays an essential role in the maintenance of the homeostasis of glucose in the blood (Wu et al., 2001). In this study, glucose and fatty acid metabolism DEGs were mainly affected by the inclusion level of SBM and/or the absence of animal proteins. Specific DEGs in each of these test comparisons using the control diet (Diet 1) as baseline will be discussed in more detail in the following sections.
Diet 1 (30%FM-control) vs. Diet 2 (15%FM + SBM)
In this test comparison, several DEGs involved in metabolism were identified in fish fed Diet 2 relative to fish fed Diet 1. The most striking result is the decrease in insulin (INS; -57.33-fold) in fish fed Diet 2 relative to fish fed Diet 1, with this being the largest DEG by magnitude in this comparison. This could demonstrate higher circulating levels of glucose in fish fed Diet 2, necessitating the higher insulin demand. Diet 2 contained the highest amount of SBM inclusion (Table 1) of any of the test diets, and clearly replacement of FM with SBM decreased plasma insulin, as has been observed in other fish species (Zhang et al., 2019). Additionally, 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 4 (PFKFB4; −3.12) can reduce hepatic glucose production by increasing glycolysis and inhibiting gluconeogenesis, thereby lowering blood glucose, indicating possible lower glucose demand in fish fed Diet 2 (Wu et al., 2001). Interestingly, this gene is also associated with the adaptation of carnivorous fish to diets including plant meals (Metón et al., 2000). Glucose-6-phosphate 1-dehydrogenase (G6PD; 1.68-fold) is the rate-limiting enzyme of the oxidative pentose phosphate pathway (PPP) (Ge et al., 2020). G6PD gene expression, which plays a role in glucose oxidation and increased gluconeogenesis, has been previously shown to be increased in fish after high-carbohydrate feeding (Jiang et al., 2019). Overall, fish fed Diet 2 had significant differential expression of genes involved in metabolic activity, particularly glucose metabolism (Supplementary Table S2).
Diet 1 (30%FM-control) vs. Diet 3 (15%FM + PBM + SBM)
In this test comparison, only 11 genes were differentially expressed in fish fed Diet 3 relative to fish fed Diet 1. Only 2 DEGs are associated with metabolism, including dipeptidyl peptidase 4 (DPP4; 3.25-fold) and trehalase (TREH; 29.33-fold). DPP4 is a well-characterized gene which regulates peptide hormones and regulates the endocrine pathway, and when inhibited enhances postprandial metabolism and can be used to treat hyperglycemia (Trzaskalski et al., 2020).
Diet 1 (30%FM-control) vs. Diet 4 (PBM + SBM)
In this test comparison, only 52 genes were differentially expressed in the FM diet relative to FM-free diet, including insulin receptor substrate 2 (IRS2; −1.56-fold) and glucose-6-phosphate 1-dehydrogenase (G6PD; 1.80-fold). These genes are involved in insulin regulation and glucose metabolism. IRS2 has been shown to control adipose tissue function and glucose metabolism as well as mediating insulin’s actions (Turaihi et al., 2018). With IRS2 being downregulated in fish fed Diet 4 relative to fish fed Diet 1, this may indicate that there could be adipose tissue activity in Diet 4-fed fish. As stated earlier, G6PD plays a role in glucose oxidation and gluconeogenesis (Metón et al., 2000), which given the upregulation in fish fed Diet 4 could indicate higher metabolic activity. There was also a higher expression of olfactory marker protein (O2; 3.57-fold) in this study which indicates the acceptability of the SPC, CM and commercial protein blend. Previous studies related to plant diet indicated the higher expression of olfactory proteins raises the potential that various extrinsic variables involving olfactory sensing such as locating food, hunger pains, and even the need to avoid predators, play a substantial part in the control of eating behavior in rainbow trout (Heraud et al., 2022).
Diet 1 (30%FM-control) vs. Diet 5 (all-plant)
In this test comparison, several DEGs involved in metabolism were identified in fish fed Diet 5 relative to fish fed Diet 1, including acyl-CoA synthetase family member 2 (ACSS2; -2.77-fold), insulin-like growth factor II (IGFII; −1.84-fold), ATP synthase F(0) complex subunit C3 (ATP5G3; −1.59-fold), acyl-CoA desaturase (FADS2; 6.98-fold), ATP-citrate synthase (ACLY; 5.70-fold), insulin-like growth factor-binding protein 3 (IGFBP3; 2.98-fold), and glucose-6-phosphate 1-dehydrogenase (G6PD; 2.34-fold). ACSS2 plays a key role in lipogenesis (Xu et al., 2018). IGFII plays a key role in cell regulation and growth, and interestingly, while there was a clear difference in expression in the current study between fish fed Diet 1 and fish fed Diet 5 (all-plant), previous studies in tilapia have demonstrated that IGFII was not influenced by diet (Herath et al., 2018). It should be noted, however, that Herath et al. (2018) formulated their FM-free diet to contain PBM; thus, it was not an all-plant diet as in the present study. FADS2 is the key enzyme in long-chain polyunsaturated fatty acid (LC-PUFA) biosynthesis. Importantly to fish fatty acid metabolism, ACLY is involved in the production of acetyl-CoA, before FADS2 converts dietary PUFAs, such as linoleic acid and alpha-linolenic acid into HUFAs, such as DHA.
Diet 1 (30%FM-control) vs. Diet 6 (Pro-Cision™ + SBM)
In this test comparison, several DEGs involved in metabolism were identified in fish fed Diet 6 relative to fish fed Diet 1, including fructose-1,6-bisphosphatase 1 (FBP1; -2.34-fold), several cytochrome c oxidases (CCO subunit 1; -1.79-fold; subunit 4, -1.75-; subunit 2, -1.63-fold), several mitochondrial NADH dehydrogenases (subunit 2; -1.82-fold; subunit 5; -1,79-fold), ATP synthase F(0) complex subunit C2 (ATP5MC2; -1.56-fold), glucose-6-phosphate 1-dehydrogenase (G6PD; 2.18-fold), and electron transfer flavoprotein-ubiquinone oxidoreductase (ETFQO; 1.7-fold). CCOs, NADH dehydrogenases, and ATP synthases are involved in the respiratory electron transport chain, while G6PD is involved in the PPP and FBP1 is a key protein in gluconeogenesis via glycolysis. Fish fed Diet 6 had much lower metabolic activity relative to fish fed Diet 1, but exhibited higher glycolytic activity, while the latter fish had higher PPP and ETFQO activity. ETFQO is a key component of the electron transport chain, and it acts as a key linkage between the oxidation of fatty acids and some amino acids to the mitochondrial respiratory chain (Zhang et al., 2006).
Oxidative stress, immune function, and apoptosis
Multiple links have been made between diet, growth performance, and immune response capability (Burrells et al., 2001; Hossain et al., 2016; Kader et al., 2018). It is in this context that we take a closer look at DEGs within each of these test diet comparisons to the control diet (Diet 1).
One of the initial lines of defense used by fish against infections consisted of small peptides known as antimicrobial peptides (AMPs). Fish hepcidins have been discovered in at least 37 fish species since they were initially discovered and isolated for the first time from the hybrid striped bass (Hirono et al., 2005). The liver and kidneys have high levels of fish hepcidin expression, while the heart, brain, skin, gut, blood cells, spleen, stomach, and muscle have low levels of expression. These hepcidin are short, cysteine-rich peptide discovered in fish in two forms, HAMP1 and HAMP2 (Neves et al., 2017). Hepcidin was found to be upregulated in almost all the diet groups including poultry by-product meal (PBM), soy protein concentrate (SPC), canola meal (CM) and a commercial protein blend (Pro-Cision™), except diet 2 soybean meal (SBM) (Table 6). Previous study has shown that these bass hepcidin can significantly reduce bacteria, fungi and yeast (Lauth et al., 2005; Shen et al., 2019). Further, numerous physiological and pathological processes have been linked to ribosomal proteins, and many of them act as antimicrobial proteins in the body’s response to infections. Ribosomal proteins can function as antimicrobial proteins in eukaryotic innate immune systems. The 40S19 was shown to be significantly downregulated in all the test diets.
Diet 1 (30%FM-control) vs. Diet 2 (15%FM + SBM)
In this test comparison, several DEGs involved in immune function and apoptosis were identified in fish fed Diet 2 relative to fish fed Diet 1, with a few of the key DEGs discussed below. The complement regulatory protein DAF blocks complement deposition on host cells, which also acts as a regulator of phagocytosis of cells (Dernstedt et al., 2021). Heat shock protein 70 (HSP70; -2.20-fold) and hypoxia up-regulated protein 1 (HYOU1; -1.84-fold) are key regulators of cellular stress. HSP70 has been well studied in many animal models and plays a key role in protein folding and has an important role in the response to oxidative stress (Turner et al., 1999), and its overexpression usually indicates that cells are suffering from oxidative stress. HYOU1 can induce an unfolded protein response in the endoplasmic reticulum as a response to cellular stress (Rao et al., 2021). Apoptosis-inducing factor 3 (AIFM3; 1.88-fold) is part of a cascade that contributes to mitochondrial cell death and participates in the assembly of the respiratory chain as part of normal mitochondrial bioenergetics (Bano and Prehn, 2018). However, long term exposure to oxidative stress has also been shown to induce apoptosis (Dastor and Dreyer, 2001).
Diet 1 (30%FM-control) vs. Diet 4 (PBM + SBM)
In this test comparison, several DEGs involved in immune function and apoptosis were identified in fish fed Diet 4 relative to fish fed Diet 1, a few of which will be discussed further below. Spondin-2 (SPON2; -11.13-fold) expression identifies pathogens and triggers an immune response (Klar et al., 1992). Von Willebrand factor A domain-containing protein 7 (VWA7; -9.84-fold) is part of the major histocompatibility complex class III region, which has traditionally been associated with the inflammation response region of the MHC (Deakin et al., 2006). MPV17-like (MPV17L; -5.94-fold) is located in the membrane of the peroxisome and is involved in the regulation of several reactive oxygen species, including glutathione peroxidase and several catalase genes (Iida et al., 2006). Beta-galactoside-binding lectin-like (BGBP; -3.43-fold) is a negative regulator of the cell cycle and plays a role in T cell immune response (Allione et al., 1998). Each of these immune and cellular stress genes is more upregulated in fish fed Diet 4 than in fish fed Diet 1, indicating that these fish were experiencing cellular stress, possibly due to the presence of dietary factors such as ANFs. Immune response to ANFs has been demonstrated in numerous studies involving alternative formulations, including soy and alternative animal protein diets. Caballero-Solaros et al. (2018) demonstrated changes in numerous immune response genes in Atlantic salmon fed FM, animal byproduct meal, or soy protein concentrate, with more genes differentially expressed in fish fed the plant protein diet. Likewise, Juárez et al. (2021) found higher activation of immune response genes in juvenile Pacific yellowtail (Seriola lalandi) fed a 25% SBM feed compared to a FM feed.
Diet 1 (30%FM-control) vs. Diet 5 (all-plant)
In this test comparison, several DEGs involved in immune function and apoptosis were identified in fish fed Diet 5 relative to fish fed Diet 1, a few of which will be discussed further below. Tumor necrosis factor ligand superfamily member 12 (TNFSF12; -2.35-fold) works with other proapoptotic TNFSF ligands to facilitate cytotoxicity in many cell types, including activated monocytes, dendritic cells, NK cells, and T cells (Wang et al., 2013). Death-associated protein 1 (DAP1; -1.71-fold) is known to function both in apoptosis and autophagy (Xia et al., 2017). Oxidation resistance protein 1 (OXR1; -1.61-fold) is a gene which has previously been shown to be induced under oxidative stress and is vital for defense against reactive oxygen species (Oliver et al., 2011). Complement C4 (C4; 1.59-fold) is one of the chief constituents of innate immunity which participates in a pathway that is responsible for the immediate recognition and elimination of invading microbes and plays an essential role in both the classical and lectin complement pathways (Wang and Liu, 2021). Macrophage-expressed gene 1 protein (MPEG1; 1.68-fold) is responsible for encoding a protein responsible for host defense against invading pathogens (Benard et al., 2015). Caspase recruitment domain-containing protein 19 (CARD19; 1.69-fold) appears to be involved in the regulation of mitochondrial function which affects pro-inflammatory innate immune responses (Rios et al., 2018). Toll-like receptor 8 (TLR8; 1.82-fold) is one of many toll-like receptors that function by sensing distinct pathogen associated molecular patterns and initiate inflammatory reactions important for innate and adaptive immune responses (Moen et al., 2019). The major histocompatibility complex (MHC) class II molecules (MHC class II; 2.43- fold) involved in the presenting of the processed antigens to the CD4(+) T-lymphocytes for the commencement of the antigen-specific immune response (Deakin et al., 2006) (Table 5 and Supplementary Table S7).
In the current study, there was significantly higher expression of immune genes with the GO analyzed by GSEA. The increased expression of immune related genes in the present study indicates the activation of innate and adaptive immune system, as indicated by significant mRNA expression of genes related to inflammatory responses, including proinflammatory cytokines of TNF-α, IL1β, IL6, and IL8 as these act as biomarkers for the activation of inflammatory responses (Han et al., 2022). Further, the two highest upregulated genes in WB fed all plant proteins were claudin-15 (∼26-fold) and matrix metalloprotease 16 (∼45-fold), which are responsible for tight junction permeability and glucose metabolism, as well as tissue remodeling and immunity. Over our results indicate that an all-plant protein diet greatly perturbs the white bass immune system and inflammasome.
Diet 1 (30%FM-control) vs. Diet 6 (Pro-Cision™ + SBM)
In this test comparison, several DEGs involved in immune function and apoptosis were identified in fish fed Diet 6 relative to fish fed Diet 1, some of which will be discussed below. As mentioned earlier, Complement C4 (C4; -1.72-fold) is one of the chief constituents of the complement immune response and apoptosis-inducing factor 3 (AIF3; 1.86-fold) is part of a cascade that contributes to mitochondrial cell death. Complement C5 (C5; -1.65-fold) plays an important role in the inflammatory and cell killing processes of the complement system. Complement C1 is a large complex within the complement that triggers the destruction of invading pathogens via lysis or by stimulation of innate and adaptive immune processes, with C1r-A subcomponent (C1r; -1.73-fold) being the auto-activation subunit following pathogen binding (Almitairi et al., 2018). All of these genes were upregulated in fish fed Diet 6 except AIF3. As mentioned previously, numerous other studies demonstrating fish response to ANFs have been published, including those mentioned earlier in the Diet1:Diet5 comparison.
Conclusion
Transcriptome analysis of the liver of WB revealed differentially expressed genes with unique gene ontology when fish were fed various alternative diets as compared to a standard, high FM diet. Juvenile WB can be fed diets without fish meal without negative effects on growth, survival, and body composition. However, WB fed an all-plant diet (Diet 5) or a protein blended diet (Diet 6) to replace FM had reduced growth, and diet greatly affected the immune system of fish fed these diets. Partial replacement diets and a diet formulated similar to a commercial diet showed an increased biological process involving metabolic and cellular process related to nutrient utilization. Olfactory genes may be responsive to diet intake and in turn growth differences observed, of which further research on palatability for white bass for boosting feed intake is warranted. These factors, and others reported in the present study, indicate possible reasons for the reduced growth observed in fish fed an all-plant diet and a diet with a commercial protein blended product.
Data availability statement
The datasets generated for this study can be found in the NCBI Gene Expression Omnibus (GEO) repository and can be accessed under the accession number GSE220235. All other data that support the findings of this study have been included in the manuscript and Supplementary Material.
Ethics statement
The animal studies were approved by USDA-ARS-HKDSNARC institutional animal care and use committee (IACUC). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
SF: Formal Analysis, Supervision, Validation, Writing–review and editing. JA: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Supervision, Validation, Writing–original draft, Writing–review and editing. NS: Formal Analysis, Investigation, Validation, Writing–original draft, Writing–review and editing. BB: Funding acquisition, Project administration, Resources, Supervision, Writing–review and editing. SR: Conceptualization, Investigation, Methodology, Resources, Writing–review and editing. BG: Investigation, Resources, Writing–review and editing. KR: Investigation, Methodology, Writing–review and editing. MM: Formal Analysis, Investigation, Methodology, Project administration, Supervision, Writing–review and editing. GH: Conceptualization, Formal Analysis, Investigation, Methodology, Writing–review and editing. CW: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This project was supported by funds appropriated to the United States Department of Agriculture (ARS Research Projects #6028-31630-009-000-D and #6010-32000-027-000-D). This research was supported, in part, by an appointment (Research Fellowship to NS) to the Agricultural Research Service (ARS) Research Participation Program administered by the Oak Ridge Institute for Science and education (ORISE) through an interagency agreement between the U.S. Department of Energy (DOE) and the U.S. Department of Agriculture (USDA). ORISE is managed by Oak Ridge Associated Universities (ORAU) under DOE contract number DE-SC0014664.
Acknowledgments
The authors thank Sarah Bennett and Rebecca Furbeck for assistance in processing the feeds and H.J. Baker Bros., Tuscola, TX, United States for providing the Pro-Cision™ used in this research. The authors are grateful for the technical staff at HKDSNARC, especially Matt Barnett, who assisted in all aspects of sample collection, and sample preparation.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
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.
Author disclaimer
All opinions expressed in this paper are the author’s and do not necessarily reflect the policies and views of USDA, DOW, or ORAU/ORISE. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture. The USDA is an equal opportunity provider and employer.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1308690/full#supplementary-material
References
Akbari A., Wu J. (2015). An integrated method of isolating napin and cruciferin from defatted canola meal. LWT-Food Sci. Technol. 64 (1), 308–315. doi:10.1016/j.lwt.2015.05.046
Akiyama D. M., Coelho S. R., Lawrence A. L., Robinson E. H. (1989). Apparent digestibility of feedstuffs by the marine shrimp Penaeus vannamei BOONE. J. Jpn. Soc. Fish. Sci. 55 (1), 91–98. doi:10.2331/suisan.55.91
Allione A., Wells V., Forni G., Mallucci L., Novelli F. (1998). β-Galactoside-Binding protein (βGBP) alters the cell cycle, up-regulates expression of the α- and β-chains of the IFN-γ receptor, and triggers IFN-γ-Mediated apoptosis of activated human T lymphocytes. J. Immunol. 161 (5), 2114–2119. doi:10.4049/jimmunol.161.5.2114
Almitairi J. O. M., Venkatraman Girija U., Furze C. M., Simpson-Gray X., Badakshi F., Marshall J. E., et al. (2018). Structure of the C1r–C1s interaction of the C1 complex of complement activation. Proc. Natl. Acad. Sci. 115 (4), 768–773. doi:10.1073/pnas.1718709115
Bano D., Prehn J. H. M. (2018). Apoptosis-inducing factor (AIF) in physiology and disease: the tale of a repented natural born killer. EBioMedicine 30, 29–37. doi:10.1016/j.ebiom.2018.03.016
Benard E. L., Racz P. I., Rougeot J., Nezhinsky A. E., Verbeek F. J., Spaink H. P., et al. (2015). Macrophage-expressed perforins mpeg1 and mpeg1. 2 have an anti-bacterial function in zebrafish. J. innate Immun. 7 (2), 136–152. doi:10.1159/000366103
Burrells C., Williams P. D., Forno P. F. (2001). Dietary nucleotides: a novel supplement in fish feeds. Aquaculture 199 (1–2), 159–169. doi:10.1016/S0044-8486(01)00577-4
Caballero-Solares A., Xue X., Parrish C. C., Foroutani M. B., Taylor R. G., Rise M. L. (2018). 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 genomics 19, 796–826. doi:10.1186/s12864-018-5188-6
Chávez-Calvillo G., Perez-Rueda E., Lizama G., Zúñiga Aguilar J., Gaxiola G., Cuzon G., et al. (2010). Differential gene expression in Litopenaeus vannamei shrimp in response to diet changes. Aquaculture 300 (1–4), 137–141. doi:10.1016/j.aquaculture.2009.11.027
Childress C. J., Fuller S., Rawles S., Beck B., Gaylord T., Barrows F., et al. (2016). Lysine supplementation of commercial fishmeal-free diet in hybrid striped bass Morone chrysops x M. saxatilis affects expression of growth-related genes. Aquac. Nutr. 22 (4), 738–744. doi:10.1111/anu.12300
Company R., Astola A., Pendón C., Valdivia M. M., Pérez-Sánchez J. (2001). Somatotropic regulation of fish growth and adiposity: growth hormone (GH) and somatolactin (SL) relationship. Comp. Biochem. Physiology Part C Toxicol. Pharmacol. 130 (4), 435–445. doi:10.1016/S1532-0456(01)00269-1
Cooper C. A., Handy R. D., Bury N. R. (2006). The effects of dietary iron concentration on gastrointestinal and branchial assimilation of both iron and cadmium in zebrafish (Danio rerio). Aquat. Toxicol. 79 (2), 167–175. doi:10.1016/j.aquatox.2006.06.008
Coyle S. D., Tidwell J. H., Webster C. D. (2000). Response of largemouth bass Micropterus salmoides to dietary supplementation of lysine, methionine, and highly unsaturated fatty acids. J. World Aquac. Soc. 31 (1), 89–95. doi:10.1111/j.1749-7345.2000.tb00702.x
Dairiki J. K., Dias C. T., dos S., Cyrino J. E. P. (2007). Lysine requirements of largemouth bass, Micropterus salmoides: a comparison of methods of analysis of dose-response trials data. J. Appl. Aquac. 19 (4), 1–27. doi:10.1300/J028v19n04_01
Dastoor Z., Dreyer J.-L. (2001). Potential role of nuclear translocation of glyceraldehyde-3-phosphate dehydrogenase in apoptosis and oxidative stress. J. cell Sci. 114 (9), 1643–1653. doi:10.1242/jcs.114.9.1643
Deakin J. E., Papenfuss A. T., Belov K., Cross J. G. R., Coggill P., Palmer S., et al. (2006). Evolution and comparative analysis of the MHC Class III inflammatory region. Bmc Genomics 7 (1), 281–314. doi:10.1186/1471-2164-7-281
Dernstedt A., Leidig J., Holm A., Kerkman P. F., Mjösberg J., Ahlm C., et al. (2021). Regulation of decay accelerating factor primes human germinal center B cells for phagocytosis. Front. Immunol. 11, 599647. doi:10.3389/fimmu.2020.599647
Dossou S., Koshio S., Ishikawa M., Yokoyama S., Dawood M. A., El Basuini M. F., et al. (2018). Effect of partial replacement of fish meal by fermented rapeseed meal on growth, immune response and oxidative condition of red sea bream juvenile, Pagrus major. Aquaculture 490, 228–235. doi:10.1016/j.aquaculture.2018.02.010
Fox J. M., Lawrence A. L., Smith F. (2004). Development of a low-fish meal feed formulation for commercial production of Litopenaeus vannamei. Av. En. Nutr. Acuicola.
Francis G., Makkar H. P. S., Becker K. (2001). Antinutritional factors present in plant-derived alternate fish feed ingredients and their effects in fish. Aquaculture 199 (3–4), 197–227. doi:10.1016/S0044-8486(01)00526-9
Fu L., Niu B., Zhu Z., Wu S., Li W. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28 (23), 3150–3152. doi:10.1093/bioinformatics/bts565
Gatlin D. M., Barrows F. T., Brown P., Dabrowski K., Gaylord T. G., Hardy R. W., et al. (2007). Expanding the utilization of sustainable plant products in aquafeeds: a review. Aquac. Res. 38 (6), 551–579. doi:10.1111/j.1365-2109.2007.01704.x
Gaylord T. G., Rawles S. D., Gatlin D. M. (2004). Amino acid availability from animal, blended, and plant feedstuffs for hybrid striped bass (Morone chrysops× M. saxatilis). Aquac. Nutr. 10 (5), 345–352. doi:10.1111/j.1365-2095.2004.00310.x
Ge T., Yang J., Zhou S., Wang Y., Li Y., Tong X. (2020). The role of the pentose phosphate pathway in diabetes and cancer. Front. Endocrinol. 11, 365. doi:10.3389/fendo.2020.00365
Geay F., Ferraresso S., Zambonino-Infante J. L., Bargelloni L., Quentel C., Vandeputte M., et al. (2011). Effects of the total replacement of fish-based diet with plant-based diet on the hepatic transcriptome of two European sea bass (Dicentrarchus labrax) half-sibfamilies showing different growth rates with the plant-based diet. BMC genomics 12 (1), 522–618. doi:10.1186/1471-2164-12-522
Gilani G. S., Cockell K. A., Sepehr E. (2005). Effects of antinutritional factors on protein digestibility and amino acid availability in foods. J. AOAC Int. 88 (3), 967–987. doi:10.1093/jaoac/88.3.967
Glencross B., Hawkins W., Curnow J. (2004). Nutritional assessment of Australian canola meals. I. Evaluation of canola oil extraction method and meal processing conditions on the digestible value of canola meals fed to the red seabream (Pagrus auratus, Paulin). Aquac. Res. 35 (1), 15–24. doi:10.1111/j.1365-2109.2004.00974.x
Gómez-Requeni P., Mingarro M., Calduch-Giner J., Médale F., Martin S., Houlihan D., et al. (2004). Protein growth performance, amino acid utilisation and somatotropic axis responsiveness to fish meal replacement by plant protein sources in gilthead sea bream (Sparus aurata). Aquaculture 232 (1–4), 493–510. doi:10.1016/S0044-8486(03)00532-5
Han Y.-K., Xu Y. C., Luo Z., Zhao T., Zheng H., Tan X. Y. (2022). Fish meal replacement by mixed plant protein in the diets for juvenile yellow catfish Pelteobagrus fulvidraco: effects on growth performance and health status. Aquac. Nutr. 2022, 2677885. doi:10.1155/2022/2677885
Herath S. S., Haga Y., Satoh S. (2018). Interactive effects of salinity and complete fishmeal replacement on growth, food consumption, and gene expression of hepatic IGF-I, IGF-II and growth hormone receptors in Nile tilapia, Oreochromis niloticus (L). Aquac. Res. 49 (6), 2128–2139. doi:10.1111/are.13667
Heraud C., Hirschinger T., Baranek E., Larroquet L., Surget A., Sandres F., et al. (2022). Detection and modulation of olfactory sensing receptors in carnivorous rainbow trout (Oncorhynchus mykiss) fed from first feeding with plant-based diet. Int. J. Mol. Sci. 23, 2123. doi:10.3390/ijms23042123
Higgs D. A., McBride J. R., Markert J. R., Dosanjh B. S., Plotnikoff M., Clarke W. (1982). Evaluation of Tower and Candle rapeseed (canola) meal and Bronowski rapeseed protein concentrate as protein supplements in practical dry diets for juvenile chinook salmon (Oncorhynchus tshawytscha). Aquaculture 29 (1–2), 1–31. doi:10.1016/0044-8486(82)90030-8
Hirono I., Hwang J. Y., Ono Y., Kurobe T., Ohira T., Nozaki R., et al. (2005). Two different types of hepcidins from the Japanese flounder Paralichthys olivaceus. FEBS J. 272 (20), 5257–5264. doi:10.1111/j.1742-4658.2005.04922.x
Hossain M. S., Koshio S., Ishikawa M., Yokoyama S., Sony N. M. (2016). Dietary nucleotide administration influences growth, immune responses and oxidative stress resistance of juvenile red sea bream (Pagrus major). Aquaculture 455, 41–49. doi:10.1016/j.aquaculture.2016.01.008
Iida R., Yasuda T., Tsubota E., Takatsuka H., Matsuki T., Kishi K. (2006). Human Mpv17-like protein is localized in peroxisomes and regulates expression of antioxidant enzymes. Biochem. biophysical Res. Commun. 344 (3), 948–954. doi:10.1016/j.bbrc.2006.04.008
Jiang G. Z., Shi H. J., Xu C., Zhang D. D., Liu W. B., Li X. F. (2019). Glucose-6-phosphate dehydrogenase in blunt snout bream Megalobrama amblycephala: molecular characterization, tissue distribution, and the responsiveness to dietary carbohydrate levels. Fish. Physiol. Biochem. Feb 45 (1), 401–415. doi:10.1007/s10695-018-0572-3
Joana G. da S., Oliva-Teles A., Oliva-Teles A. (1998). Apparent digestibility coefficients of feedstuffs in seabass (Dicentrarchus labrax) juveniles. Aquat. Living Resour. 11 (3), 187–191. doi:10.1016/S0990-7440(98)80115-0
Juárez O. E., Lafarga-De la Cruz F., Lazo J. P., Delgado-Vega R., Chávez-García D., López-Landavery E., et al. (2021). Transcriptomic assessment of dietary fishmeal partial replacement by soybean meal and prebiotics inclusion in the liver of juvenile Pacific yellowtail (Seriola lalandi). Mol. Biol. Rep. 48, 7127–7140. doi:10.1007/s11033-021-06703-4
Kader M. A., Bulbul M., Abol-Munafi A. B., Asaduzzaman M., Mian S., Noordin N. B. M., et al. (2018). Modulation of growth performance, immunological responses and disease resistance of juvenile Nile tilapia (Oreochromis niloticus)(Linnaeus, 1758) by supplementing dietary inosine monophosphate. Aquac. Rep. 10, 23–31. doi:10.1016/j.aqrep.2018.03.003
Kaushik S. J. (2002). “European sea bass, Dicentrarchus labrax,” in Nutrient requirements and feeding of finfish for aquaculture (Wallingford, United Kingdom: CABI Publishing), 28–39. doi:10.1079/9780851995199.0028
Kent W. J. (2002). BLAT--the BLAST-like alignment tool. Genome Res. 12 (4), 656–664. doi:10.1101/gr.229202
Klar A., Baldassare M., Jessell T. M. (1992). F-spondin: a gene expressed at high levels in the floor plate encodes a secreted protein that promotes neural cell adhesion and neurite extension. Cell 69 (1), 95–110. doi:10.1016/0092-8674(92)90121-R
Krueger F., James F., Ewels P., Afyounian E., Schuster-Boeckler B. (2021) 'FelixKrueger/TrimGalore: v0.6.7' - via Zenodo (0.6.7). doi:10.5281/zenodo.5127899
Kumar V., Khalil W. K. B., Weiler U., Becker K. (2011). Influences of incorporating detoxified Jatropha curcas kernel meal in common carp (Cyprinus carpio L.) diet on the expression of growth hormone-and insulin-like growth factor-1-encoding genes. J. animal physiology animal Nutr. 97 (1), 97–108. doi:10.1111/j.1439-0396.2011.01247.x
Lange M. D., Abernathy J., Farmer B. D. (2019). Evaluation of a recombinant Flavobacterium columnare DnaK protein vaccine as a means of protection against columnaris disease in channel catfish (Ictalurus punctatus). Front. Immunol. 10, 1175. doi:10.3389/fimmu.2019.01175
Lange M. D., Farmer B. D., Abernathy J. (2018). Catfish mucus alters the Flavobacterium columnare transcriptome. FEMS Microbiol. Lett. 365 (22), fny244. doi:10.1093/femsle/fny244
Langmead B., Salzberg S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9 (4), 357–359. doi:10.1038/nmeth.1923
Lauth X., Babon J. J., Stannard J. A., Singh S., Nizet V., Carlberg J. M., et al. (2005). Bass hepcidin synthesis, solution structure, antimicrobial activities and synergism, and in vivo hepatic response to bacterial infections. J. Biol. Chem. 280 (10), 9272–9282. doi:10.1074/jbc.M411154200
Lemos D., Lawrence A. L., Siccardi A. J. (2009). Prediction of apparent protein digestibility of ingredients and diets by in vitro pH-stat degree of protein hydrolysis with species-specific enzymes for juvenile Pacific white shrimp Litopenaeus vannamei. Aquaculture 295 (1–2), 89–98. doi:10.1016/j.aquaculture.2009.06.011
Li C., Beck B. H., Fuller S. A., Peatman E. (2014). Transcriptome annotation and marker discovery in white bass (Morone chrysops) and striped bass (Morone saxatilis). Anim. Genet. 45 (6), 885–887. doi:10.1111/age.12211
Lim C., Klesius P. H., Higgs D. A. (1998). Substitution of canola meal for soybean meal in diets for channel catfish Ictalurus punctatus. J. World Aquac. Soc. 29 (2), 161–168. doi:10.1111/j.1749-7345.1998.tb00975.x
Lorenz E. K., Sabioni R. E., Volkoff H., Cyrino J. E. P. (2022). Growth performance, health, and gene expression of appetite-regulating hormones in Dourado Salminus brasiliensis, fed vegetable-based diets supplemented with swine liver hydrolysate. Aquaculture 548, 737640. doi:10.1016/j.aquaculture.2021.737640
Love M. I., Huber W., Anders S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi:10.1186/s13059-014-0550-8
Mazumder A., Dwivedi A., Du Plessis J. (2016). Sinigrin and its therapeutic benefits. Molecules 21 (4), 416. doi:10.3390/molecules21040416
McCann K. M., Rawles S. D., Lochmann R. T., McEntire M. E., Sealey W. M., Gaylord T. G., et al. (2021). Dietary replacement of fishmeal with commercial protein blends designed for aquafeeds in hybrid striped bass (Morone chrysops♂× Morone saxatilis♀): digestibility, growth, body composition, and nutrient retention. Aquac. Rep. 21, 100903. doi:10.1016/j.aqrep.2021.100903
Meiler K. A., Kumar V. (2021). Organic and inorganic zinc in the diet of a commercial strain of diploid and triploid rainbow trout (Oncorhynchus mykiss): effects on performance and mineral retention. Aquaculture 545, 737126. doi:10.1016/j.aquaculture.2021.737126
Metón I., Caseras A., Fernández F., Baanante I. V. (2000). 6-Phosphofructo-2-kinase/fructose-2, 6-bisphosphatase gene expression is regulated by diet composition and ration size in liver of gilthead sea bream, Sparus aurata. Biochimica Biophysica Acta (BBA)-Gene Struct. Expr. 1491 (1–3), 220–228. doi:10.1016/S0167-4781(00)00040-3
Moen S. H., Ehrnström B., Kojen J. F., Yurchenko M., Beckwith K. S., Afset J. E., et al. (2019). Human Toll-like receptor 8 (TLR8) is an important sensor of pyogenic bacteria, and is attenuated by cell surface TLR signaling. Front. Immunol. 10, 1209. doi:10.3389/fimmu.2019.01209
Nengas I., Alexis M. N., Davies S. J. (1999). High inclusion levels of poultry meals and related byproducts in diets for gilthead seabream Sparus aurata L. Aquaculture 179 (1–4), 13–23. doi:10.1016/S0044-8486(99)00148-9
Neves J. V., Ramos M. F., Moreira A. C., Silva T., Gomes M. S., Rodrigues P. N. S. (2017). Hamp1 but not Hamp2 regulates ferroportin in fish with two functionally distinct hepcidin types. Sci. Rep. 7 (1), 14793. doi:10.1038/s41598-017-14933-5
Ngo D. T., Wade N. M., Pirozzi I., Glencross B. D. (2016). Effects of canola meal on growth, feed utilisation, plasma biochemistry, histology of digestive organs and hepatic gene expression of barramundi (Asian seabass; Lates calcarifer). Aquaculture 464, 95–105. doi:10.1016/j.aquaculture.2016.06.020
NRC (National Research Council) (2011). Nutrient requirements of fish and shrimp. Washington, DC, United States: National Academies Press.
Oliver P. L., Finelli M. J., Edwards B., Bitoun E., Butts D. L., Becker E. B. E., et al. (2011). Oxr1 is essential for protection against oxidative stress-induced neurodegeneration. PLoS Genet. 7 (10), e1002338. doi:10.1371/journal.pgen.1002338
OmicsBox (2019). Bioinformatics made easy, BioBam bioinformatics. Available at: https://www.biobam.com/omicsbox.
Panserat S., Hortopan G., Plagnes-Juan E., Kolditz C., Lansard M., Skiba-Cassy S., et al. (2009). Differential gene expression after total replacement of dietary fish meal and fish oil by plant products in rainbow trout (Oncorhynchus mykiss) liver. Aquaculture 294 (1–2), 123–131. doi:10.1016/j.aquaculture.2009.05.013
Panserat S., Kolditz C., Richard N., Plagnes-Juan E., Piumi F., Esquerré D., et al. (2008). Hepatic gene expression profiles in juvenile rainbow trout (Oncorhynchus mykiss) fed fishmeal or fish oil-free diets. Br. J. Nutr. 100 (5), 953–967. doi:10.1017/S0007114508981411
Pfaffl M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic acids Res. 29 (9), e45. doi:10.1093/nar/29.9.e45
Portz L., Cyrino J. E. P. (2003). Comparison of the amino acid contents of roe, whole body and muscle tissue and their A/E ratios for largemouth bass Micropterus salmoides (Lacepéde, 1802). Aquac. Res. 34 (8), 585–592. doi:10.1046/j.1365-2109.2003.00834.x
Portz L., Cyrino J. E. P. (2004). Digestibility of nutrients and amino acids of different protein sources in practical diets by largemouth bass Micropterus salmoides (Lacepéde, 1802). Aquac. Res. 35 (4), 312–320. doi:10.1111/j.1365-2109.2004.00984.x
Portz L., Cyrino J. E. P., Martino R. C. (2001). Growth and body composition of juvenile largemouth bass Micropterus salmoides in response to dietary protein and energy levels. Aquac. Nutr. 7 (4), 247–254. doi:10.1046/j.1365-2095.2001.00182.x
Rahman M. M., Li X., Sharifuzzaman S., He M., Poolsawat L., Yang H., et al. (2021). Dietary threonine requirement of juvenile largemouth bass, Micropterus salmoides. Aquaculture 543, 736884. doi:10.1016/j.aquaculture.2021.736884
Rao S., Oyang L., Liang J., Yi P., Han Y., Luo X., et al. (2021). Biological function of HYOU1 in tumors and other diseases. Onco Targets Ther. 14, 1727–1735. doi:10.2147/OTT.S297332
Rawles S. D., Fuller A., Green B. W., Abernathy J. W., Straus D. L., Deshotel M. B., et al. (2022). Growth, body composition, and survival of juvenile white bass (Morone chrysops) when dietary fish meal is partially or totally replaced by soybean meal, poultry by-product meal, an all-plant protein blend or a commercial plant-animal protein blend. Aquac. Rep. 26, 101307. doi:10.1016/j.aqrep.2022.101307
Rios K. E., Liem A., Kashyap A., Maynard S., Paul S., Schaefer B. C. (2018). CARD19: a mitochondrial caspase recruitment domain protein with a role in regulating pro-inflammatory innate immune responses. J. Immunol. 200 (1), 169.18. doi:10.4049/jimmunol.200.Supp.169.18
Roberts A., Pachter L. (2013). Streaming fragment assignment for real-time analysis of sequencing experiments. Nat. methods 10 (1), 71–73. doi:10.1038/nmeth.2251
Schroeter J. C., Peterson B. C., Bledsoe J., Li M., Small B. C. (2018). Targeted gene expression panels and microbiota analysis provide insight into the effects of alternative production diet formulations on channel catfish nutritional physiology. Aquaculture 489, 46–55. doi:10.1016/j.aquaculture.2018.02.004
Shen Y., Zhao Z., Zhao J., Chen X., Cao M., Wu M. (2019). Expression and functional analysis of hepcidin from Mandarin fish (Siniperca chuatsi). Int. J. Mol. Sci. 20, 5602. doi:10.3390/ijms20225602
Shiel B. P., Hall N. E., Cooke I. R., Robinson N. A., Strugnell J. M. (2017). Epipodial tentacle gene expression and predetermined resilience to summer mortality in the commercially important greenlip abalone, Haliotis laevigata. Mar. Biotechnol. 19, 191–205. doi:10.1007/s10126-017-9742-z
Simão F. A., Waterhouse R. M., Ioannidis P., Kriventseva E. V., Zdobnov E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinforma. Oxf. Engl. 31 (19), 3210–3212. doi:10.1093/bioinformatics/btv351
Tacchi L., Secombes C. J., Bickerdike R., Adler M. A., Venegas C., Takle H., et al. (2012). Transcriptomic and physiological responses to fishmeal substitution with plant proteins in formulated feed in farmed Atlantic salmon (Salmo salar). BMC genomics 13 (1), 363–421. doi:10.1186/1471-2164-13-363
Tidwell J. H., Coyle S. D., Bright L. A., Yasharian D. (2005). Evaluation of plant and animal source proteins for replacement of fish meal in practical diets for the largemouth bass Micropterus salmoides. J. World Aquac. Soc. 36 (4), 454–463. doi:10.1111/j.1749-7345.2005.tb00393.x
Tidwell J. H., Webster C. D., Coyle S. D. (1996). Effects of dietary protein level on second year growth and water quality for largemouth bass (Micropterus salmoides) raised in ponds. Aquaculture 145 (1–4), 213–223. doi:10.1016/S0044-8486(96)01356-7
Trzaskalski N. A., Fadzeyeva E., Mulvihill E. E. (2020). Dipeptidyl peptidase-4 at the interface between inflammation and metabolism. Clin. Med. Insights Endocrinol. Diabetes 13, 1179551420912972. doi:10.1177/1179551420912972
Turaihi A. H., Bakker W., van Hinsbergh V. W. M., Serné E. H., Smulders Y. M., Niessen H. W. M., et al. (2018). Insulin receptor substrate 2 controls insulin-mediated vasoreactivity and perivascular adipose tissue function in muscle. Front. physiology 9, 245. doi:10.3389/fphys.2018.00245
Turner C. P., Panter S. S., Sharp F. R. (1999). Anti-oxidants prevent focal rat brain injury as assessed by induction of heat shock proteins (HSP70, HO-1/HSP32, HSP47) following subarachnoid injections of lysed blood. Mol. brain Res. 65 (1), 87–102. doi:10.1016/S0169-328X(98)00340-4
Underwood W., Anthony R. (2020) AVMA guidelines for the euthanasia of animals: 2020 edition, Retrieved on March, 2013(30), pp. 2020–2021.
Untergasser A., Cutcutache I., Koressaar T., Ye J., Faircloth B. C., Remm M., et al. (2012). Primer3--new capabilities and interfaces. Nucleic Acids Res. 1 (15), e115. 40. doi:10.1093/nar/gks596
Wang H., Liu M. (2021). Complement C4, infections, and autoimmune diseases. Front. Immunol. 12, 694928. doi:10.3389/fimmu.2021.694928
Wang H.-Y., Ma C. A., Zhao Y., Fan X., Zhou Q., Edmonds P., et al. (2013). Antibody deficiency associated with an inherited autosomal dominant mutation in TWEAK. Proc. Natl. Acad. Sci. 110 (13), 5127–5132. doi:10.1073/pnas.1221211110
Wang Y., Wang F., Ji W. X., Han H., Li P. (2015). Optimizing dietary protein sources for Japanese sea bass (Lateolabrax japonicus) with an emphasis on using poultry by-product meal to substitute fish meal. Aquac. Res. 46 (4), 874–883. doi:10.1111/are.12242
Webster C. D., Thompson K. R., Morgan A. M., Grisby E. J., Gannam A. L. (2000). Use of hempseed meal, poultry by-product meal, and canola meal in practical diets without fish meal for sunshine bass (Morone chrysops× M. saxatilis). Aquaculture 188 (3–4), 299–309. doi:10.1016/S0044-8486(00)00338-0
Webster C. D., Tidwell J. H., Goodgame L. S., Yancey D. H., Mackey L. (1992). Use of soybean meal and distillers grains with solubles as partial or total replacement of fish meal in diets for channel catfish, Ictalurus punctatus. Aquaculture 106 (3–4), 301–309. doi:10.1016/0044-8486(92)90262-J
Webster C. D., Tiu L. G., Morgan A. M., Gannam A. (1999). Effect of partial and total replacement of fish meal on growth and body composition of sunshine bass Morone chrysops× M. saxatilis fed practical diets. J. World Aquac. Soc. 30 (4), 443–453. doi:10.1111/j.1749-7345.1999.tb00992.x
Wu C., Okar D. A., Newgard C. B., Lange A. J. (2001). Overexpression of 6-phosphofructo-2-kinase/fructose-2, 6-bisphosphatase in mouse liver lowers blood glucose by suppressing hepatic glucose production. J. Clin. investigation 107 (1), 91–98. doi:10.1172/JCI11103
Wu Y. B., Ren X., Chai X., Li P., Wang Y. (2018). Replacing fish meal with a blend of poultry by-product meal and feather meal in diets for giant croaker (Nibea japonica). Aquac. Nutr. 24 (3), 1085–1091. doi:10.1111/anu.12647
Xia W.-L., Kang L. H., Liu C. B., Kang C. J. (2017). Death associated protein 1 (DAP 1) positively regulates virus replication and apoptosis of hemocytes in shrimp Marsupenaeus japonicus. Fish Shellfish Immunol. 63, 304–313. doi:10.1016/j.fsi.2017.02.014
Xu H., Luo J., Ma G., Zhang X., Yao D., Li M., et al. (2018). Acyl-CoA synthetase short-chain family member 2 (ACSS2) is regulated by SREBP-1 and plays a role in fatty acid synthesis in caprine mammary epithelial cells. J. Cell. physiology 233 (2), 1005–1016. doi:10.1002/jcp.25954
Zhang J., Frerman F. E., Kim J.-J. P. (2006). Structure of electron transfer flavoprotein-ubiquinone oxidoreductase and electron transfer to the mitochondrial ubiquinone pool. Proc. Natl. Acad. Sci. 103 (44), 16212–16217. doi:10.1073/pnas.0604567103
Zhang X.-D., Zhang J.-W., Wang H.-Z., Lin B., Chen L., Li G., et al. (2019). Evaluation of soybean meal as alternative to fish meal in diet for juvenile Asian red-tailed catfish (Hemibagrus wyckioides). Aquacult Nutr. 25, 1036–1049. doi:10.1111/anu.12921
Keywords: alternative diets, carnivorous fish, fish meal, plant meal, transcriptomics, RNAseq, white bass, temperate basses
Citation: Fuller SA, Abernathy JW, Sankappa NM, Beck BH, Rawles SD, Green BW, Rosentrater KA, McEntire ME, Huskey G and Webster CD (2024) Hepatic transcriptome analyses of juvenile white bass (Morone chrysops) when fed diets where fish meal is partially or totally replaced by alternative protein sources. Front. Physiol. 14:1308690. doi: 10.3389/fphys.2023.1308690
Received: 06 October 2023; Accepted: 26 December 2023;
Published: 15 January 2024.
Edited by:
Vikash Kumar, Central Inland Fisheries Research Institute (ICAR), IndiaReviewed by:
Basanta Kumar Das, Central Inland Fisheries Research Institute (ICAR), IndiaFarhana Hoque, Central Institute of Freshwater Aquaculture (ICAR), India
Copyright © 2024 Fuller, Abernathy, Sankappa, Beck, Rawles, Green, Rosentrater, McEntire, Huskey and Webster. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: S. Adam Fuller, YWRhbS5mdWxsZXJAdXNkYS5nb3Y=; Jason W. Abernathy, amFzb24uYWJlcm5hdGh5QHVzZGEuZ292; Carl D. Webster, Y2FybC53ZWJzdGVyQHVzZGEuZ292
†These authors share first authorship