Abstract

Sulfate-based alginate polysaccharide (SAP) is a novelty marine prebiotic. To investigate the beneficial effects of SAP in Pacific white shrimp (Litopenaeus vannamei) fed low-fishmeal diets, six diets (whole fishmeal group: FM; fishmeal replacement with cottonseed protein concentrate and SAP supplementary groups: SAP0, SAP1, SAP2, SAP3, and SAP4) were formulated and fed shrimp for 56 days. The results showed that SAP2 and SAP3 groups reached the best weight gain and specific growth rate (). In serum, the activities of lysozyme and acid phosphatase showed the trend of firstly increased and then decreased (). In the gut, the highest activities of trypsin, lipase, and amylase were found in SAP2 and SAP3 groups (); the histological indexes gradually improved with SAP levels increased (); sequencing of microbiota revealed that the composition and structure of microbiota have been improved, especially in the decreasing abundance of Vibrio, Pseudoalteromonas, and Candidatus Bacilloplasma at genus level. Besides, transcriptomics revealed a degree of similarity in differential gene expression patterns in shrimp; the comparison of RNA-Seq and qRT-PCR verified that SAP improved antioxidant and immunity in shrimp. The challenge revealed that SAP strengthened the resistance against Vibrio parahaemolyticus. Totally, the supplementary SAP to low-fishmeal diets improved growth, intestinal health, and disease resistance in shrimp. Based on the polynomial curve analysis of specific growth rate among SAPs groups, the optimum SAP supplementary level was 1.91%.

1. Introduction

Pacific white shrimp (Litopenaeus vannamei) is a widely distributed omnivorous crustacean. The shrimp performs excellent traits, such as strong adaptability, rapid growth rate, and abundant nutritional value [1]. In aquatic commercial market, the shrimp obtains high economic value, which is well popular with the consumer [2, 3].

In recent years, aquaculture has made great progress, providing high-quality food protein for the public. However, the global fishmeal production resource is in short supply, due to the El Niño phenomenon [4], overfishing, and marine pollution [5]. According to the statistics conducted by FAO (Food and Agriculture Organization) [6], the global demand for fishmeal has increased year by year, yet the fishmeal production cannot meet the aquatic industry development [7]. The plant protein possesses the advantages of wide sources, multiple types, large yields, and low prices, comparing with protein sources obtained from other channels. Plant protein has broad application prospect and promotion value. In brief, utilizing plant protein to replace fishmeal is one of the solutions for sustainable aquaculture [8, 9].

However, the aquatic animals generally possess low tolerance to dietary plant proteins, due to the negative factors such as antinutritional factors, imbalanced essential amino acids, unreasonable fatty acid composition, and poor palatability [10]. This has led to restriction on the application of plant protein in aquafeeds. According to related studies conducted on high-level plant protein feed, Pacific white shrimp would perform reduced growth performance [11], body composition change [12], lowered digestion [13], and antioxidant obstacle [14]. In order to improve the utilization efficiency of plant protein, it is a question worth exploring that relieves the negative impact of fishmeal replacement [15]. Cottonseed protein concentrate (CPC) is a new type of dietary raw material made from cottonseed, which is made through leaching, extraction, and concentration. It has the characteristics of high protein ingredient and low antinutritional content [16]. Cottonseed protein concentrate is considered to be a protein source with application prospects in the aquaculture field. A related study found that CPC can be used as a partial alternative protein source for hybrid grouper (♀Epinephelus fuscoguttatus×♂E. lanceolatus) [17] and juvenile golden pompano (Trachinotus ovatus) [18]. However, Wan et al.[11] conducted a study in L. vannamei and found that high level of CPC replacement fishmeal would slow down the shrimp growth and disturb indexes of hemolymph and hematology.

To eliminate the negative effects of reducing fishmeal, algal polysaccharide is the potentially ameliorative modifier. With antibiotics not being recommended in aquaculture, algal polysaccharide is also an environmentally friendly feed supplementation. The seaweed polysaccharides come from enzymatically hydrolyzed marine algae using special craftsmanship. In recent years, related researchers have devoted themselves to the pharmacological study of seaweed polysaccharides, which has laid theoretical foundation for seaweed polysaccharide application in aquaculture. Moreover, the global field is rich in seaweed resources, which have a broad prospect of exploitation and application [19, 20]. Related studies have shown that algal polysaccharide possesses a variety of physiologically active functions in maintaining organism health. Chotigeat et al. [21] proved that seaweed polysaccharides can strengthen disease resistance in black tiger shrimp (Penaeus monodon). Besides, Supamattaya et al. [22] confirmed that the seaweed extract could improve immune response and health condition in Pacific white shrimp. Another study has found that three different seaweed polysaccharides both can maintain intestinal microbiota homeostasis, promoting probiotic proliferation and production of short-chain fatty acids, yet inhibiting conditional pathogens [23]. Wcl et al. [24] conducted a study in banana shrimp (Fenneropenaeus merguiensis) and found that dietary seaweed (Enteromorpha) polysaccharide improves growth performance involved in regulation of immune responses, intestinal histology, and microbial community. As a prebiotic and potential immunomodulator, the research and application of seaweed polysaccharides in aquatic animals are still in progress [25, 26]. It is worth noting that the sulfate-based alginate polysaccharide (SAP) is a novel seaweed allosteric extract that fermented from brown algae, red algae, and green algae, whose sulfate anion is abundant in molecular structure. Due to the multiple biological activities, the sulfate-based alginate polysaccharide has potentially attempted prospects that needed further study in aquaculture.

In this context, it is an interesting issue to investigate SAP supplementary effects in shrimp fed reducing dietary fishmeal. The present study focused on shrimp growth performance, serum biochemical indexes, intestinal morphology, intestinal microbiota, and gut-related gene expression by virtue of transcriptome analysis. The present study provides deeper insight for dietary supplementation in low fishmeal feed, promoting the environmentally friendly development of aquaculture.

2. Materials and Methods

2.1. Experimental Diets

The experimental sulfate-based alginate polysaccharide (SAP) was provided by Haixingyuan Biological Technology Co., Ltd. (Qingdao, China), which was prepared by specially enzymatic hydrolysis of brown algae, red algae, and green algae. The sugar hydroxyl group of SAP contains abundant sulfuric acid groups, which is the active site of SAP. Figure 1 shows the absolute molecular configuration of SAP.

To prepare diets, fish meal, cottonseed protein concentrate, shrimp shell meal, soybean meal, corn protein flour, peanut meal, and wheat flour were used as protein sources, while fish oil and soy lecithin were used as lipid sources. Firstly, the fishmeal level of the reference group (FM) was set at 25%. Based on that, 30% fishmeal protein was replaced by cottonseed protein concentrate to establish five SAP experimental groups. In the SAP experimental groups, sulfate-based alginate polysaccharide (SAP) supplementation was added at 0%, 1%, 2%, 3%, and 4% levels, respectively. The ingredient compositions and nutrient contents of test diets are shown in Table 1; besides, the essential amino acid (EAA) profiles of test diets are shown in Table 2.

In the process of feed production, various feed materials were crushed and passed through an 80-mesh sieve, accurately weighed according to the feed formulation. The microcomponents were added by the step-by-step expansion principle and then being mixed with the bulk raw materials evenly. They were extruded with a twin-screw granulator (Huagong Electromechanical Technology Co., Ltd., Guangzhou, China) into two kinds of pellet feeds of 1.0 and 1.5 mm. To improve the digestibility and stability of feed, the pellets were cured at 60°C for 30 min, then air-dried for 48 hours, and finally stored at -20°C until use.

2.2. Feeding Management

The culture experiment was approved by the Animal Ethics Committee of Guangdong Ocean University. The culture trial was conducted in 0.3 m3 tanks of indoor culture system at the Marine Biological Research Base of Guangdong Ocean University (Donghai Island, Zhanjiang). The experimental shrimp larvae were obtained from the Southern Base of Marine-culture Seed Project of National 863 Program (Zhanjiang, China) and temporarily reared in outdoor concrete pond for 4 weeks. Before the formal experiment, 720 shrimp larvae of uniform size, healthy physique, and active vitality were randomly selected, whose initial weight was  g. They were randomly divided into 6 groups with 3 replicates, and 40 shrimps were farmed in each tank. The culture experiment lasted for 56 days.

During culture, the feeding amount was calculated as 8%-10% of body weight, which also would be adjusted according to weather, water quality, and shrimp vigor. Besides, the feeding was conducted four times in each day at 07:00, 11:00, 17:00, and 21:00, respectively. The shrimp food intake was observed after feeding for 30 to 60 minutes. The water was changed every other day at the initial stage of experiment and every day at the last stage. The water change was ranging from 1/3 to 1/2 of the tank. Oxygenation was continued during the culture period; thus, the dissolved oxygen concentration was higher than 6.7 mg/L. The water temperature ranged from 28.4°C to 31.2°C. The total alkalinity ranged from 100 to 110 ppm. The salinity fluctuation was 26~28. The pH ranged from 7.8 to 8.2. Besides, the ammonia nitrogen concentration was less than 0.03 mg/L.

At the end of the feeding trial, 30 shrimps were randomly selected from each treatment for the challenge experience against Vibrio parahaemolyticus. The resource of V. parahaemolyticus was obtained from the Key Laboratory of Aquatic Economic Animal Disease Control in Guangdong Higher Education Institution (Zhanjiang, Guangdong). At the beginning of the challenge experiment, 1 mL syringes were used and 200 μL of  V. parahaemolyticus live suspension was injected into the shrimp intramuscularly. The mortality rate was counted every 24 hours postinjection. Subsequently, the cumulative mortality of shrimp after injection to the stabilization period was added up to 7 days.

2.3. Sample Collection

Before sampling, shrimp were fasted for 24 h, then counted and weighed in repetition units. Five shrimps were randomly selected in each tank, blotted with filter paper, and stored at -20°C for routine nutrient analysis of the whole body. In each tank, another 8 shrimps were randomly selected, whose blood was drawn from the pericardial cavity with sterile syringes. The blood was rested at 4°C for 12 hours and centrifuged at 4000 rpm for 10 min (Thermo Fisher, USA). The supernatant was aspirated, and then, the serum was stored at -80°C for later serum index testing. Another 10 shrimps in each replicate were randomly selected and sacrificed, whose intestine tissue was separated on aseptic workbench, washed with 0.65% normal saline, and immediately stored at -80°C for further test of enzyme activity, gene expression, intestinal microbiota, and transcriptome analysis.

Three additional shrimps were selected from each replicate, whose intestines were taken and preserved in 4% formaldehyde solution. The intestinal Alcian blue-periodic acid Schiff sections were made by Sevier Biotechnology Co., Ltd. (Guangzhou, China); the muscle and fold were later observed using a light microscope (Leica, DM6000), then measured using cellSens Standard 1.8 software, and photographed using LAS 3.8 software. Ten viewing fields of each slice were randomly selected for measurement.

2.4. Growth Parameters and Biochemical Analysis

Growth parameters were calculated according to the following formulas:

The approximate-analysis detection of moisture, crude protein, and crude lipid was performed following the standard procedures of the Association of Official Analytical Chemists [27].

2.5. Determination of Biochemical Parameters and Enzyme Activity

After absorbing the surface water with filter paper, the intestinal tissues were accurately weighed. According to the ratio of , normal saline was added to mix with gut tissues. The intestinal tissues were water-bathed at 0°C, shred, and then homogenized and later centrifuged at 2500-3000 rpm for 10 min (Thermo Fisher, USA). Finally, the supernatant was taken and stored at -80°C for enzyme activity assay.

The serum biochemical parameters were determined by commercial kits obtained from Jiancheng Biological Engineering Institute (96T, Nanjing, China), including Malondialdehyde (MDA), superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GSH-Px), and phenol oxidase (PO), which adopted the methods described by Ray et al. [28]. In addition, the following indexes in serum was detected using the methods described by Ayiku et al. [29], including lysozyme (LZM), glutamic pyruvic transaminase (GPT), glutamic oxalacetic transaminase (GOT), alkaline phosphatase (AKP), and acid phosphatase (ACP). Besides, the activity of intestinal trypsin, lipase, and amylase was detected using the 96T Elisa kits obtained from Meilian Biotechnology Co., Ltd. (Shanghai, China). Totally, the experimental operations were carried out strictly according to the kit instructions. Each sample was tested in triplicate.

2.6. Intestinal Microbiome Analysis

The intestinal microbiota (16S amplified region) of L. vannamei was detected using microbiome 16S rDNA sequencing technology, with the assistance of Gene Denovo Biotechnology Co., Ltd. (Guangzhou, China). Samples for microbiome analysis were three replicates per treatment. The soil DNA kit (Magen Co., Ltd., China) was used to extract microbiota DNA. PCR (polymerase chain reaction) was used to amplify the 16S rDNA region of the chromosomal DNA. In addition, the PCR was conducted according to the experimental steps: 95°C for 2 min, after that 27 cycles at 98°C for 10 s, then 62°C for 30 s and 68°C for 30 s, and finally 68°C for 10 min. Besides, the adopted primer sequences were 5-CCTACGGGNGGCWGCAG-3 (341F) and 5-GGACTACHVGGGTATCTAAT-3 (806R), respectively. Amplification products were performed separation using 2% agarose gels and then operated purification using the DNA Gel Extraction Kit obtained from Axygen Biosciences (Silicon Valley, USA). Moreover, the amplification products were quantified performing the ABI SteponePlus real-time fluorescence quantitative PCR system (Life Technologies, USA). Subsequently, the purified amplification products were pooled equimolarly and performed pair-end sequencing () on the Illumina Hiseq2500 platform (San Diego, USA) following the standard protocol. Ultimately, all raw data generated by the present microbiome were deposited in the sequence read archive (SRA) of the National Center for Biotechnology Information (NCBI) under the registration number PRJNA741699.

After performing the noisy sequence filtration and chimerism checking, the clustering was performed using 97% similarity as the threshold (VSEARCH, version 10.00)[30]. Based on the above, the clustered data was divided into identical operational taxonomic units (OTUs)[31]. The number of OTUs was plotted against the total number of sequences drawing sparse curves, to evaluate the sequencing depth. Besides, the intergroup Venn analysis was performed using the microbial indication analysis function of the Omicsmart platform (https://www.omicshare.com/tools/), to identify unique and shared OTUs. The alpha indexes were calculated operating the Mothur software (version v.1.30, http://www.mothur.org/). In addition, the overall bacterial composition was compared using principal coordinate analysis (PCoA) function on the Omicsmart platform (https://www.omicshare.com/tools/). Microbial species annotation was performed using RDP Classifier (version 2.2, http://sourceforge.net/projects/rdpclassifier/), referring to the SILVA database (https://www. arbsilva.de/).

2.7. Intestinal Transcriptome Analysis
2.7.1. RNA Extraction, cDNA Library Construction, and Illumina Sequencing

Total RNA from Pacific white shrimp intestine tissue was extracted using TRIzol Kit obtained from Invitrogen (California, USA), strictly following the instructions. Then, the total RNA was treated with RNase-free DNase I obtained from Tiangen Biochemical Technology Co., Ltd. (Beijing, China). The quality of RNA was tested by using Agilent 2100 Bioanalyzer (Agilent Technologies). Besides, the RNA quantity was measured by a spectrophotometer (NanoDrop Technologies). The RNA integrity was evaluated by using 1.5% agarose-gel electrophoresis, and later examined by absorption spectroscopy. After that, 5 μg high-quality RNA of each group were pooled and sequenced on a HiSeq2000 platform (Illumina, USA), with the operating assistance of Gene Denovo Biotechnology Co., Ltd. (Guangzhou, China). In order to deeply comprehend the profile of gene expression, the entire transcriptome de novo assembly was sequenced using a 100 bp paired-end Solexa strategy, including the first- and second-strand cDNA synthesis. The sequencing data generated in this transcriptome study has been uploaded in the Genome Sequence Archive of the National Genomics Data Center (https://ngdc.cncb.ac.cn/). Besides, the assigned accession of the submission is CRA004649.

2.7.2. De Novo Assembly and Unigene Annotation

To prepare for assembly, raw reads produced by Illumina Hiseq 2000 sequencing platform from seven digital libraries were cleaned by using a software tool in SeqPrep (https://github.com/jstjohn/SeqPrep) and Sickle (https://github.com/najoshi/sickle), according to filter adaptor reads, including the reads (>10% ) which were below 20 bases () and low-quality reads with unclear bases. After that, de novo assembly was performed using short read assembling program, Trinity, whose version was numbertrinityrnaseq-r2014-04-13, as described by Zhao et al. [32]. The transcript clustering into unigenes was conducted by virtue of TIGR gene indices clustering tools in accordance with sequence similarity and length, which adopted the method of Geo et al. [33]. The assembled unigenes were conducted for gene similarity by virtue of BLASTx tool (http://www.ncbi.nlm.nih.gov/BLAST/) with an -value approach of 1−5 to the following databases in National Coalition Building Institute, including nonredundant protein database (http://www.ncbi.nlm.nih.gov), Swiss-Prot protein database (http://www.expasy.ch/sprot), COG/KOG database (http://www.ncbi.nlm.nih.gov/COG), Kyoto Encyclopedia of Genes and Genomes database (http://www.genome.jp/kegg), and GO used Blast2GO, as described by Stefan et al. [34]. And then, protein functional annotations may be received based on the best alignments.

2.7.3. Analysis of Differentially Expressed Genes

As described by Langmead and Salzberg [35], the expression levels of genes were evaluated by reflecting clean reads to unigenes using Bowtie2; after that, they computed for each sample with RSEM using the method described by Dewey and Li [36]. Based on the suggestion conducted by Mortazavi et al. [37], all gene abundance was normalized and calculated using uniquely mapped reads, according to read per kilobase of exon model in per million mapped reads. Based on the method described by Smyth [38], the identification of different gene expressions between treatment groups was conducted by virtue of using R package. In multiple tests, the significance of gene expression was assessed by the value of false discovery rate (FDR) below 0.001. Adopting the following judgment principle, if the FDR was below 0.001, besides a higher than twofold change in expression across libraries (log2 ratio value above 1) was observed, the gene would be considered as significant difference. When the value of log2 fold was more than 1, the transcripts would be considered to upregulate. On the contrary, if the value of log2 fold was below 1, the downregulated transcripts were observed. What is more, the coexpression gene of different gene expressions between treatment groups was deeply understood by KEGG pathway enrichment, which was based on the Kyoto Encyclopedia of Genes and Genomes database (https://www.kegg.jp/kegg/).

2.7.4. Validation of Real-Time Quantitative PCR Analysis

To verify the RNA-Seq data of transcriptome in Pacific white shrimp, the extracted RNA which were used before for transcriptome were detected by (quantitative real-time) qRT-PCR analysis. In this context, nine significantly differentially expressed genes related to antioxidant and innate immunity were opted after comprehensive comparison. The specific primers used in this study are shown in Table 3, which were synthesized in Shenggong Biological Engineering Co., Ltd. (Shanghai, China).

The RT-qPCR reaction was performed using Applied Biosystems 7500 Real-Time PCR System obtained from Life Technologies (California, USA), whose experimental steps adopt SYBR @ Premix Ex Taq™ kit (Takara, Kusatsu). All experimental operations were carried out in strict accordance with the instructions. The thermal cycling conditions were conducted following 95°C for 2 min, after that 40 cycles of 95°C for 10 s, and then 72°C for 20 s. Besides, all reactions were operated in triplicate. According to the test method described by Vieira-Girão et al. [39], β-actin was set as the reference gene in shrimp. The gene expression was calculated using the 2ΔΔCT method, described by Du et al. [40].

2.8. Statistical Analysis

The statistical analysis was conducted using SPSS software 19.0 (IBM, USA). All data are expressed as (SEM). The data normality and variance homogeneity were tested adopting Shapiro-Wilk and Levene’s test, respectively. Two following analysis methods were conducted in the present study: (1) the FM group and SAP0 group were subjected to -test; (2) one-way analysis of variance (ANOVA) was conducted to analyze SAP groups (including SAP0, SAP1, SAP2, SAP3, and SAP4 groups), and then, multiple comparison between groups was tested using the Tukey method. Besides, was considered as the statistical significance.

2.9. Ethics Statement

The experimental procedure has been approved by the animal care and ethics committee of Guangdong Ocean University (GDOU).

3. Results

3.1. Growth Performance

As shown in Table 4, there were no significant differences among groups in initial weight (IW) and SR. The final weight (FW) and FCR of the FM group were significantly different from those of the SAP0 group (). The maximum and minimum of FW appeared in the SAP2 group and SAP4 group, respectively (). The WG and SGR in SAP2 and SAP3 groups were significantly higher than those in other groups (). Besides, the SAP4 group was the minimum of WG and SGR (). On the contrary, the FCR in the SAP4 group was the maximum between groups, while FCR in SAP1, SAP2, and SAP3 groups was significantly lower than that in other groups (). Based on the polynomial curve analysis of SGR among SAP groups, Figure 2 suggests that the optimal level of SAP supplementation for shrimp was 1.91%.

3.2. Whole-Body Proximate Composition

As shown in Table 5, the dry matter of shrimps in the FM group was significantly different from that in the SAP0 group (), but there was no significant difference among the SAP treatments (). The crude protein of the FM group was significantly higher than that of the SAP0 group; and crude protein of SAP0 and SAP1 groups was significantly lower than that of the other groups (); besides, the crude protein of SAP2, SAP3, and SAP4 groups had no difference between groups (). The crude lipid of the FM group was significantly higher than that of the SAP0 group (); the crude lipid content of the SAP group showed a trend of firstly decreasing and then increasing as follows: the SAP0 and SAP4 groups were significantly higher than the other groups, while the SAP2 and SAP3 groups were the least significant values (). The ash content of the FM group was significantly lower than that of the SAP0 group (); yet, there was no significant difference in ash content among SAP groups ().

3.3. Serum Biochemical Parameters

Table 6 shows the serum biochemical parameters in shrimp. Among the 10 parameters that have been investigated, the values in the FM group were significantly different from those in the SAP0 group (). The MDA content of SAP0 and SAP1 groups was significantly higher than that of the SAP2, SAP3, and SAP4 group (). SOD activity showed an upward trend among SAP groups; besides, the SAP0 group was significantly lower than the other groups (). The CAT activity of SAP0 and SAP1 groups was significantly lower than that of SAP2, SAP3, and SAP4 groups (). The GSH-Px activity in the SAP0 group was the minimum, while that in the SAP4 group was the maximum value (). There was no significant difference in GSH-Px activity among SAP1, SAP2, and SAP3 groups (). The PO activity of SAP3 and SAP4 was significantly higher than that of SAP0 and SAP2 (), while SAP0, SAP1, and SAP2 had no significant difference in PO activity (). The LZM activity showed an upward trend among SAP groups, and the AP0 group was the minimum; besides, the SAP4 group has the maximum value (). The activity change of GPT and GOT had similar downward trend as follows: the SAP0 group has the significant maximum values, while the minimum was the SAP3 or SAP4 group (). There was no significant difference in AKP activity among SAP groups (). The ACP activity in SAP2 and SAP3 groups was significantly higher than that in other groups, yet the minimum was in the SAP0 group ().

3.4. Intestine Digestive Enzymes

Table 7 shows the change in digestive enzyme activities of shrimp intestine. The activities of trypsin, lipase, and amylase in the FM group were significantly higher than those in the SAP0 group (). The trypsin activity was significantly lowest in the SAP0 group (), while there was no significant difference among SAP1, SAP2, SAP3, and SAP4 groups (). The lipase activity of the SAP2 group was the significant maximum, while the minimum value appeared in the SAP0 group (). The maximum of amylase activity appeared in the SAP3 group, while the SAP0 and SAP1 groups were significantly lower than in other groups ().

3.5. Intestinal Morphology

Figure 3 shows the intestinal histology of Litopenaeus vannamei. After feeding with SAP supplementation for 8 weeks, the shrimps showed relatively favorable intestinal morphology. As shown in Figure 3(b), the intestinal histology indexes of the FM group were significantly different from those of the SAP0 group, including muscle thickness, fold height, and fold width (). The muscle thickness among SAP treatments showed a trend of increasing firstly and then decreasing; the SAP2 group was the maximum, while the SAP0 group was the minimum significantly (). The villus height of SAP2 and SAP3 was significantly higher than that of SAP0, SAP3, and SAP4 groups (); yet, there was no difference in villus height between SAP0, SAP3, and SAP4 groups (). The villus width of SAP2 and SAP3 groups was significantly higher than that of the SAP0 group (), while there was no significance in villus width among SAP0, SAP1, and SAP4 groups ().

3.6. Intestinal Microbiota Analysis
3.6.1. Richness and Diversity Analysis

Figure 4 shows the alpha-index analysis of shrimp microbiota. The core OTU (operational taxonomic unit) accounted for 181. Clearly, the proportions of shared OTUs within each group were 84.19%, 79.39%, 90.05%, 88.73%, 71.83%, and 74.49%, respectively. This shows that there was a stable core of dominant bacteria in intestinal microbiota. Besides, the microbiota coverage of each group was all close to 99%, which suggests that the sequencing depth of this intestinal microbiota analysis was sufficient. Thus, the experimental results were representative. The ACE and Chao1 indexes of the FM group were significantly higher than those of the SAP0 group (). Moreover, the ACE and Chao1 indexes had similar trend as follows: SAP0 and SAP2 were significantly lower than SAP1, SAP3, and SAP4 (). This suggests that the microbiota richness of SAP0 and SAP2 groups was lower than that of the other groups. On the other hand, the Shannon and Simpson indexes of the FM group were not significantly different from those of the SAP0 group (). In addition, Shannon and Simpson indexes decreased firstly and then increased with the increasing SAP supplementation levels. Besides, the lowest value appeared in the SAP2 group (). The two indexes have inverse proportion relationship with the microbiota diversity. Thus, the microbiota diversity in shrimp gut firstly increased and then decreased as the SAP supplementation increased. And the highest diversity of microbiota occurred in the SAP2 group.

3.6.2. Microbial Community Analysis

Figure 5 reveals the comparison of the overall bacterial composition between treatments. On the one hand, the SAP-containing groups were close to each other (in blue frame), indicating that the bacterial composition was closer inside the SAP-containing groups. Moreover, the layout of the FM and SAP0 groups was clearly separated from the SAP-containing groups, reflecting shrimp fed with or without supplementing SAP that caused structural difference on the bacterial composition of shrimp intestinal microbiota.

At the phylum level of intestinal microbiota (Figure 6(a)), the relative abundance of Verrucomicrobia in the FM group was significantly higher than that in the SAP0 group (). Yet, the relative abundance of Tenericutes and Cyanobacteria in the FM group was significantly lower than that in the SAP0 group (). As the SAP supplements increased, the relative abundances of Bacteroidetes and Verrucomicrobia showed an upward trend; the minimum values of their abundance appeared in the SAP0 group, while the maximums appeared in SAP2 or SAP3 groups (). The relative abundances of Proteobacteria and Tenericutes showed a downward trend; the maximum occurred in the SAP0 group, while the SAP4 or SAP3 group was the significant minimum (). There was no significant difference in the relative abundance of Actinobacteria and Planctomycetes among groups (). The abundance of Firmicutes in the SAP1 group was significantly lower than that of other groups (). The abundance of Cyanobacteria in the SAP0 group was significantly higher than that of other groups, while the abundance in SAP1 and SAP2 groups was significantly lower than that in SAP3 and SAP4 groups ().

As shown in Figure 6(b), the relative abundance of Vibrio, Hoppeia, Haloferula, Motilimonas, and Actibacter in the FM group was significantly higher than that in the SAP0 group (). Moreover, the relative abundance of Halocynthiibacter, Pseudoalteromonas, and Candidatus Bacilloplasma in the FM group was significantly lower than that in the FM group (). The relative abundance of Vibrio and Candidatus Bacilloplasma showed a downward trend between SAP groups; besides, the SAP0 and SAP4 groups were the maximum and minimum, respectively (). The abundance of Hoppeia and Haloferula in SAP0 and SAP2 groups was significantly lower than that in other groups (). The relative abundances of Haloferula, Motilimonas, and Tamlana increased firstly and then decreased with the increased SAP supplement levels; besides, the maximum and minimum were SAP3 and SAP0 groups, respectively (). The relative abundance of Halocynthiibacter in the SAP0 group was significantly higher than that in other groups, while the abundance of Halocynthiibacter in the SAP1 group was the minimum value (). The abundance of Pseudoalteromonas in the SAP0 group was significantly higher than that in other groups, while the abundance in SAP1 and SAP2 groups was lower than that in SAP3 and SAP4 groups (). The abundance of Tenacibaculum firstly decreased and then increased; besides, the relative abundance of SAP1, SAP2, and SAP3 groups was significantly lower than that of SAP0 and SAP4 groups (). The relative abundance of Actibacter in the SAP1 group was significantly higher than that in other SAP treatments ().

3.7. Intestinal Transcriptome Analysis
3.7.1. Illumina Sequencing and De Novo Assembly

This experiment was aimed at testing the intestinal transcriptome of Pacific white shrimp using Illumina HiSeq™ 2000, which is compared with a reference transcriptome detection method. As shown in Table 8, the raw data (bp) produced by each treatment group was above 6 billion. Among them, the proportion of clean reads Q20 was over 97%. Moreover, the proportion of clean reads Q30 was above 92%. After splicing and redundancy removal, all clean reading strips were further assembled, and the range of N50 fluctuated from 47892 bp to 60239 bp. Besides, the GC percentage was greater than 43%. In brief, the results show off the good quality control of the transcriptome data.

3.7.2. Function Annotation, Pathway Assignment, and Differentially Expressed Gene Analysis

In order to express the characteristics of gene products, the GO (gene ontology) enrichment analysis has been divided into three major classification categories (Figure 7A), including biological process, molecular function, and cellular component. Regarding the biological process, each subcategory displayed the downward trend, among which the subcategories with higher values were cellular process, single-organism process, metabolic process, biological regulation, and regulation of biological process. In terms of molecular function, the two more prominent categories were binding and catalytic activity. Moreover, the prominent subcategories with higher values were cell, cell part, organelle, membrane, and organelle part in the cellular component category.

In addition, the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis (Figure 7(b)) was mainly divided into six aspects, namely, metabolism, human diseases, organismal systems, genetic information processing, cellular processes, and environmental information processing. Regarding the metabolism, global and overview maps occupied the highest and most prominent data. Moreover, the infectious diseases and cancers had relatively high values in the classification of human diseases. Additionally, the endocrine system, immune system, digestive system, and nervous system occupied the main front row position in organic systems. On the other hand, the four subdirectories accounted for relatively close proportions in genetic information processing, namely, translation; folding, sorting, and degradation; transcription; and replication and repair. In cellular processes, they were located in numerical order as follows: transport and catabolism, cell growth and death, cellular community eukaryotes, and cell motility. Finally, the signal transduction was the most important approach in the environmental information processing.

The representative comparisons (FM-vs-SAP0 and SAP0-vs-SAP2) were selected for analysis, performing the differentially expressed gene analysis in Figure 8. The results showed that there were more upregulated unigenes in the comparison of FM-vs-SAP0; besides, the extra amount was 334. What is more, the downregulation of unigenes is 406 more than the upregulation of unigenes in the comparison of SAP0-vs-SAP2. Taken together, the Venn diagram analysis found that there were 980 unigenes with similar expression patterns, the proportion of which was 41.58%. At the same time, the heat map analysis suggested that the differential gene expression between FM/SAP0 and SAP2/SAP3 was highly related (green frames in Figure 8(d)). These results recommended that the differential gene expression patterns of L. vannamei fed different diets had a certain degree of similarity, which implemented the conservative strategy.

3.7.3. Validation of Transcriptome by Quantitative Real-Time PCR

In order to identify and verify the transcriptome expression profile, three treatments (including the FM, SAP0, and SAP2) had been selected, to compare the expression of significantly differential genes involved in nonspecific immunity and antioxidant capacity after comprehensive comparison.

Reflecting from a macro point of view, the comparison between RNA-Seq and qRT-PCR showed the degree of similarity, which indicates that the transcriptome data had a certain degree of accuracy and reliability (Figure 9). Specifically, the expression of signal molecules (Toll-1 and Akirin) and antibacterial genes (ALF, Crustin, and LZM) in the SAP2 group was significantly upregulated (), comparing with the SAP0 group. Moreover, the expression levels of genes related to antioxidation (HSP70 and CAT) were significantly upregulated in the FM and SAP2 groups ().

3.8. Challenge Test

The cumulative mortality of challenge test is demonstrated in Figure 10. The cumulative mortality in the FM group was significantly lower than that in the SAP0 group (). The cumulative mortality displayed a trend that decreased firstly and then increased among SAP groups: the SAP0 group was significantly higher than the other SAP groups, while the SAP3 group was the least significant in terms of cumulative mortality (). There was no significant difference between SAP1, SAP2, and SAP4 groups ().

4. Discussion

The Pacific white shrimp is a fast-growing marine organism. During its growing larval stage, amino acids play important roles in the protein synthesis. Therefore, amino acid nutritional intake is a key prerequisite for meeting growth requirements (Hiroki & [41]). In addition, relevant studies found that dietary plant proteins generally contain antinutritional factors and unbalanced amino acid composition, which would cause inhibitory response on the growth of aquatic animals, such as hybrid grouper (Epinephelus fuscoguttatus♀×E. lanceolatus♂) [16] and golden pompano (Trachinotus ovatus) [18]. In the present experiment, the SAP0 treatment was simultaneously lacking supplementation in the context of replacing fishmeal. Correspondingly, shrimp in the SAP0 group displayed worse growth performance. On the other hand, the WG and SGR of shrimp were improved with the increase of SAP supplementation. This suggests that the growth performance of shrimp benefited from SAP supplementation. This may be caused by the fact that fermented seaweed is rich in various mineral elements, organic acids, polysaccharides, and natural antibiotics. These rich active substances help to enhance protein synthesis in shrimp and thus promote growth [42, 43]. This is also verified by the higher crude protein content in the whole-body component of SAP2, SAP3, and SAP4. This result is similar to the conclusion reported by Islam et al. [44], who found that supplementary fermented seaweed could improve organism growth and protein synthesis. On the other hand, the results of this experiment also found that only moderate SAP supplements were beneficial. Excessive seaweed polysaccharide may cause excessive fermentation of intestinal microorganisms and imbalance of intestinal microenvironment. At the same time, the relatively excessive polysaccharide would reduce the viscosity and hygroscopicity of the feed and lower the palatability of the feed [45]. The factors mentioned above are not conducive to the growth of shrimps. This may explain why the growth-promoting effect is not excellent in high-concentration SAP treatments.

Compared to marine bony fishes, the Pacific white shrimp is a relatively simple species in terms of body structure. The shrimp has an open blood system where blood flows throughout the body. Therefore, serum is an important dynamic response indicator of shrimp body functions [28]. In this experiment, we focus on the antioxidant and nonspecific immunity indicators of shrimp serum. In one respect, antioxidant capacity is an important means to remove metabolic waste and maintain internal environment balance [16, 46]. The SAP0 group contained the highest MDA content. Moreover, the enzyme activities related to SOD, CAT, GSH-Px, and PO increased in FM and SAP treatments. Specifically, SOD is an antioxidant metalloenzyme that catalyzes the dismutation of superoxide anion radicals, which plays a vital role in the oxidative and antioxidant balance of the organism [47]. In addition, CAT is a marker enzyme of peroxisome, and its enzymatic activity provides the antioxidant defense mechanism for the organism [48]. What is more, GSH-Px is an important peroxide-decomposing enzyme widely present in the body, which reduces toxic peroxides to nontoxic hydroxyl compounds and promotes the decomposition of H2O2, thus protecting the structure and function of cell membranes from interference and damage by peroxides [15]. Besides, the PO plays an important role in the process of resisting pathogenic microbial infestation, which mainly includes some resistance-related antioxidant enzyme reactions in the phenolic metabolic system [49]. The elevated serum antioxidant capacity of the shrimp organism may be related to dietary SAP supplement. According to relevant study, seaweed extracts are rich in phenolic compounds, the pigments of which are similar to carotenoids [50]. Besides the active substances have the effect of scavenging free radicals and enhancing the antioxidant capacity [51]. Furthermore, the special feature of SAP is closely related to its structural complexity. The molecular structure of SAP is branched, and then, some branching structures can be sulfated. SAP has highly variable functions and structures, which depend on the repetition of oligosaccharide motifs in the polysaccharide, the nature of the oligosaccharide chains, the degree of branching, the molecular weight, and the degree of sulfation [52]. Totally, these highly variable molecular structures confer antioxidant activity to SAP.

As the invertebrates, nonspecific immunity of Pacific white shrimp plays an important role in its maintenance of organism resistance and health. Lysozyme can effectively hydrolyze the peptidoglycan of bacterial cell wall, which can be used as a natural anti-infective substance with bactericidal effect [53]. Moreover, AKP and ACP are major indicators of growth and immune function in aquatic animals, which play important roles in cell regulation and nutrient transport [28]. Thus, AKP, ACP, and lysozyme are closely related in the nonspecific immunity in shrimp organism. The present result showed that SAP0 contained lower activities of LZM, AKP, and ACP, compared to the FM group. Related researches have suggested that dietary fishmeal replacement with plant protein would decrease the nonspecific immunity in aquatic animals, which was verified in hybrid grouper [54] and Pacific white shrimp [55]. From another perspective, the activities of LZM, AKP, and ACP showed an elevated trend in the SAP-treated groups. This suggests that SAP supplementation enhanced the nonspecific immunity activity in shrimp serum. Similarly, the efficacy of seaweed extract analogues in enhancing organism-specific immunity is proved in the results studied on Labeo rohita [56]. On the other hand, the results of this experiment showed that the GPT and GOT activities were significantly highest in the SAP0 group. GPT and GOT play major roles in the hepatopancreas of shrimp. Besides, the GOT and GPT activities in serum are pivotal indexes of the functional stability and organism health [29]. In this study, the GOT and GPT activities increased in the SAP0 group. This indicates to some extent that the fishmeal replacement diet may cause impairment of hepatopancreas function in shrimp. In turn, the activities of GOT and GPT decreased as the levels of SAP supplementation increased. This responded that SAP supplementation improved the hepatopancreas health of shrimp while enhancing nonspecific immunity, again confirming the prebiotic role of seaweed extract analogues [57].

The digestive function relies on the mechanical gastrointestinal motility and chemical action of digestive enzymes. In general, aquatic animals have relatively simple digestive tract and weaker digestive capacity. Besides, the digestive enzymes play crucial roles in food digestion and nutrient absorption through hydrolysis [58]. This experiment revealed that the digestive enzyme activities were significantly lower in the SAP0 group than in the FM group, including the trypsin, lipase, and amylase. Generally, dietary plant protein property would cause decreased activity of digestive enzymes in aquatic animals, which was corroborated by the results in trout (Oncorhynchus mykiss) and sea bream (Sparus aurata) [59]. A related study suggested that the inhibitory effect of plant protein diet on digestive enzyme activity may be related to the presence of antinutritional factors [60]. In the present experiment, the improvement in digestive enzyme activities was obtained with the addition of SAP supplementation. The improved digestive capacity of the shrimp may be related to the role of acidifiers belonging to SAP supplementation. The polysaccharides contained in SAP cannot be digested and absorbed by shrimp but can be broken down by microorganisms in the intestinal tract. Above all, microbial fermentation can hydrolyze polysaccharides and oligosaccharides to produce some beneficial short-chain fatty acids, such as acetic acid and butyric acid. In this context, microbial fermentation causes a decrease of pH in the intestine, which promotes digestive tract enzymatic activation and nutrient digestion [61].

The intestine is the main organ for the digestion and nutrient absorption in shrimp. A normal histology in the gut maximizes the surface area, which is a good guarantee for the active function of the gut [58]. In the present experiment, we focused on the thickness of muscle layer, villus height, and villus width of the shrimp intestine. Relative to the FM group, the intestinal histology indexes were relatively poor in the SAP0 group. In contrast to the SAP0 group, the SAP2 and SAP3 treatments were again improved, especially in terms of villus height and villus width. Microvilli are numerous tiny protrusions in the intestine that are efficient means of increasing the contact area between the digestive tract and food [28]. The improved intestinal morphology facilitates the shrimp to enhance digestion and absorption, which corresponds to the enhanced digestive capability in this experiment. Similarly, it was found that dietary seaweed polysaccharides also alleviated intestinal mucosal inflammation and improved intestinal histology in mice (Mus musculus) [62]. This is because algae oligosaccharide dose activates the mannose receptor signaling pathway by eliminating genetic interference and ultimately promotes intestinal epithelial cell recovery [63].

Intestinal microbiota is the general term for the microbial society in the organism’s gut. In recent years, intestinal microbiota sequencing has received a lot of attention. This is because studies have found that the intestinal microbiota is closely related to the digestion, immunity, metabolism, and health of the host [64, 65]. At the same time, an ecologically balanced intestinal microflora is a positive guarantee of intestinal function and a window into intestinal health [66]. In this experiment, Venn analysis elaborated that there was a core proportion in gut microbial structure, even though shrimp from different treatments received different diets. The core of the bacterial composition may be related to the indigenous bacteria. A related study has found that the microflora in the intestine consists of mobile and indigenous bacteria. The indigenous bacteria mainly colonize the organism gut when the host larvae are incubated. Besides, the composition of indigenous bacteria is influenced by the local environment and is somewhat stable [67, 68]. The present results explained the dynamics of the gut microflora in the context of shrimp receiving different diets. In other words, shrimp fed SAP supplementation obtained higher similarity in the social composition of the microflora. Specifically, the intestinal microbiota is the main point of focus for the probiotic effect of polysaccharide prebiotics [23]. Consistent with the above, polysaccharides and oligosaccharides that cannot be digested by the host but can be fermented by beneficial bacteria into short-chain fatty acids thus improve the intestinal microecology. On the one hand, the proliferation of beneficial bacteria would preferentially occupy the colonial ecological niche of the intestine. In addition, the acidic microecological environment caused by fermentation has an inhibitory effect on conditionally pathogenic bacteria, mainly gram-negative bacteria. Finally, prebiotics themselves can effectively reduce the adhesion protein expression of harmful bacteria, which can also reduce the colonization and further multiplication of harmful bacteria [69]. The present results found that the relative abundance of beneficial bacteria (Bacteroidetes) increased at the phylum level, while the relative abundance of conditional pathogens showed the decreasing trend, such as Proteobacteria and Cyanobacteria. Moreover, the results mainly revealed decreasing relative abundance of harmful bacteria at the genus level, such as Vibrio, Pseudoalteromonas, and Candidatus Bacilloplasma. A study conducted by Zhang et al. [62] found an increase in beneficial microorganisms (Ceratobacteraceae) and a decrease in harmful bacteria (Desulfovibrio) in fucoidan oligosaccharide-dosed mice (Mus musculus). Besides, a Pacific white shrimp study came to a similar conclusion that different types of seaweed extracts both improved the gut microflora, specifically increasing the relative abundance of beneficial bacteria (Bacteroidetes, Firmicutes, and Nontuberculous) and decreasing the relative abundance of harmful bacteria (γ-Proteobacteria and Vibrioceae) [70].

Transcriptome sequencing is the sum of all RNAs that can be transcribed in a given functional state, mainly including mRNAs and noncoding RNAs [71, 72]. The transcriptome is also the most important regulatory modality of organism. In recent years, transcriptomics has been increasingly studied in aquaculture yield, such as Pacific white shrimp [73] and Sparus aurata [74]. In this experiment, transcriptomic sequencing of the shrimp gut was performed, whose experimental data showed good quality control. The transcriptomic analysis revealed that the regulation of shrimp fed different diets was mainly concentrated in the biological process during the gene ontology enrichment analysis. This is because of more unigenes enriched in the biological process, in comparison to other collection. On the other hand, it was found that the pathways affected by SAP were mainly concentrated in the metabolism, diseases, and organismal system collections, based on the results of KEGG pathway enrichment. Notably, signal transduction was found to be enriched to a relatively high number of unigenes, when further subdivided into the next subset. Whether this is suggesting that the regulatory effect of dietary SAP supplementation was based on the signal transduction is an interesting and worthwhile question to explore. A study conducted by Zhao et al. [63] found that marine fucoidan oligosaccharides alleviated intestinal mucosa in mice. The mechanism was exactly based on the role of fucoidan oligosaccharides in signaling to intestinal cells through the mannose receptor signaling pathway. Further, the present experiment investigated significantly differential gene expression and revealed a degree of similarity in the differential gene expression patterns of Pacific white shrimp fed different diets, which performed a conservative metabolic strategy. The reason for this may be that shrimp is relatively simple and stress-prone invertebrates, thus adopting the cautious metabolic strategy which is more conducive to maintain the homeostasis of their internal environment and survival [75].

The results of RNA-Seq and qRT-PCR were compared in this experiment and showed generally consistent trends, indicating a degree of accuracy and reliability of the transcriptome data [76]. Consistent with the above, in order to explore the signaling role of SAP supplementation, this experiment revealed that the gene expression levels of signaling molecules (Toll-1, Akirin, and TRAF6) were upregulated in the SAP2 group, relative to the SAP0 group. The nonspecific immune system is an important line of defense. Toll-1 was found to be actively involved in intrinsic immune signaling [77]. In addition, Akirin is highly conserved evolutionarily and actively involved in regulating the transcription of immune-related target genes in the TLR (Toll-like receptors) pathway downstream [78]. Moreover, TRAF6 is an important component of the TLR superfamily, which can activate the promoters of antimicrobial peptide genes during signaling [79]. After that, the present results suggested a relative upregulation in the expression of antimicrobial peptide-related genes, including ALF, Crustin, and LZM. On this basis, this can be deduced that SAP supplementation may be used to upregulate the expression of antimicrobial peptide genes by regulating nonspecific immune signaling and thus achieving enhanced nonspecific immunity [80]. On the other hand, the FM and SAP2 groups showed increased expression of genes related to antioxidant capacity, including HSP70, SOD, and CAT. HSP70 is a conserved protein used to maintain the homeostasis, which involves in antistress and synergistic immune effects [81]. In other words, dietary fishmeal substitution would decrease the antioxidant capacity of aquatic animals [11]. However, the upregulation of antioxidant-related gene levels in the SAP2 group was found, which was also corroborated by elevated antioxidant enzyme activities in serum. Specifically, SAP supplementation could improve the antioxidant capacity by increasing the antioxidant-related gene expression and enzyme activity in shrimp. This is similar to the suggestion reported on grass carp (Ctenopharyngodon idellus) [57].

The data collected from the challenging experiment postfeeding trial was also a solid argument for the above inference. In this experiment, the cumulative mortality rate of Vibrio parahaemolyticus challenging in the SAP group was all significantly lower than that in the SAP0 group. This indicates that the disease resistance of SAP treatments was significantly improved by the SAP supplementation. In other words, the increased disease resistance of treatments receiving SAP dietary doses may be associated with the upregulated antimicrobial peptide gene expression and the increased antioxidant capacity [21, 82]. Similarly, related experiments on Pacific white shrimp found that macroalgae extractions improved the immune response and gut health under the WSSV (white spot syndrome virus) challenge condition [70].

5. Conclusion

In summary, the results of this experiment demonstrate that SAP supplementation improved growth, serum antioxidant and nonspecific immunity, intestinal health, and disease resistance in Pacific white shrimp fed fishmeal replacement with cottonseed protein concentrate. Moreover, the optimal level of dietary SAP supplementation for Pacific white shrimp could be inferred at 1.91%. In this context, SAP can be considered as the potential improver for shrimp’s low fishmeal feed.

Data Availability

All datasets generated for this study are included in the article/supplementary material. All raw data generated by the present gut microbiome were deposited in the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) under the registration number PRJNA741699. The sequencing data generated in this transcriptome study has been uploaded in the Genome Sequence Archive of the National Genomics Data Center (https://ngdc.cncb.ac.cn/), with the assigned accession number of CRA004649.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

Hongyu Liu, Beiping Tan, Xiaohui Dong, Qihui Yang, and Shuyan Chi were involved in conceptualization; Guofeng Chen, Hongyu Liu, Shuang Zhang, and Min Yao designed the experiments; Baoqiang Liu and Jiajun Chen provided necessary experience assistance; Guofeng Chen was involved in data analysis and writing original draft; Guofeng Chen and Hongyu Liu reviewed and edited the manuscript. All authors approved the final version of the manuscript for submission.

Acknowledgments

This work was economically supported by the National Key R&D Program of China (2019YFD0900200), National Natural Science Foundation of China (No. 31772864), Natural Science Foundation of Guangdong Province (2018A030313154 and 2020A1515011129), and Guangdong Enterprise Science and Technology Commissioner Project (GDKTP2020017800).