# Evaluation of Cottonseed Meal as an Alternative to Fish Meal in Diet for Juvenile Asian Red-Tailed CatfishHemibagrus wyckioides

Authors: Biwei Li; Linhai Su; Yan Sun; Han Huang; Junming Deng; Zhiyong Cao
Journal: Aquaculture Nutrition
 (2023)

---

## Abstract

A 10-week trial was performed to investigate the effects of replacing fishmeal with cottonseed meal (CSM) on the growth rate, protein metabolism, and antioxidant response of Asian red-tailed catfishHemibagrus wyckioides. Five isonitrogenous and isocaloric diets (C0, C8.5, C17.2, C25.7, and C34.4) were prepared to contain 0%, 8.5%, 17.2%, 25.7%, and 34.4% CSM replacing fishmeal, respectively. The weight gain, daily growth coefficient, pepsin, and intestinal amylase activities initially increased and then decreased with the raising dietary CSM levels; the highest values were observed in the C17.2 group (P<0.05). However, feed cost exhibited the opposite trend. With the increasing dietary CSM levels, the protein efficiency ratio and intestinal trypsin activity decreased but feed conversion rate increased gradually; while no differences were observed among the C0, C8.5, and C17.2 groups (P>0.05). Dietary CSM inclusion regardless of levels increased the plasma growth hormone level as well as hepatic aspartate aminotransferase (AST) and γ-glutamyl transpeptidase activities but decreased the plasma glutamate dehydrogenase and AST activities (P<0.05). With the increasing dietary CSM levels, the plasma alkaline phosphatase (AKP) and hepatic superoxide dismutase activities decreased but malondialdehyde content increased gradually, while no differences were observed among the C0, C8.5, and C17.2 groups (P>0.05). The plasma immunoglobulin M content and hepatic glutathione reductase activity initially increased but then decreased with the raising dietary CSM levels; the highest values were found in the C17.2 group. These results indicated that dietary CSM inclusion level up to 17.2% improved the growth rate, feed cost, digestive enzyme activity, and protein metabolism without compromising antioxidant capacity of H. wyckioide, whereas these parameters were depressed by further inclusion of CSM. CSM is a potentially cost-effective alternative plant protein source in diet of H. wyckioide.


---

## Body

## 1. Introduction

As the most important protein source in aquafeed, fishmeal (FM) possesses not only balance amino acids and active substances (e.g. taurine, hydroxyproline, cholesterol) but also incomparable palatability [1]. However, FM can no longer meet the needs of aquaculture because of decreasing fishery resources and increasing FM demand; and which become an enormous barrier for aquaculture expansion. Thus, the pursuit for alternatives to FM has always been hotspot of aquaculture [2]. Among alternative protein sources for FM, plant protein sources appear to be the most appropriate because of their availability, eco-friendly, sustainability, and price [3]. Cottonseed meal (CSM) is a by-product obtained from the process of extracting oil from cottonseed. According to the China Statistical Yearbook 2021, the planting area of cottonseed in China reached 3.03 million hectares, and the output reached 10.32 million tons. Behind soybean meal and rapeseed meal, CSM ranks third in terms of tonnage among the commercially available plant protein sources. CSM is considered a valuable plant protein source as an FM substitute in aquafeed because of its relatively favorable amino acid profiles [4]. Previous studies showed that CSM could replace 15–50% FM without compromising growth rate of Ussuri catfish Pseudobagrus ussuriensis [5], blunt snout bream Megalobrama amblycephala [6], grass carp Ctenopharyngodon idellus [7], silver sillago Sillago sihama [8], red drum Sciaenops ocellatus [9], and turbot Scophthalmus maximus [10, 11]. However, the further substitution were limited by several factors, including relatively low contents of lysine and methionine as well as high level of free gossypol [12].Asian red-tailed catfishHemibagrus wyckioides, an omnivorous freshwater fish, is widely distributed in the Mekong River Basin [13]. This species is very popular with farmers and consumers for its fast growth rate, large size, extensive adaptability, and superior resistance to disease [14, 15]. In addition, this catfish do not exhibit cannibalism phenomenon unlike other catfish [16]. Although this species is widely cultured in ponds and cages throughout the Lancang-Mekong River Area, there is no special compound feed for H. wyckiodies due to limited information on the requirements of nutrients [13, 17–19] and utilization of feed ingredients [20, 21]. Hence, an assessment of the utilization efficiency of various plant protein sources by H. wyckioides is urgently needed. To date, there is no data on the influence of CSM inclusion in diet for H. wyckioides. The aim of the present study was to assess the effects of replacing FM with CSM on growth performance, protein metabolism, and antioxidant response of H. wyckioides.


## 2. Materials and Methods

### 2.1. Experimental Diets

Using FM (crude protein, 73.6%) and CSM (crude protein, 49.6%; free gossypol, 1.12 mg/g) as the primary protein sources, five isonitrogenous (41% crude protein) and isocaloric (21 kJ/g gross energy) diets (C0, C8.5, C17.2, C25.7, and C34.4) were prepared to replace 0%, 12%, 24%, 36%, and 48% FM protein with 0%, 8.5%, 17.2%, 25.7%, and 34.4% CSM, respectively. DL-methionine, L-lysine, and L-leucine were also supplemented to each diet to provide equal amounts of methionine, lysine, and leucine in all experimental diets. The feed ingredients and proximate composition of experimental diets are listed in Table1, and the corresponding amino acid profiles are shown in Table 2.Table 1
Ingredients and proximate composition (% dry matter) of five experimental diets.
Price6Dietary inclusion levels of cottonseed mealUS$/kg0%8.5%17.2%25.7%34.4%IngredientsFish meal11.9052.0045.8039.5033.3027.00Cottonseed meal10.890.008.5017.2025.7034.40Fish oil2.410.901.401.902.402.90Soybean oil1.474.804.754.704.654.60DL-Methionine32.930.000.090.190.280.38L-Lysine31.360.000.170.340.510.68L-leucine9.140.000.080.160.250.33Wheat middling0.4925.5524.4623.2621.1619.96α-Starch0.9313.0011.009.008.006.00Soybean lecithin0.960.500.500.500.500.50Vitamin C22.290.020.020.020.020.02Ca(H2PO4)20.611.201.201.201.201.20NaCl0.090.200.200.200.200.20Ethoxyquin (30%)4.570.030.030.030.030.03Choline chloride (40%)1.040.300.300.300.300.30Vitamin mixture44.290.500.500.500.500.50Mineral mixture51.860.500.500.500.500.50Proximate compositionDry matter90.7988.0989.8589.0488.61Crude protein40.7241.4141.4541.2741.12Crude lipid10.8810.4310.0610.6910.62Ash11.0610.8010.319.749.28Cross energy (KJ/g)20.5920.4020.6120.7120.72Free gossypol (mg/g)−0.100.200.290.39Feed price (US$/kg)1.401.351.311.271.231Supplied by Kunming Tianyuan Feed Co., Ltd. (Yunnan, China); fish meal, 73.6% crude protein, 9.0% crude lipid; cottonseed meal, 49.6% crude protein, 1.3% crude lipid; wheat middling, 13.5% crude protein, 0.86% crude lipid. 2L-Ascorbate-2-polyphosphate (35%), supplied by Galaxy Chemicals Co., Ltd. (Hubei, China). 3Supplied by Shanghai Hanhong Chemical Co., Ltd. (Shanghai, China). 4Vitamin premix (g/kg of mixture): retinyl acetate (2800000 IU/g), 2; cholecalciferol, 0.03; DL-α-tocopheryl acetate, 30; menadione, 3; thiamine hydrochloride, 8; riboflavin, 11; pyridoxine hydrochloride, 8; vitamin B12, 0.02; ascorbic acid, 50; folic acid, 1; biotin, 0.1; niacin, 30; calcium D-pantothenate, 32; and inositol, 25. 5Mineral premix (g/kg of mixture): MgSO4•7H2O, 180; KI, 1; FeSO4•H2O, 260; ZnSO4•H2O, 180; CuSO4• 5H2O, 25; Na2Se2O3, 0.01; MnSO4•H2O, 180; and CoCl2•6H2O, 0.75. 6Price of feed ingredients is calculated with the exchange rate of US$ to RMB at 7.00.Table 2
Amino acid composition of the experiment diets (% dry matter).
Dietary inclusion levels of cottonseed meal0%8.5%17.2%25.7%34.4%Essential amino acidsArginine2.162.422.692.943.21Histidine1.361.321.271.231.18Isoleucine1.401.351.301.241.19Leucine2.542.552.562.542.57Lysine2.862.892.872.832.86Methionine1.011.041.021.051.03Phenylalanine1.521.551.591.621.66Threonine1.401.351.301.251.20Valine1.961.921.881.831.79Tryptophan10.440.430.430.420.41Non-essential amino acidsAlanine2.152.041.921.801.69Aspartic acid3.403.383.363.333.31Cysteine0.360.370.380.390.39Glutamic acid6.426.646.877.037.26Glycine2.392.282.172.061.95Proline1.841.801.761.701.66Serine1.371.371.381.371.38Tyrosine1.061.041.020.990.971Tryptophan content of experimental diets was estimated according to the ingredient amino acid composition. The content of tryptophan in fish meal, cottonseed meal, and wheat middling was 0.76, 0.51, and 0.17% of dry matter, respectively.All feed ingredients except lipid sources were crushed and passed through a sieve with a diameter of 320μm. Soybean lecithin was preblended in the mixture of soybean oil and fish oil. All ingredients were thoroughly mixed; the lipid mixture was added, and then thoroughly mixed again. Suitable water was added to make a dough, and the dough was squeezed into 1.5-mm pellet by a pellet feed maker (KS-180; Jiangsu Jingu Rice Mill Co., Ltd., Zhenjiang, China). The wet pellets were dried at 40°C and then stored at −20°C until use.


### 2.2. Fish and Experimental Procedure

Before the initiation of feeding trial, juvenileH. wyckioides originated from the same batch were fed commercial rainbow trout diet (TR-2242, Chile) for 14 days to acclimate the experimental conditions. After 24-h fast, healthy and uniform juveniles (initial body weight 3.21±0.02 g) were selected and randomly assigned into five dietary groups with triplicates of 35 juveniles per tank (1.0m×0.6m×0.5m). During the feeding trial, the experimental diets were manually fed twice a day (07 : 00, 17 : 00) to apparent satiation for 10 weeks. The cultivation water was recirculated through a filtration system composed of filter sponge, coral stone and active carbon to remove particulate matter and hazardous substances, and water temperature was kept at 27±1°C by a heating rod. All tanks were provided with continuous aeration (dissolved oxygen ≥6 mg/L) and natural photoperiod (13.5 h light/10.5 h dark at the beginning of July and 12 h light/12 h dark at the end of September).


### 2.3. Sample Collection

After the feeding trail, all fish were anesthetized with eugenol (1 : 12000; Shanghai Reagent Corporation, Shanghai, China) after 24-h fast. All fish from each tank were weighed and counted to calculate the growth rate and feed utilization. Five fish per tank were randomly collected and stored at –20°C for body composition analysis. Blood samples were collected from the caudal vein of ten fish per tank using a sterile 1-ml syringe. One part of blood sample was transferred into an Eppendorf tube without anticoagulant, and the other part was collected in a heparinized tube. Serum/plasma samples were allowed to clot for 4 h at 4°C, then centrifuged at 4000 g for 10 min at 4°C, the supernatants were collected and stored at –80°C. Liver, stomach, and midgut (the part from anterior valvula intestine to anterior rectum) were dissected from five fish per tank and stored at –80°C until analysis. Three additional samples of liver and dorsal muscle were collected from each tank and pooled in a 1.5-ml Sterile cryovials (Axygen®), and stored at –80°C for analysis of relative expression of protein metabolism-related genes (mammalian target of rapamycin (mTOR), adenosine monophosphate deaminase 1 (AMPD1), glutamate dehydrogenase (GDH), and insulin-like growth factor 1 (IGF-1).


### 2.4. Chemical Analysis

#### 2.4.1. Proximate Composition

The proximate composition of feed ingredients, experimental diets, and whole-body samples were measured following the method of AOAC [22]. Dry matter was determined by drying the sample to a constant weight at 105°C; crude protein content was examined by the Kjeldahl method (N×6.25); crude lipid content was determined using the Soxhlet method with ether extraction; crude ash content was determined after burning in a muffle furnace at 550°C for 6 h. A bomb calorimeter (Parr 1351; Parr Instrument Co., Moline, IL, USA) was used to measure gross energy. Free gossypol content in diets was determined by high-performance liquid chromatography [23]. The amino acid composition of diets was determined using the method described by Mai et al. [24].


#### 2.4.2. Digestibility Enzyme Activity

To obtain an adequate crude enzyme extract solution, the amount of physiological saline solution (0.9% NaCl) needed to add to the wet stomach/midgut sample was determined by a preliminary study. Approximately 0.2 g wet stomach/midgut plus 1.8 mL cold 0.9% NaCl solution were transferred to a 5-mL centrifugal tube and then homogenized. The homogenate was centrifuged at 6,000 g for 15 min at 4°C (MX-160; Tomy Industry, Tokyo, Japan). The supernatant was transferred to a new test tube for analysis of digestive enzyme activity. Pepsin, trypsin, lipase, and amylase activities were assayed by commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). One unit of pepsin/trypsin activity was defined as 1 mmol of tyrosine equivalent released per min at 37°C. One unit of lipase activity was defined as the amount of enzyme required to produce 1μmol of p-nitrophenol at 37°C for 1 minute. One unit of amylase activity was defined as the amount of enzyme hydrolyzing 0.1 mg of starch in 30 minutes at 37°C. Specific enzyme activity was expressed as enzyme activity per gram protein. The amount of soluble protein in supernatant was determined following the approach of Marion and Bradford [25].


#### 2.4.3. Protein Metabolism-Related Parameters

Blood urea nitrogen (BUN) and albumin contents as well as alanine aminotransferase (ALT) and aspartate aminotransferase (AST) activities in plasma were measured using commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Plasma insulin, growth hormone (GH) and IGF-1 levels, and GDH and AMPD activities were measured by enzyme-linked immunosorbent assay (ELISA) kits (R&D Systems Inc., Minneapolis, MN). Hepatic GDH, AMPD, AST, and ALT activities were determined using the same kits as for plasma and expressed as enzyme activity per gram protein. Hepaticγ-glutamyl transferase (γ-GT) activity was measured by ELISA kit (R&D Systems Inc., Minneapolis, MN) and expressed as enzyme activity per gram protein.


#### 2.4.4. Antioxidant-Related Parameters

Immunoglobulin M (IgM) and malondialdehyde (MDA) contents, alkaline phosphatase (AKP) activity in plasma; and superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), glutathione peroxidase (GSH-Px), glutathione reductase (GR), AKP, and lactate dehydrogenase (LDH) activities in liver were analyzed using the commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Nitric oxide (NO) content in liver was measured using a nitrate reductase method [26].


#### 2.4.5. RNA Extraction, cDNA Synthesis, and Real-Time PCR Analysis

The procedures of RNA isolation, reverse transcription, and quantitative real-time PCR were performed as described by Zhang et al. [20]. Briefly, total tissue RNA was extracted with an RNAiso Plus Kit (Takara Bio, Inc., Dalian, China) and then quantified by a Nanodrop 2000 spectrophotometer (Thermo Scientific, Shanghai, China). Agarose gel (1%) electrophoresis was used to confirm the completeness of RNA sequences and then total RNA was immediately reverse-transcribed to obtain cDNA (Thermo Scientific Revert Aid Firs Strand cDNA Synthesis Kit), which was stored at –20°C until use. The mRNA levels of target genes were assessed using real-time PCR. Primers of target genes were devised by using Primer Premier 5.0 software in accordance with the nucleotide sequences of similar species in GenBank. Gene sequences that could not be acquired in GenBank were cloned, amplified, and sequenced. The new gene sequences were used to design specific real-time PCR primers. All primers were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). The template amount was normalized by using the β-actin gene as an internal reference (Table 3). The 2–ΔΔCt method was used to analyze the relative expression level of the above genes [27].Table 3
Sequence of primers used for real-time quantitative PCR.
NamePrimer sequence (5′−3′)Length(bp)Tm(°C)AMPD1ForwardCCACTATTGACCCACAGTCATACC2456.7ReverseTATGCTTGGATTCATCGTCAACAC24GDHForwardTCAAAATCAACCCCAAAAACTTCT2459.6ReverseATCAGGGGCAGGGACATCAATA22IGF-1ForwardGGGGACCGGGGCTTTTATT1960.1ReverseGTGTGCCGTTGCTCTCGTA19mTORForwardAAGCCGCGTCACATCACACC2059.1ReverseATCAAAGCGCTCCTCCATCAG21β-ActinForwardGGCCGTGACCTGACTGAATACCTC2461.3ReverseAATGCCCATCTCCTGCTCAAAGTC24AMPD1: adenosine monophosphate deaminase 1; GDH: glutamate dehydrogenase; IGF-1: insulin-like growth factors-1; mTOR: mammalian target of rapamycin.


#### 2.4.6. Calculation and Statistical Analysis

The following formulas were used to calculate various parameters and indices:(1)MeanmetabolicbodyweightMBW=Wi/10000.75+Wf/10000.752,Feedintakeg/kgMBWperday=DI/MBWd,WeightgainWG=Wf–WiWi,DailygrowthcoefficientDGC,%/d=Wf1/3–Wi1/3d×100,FeedconversionrateFCR=DIWf–Wi,ProteinefficiencyratioPER=Wf–Wiproteinintake,FeedcostUS$/kgfishgain=FP×FCR.where Wi and Wf are the initial and final body weights (g); d is the feeding days; and DI is the dry diet intake per fish (g DM/fish); and FP (feed price, US$/kg feed) is calculated as feedingredientscost+manufacturecost.All data were analyzed using one-way analysis of variance (ANOVA) in SPSS 17 (SPSS Inc., Chicago, IL, USA) for Windows. If the difference was significant (P<0.05), Tukey’s test was used for multiple comparisons. All results are expressed as mean±standarderror (n=3). Regression analysis was also performed to determine the degree of relationship between various parameters and replacement level of FM by CSM.


## 2.1. Experimental Diets

Using FM (crude protein, 73.6%) and CSM (crude protein, 49.6%; free gossypol, 1.12 mg/g) as the primary protein sources, five isonitrogenous (41% crude protein) and isocaloric (21 kJ/g gross energy) diets (C0, C8.5, C17.2, C25.7, and C34.4) were prepared to replace 0%, 12%, 24%, 36%, and 48% FM protein with 0%, 8.5%, 17.2%, 25.7%, and 34.4% CSM, respectively. DL-methionine, L-lysine, and L-leucine were also supplemented to each diet to provide equal amounts of methionine, lysine, and leucine in all experimental diets. The feed ingredients and proximate composition of experimental diets are listed in Table1, and the corresponding amino acid profiles are shown in Table 2.Table 1
Ingredients and proximate composition (% dry matter) of five experimental diets.
Price6Dietary inclusion levels of cottonseed mealUS$/kg0%8.5%17.2%25.7%34.4%IngredientsFish meal11.9052.0045.8039.5033.3027.00Cottonseed meal10.890.008.5017.2025.7034.40Fish oil2.410.901.401.902.402.90Soybean oil1.474.804.754.704.654.60DL-Methionine32.930.000.090.190.280.38L-Lysine31.360.000.170.340.510.68L-leucine9.140.000.080.160.250.33Wheat middling0.4925.5524.4623.2621.1619.96α-Starch0.9313.0011.009.008.006.00Soybean lecithin0.960.500.500.500.500.50Vitamin C22.290.020.020.020.020.02Ca(H2PO4)20.611.201.201.201.201.20NaCl0.090.200.200.200.200.20Ethoxyquin (30%)4.570.030.030.030.030.03Choline chloride (40%)1.040.300.300.300.300.30Vitamin mixture44.290.500.500.500.500.50Mineral mixture51.860.500.500.500.500.50Proximate compositionDry matter90.7988.0989.8589.0488.61Crude protein40.7241.4141.4541.2741.12Crude lipid10.8810.4310.0610.6910.62Ash11.0610.8010.319.749.28Cross energy (KJ/g)20.5920.4020.6120.7120.72Free gossypol (mg/g)−0.100.200.290.39Feed price (US$/kg)1.401.351.311.271.231Supplied by Kunming Tianyuan Feed Co., Ltd. (Yunnan, China); fish meal, 73.6% crude protein, 9.0% crude lipid; cottonseed meal, 49.6% crude protein, 1.3% crude lipid; wheat middling, 13.5% crude protein, 0.86% crude lipid. 2L-Ascorbate-2-polyphosphate (35%), supplied by Galaxy Chemicals Co., Ltd. (Hubei, China). 3Supplied by Shanghai Hanhong Chemical Co., Ltd. (Shanghai, China). 4Vitamin premix (g/kg of mixture): retinyl acetate (2800000 IU/g), 2; cholecalciferol, 0.03; DL-α-tocopheryl acetate, 30; menadione, 3; thiamine hydrochloride, 8; riboflavin, 11; pyridoxine hydrochloride, 8; vitamin B12, 0.02; ascorbic acid, 50; folic acid, 1; biotin, 0.1; niacin, 30; calcium D-pantothenate, 32; and inositol, 25. 5Mineral premix (g/kg of mixture): MgSO4•7H2O, 180; KI, 1; FeSO4•H2O, 260; ZnSO4•H2O, 180; CuSO4• 5H2O, 25; Na2Se2O3, 0.01; MnSO4•H2O, 180; and CoCl2•6H2O, 0.75. 6Price of feed ingredients is calculated with the exchange rate of US$ to RMB at 7.00.Table 2
Amino acid composition of the experiment diets (% dry matter).
Dietary inclusion levels of cottonseed meal0%8.5%17.2%25.7%34.4%Essential amino acidsArginine2.162.422.692.943.21Histidine1.361.321.271.231.18Isoleucine1.401.351.301.241.19Leucine2.542.552.562.542.57Lysine2.862.892.872.832.86Methionine1.011.041.021.051.03Phenylalanine1.521.551.591.621.66Threonine1.401.351.301.251.20Valine1.961.921.881.831.79Tryptophan10.440.430.430.420.41Non-essential amino acidsAlanine2.152.041.921.801.69Aspartic acid3.403.383.363.333.31Cysteine0.360.370.380.390.39Glutamic acid6.426.646.877.037.26Glycine2.392.282.172.061.95Proline1.841.801.761.701.66Serine1.371.371.381.371.38Tyrosine1.061.041.020.990.971Tryptophan content of experimental diets was estimated according to the ingredient amino acid composition. The content of tryptophan in fish meal, cottonseed meal, and wheat middling was 0.76, 0.51, and 0.17% of dry matter, respectively.All feed ingredients except lipid sources were crushed and passed through a sieve with a diameter of 320μm. Soybean lecithin was preblended in the mixture of soybean oil and fish oil. All ingredients were thoroughly mixed; the lipid mixture was added, and then thoroughly mixed again. Suitable water was added to make a dough, and the dough was squeezed into 1.5-mm pellet by a pellet feed maker (KS-180; Jiangsu Jingu Rice Mill Co., Ltd., Zhenjiang, China). The wet pellets were dried at 40°C and then stored at −20°C until use.


## 2.2. Fish and Experimental Procedure

Before the initiation of feeding trial, juvenileH. wyckioides originated from the same batch were fed commercial rainbow trout diet (TR-2242, Chile) for 14 days to acclimate the experimental conditions. After 24-h fast, healthy and uniform juveniles (initial body weight 3.21±0.02 g) were selected and randomly assigned into five dietary groups with triplicates of 35 juveniles per tank (1.0m×0.6m×0.5m). During the feeding trial, the experimental diets were manually fed twice a day (07 : 00, 17 : 00) to apparent satiation for 10 weeks. The cultivation water was recirculated through a filtration system composed of filter sponge, coral stone and active carbon to remove particulate matter and hazardous substances, and water temperature was kept at 27±1°C by a heating rod. All tanks were provided with continuous aeration (dissolved oxygen ≥6 mg/L) and natural photoperiod (13.5 h light/10.5 h dark at the beginning of July and 12 h light/12 h dark at the end of September).


## 2.3. Sample Collection

After the feeding trail, all fish were anesthetized with eugenol (1 : 12000; Shanghai Reagent Corporation, Shanghai, China) after 24-h fast. All fish from each tank were weighed and counted to calculate the growth rate and feed utilization. Five fish per tank were randomly collected and stored at –20°C for body composition analysis. Blood samples were collected from the caudal vein of ten fish per tank using a sterile 1-ml syringe. One part of blood sample was transferred into an Eppendorf tube without anticoagulant, and the other part was collected in a heparinized tube. Serum/plasma samples were allowed to clot for 4 h at 4°C, then centrifuged at 4000 g for 10 min at 4°C, the supernatants were collected and stored at –80°C. Liver, stomach, and midgut (the part from anterior valvula intestine to anterior rectum) were dissected from five fish per tank and stored at –80°C until analysis. Three additional samples of liver and dorsal muscle were collected from each tank and pooled in a 1.5-ml Sterile cryovials (Axygen®), and stored at –80°C for analysis of relative expression of protein metabolism-related genes (mammalian target of rapamycin (mTOR), adenosine monophosphate deaminase 1 (AMPD1), glutamate dehydrogenase (GDH), and insulin-like growth factor 1 (IGF-1).


## 2.4. Chemical Analysis

### 2.4.1. Proximate Composition

The proximate composition of feed ingredients, experimental diets, and whole-body samples were measured following the method of AOAC [22]. Dry matter was determined by drying the sample to a constant weight at 105°C; crude protein content was examined by the Kjeldahl method (N×6.25); crude lipid content was determined using the Soxhlet method with ether extraction; crude ash content was determined after burning in a muffle furnace at 550°C for 6 h. A bomb calorimeter (Parr 1351; Parr Instrument Co., Moline, IL, USA) was used to measure gross energy. Free gossypol content in diets was determined by high-performance liquid chromatography [23]. The amino acid composition of diets was determined using the method described by Mai et al. [24].


### 2.4.2. Digestibility Enzyme Activity

To obtain an adequate crude enzyme extract solution, the amount of physiological saline solution (0.9% NaCl) needed to add to the wet stomach/midgut sample was determined by a preliminary study. Approximately 0.2 g wet stomach/midgut plus 1.8 mL cold 0.9% NaCl solution were transferred to a 5-mL centrifugal tube and then homogenized. The homogenate was centrifuged at 6,000 g for 15 min at 4°C (MX-160; Tomy Industry, Tokyo, Japan). The supernatant was transferred to a new test tube for analysis of digestive enzyme activity. Pepsin, trypsin, lipase, and amylase activities were assayed by commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). One unit of pepsin/trypsin activity was defined as 1 mmol of tyrosine equivalent released per min at 37°C. One unit of lipase activity was defined as the amount of enzyme required to produce 1μmol of p-nitrophenol at 37°C for 1 minute. One unit of amylase activity was defined as the amount of enzyme hydrolyzing 0.1 mg of starch in 30 minutes at 37°C. Specific enzyme activity was expressed as enzyme activity per gram protein. The amount of soluble protein in supernatant was determined following the approach of Marion and Bradford [25].


### 2.4.3. Protein Metabolism-Related Parameters

Blood urea nitrogen (BUN) and albumin contents as well as alanine aminotransferase (ALT) and aspartate aminotransferase (AST) activities in plasma were measured using commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Plasma insulin, growth hormone (GH) and IGF-1 levels, and GDH and AMPD activities were measured by enzyme-linked immunosorbent assay (ELISA) kits (R&D Systems Inc., Minneapolis, MN). Hepatic GDH, AMPD, AST, and ALT activities were determined using the same kits as for plasma and expressed as enzyme activity per gram protein. Hepaticγ-glutamyl transferase (γ-GT) activity was measured by ELISA kit (R&D Systems Inc., Minneapolis, MN) and expressed as enzyme activity per gram protein.


### 2.4.4. Antioxidant-Related Parameters

Immunoglobulin M (IgM) and malondialdehyde (MDA) contents, alkaline phosphatase (AKP) activity in plasma; and superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), glutathione peroxidase (GSH-Px), glutathione reductase (GR), AKP, and lactate dehydrogenase (LDH) activities in liver were analyzed using the commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Nitric oxide (NO) content in liver was measured using a nitrate reductase method [26].


### 2.4.5. RNA Extraction, cDNA Synthesis, and Real-Time PCR Analysis

The procedures of RNA isolation, reverse transcription, and quantitative real-time PCR were performed as described by Zhang et al. [20]. Briefly, total tissue RNA was extracted with an RNAiso Plus Kit (Takara Bio, Inc., Dalian, China) and then quantified by a Nanodrop 2000 spectrophotometer (Thermo Scientific, Shanghai, China). Agarose gel (1%) electrophoresis was used to confirm the completeness of RNA sequences and then total RNA was immediately reverse-transcribed to obtain cDNA (Thermo Scientific Revert Aid Firs Strand cDNA Synthesis Kit), which was stored at –20°C until use. The mRNA levels of target genes were assessed using real-time PCR. Primers of target genes were devised by using Primer Premier 5.0 software in accordance with the nucleotide sequences of similar species in GenBank. Gene sequences that could not be acquired in GenBank were cloned, amplified, and sequenced. The new gene sequences were used to design specific real-time PCR primers. All primers were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). The template amount was normalized by using the β-actin gene as an internal reference (Table 3). The 2–ΔΔCt method was used to analyze the relative expression level of the above genes [27].Table 3
Sequence of primers used for real-time quantitative PCR.
NamePrimer sequence (5′−3′)Length(bp)Tm(°C)AMPD1ForwardCCACTATTGACCCACAGTCATACC2456.7ReverseTATGCTTGGATTCATCGTCAACAC24GDHForwardTCAAAATCAACCCCAAAAACTTCT2459.6ReverseATCAGGGGCAGGGACATCAATA22IGF-1ForwardGGGGACCGGGGCTTTTATT1960.1ReverseGTGTGCCGTTGCTCTCGTA19mTORForwardAAGCCGCGTCACATCACACC2059.1ReverseATCAAAGCGCTCCTCCATCAG21β-ActinForwardGGCCGTGACCTGACTGAATACCTC2461.3ReverseAATGCCCATCTCCTGCTCAAAGTC24AMPD1: adenosine monophosphate deaminase 1; GDH: glutamate dehydrogenase; IGF-1: insulin-like growth factors-1; mTOR: mammalian target of rapamycin.


### 2.4.6. Calculation and Statistical Analysis

The following formulas were used to calculate various parameters and indices:(1)MeanmetabolicbodyweightMBW=Wi/10000.75+Wf/10000.752,Feedintakeg/kgMBWperday=DI/MBWd,WeightgainWG=Wf–WiWi,DailygrowthcoefficientDGC,%/d=Wf1/3–Wi1/3d×100,FeedconversionrateFCR=DIWf–Wi,ProteinefficiencyratioPER=Wf–Wiproteinintake,FeedcostUS$/kgfishgain=FP×FCR.where Wi and Wf are the initial and final body weights (g); d is the feeding days; and DI is the dry diet intake per fish (g DM/fish); and FP (feed price, US$/kg feed) is calculated as feedingredientscost+manufacturecost.All data were analyzed using one-way analysis of variance (ANOVA) in SPSS 17 (SPSS Inc., Chicago, IL, USA) for Windows. If the difference was significant (P<0.05), Tukey’s test was used for multiple comparisons. All results are expressed as mean±standarderror (n=3). Regression analysis was also performed to determine the degree of relationship between various parameters and replacement level of FM by CSM.


## 2.4.1. Proximate Composition

The proximate composition of feed ingredients, experimental diets, and whole-body samples were measured following the method of AOAC [22]. Dry matter was determined by drying the sample to a constant weight at 105°C; crude protein content was examined by the Kjeldahl method (N×6.25); crude lipid content was determined using the Soxhlet method with ether extraction; crude ash content was determined after burning in a muffle furnace at 550°C for 6 h. A bomb calorimeter (Parr 1351; Parr Instrument Co., Moline, IL, USA) was used to measure gross energy. Free gossypol content in diets was determined by high-performance liquid chromatography [23]. The amino acid composition of diets was determined using the method described by Mai et al. [24].


## 2.4.2. Digestibility Enzyme Activity

To obtain an adequate crude enzyme extract solution, the amount of physiological saline solution (0.9% NaCl) needed to add to the wet stomach/midgut sample was determined by a preliminary study. Approximately 0.2 g wet stomach/midgut plus 1.8 mL cold 0.9% NaCl solution were transferred to a 5-mL centrifugal tube and then homogenized. The homogenate was centrifuged at 6,000 g for 15 min at 4°C (MX-160; Tomy Industry, Tokyo, Japan). The supernatant was transferred to a new test tube for analysis of digestive enzyme activity. Pepsin, trypsin, lipase, and amylase activities were assayed by commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). One unit of pepsin/trypsin activity was defined as 1 mmol of tyrosine equivalent released per min at 37°C. One unit of lipase activity was defined as the amount of enzyme required to produce 1μmol of p-nitrophenol at 37°C for 1 minute. One unit of amylase activity was defined as the amount of enzyme hydrolyzing 0.1 mg of starch in 30 minutes at 37°C. Specific enzyme activity was expressed as enzyme activity per gram protein. The amount of soluble protein in supernatant was determined following the approach of Marion and Bradford [25].


## 2.4.3. Protein Metabolism-Related Parameters

Blood urea nitrogen (BUN) and albumin contents as well as alanine aminotransferase (ALT) and aspartate aminotransferase (AST) activities in plasma were measured using commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Plasma insulin, growth hormone (GH) and IGF-1 levels, and GDH and AMPD activities were measured by enzyme-linked immunosorbent assay (ELISA) kits (R&D Systems Inc., Minneapolis, MN). Hepatic GDH, AMPD, AST, and ALT activities were determined using the same kits as for plasma and expressed as enzyme activity per gram protein. Hepaticγ-glutamyl transferase (γ-GT) activity was measured by ELISA kit (R&D Systems Inc., Minneapolis, MN) and expressed as enzyme activity per gram protein.


## 2.4.4. Antioxidant-Related Parameters

Immunoglobulin M (IgM) and malondialdehyde (MDA) contents, alkaline phosphatase (AKP) activity in plasma; and superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), glutathione peroxidase (GSH-Px), glutathione reductase (GR), AKP, and lactate dehydrogenase (LDH) activities in liver were analyzed using the commercial kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). Nitric oxide (NO) content in liver was measured using a nitrate reductase method [26].


## 2.4.5. RNA Extraction, cDNA Synthesis, and Real-Time PCR Analysis

The procedures of RNA isolation, reverse transcription, and quantitative real-time PCR were performed as described by Zhang et al. [20]. Briefly, total tissue RNA was extracted with an RNAiso Plus Kit (Takara Bio, Inc., Dalian, China) and then quantified by a Nanodrop 2000 spectrophotometer (Thermo Scientific, Shanghai, China). Agarose gel (1%) electrophoresis was used to confirm the completeness of RNA sequences and then total RNA was immediately reverse-transcribed to obtain cDNA (Thermo Scientific Revert Aid Firs Strand cDNA Synthesis Kit), which was stored at –20°C until use. The mRNA levels of target genes were assessed using real-time PCR. Primers of target genes were devised by using Primer Premier 5.0 software in accordance with the nucleotide sequences of similar species in GenBank. Gene sequences that could not be acquired in GenBank were cloned, amplified, and sequenced. The new gene sequences were used to design specific real-time PCR primers. All primers were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China). The template amount was normalized by using the β-actin gene as an internal reference (Table 3). The 2–ΔΔCt method was used to analyze the relative expression level of the above genes [27].Table 3
Sequence of primers used for real-time quantitative PCR.
NamePrimer sequence (5′−3′)Length(bp)Tm(°C)AMPD1ForwardCCACTATTGACCCACAGTCATACC2456.7ReverseTATGCTTGGATTCATCGTCAACAC24GDHForwardTCAAAATCAACCCCAAAAACTTCT2459.6ReverseATCAGGGGCAGGGACATCAATA22IGF-1ForwardGGGGACCGGGGCTTTTATT1960.1ReverseGTGTGCCGTTGCTCTCGTA19mTORForwardAAGCCGCGTCACATCACACC2059.1ReverseATCAAAGCGCTCCTCCATCAG21β-ActinForwardGGCCGTGACCTGACTGAATACCTC2461.3ReverseAATGCCCATCTCCTGCTCAAAGTC24AMPD1: adenosine monophosphate deaminase 1; GDH: glutamate dehydrogenase; IGF-1: insulin-like growth factors-1; mTOR: mammalian target of rapamycin.


## 2.4.6. Calculation and Statistical Analysis

The following formulas were used to calculate various parameters and indices:(1)MeanmetabolicbodyweightMBW=Wi/10000.75+Wf/10000.752,Feedintakeg/kgMBWperday=DI/MBWd,WeightgainWG=Wf–WiWi,DailygrowthcoefficientDGC,%/d=Wf1/3–Wi1/3d×100,FeedconversionrateFCR=DIWf–Wi,ProteinefficiencyratioPER=Wf–Wiproteinintake,FeedcostUS$/kgfishgain=FP×FCR.where Wi and Wf are the initial and final body weights (g); d is the feeding days; and DI is the dry diet intake per fish (g DM/fish); and FP (feed price, US$/kg feed) is calculated as feedingredientscost+manufacturecost.All data were analyzed using one-way analysis of variance (ANOVA) in SPSS 17 (SPSS Inc., Chicago, IL, USA) for Windows. If the difference was significant (P<0.05), Tukey’s test was used for multiple comparisons. All results are expressed as mean±standarderror (n=3). Regression analysis was also performed to determine the degree of relationship between various parameters and replacement level of FM by CSM.


## 3. Results

### 3.1. Growth Performance

After the feeding trial, the survival rate of fish ranged from 96.67% to 98.89%, no statistical difference was found among the dietary groups (P>0.05; Table 4). The final body weight (y=20.240+14.352x–31.630x2, P=0.278), WG (y=5.353+4.512x–9.673x2, P=0.318), and DGC (y=2.240+1.092x–2.397x2, P=0.313) initially increased and then declined with the raising dietary CSM content; the highest growth rate was observed in the C17.2 group. Conversely, the feed cost firstly decreased and then increased slightly with the raising dietary CSM content; the lowest value was found in the C17.2 group (y=1.251–0.382x+0.493x2, P=0.005). The PER (y=2.747–0.414x, P=0.002) decreased but FCR (y=0.885+0.133x, P=0.006) increased with the raising dietary CSM content, but no significant differences were found among the C0, C8.5, and C17.2 groups. However, no statistical difference was found in the feed intake (g/kg MBW per day) among the dietary groups (P>0.05).Table 4
Growth performance ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueFinal body weight (g)20.64±0.62ab20.25±0.77ab23.23±0.88b20.74±1.24ab19.90±1.03ay=20.240+14.352x–31.630x20.1920.278Weight gain5.50±0.22ab5.32±0.18ab6.32±0.26b5.57±0.41ab5.29±0.34ay=5.353+4.512x–9.673x20.1740.318DGC (%/day)2.27±0.05ab2.23±0.05ab2.47±0.06b2.28±0.10ab2.21±0.09ay=2.240+1.092x–2.397x20.1760.313Feed intake (g/kg MBW/d)8.21±0.228.21±0.148.81±0.218.57±0.308.66±0.33y=8.244+1.039x0.1720.124FCR0.89±0.01a0.90±0.01a0.90±0.01a0.93±0.02ab0.96±0.02by=0.885+0.133x0.4590.006PER2.75±0.03c2.66±0.04abc2.69±0.02bc2.61±0.05ab2.53±0.06ay=2.747–0.414x0.5430.002Feed cost (US$/kg fish gain)1.25±0.01b1.22±0.02ab1.17±0.01a1.18±0.02ab1.18±0.01aby=1.251–0.382x+0.493x20.5880.005Survival rate (%)98.89±1.1196.67±1.9398.89±1.1196.67±1.9396.67±1.93y=98.445–3.706x0.0600.380Values are presented asmeans±standarderror (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). DGC: daily growth coefficient; MBW: mean metabolic body weight; FCR: feed conversion rate; PER: protein efficiency ratio; feed cost is converted from RMB into US$ with the exchange rate of US$ to RMB at 7.00.


### 3.2. Digestive Enzyme Activity

The activities of pepsin (y=15.430+23.471x–50.761x2, P=0.157) and intestinal amylase (y=2.792+2.198x–6.101x2, P=0.334) initially raised and then decreased with the increasing dietary CSM content (Table 5); the highest activities were observed in the C17.2 group. The intestinal trypsin activity (y=1.309–0.683x, P=0.007) decreased gradually with the raising dietary CSM level, and that was significantly higher in the C0 group compared to the C25.7 and C34.4 groups (P<0.05). However, no statistical difference was found in the intestinal lipase activity among the dietary groups (P>0.05).Table 5
Digestive enzymes activities in the gastrointestinal tract ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueStomachPepsin (U/mg protein)15.67±0.35a16.33 ± 1.90a20.25±1.22b15.67±0.85a15.46±0.39ay=15.430+23.471x–50.761x20.2660.157IntestineTrypsin (U/μg protein)1.38±0.06b1.14±0.03ab1.19±0.12ab0.98±0.11a1.04±0.02ay=1.309–0.683x0.4450.007Lipase (U/g protein)39.46±3.9338.00±3.7833.40±2.1346.29±6.9936.68±3.17y=38.217+2.292x0.0030.853Amylase (U/mg protein)2.70±0.27ab3.07±0.29ab3.25±0.07b2.32±0.22a2.63±0.27aby=2.792+2.198x–6.101x20.1670.334Values are presented as means ± standard error (n =3). Mean values in the same row with different superscripts are significant different (P <0.05).


### 3.3. Protein Metabolism-Related Parameters

The plasma GH content (y=17.758+56.866x–103.737x2, P<0.001) and AMPD activity (y=14.234+57.976x–97.520x2, P=0.019) as well as hepatic AST activity (y=13.362+111.039x–204.613x2, P=0.005) initially increased and then decreased with the increasing dietary CSM content; and these values were lower in the C0 group compared to the other groups (P<0.05, Table 6). The plasma GDH content in the C0 group or hepatic AMPD activity in the C8.5 group were higher than that in the other groups; in contrast, the hepatic γ-GT activity was lower in the C0 group as compared to the other groups (P<0.05). The plasma AST activity (y=29.249–44.236x, P<0.001) gradually declined with the raising dietary CSM content. However, no statistical differences were found in plasma albumin, BUN, insulin, IGF-1 contents, and ALT activity as well as hepatic GDH and ALT activities among the dietary groups (P>0.05).Table 6
Protein metabolism-related parameters in plasma and liver ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valuePlasmaAlbumin (g/L)15.34±2.1916.09±1.0813.47±1.5011.97±0.9912.72±1.14y=15.791–7.797x0.2640.050BUN (mmol/L)9.86±1.0110.96±1.2613.56±1.5410.27±0.6312.51±1.07y=10.511+3.844x0.0930.268Insulin (mU/L)8.11±0.458.25±1.207.80±0.408.52±0.818.35±0.59y=8.057+0.628x0.0100.728GH (μg/L)17.60±0.31a23.61±0.43c24.81±1.30c25.09±0.83c21.10±0.56by=17.758+56.866x–103.737x20.851<0.001IGF-1 (μg/L)19.90±1.2822.33±0.2321.08±2.4920.19±1.1418.89±0.92y=20.294+13.438x–35.202x20.1990.264GDH (U/L)9.63±0.42b6.04±0.43a5.95±0.61a7.06±0.33a6.38±0.50ay=9.115–21.297x+34.838x20.5680.007AMPD (U/L)14.05±0.23a19.53±0.98ab24.39±0.32b20.25±2.84ab20.39±3.25aby=14.234+57.976x–97.520x20.4840.019AST (U/L)29.21±0.19c22.77±1.67c20.89±1.46b12.47±0.80b7.82±0.73ay=29.249–44.236x0.935<0.001ALT (U/L)12.26±0.7813.93±0.6210.10±2.1410.77±0.6910.11±1.02y=12.925–6.214x0.2220.076LiverGDH (U/g protein)0.17±0.000.15±0.010.15±0.010.16±0.010.15±0.01y=0.161–0.019x0.0770.318AMPD (U/g protein)0.71±0.02a4.15±0.20b0.99±0.07a0.74±0.06a0.73±0.03ay=1.565+6.682x–19.742x20.2460.183γ-GT (U/g protein)0.17±0.01a0.39±0.04b0.36±0.04b0.36±0.04b0.42±0.03by=0.202+1.093x–1.438x20.5710.006AST (U/g protein)11.00±1.20a28.61±1.12b27.82±1.89b22.49±3.87b21.75±1.49by=13.362+111.039x–204.613x20.5920.005ALT (U/g protein)12.46±2.2616.62±2.1815.07±1.8418.65±3.6614.41±2.26y=12.586+32.723x–57.854x20.1550.363Values are presented asmeans±standard error (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). BUN: blood urea nitrogen; GH: growth hormone; IGF-1: insulin-like growth factors-1; GDH: glutamate dehydrogenase; AMPD: adenosine monophosphate deaminase; AST: aspartate aminotransferase; ALT: alanine aminotransferase; γ-GT: γ-glutamyl transferase.The hepatic IGF-1 expression level was higher in the C0 group compared to the C8.5 and C17.2 groups (P<0.05, Table 7). Conversely, the hepatic mTOR expression level was higher in the C8.5 group compared to the C0, C25.7, and C34.4 groups (P<0.05). The relative expression levels of hepatic AMPD and muscle GDH were higher in the C0 group compared to the other groups (P<0.05). The hepatic GDH expression level was higher in the C0, C17.2, and C34.4 groups compared to the C8.5 and C25.7 groups (P<0.05). The muscle AMPD expression level was lower in the C8.5 group compared to the C17.2 and C34.4 groups (P<0.05).Table 7
Relative mRNA expression of protein metabolism-related enzymes inHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueLiverIGF-11.04±0.29b0.46±0.11a0.14±0.01a0.63±0.07ab0.52±0.04aby=0.981–4.908x+8.705x20.5590.057mTOR1.01±0.10a3.92±0.48b2.73±0.83ab0.93±0.04a1.49±0.05ay=1.813+11.421x–28.286x20.3100.188AMPD11.01±0.12b0.05±0.02a0.08±0.03a0.02±0.01a0.05±0.03ay=0.856–5.769x+8.799x20.792<0.001GDH1.00±0.01b0.42±0.03a0.73±0.06b0.39±0.12a0.76±0.16by=0.912–3.167x+5.761x20.3440.121MuscleAMPD1.02±0.20ab0.41±0.03a1.43±0.42b0.76±0.02ab1.56±0.07by=0.749+1.187x0.1810.220GDH1.00±0.01b0.41±0.01a0.40±0.02a0.29±0.06a0.25±0.01ay=0.933–3.645x+4.828x20.872<0.001Values are presented asmeans±standarderror (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). IGF-1: insulin-like growth factors-1; mTOR: mammalian target of rapamycin; AMPD1, adenosine monophosphate deaminase 1; GDH, glutamate dehydrogenase.


### 3.4. Immune and Antioxidant-Related Indexes

The plasma IgM content (y=0.061+0.743x–1.455x2, P=0.002) and hepatic GR activity (y=8.452+14.149x–37.782x2, P=0.049) initially increased and then decreased with the increasing dietary CSM content; and the highest values were found in the C17.2 group (P<0.05; Table 8). The plasma MDA content (y=4.497+28.814x, P<0.001) increased gradually with the raising dietary CSM content, whereas no statistical difference was observed among the C0, C8.5, C17.2, and C25.7 groups (P>0.05). Conversely, the plasma AKP (y=16.747–31.908x, P<0.001) and hepatic SOD activities (y=33.946–9.580x, P=0.047) declined generally with the raising dietary CSM content. The hepatic NO content was higher in the C0 group compared to the C8.5 group (P<0.05). However, no statistical differences were observed in the hepatic CAT, POD, GSH-Px, AKP, and LDH activities among the dietary groups (P>0.05).Table 8
Antioxidant-related parameters in plasma and liver ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valuePlasmaIgM (g/L)0.05±0.01a0.15±0.01c0.15±0.03c0.13±0.01bc0.09±0.01aby=0.061+0.743x–1.455x20.6570.002MDA (nmol/ml)7.21±0.82a7.52±0.43a8.37±1.05a11.40±0.14a22.56±2.55by=4.497+28.814x0.652<0.001AKP (U/dl)15.18±2.41b17.66±2.56b6.38±1.19a2.74±0.72a3.49±0.74ay=16.747–31.908x0.669<0.001LiverSOD (U/mg protein)34.46±0.73b32.07±2.01ab29.94±1.78ab33.65±2.95ab27.72±1.30ay=33.946–9.580x0.2310.047CAT (U/mg protein)10.46±0.3010.81±0.3510.97±0.2010.16±0.5110.07±0.33y=10.531+2.677x–7.365x20.3030.115POD (U/mg protein)1.17±0.141.08±0.081.13±0.051.08±0.030.97±0.03y=1.167–0.328x0.2160.081GSH-Px (U/μg protein)0.14±0.000.11±0.010.11±0.010.12±0.010.12±0.01y=126.150–23.065x0.1150.216GR (U/g protein)7.99±0.85ab9.21±1.09b12.78±1.28c6.38±0.12ab5.50±0.34ay=8.452+14.149x–37.782x20.3510.049AKP (U/g protein)12.20±1.5715.77±2.8614.86±2.9310.93±1.0611.19±1.55y=14.262–5.136x0.0680.347LDH (U/mg protein)4.07±0.124.41±0.243.93±0.314.39±0.193.64±0.38y=4.088+1.969x–5.142x20.2190.227NO (μmol/g protein)6.64±0.40b4.54±0.42a4.71±0.81ab5.85±0.42ab5.25±0.71aby=6.080–6.609x+10.073x20.1240.453Values are presented asmeans±standard error (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). IgM: immunoglobulin M; MDA: malondialdehyde; AKP: alkaline phosphatase; SOD: superoxide dismutase; CAT: catalase; POD: peroxidase; GSH-Px: glutathione peroxidase; GR: glutathione reductase; LDH: lactate dehydrogenase; NO: nitric oxide.


### 3.5. Whole-Body Composition

The whole-body crude lipid content (y=10.746–2.125x, P=0.004) declined generally with the increasing dietary CSM content; and that was lower in the C34.4 group compared to the C0 and C17.2 groups (P<0.05, Table 9). Conversely, the whole-body moisture content (y=69.971+2.244x, P=0.035) increased gradually with the raising dietary CSM content, and that was higher in the C34.4 group compared to the C0 group (P<0.05). However, no statistical differences in the whole-body crude protein and ash contents were observed among the dietary groups (P>0.05).Table 9
The whole-body composition ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
InitialDietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueMoisture (%)85.4769.76±0.20a70.77±0.11ab70.28±0.44ab70.47±0.49ab71.26±0.32by=69.971+2.244x0.3000.035Crude protein (%)11.5814.90±0.0915.13±0.0414.90±0.0615.22±0.3314.88±0.20y=14.917+1.393x–2.827x20.0560.707Crude lipid (%)1.4710.75±0.06b10.26±0.14ab10.51±0.07b10.09±0.35ab9.57±0.33ay=10.746–2.125x0.4910.004Ash (%)1.723.02±0.033.11±0.063.07±0.063.04±0.113.02±0.01y=3.036+0.452x–1.075x20.0870.580Values are presented asmeans±standard error (n=3). Mean values in the same row with different superscripts are significant different (P<0.05).


## 3.1. Growth Performance

After the feeding trial, the survival rate of fish ranged from 96.67% to 98.89%, no statistical difference was found among the dietary groups (P>0.05; Table 4). The final body weight (y=20.240+14.352x–31.630x2, P=0.278), WG (y=5.353+4.512x–9.673x2, P=0.318), and DGC (y=2.240+1.092x–2.397x2, P=0.313) initially increased and then declined with the raising dietary CSM content; the highest growth rate was observed in the C17.2 group. Conversely, the feed cost firstly decreased and then increased slightly with the raising dietary CSM content; the lowest value was found in the C17.2 group (y=1.251–0.382x+0.493x2, P=0.005). The PER (y=2.747–0.414x, P=0.002) decreased but FCR (y=0.885+0.133x, P=0.006) increased with the raising dietary CSM content, but no significant differences were found among the C0, C8.5, and C17.2 groups. However, no statistical difference was found in the feed intake (g/kg MBW per day) among the dietary groups (P>0.05).Table 4
Growth performance ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueFinal body weight (g)20.64±0.62ab20.25±0.77ab23.23±0.88b20.74±1.24ab19.90±1.03ay=20.240+14.352x–31.630x20.1920.278Weight gain5.50±0.22ab5.32±0.18ab6.32±0.26b5.57±0.41ab5.29±0.34ay=5.353+4.512x–9.673x20.1740.318DGC (%/day)2.27±0.05ab2.23±0.05ab2.47±0.06b2.28±0.10ab2.21±0.09ay=2.240+1.092x–2.397x20.1760.313Feed intake (g/kg MBW/d)8.21±0.228.21±0.148.81±0.218.57±0.308.66±0.33y=8.244+1.039x0.1720.124FCR0.89±0.01a0.90±0.01a0.90±0.01a0.93±0.02ab0.96±0.02by=0.885+0.133x0.4590.006PER2.75±0.03c2.66±0.04abc2.69±0.02bc2.61±0.05ab2.53±0.06ay=2.747–0.414x0.5430.002Feed cost (US$/kg fish gain)1.25±0.01b1.22±0.02ab1.17±0.01a1.18±0.02ab1.18±0.01aby=1.251–0.382x+0.493x20.5880.005Survival rate (%)98.89±1.1196.67±1.9398.89±1.1196.67±1.9396.67±1.93y=98.445–3.706x0.0600.380Values are presented asmeans±standarderror (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). DGC: daily growth coefficient; MBW: mean metabolic body weight; FCR: feed conversion rate; PER: protein efficiency ratio; feed cost is converted from RMB into US$ with the exchange rate of US$ to RMB at 7.00.


## 3.2. Digestive Enzyme Activity

The activities of pepsin (y=15.430+23.471x–50.761x2, P=0.157) and intestinal amylase (y=2.792+2.198x–6.101x2, P=0.334) initially raised and then decreased with the increasing dietary CSM content (Table 5); the highest activities were observed in the C17.2 group. The intestinal trypsin activity (y=1.309–0.683x, P=0.007) decreased gradually with the raising dietary CSM level, and that was significantly higher in the C0 group compared to the C25.7 and C34.4 groups (P<0.05). However, no statistical difference was found in the intestinal lipase activity among the dietary groups (P>0.05).Table 5
Digestive enzymes activities in the gastrointestinal tract ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueStomachPepsin (U/mg protein)15.67±0.35a16.33 ± 1.90a20.25±1.22b15.67±0.85a15.46±0.39ay=15.430+23.471x–50.761x20.2660.157IntestineTrypsin (U/μg protein)1.38±0.06b1.14±0.03ab1.19±0.12ab0.98±0.11a1.04±0.02ay=1.309–0.683x0.4450.007Lipase (U/g protein)39.46±3.9338.00±3.7833.40±2.1346.29±6.9936.68±3.17y=38.217+2.292x0.0030.853Amylase (U/mg protein)2.70±0.27ab3.07±0.29ab3.25±0.07b2.32±0.22a2.63±0.27aby=2.792+2.198x–6.101x20.1670.334Values are presented as means ± standard error (n =3). Mean values in the same row with different superscripts are significant different (P <0.05).


## 3.3. Protein Metabolism-Related Parameters

The plasma GH content (y=17.758+56.866x–103.737x2, P<0.001) and AMPD activity (y=14.234+57.976x–97.520x2, P=0.019) as well as hepatic AST activity (y=13.362+111.039x–204.613x2, P=0.005) initially increased and then decreased with the increasing dietary CSM content; and these values were lower in the C0 group compared to the other groups (P<0.05, Table 6). The plasma GDH content in the C0 group or hepatic AMPD activity in the C8.5 group were higher than that in the other groups; in contrast, the hepatic γ-GT activity was lower in the C0 group as compared to the other groups (P<0.05). The plasma AST activity (y=29.249–44.236x, P<0.001) gradually declined with the raising dietary CSM content. However, no statistical differences were found in plasma albumin, BUN, insulin, IGF-1 contents, and ALT activity as well as hepatic GDH and ALT activities among the dietary groups (P>0.05).Table 6
Protein metabolism-related parameters in plasma and liver ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valuePlasmaAlbumin (g/L)15.34±2.1916.09±1.0813.47±1.5011.97±0.9912.72±1.14y=15.791–7.797x0.2640.050BUN (mmol/L)9.86±1.0110.96±1.2613.56±1.5410.27±0.6312.51±1.07y=10.511+3.844x0.0930.268Insulin (mU/L)8.11±0.458.25±1.207.80±0.408.52±0.818.35±0.59y=8.057+0.628x0.0100.728GH (μg/L)17.60±0.31a23.61±0.43c24.81±1.30c25.09±0.83c21.10±0.56by=17.758+56.866x–103.737x20.851<0.001IGF-1 (μg/L)19.90±1.2822.33±0.2321.08±2.4920.19±1.1418.89±0.92y=20.294+13.438x–35.202x20.1990.264GDH (U/L)9.63±0.42b6.04±0.43a5.95±0.61a7.06±0.33a6.38±0.50ay=9.115–21.297x+34.838x20.5680.007AMPD (U/L)14.05±0.23a19.53±0.98ab24.39±0.32b20.25±2.84ab20.39±3.25aby=14.234+57.976x–97.520x20.4840.019AST (U/L)29.21±0.19c22.77±1.67c20.89±1.46b12.47±0.80b7.82±0.73ay=29.249–44.236x0.935<0.001ALT (U/L)12.26±0.7813.93±0.6210.10±2.1410.77±0.6910.11±1.02y=12.925–6.214x0.2220.076LiverGDH (U/g protein)0.17±0.000.15±0.010.15±0.010.16±0.010.15±0.01y=0.161–0.019x0.0770.318AMPD (U/g protein)0.71±0.02a4.15±0.20b0.99±0.07a0.74±0.06a0.73±0.03ay=1.565+6.682x–19.742x20.2460.183γ-GT (U/g protein)0.17±0.01a0.39±0.04b0.36±0.04b0.36±0.04b0.42±0.03by=0.202+1.093x–1.438x20.5710.006AST (U/g protein)11.00±1.20a28.61±1.12b27.82±1.89b22.49±3.87b21.75±1.49by=13.362+111.039x–204.613x20.5920.005ALT (U/g protein)12.46±2.2616.62±2.1815.07±1.8418.65±3.6614.41±2.26y=12.586+32.723x–57.854x20.1550.363Values are presented asmeans±standard error (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). BUN: blood urea nitrogen; GH: growth hormone; IGF-1: insulin-like growth factors-1; GDH: glutamate dehydrogenase; AMPD: adenosine monophosphate deaminase; AST: aspartate aminotransferase; ALT: alanine aminotransferase; γ-GT: γ-glutamyl transferase.The hepatic IGF-1 expression level was higher in the C0 group compared to the C8.5 and C17.2 groups (P<0.05, Table 7). Conversely, the hepatic mTOR expression level was higher in the C8.5 group compared to the C0, C25.7, and C34.4 groups (P<0.05). The relative expression levels of hepatic AMPD and muscle GDH were higher in the C0 group compared to the other groups (P<0.05). The hepatic GDH expression level was higher in the C0, C17.2, and C34.4 groups compared to the C8.5 and C25.7 groups (P<0.05). The muscle AMPD expression level was lower in the C8.5 group compared to the C17.2 and C34.4 groups (P<0.05).Table 7
Relative mRNA expression of protein metabolism-related enzymes inHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueLiverIGF-11.04±0.29b0.46±0.11a0.14±0.01a0.63±0.07ab0.52±0.04aby=0.981–4.908x+8.705x20.5590.057mTOR1.01±0.10a3.92±0.48b2.73±0.83ab0.93±0.04a1.49±0.05ay=1.813+11.421x–28.286x20.3100.188AMPD11.01±0.12b0.05±0.02a0.08±0.03a0.02±0.01a0.05±0.03ay=0.856–5.769x+8.799x20.792<0.001GDH1.00±0.01b0.42±0.03a0.73±0.06b0.39±0.12a0.76±0.16by=0.912–3.167x+5.761x20.3440.121MuscleAMPD1.02±0.20ab0.41±0.03a1.43±0.42b0.76±0.02ab1.56±0.07by=0.749+1.187x0.1810.220GDH1.00±0.01b0.41±0.01a0.40±0.02a0.29±0.06a0.25±0.01ay=0.933–3.645x+4.828x20.872<0.001Values are presented asmeans±standarderror (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). IGF-1: insulin-like growth factors-1; mTOR: mammalian target of rapamycin; AMPD1, adenosine monophosphate deaminase 1; GDH, glutamate dehydrogenase.


## 3.4. Immune and Antioxidant-Related Indexes

The plasma IgM content (y=0.061+0.743x–1.455x2, P=0.002) and hepatic GR activity (y=8.452+14.149x–37.782x2, P=0.049) initially increased and then decreased with the increasing dietary CSM content; and the highest values were found in the C17.2 group (P<0.05; Table 8). The plasma MDA content (y=4.497+28.814x, P<0.001) increased gradually with the raising dietary CSM content, whereas no statistical difference was observed among the C0, C8.5, C17.2, and C25.7 groups (P>0.05). Conversely, the plasma AKP (y=16.747–31.908x, P<0.001) and hepatic SOD activities (y=33.946–9.580x, P=0.047) declined generally with the raising dietary CSM content. The hepatic NO content was higher in the C0 group compared to the C8.5 group (P<0.05). However, no statistical differences were observed in the hepatic CAT, POD, GSH-Px, AKP, and LDH activities among the dietary groups (P>0.05).Table 8
Antioxidant-related parameters in plasma and liver ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
Dietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valuePlasmaIgM (g/L)0.05±0.01a0.15±0.01c0.15±0.03c0.13±0.01bc0.09±0.01aby=0.061+0.743x–1.455x20.6570.002MDA (nmol/ml)7.21±0.82a7.52±0.43a8.37±1.05a11.40±0.14a22.56±2.55by=4.497+28.814x0.652<0.001AKP (U/dl)15.18±2.41b17.66±2.56b6.38±1.19a2.74±0.72a3.49±0.74ay=16.747–31.908x0.669<0.001LiverSOD (U/mg protein)34.46±0.73b32.07±2.01ab29.94±1.78ab33.65±2.95ab27.72±1.30ay=33.946–9.580x0.2310.047CAT (U/mg protein)10.46±0.3010.81±0.3510.97±0.2010.16±0.5110.07±0.33y=10.531+2.677x–7.365x20.3030.115POD (U/mg protein)1.17±0.141.08±0.081.13±0.051.08±0.030.97±0.03y=1.167–0.328x0.2160.081GSH-Px (U/μg protein)0.14±0.000.11±0.010.11±0.010.12±0.010.12±0.01y=126.150–23.065x0.1150.216GR (U/g protein)7.99±0.85ab9.21±1.09b12.78±1.28c6.38±0.12ab5.50±0.34ay=8.452+14.149x–37.782x20.3510.049AKP (U/g protein)12.20±1.5715.77±2.8614.86±2.9310.93±1.0611.19±1.55y=14.262–5.136x0.0680.347LDH (U/mg protein)4.07±0.124.41±0.243.93±0.314.39±0.193.64±0.38y=4.088+1.969x–5.142x20.2190.227NO (μmol/g protein)6.64±0.40b4.54±0.42a4.71±0.81ab5.85±0.42ab5.25±0.71aby=6.080–6.609x+10.073x20.1240.453Values are presented asmeans±standard error (n=3). Mean values in the same row with different superscripts are significant different (P<0.05). IgM: immunoglobulin M; MDA: malondialdehyde; AKP: alkaline phosphatase; SOD: superoxide dismutase; CAT: catalase; POD: peroxidase; GSH-Px: glutathione peroxidase; GR: glutathione reductase; LDH: lactate dehydrogenase; NO: nitric oxide.


## 3.5. Whole-Body Composition

The whole-body crude lipid content (y=10.746–2.125x, P=0.004) declined generally with the increasing dietary CSM content; and that was lower in the C34.4 group compared to the C0 and C17.2 groups (P<0.05, Table 9). Conversely, the whole-body moisture content (y=69.971+2.244x, P=0.035) increased gradually with the raising dietary CSM content, and that was higher in the C34.4 group compared to the C0 group (P<0.05). However, no statistical differences in the whole-body crude protein and ash contents were observed among the dietary groups (P>0.05).Table 9
The whole-body composition ofHemibagrus wyckioides fed diets with various level of cottonseed meal.
InitialDietary inclusion levels of cottonseed mealRegression analysis0%8.5%17.2%25.7%34.4%Regression equationR2P valueMoisture (%)85.4769.76±0.20a70.77±0.11ab70.28±0.44ab70.47±0.49ab71.26±0.32by=69.971+2.244x0.3000.035Crude protein (%)11.5814.90±0.0915.13±0.0414.90±0.0615.22±0.3314.88±0.20y=14.917+1.393x–2.827x20.0560.707Crude lipid (%)1.4710.75±0.06b10.26±0.14ab10.51±0.07b10.09±0.35ab9.57±0.33ay=10.746–2.125x0.4910.004Ash (%)1.723.02±0.033.11±0.063.07±0.063.04±0.113.02±0.01y=3.036+0.452x–1.075x20.0870.580Values are presented asmeans±standard error (n=3). Mean values in the same row with different superscripts are significant different (P<0.05).


## 4. Discussion

CSM is a potentially cost-effective alternative plant protein source for fish, but the inclusion level of CSM in aquafeeds was limited due to the potential toxic effect of free gossypol and the deficient in essential amino acids (EAAs) such as lysine, methionine, and leucine [7, 28]. However, CSM is rich in arginine and phenylalanine [7]. In this study, lysine, methionine, and leucine were supplemented to the diets containing CSM to attain concentrations like the control diet (Table 2); arginine concentration in the experimental diets increased from 2.16% (Diet C0) to 3.21% (Diet C34.4) with the increase of dietary CSM level. Except for the above-mentioned amino acids, other EAAs were similar among the experimental diets. Additionally, free gossypol content in the diets raised from 0 g/kg (Diet C0) to 0.39 g/kg (Diet C34.4), the highest value exceeded the tolerance limit of tilapia Tilapia aurea (0.18 g/kg) [29], rainbow trout Oncorhynchus mykiss (0.25–0.29 g/kg) [30, 31], and channel catfish Ictalurus punctatus (0.30 g/kg) [32], but less than that of crucian carp Carassius auratus gibelio♀ × Cyprinus carpio♂ (0.64 g/kg) [33].The present study indicated that the growth rate ofH. wyckioides improved with the increasing dietary CSM content up to 17.2% and then decreased; 25.7% CSM could be incorporated into the diet without causing growth depression, but the inclusion level of 34.4% CSM depressed the growth of H. wyckioides. Similar results were reported in previous studies with grass carp [7], turbot [10], tilapia [34], and southern flounder Paralichthys lethostigma [35], the inclusion of suitable CSM replacing FM markedly improved the growth performance and feed utilization. Furthermore, FM can be completely replaced by CSM in diet of snubnose pompano Trachinotus blochii without adverse effect on growth, metabolism, and general health [12]. However, Wang et al. [9] suggested that dietary CSM inclusion linearly depressed the growth rate of red drum Sciaenops ocellatus. These discordant findings may be explained by the differences in fish species and size, CSM quality, diet composition, and culture conditions. The amount of CSM that can be incorporated in aquafeeds depends mainly on the balance of amino acid and the tolerance of fish to free gossypol [36]. In this study, the supplements of methionine, lysine, and leucine were added as needed to ensure that each diet contained the same concentrations of limiting amino acids. However, dietary arginine concentration gradually increased from 2.16% to 3.21% with the raising inclusion level of CSM. Thus, the growth promoting effect of 17.2% CSM (2.69% arginine) in diet was directly related to the improvement of dietary arginine level in this study. Additionally, the inclusion level of CSM was also determined by the tolerance of H. wyckioides to free gossypol. The present study found no growth depression in fish fed diets containing free gossypol levels ranging from 0.10 to 0.29 mg/g; however, the growth rate was depressed when free gossypol concentrations reached 0.39 mg/g. Thus, the current study demonstrate that 0.29 mg/g free gossypol is the tolerance threshold for H. wyckioides, the value is similar to that for rainbow trout (0.29 g/kg) [31] and channel catfish (0.30 g/kg) [32].Additionally, feed cost is the main indicator reflecting economic benefit of FM replacement. In this study, the feed cost initially decreased and then increased with the raising dietary CSM levels with the lowest value in the C17.2 group, which indicating that feeding diets containing 17.2% CSM apparently reduced the feed cost and FM reliance ofH. wyckioides farming.As is well known, the digestibility and utilization of nutrients mainly depend on digestive enzyme activity [37]. Specifically, the hydrolysis of protein firstly depends on pepsin, and the subsequent digestive process is carried out by trypsin and chymotrypsin in intestine [38]. Lipase can be detected throughout the intestine but is most active in midgut [20]. Omnivorous fish can utilize carbohydrate more efficiently than carnivorous fish because of their higher intestinal amylase activity [39]. In this study, the activities of pepsin and amylase increased with the raising dietary CSM concentrations up to 17.2% but then decreased at higher CSM levels; trypsin activity did not decrease significantly until dietary CSM content reached up 25.7%. These findings indicated that the suitable dietary arginine concentrations (2.69%) induced by dietary CSM inclusion enhanced the digestive capacity of H. wyckioides, but further inclusion of dietary CSM depressed the activities of digestive enzymes. The reduced digestive enzymes activities in fish fed high levels of CSM in this study might be attributable to the high content of free gossypol in high CSM-based diets. Previous studies reported that the free gossypol from CSM may directly inhibit certain enzymes (e.g., pepsinogen, pepsin, and trypsin) by binding the free epsilon amino groups of lysine in the gastrointestinal tract, and thereby depressed the intestinal digestive enzymes activities [7, 40].Protein metabolism is a dynamic process involving the balance between the synthesis and degradation of protein [41]. The mTOR is a serine/threonine kinase, which regulates the protein synthesis [19]. As upstream molecules, insulin/IGF-1 regulate the TOR signaling pathway and participate in a series of physiological processes including cell growth and proliferation [42]. Additionally, GH can activate the TOR signaling pathway, and thereby accelerate protein synthesis [43]. In this study, the plasma GH level and relative expression of hepatic mTOR firstly enhanced and then declined with the increasing dietary CSM content, consistent with the trend of growth performance. These results indicated that the suitable dietary CSM inclusion (8.5%–17.2%) might enhance the protein synthesis and thereby improved the growth performance of H. wyckioides, but those were restrained by the inclusion of 25.7% or more. The promoting effect of low CSM inclusion on protein synthesis may be related to the improvement of amino acids balance (increased arginine content) in diets, while the depressed protein synthesis might be explained as an effect of excessive free gossypol in diets with higher CSM concentrations. A similar observation has been found in blunt snout bream [44], turbot [45], and crab Portunus trituberculatus [4] fed diets with moderate CSM inclusion.It is well known that the activities and relative expression levels of GDH, MAPD,γ-GT, AST, and ALT are correlated with the protein or amino acid catabolism [21]. As an essential enzyme in glutamate metabolism, GDH reversibly catalyzes deamination of glutamate with the production of ammonia [46]. AMPD is an indispensable contributor to ammonia production through its catalysis of the irreversible hydrolysis of AMP to IMP and NH4+ [20]. It is well known that ALT and AST can catabolize amino acids and transfer amino groups to α-keto acid [47]; γ-GT is linked to amino acid transport via the γ-glutamyl cycle [48]. In this study, dietary CSM inclusion regardless of level depressed the plasma GDH and AST activities as well as relative expression of AMPD in liver, GDH in liver and muscle, but increased the hepatic γ-GT and AST activities. The response discrepancies of these parameters in different tissues have been observed in several earlier studies with Ussuri catfish [5] and turbot [10] fed diets with various contents of CSM. These results indicated that dietary CSM inclusions altered the amino acid catabolism in H. wyckioides, which may be related to the synchronous influences of dietary arginine and free gossypol contents. The improvement of amino acids balance (increased arginine content) in diets might inhibit the protein degradation, while free gossypol had certain negative effect on protein metabolism. Free gossypol is easy to combine with lysine, resulting in the decrease of lysine activity [4]. Thus, the high level of arginine present in CSM can promote protein synthesis but inhibit protein degradation as energy source, and thereby improved the growth performance, whereas the high concentration of free gossypol may disturb the protein metabolism of H. wyckioides.As one of the most significant immune components in teleost, IgM is the primary antibody involved in the innate immune response of fish [49]. Additionally, antioxidant enzymes including SOD, CAT, POD, GSH-Px, and GR are essential constituents of the complex immune system and key components of the enzymatic defense mechanism protecting organisms from oxidative damage [20, 21]. MDA is the end product of lipid peroxidation and an important indicator of oxidative stress, which can be used to evaluate free radical activity [21]. The biological production of NO plays a key role in nonspecific host defense; but NO overproduction can also damage animal tissue health by stimulating tumorigenesis [50]. In this study, low inclusion level (≤17.2%) of CSM generally enhanced the plasma IgM level and hepatic GR activity, while higher CSM inclusion (≥25.7%) generally increased the plasma MDA level but decreased the plasma AKP activity as well as hepatic SOD and GR activities. These findings indicated that lower CSM concentrations enhanced the immune capacity of H. wyckioides to some extent, but higher dietary CSM inclusion may depress the antioxidant capacity of H. wyckioides, which may be mainly related to the free gossypol present in CSM. The suppression of free gossypol on antioxidant capability has been verified in other aquatic animals including Ussuri catfish [5], grass carp [7], and Pacific white shrimp Litopenaeus vannamei [51].Previous studies have shown that dietary FM replacement with CSM has different effects on the proximate compositions of different fish species. The whole-body protein [8] or lipid [9, 52] contents of fish showed a decreasing trend with the increasing replacing levels of FM by CSM, while dietary CSM inclusion has no significant effect on the proximate composition of blunt snout bream [6] and grass carp [7]. In this study, the whole-body crude lipid content of H. wyckioides showed a decreasing trend with the increased dietary CSM inclusion level, which might be mainly related to the free gossypol present in CSM. The potential toxicological effects of dietary free gossypol have been shown to impair the protein and lipid deposition of fish [9, 32, 52, 53]. For example, tilapia fed CSM-based diet had a lower whole-body lipid level than fish fed CSM-based diet supplemented with iron for detoxification of gossypol [53]; the whole-body lipid content of channel catfish linearly decreased with the increasing dietary gossypol concentrations [32].


## 5. Conclusion

The suitable CSM inclusion (17.2%) enhanced the amino acids balance (increased arginine content), then improved the digestive enzymes activities and protein metabolism, and thereby promoted the growth performance ofH. wyckioides. However, excess dietary CSM inclusions (≥34.4%) overall decreased the digestive enzyme activity and antioxidant capacity of H. wyckioides, and thereby depressed the growth rate of H. wyckioides, which may be attributed to excess free gossypol in high CSM-based diets. Additionally, appropriate amount of CSM (17.2%) inclusion in diets for H. wyckioide was helpful to reduce the feed cost. CSM is a potentially cost-effective alternative plant protein source in aquafeeds.



---

*Source: 1741724-2023-01-06.xml*