Abstract The tiger grouper (Epinephelus fuscoguttatus), an important mariculture fish in Southeast Asia, faces increasing health issues in recent years. Phellodendri Cortex (PC) is a traditional Chinese herbal medicine that exhibits a variety of beneficial effects on tiger groupers. The effects of PC, however, varies with the period of dietary intervention. This study aims to investigate the long-term effects of 1% PC supplementation on tiger groupers, focusing on growth, immunity, disease resistance, and intestinal gene expression. The tiger groupers (with an initial mean weight of 27.5 ± 0.5 g) were fed with a diet of Phellodendri Cortex supplementation and a control diet for 8 weeks. Our results indicate that the long-term PC supplementation did not affect growth or Vibrio disease resistance in tiger groupers. However, the transcriptome analysis revealed potential damage to the structural and functional integrity of the groupers’ intestines. On the other hand, anti-inflammatory and cathepsin inhibition effects were also observed, offering potential benefits to fish enteritis prevention and therapy. Therefore, long-term PC supplementation in grouper culture should be applied with caution. Keywords: transcriptome, Phellodendri Cortex, Epinephelus fuscoguttatus, long-term supplementation 1. Introduction The tiger grouper (Epinephelus fuscoguttatus) is widely farmed in the subtropical sea areas of the Indian Ocean and the Pacific Ocean [[36]1]. However, the increase in the prevalence of diseases associated with bacteria, viruses, and parasites has hampered the development of grouper culture, resulting in mass economic losses in the industry [[37]2,[38]3,[39]4,[40]5]. The current therapeutic agents for the above-mentioned diseases are mainly antibiotics and chemical agents. However, the resulting adverse effects, such as drug resistance and environmental pollution, have made these disease control approaches undesirable [[41]6]. Medicinal plants and their byproducts, which are natural and environmentally friendly, have shown great potential to become alternatives to antibiotics and chemicals, since many of them have been demonstrated to be effective for disease control and prevention [[42]7,[43]8]. Phellodendri Cortex (PC, “Huangbo” in Chinese) is the dried bark of Phellodendron amurense Ruprecht. As one of the most commonly used traditional Chinese herbal medicines, PC has been used for the treatment of diarrhea, hemorrhoids, toxicity, swollen red eyes, mouth sores, and other symptoms in China for thousands of years [[44]9]. Berberine is one of its main components and its active ingredient [[45]10], and recent studies have demonstrated that PC with its berberine has anti-bacterial, anti-inflammatory, and antioxidant effects. It can also have hypoglycemic, anti-tumor, and neuroprotective effects. Because numerous studies have shown that only extremely low amounts of berberine (about 0.36%) can enter the systemic circulation [[46]11], we suspected that the majority of the regulatory effect of PC functioned through the intestine. The intestine is a multifunctional organ in fish. As the largest and most-complex secretory organ in the body of a fish, it is critical for digestion, absorption, and the regulation of the metabolism [[47]12]. The intestinal tract is also an important immune organ, since it provides biological, physical, and chemical barriers that enable fish to fight against bacterial invasion or absorption of toxic substances [[48]13,[49]14]. Intestinal health and normal functionality are critical for fish health. In recent years, to reduce costs, the use of considerable amounts of plant-derived proteins, oil, starch, and gluten in aquafeeds has become the norm in carnivorous fish farming [[50]15,[51]16,[52]17,[53]18]. However, carnivorous fish, such as Atlantic salmon and groupers, are known to develop intestinal damage and inflammation when they consume high levels of terrestrial plant components (such as soybean oils and protein). Fish that are farmed in aquaculture are susceptible to many environmental stresses induced by the poor water quality and high stocking densities that are common in many farming practices. For example, an increase in ammonia nitrogen, nitrites, and nitrate contents can cause irreversible damage to intestinal function [[54]19,[55]20]. The intestine, being a member of the gut–brain axis, exhibits abnormal molecular responses and its functions are jeopardized when exposed to environmental stresses [[56]21,[57]22,[58]23]. With the intensity in development of intensive aquaculture over the last decade, intestinal damage and inflammation, which is a non-infectious disease that leads to serious economic loss, has become a major issue linked to carnivorous fish farming [[59]24]. Fortunately, intestinal health is usually manageable to some degree. Studies have shown that microbial preparations (probiotics, prebiotics, or symbionts) can alleviate intestinal damage and inflammation through manipulation of the intestinal micro-ecosystem [[60]25]. In addition, a number of additives and drugs can also regulate the expression of tight junction-related, metabolism-related, or immune-related genes in the intestine. A previous study in our laboratory revealed that a one-week supplement of 1% PC (in relation to the basal diet) can significant downregulate the gene expressions of IL-8 and IgM in the head kidney of tiger groupers [[61]26]. PC might, therefore, be a promising short-term feed additive to reduce inflammatory signs in tiger groupers. However, the long-term effect of the PC supplement in grouper fish and how it might affect the gene expression of the intestine remain unknown. In this study combining enzyme activity analysis and transcriptome analysis, we explored the long-term (8-weeks) effects of the PC supplement on the antioxidant capacity, non-specific immunity, intestinal gene expression, and disease resistance of tiger groupers. This present research will further our understanding of the long-term effects of PC in tiger groupers and will provide insights into the molecular regulatory mechanism of PC on intestinal genes in teleost fish. 2. Materials and Methods All the experimental procedures of this study were approved by the Hainan University Animal Use and Care Committee (No. HNUAUCC-2020-00011). 2.1. Preparation of PC Extract and Experimental Diets Phellodendri Cortex (PC) was bought from a local, traditional Chinese drugstore in Hainan, China. After being washed with tap water and dried in an oven at 60 °C, the herb was ground into powder and sieved through a 40-mesh sieve (with an aperture of 0.425 mm). Then, 100 g of PC powder was suspended in 1000 mL of deionized water at room temperature for 72 h. Afterwards, the powder suspension was centrifuged at 9000× g for 5 min, and the supernatant was collected. The resulting supernatant (a water extract of PC) was stored at 4 °C for subsequent use. Commercial feed (comprising 51.16% of crude protein, 13.45% of crude lipid, and 7.41% of ash; from Tongwei^®, Tongwei Co., Ltd., Chengdu, China) was used as the basal diet. To prepare the experimental diet, the PC water extract was incorporated into the basal diet through spraying 10 mL of the PC water extract (corresponding to 1 g of PC powder) on every 100 g of the basal diet. The resulting experimental diet had a crude drug–basal diet ratio of 1% (w/w). The experimental diets were packed and stored at 4 °C until use. New batches of the experimental diets were prepared each week to ensure the freshness of the experimental diets. 2.2. Rearing In total, 240 juvenile tiger groupers, Epinephelus fuscoguttatus, were obtained from a local hatchery at Haikou City, Hainan province. These juveniles (with a mean body weight of 27.5 ± 0.5 g) were acclimated to the experimental conditions for two weeks before the feeding experiment began. During the acclimation time, fish were fed the control diets. At the end of the acclimation period, the tiger groupers were separated into two groups with 40 fish per tank in three replicates. The control group and the experimental group were fed with the control diet and the experimental diet, respectively. Fish were fed twice daily (at 10:00 and 16:00), at a feeding rate of 2% of body weight/day for 56 days (8 weeks). All tanks were maintained under natural lighting. During the entire feeding trial period, the flowing water rate was kept at 0.2 L/min and the water quality was maintained within the following ranges: a temperature of 22 ± 3 °C, a pH of 7.3 ± 0.3, a dissolved oxygen concentration of 5.8 mg/L, and a salinity of 30. 2.3. Sampling At the end of the 8-week feeding trial, blood and intestine samples were taken from five fish that were randomly chosen from each tank. Before sampling, the fish were fasted for 24 h. The fish were first anesthetized with MS222 (tricaine methanesulfonate, 100 mg/L), and then their body weight and body length were measured. Afterwards, blood was collected from the caudal artery and allowed to set in a 4 °C refrigerator for 24 h. The blood was centrifuged at 4000 rpm for 10 min, and then the supernatant was collected and stored at −80 °C for the determination of non-specific immunity parameters. After the blood sampling, the fish liver, spleen, and mid-intestine were all excised. The liver and spleen were weighed to measure the hepatosomatic index and the spleen-somatic index. The mid-intestine was frozen in liquid nitrogen for 24 h, then transferred to −80 °C until the RNA extraction. 2.4. Growth Performance, Hepatosomatic Index, and Spleen-Somatic Index At the beginning and end of the 8-week feeding trial, the total weight of all the fish in each tank was measured, and the average fish weight (the total weight/the total number of fish) of each tank was used as one replicate for the average fish weight of each group. Then the livers and spleens of these fish were excised and weighed. The weight gain rate (WGR, %), specific growth rate (SGR, %), hepatopancreas somatic index (HSI, %) and spleen-somatic index (SSI, %) were calculated as follows: WGR (%) = [100 × (final weight − initial weight)/initial weight]; SGR (% day^−1) = 100 × [ln (final weight) − ln (initial weight)]/days; HSI (%) = (liver weight/body weight) × 100; SSI (%) = (spleen weight/body weight) × 100. 2.5. Biochemical Analysis The amounts of superoxide dismutase (SOD, Cat. No. A001–3), lysozyme (LZM, Cat. No. A050-1), catalase (CAT, Cat. No. A007-1), alkaline phosphatase (ALP, Cat. No. A059-2), and total antioxidant capacity (T-AOC, Cat. No. A015-2) of the serum samples were measured using microplate methods, following the instructions of the commercial assay kits (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). These indices were measured to reflect the changes in the non-specific immunity indices. 2.6. Challenge Test A fish pathogen Vibrio harveyi strain [62]GDH11385 was isolated from a previous study [[63]27] and maintained as frozen stock in a −80 °C freezer. For the challenge test, the frozen stock was resuscitated. The bacterial suspension was adjusted to a concentration of 8.46 × 10^6 CFU/mL. On the day of the challenge test (Day 56), 30 fish from each group were intraperitoneally injected with 80 μL of the Vibrio harveyi strain [64]GDH11385 at a dose of 8.46 × 10^6 CFU/mL. At the same time, another 30 fish from each group were injected with 80 μL of PBS to serve as the control group. The applied LD[50] value of the V. harveyi strain [65]GDH11385 (8.46 × 10^6 CFU/mL × 80 μL ÷ 38 g fish body weight = 1.7 × 10^4 CFU/g fish body weight) was based on the previous work of our laboratory [[66]27], and a preliminary LD50 test was conducted before the formal challenge experiment. The mortality was observed and recorded for 7 days after the challenge tests. All dead and surviving fish in the challenge tests were examined for bacteria to verify the presence of the pathogen. 2.7. RNA Extraction, cDNA Library Construction, and Transcriptome Sequencing RNA extraction, cDNA library construction, and transcriptome sequencing were all conducted by Meiji Biotechnology Co., Ltd. (Shanghai, China). Total RNA was extracted from the middle intestines of the fish by using an RNA extraction kit (Promega, Shanghai, China) according to the manufacturer’s instructions. The purity and the concentration of the RNA were examined with a NanoDrop 2000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA). The quality of the RNA was assessed with 1% agarose gel electrophoresis. The RIN score was measured with an Agilent 2100 bioanalyzer (Agilent, Shanghai, China). It was only when the quantity of the obtained RNA was more than 1 g and of high quality (OD260/280 ≥ 1.8, OD260/230 ≥ 1.0) that it could be used to construct a single cDNA library. The eukaryotic mRNA was enriched using magnetic beads with Oligo (dT). The obtained mRNA was randomly broken-up by adding a fragmentation buffer, and then the mRNA fragments of about 300 bp were separated using magnetic beads. Afterwards, these mRNA fragments were used as the template for the first-strand cDNA synthesis using random hexamer primers and reverse transcriptase, followed with second-strand cDNA synthesis. The double-stranded cDNA was repaired by adding an End Repair Mix and adding an A tail. Finally, ten cDNA libraries (five libraries for each group) were constructed and sequenced on the Illumina Novaseq 6000 platform at Majorbio Bio-Pharm Technoloay Co., Ltd. (Shanghai, China), following the manufacturer’s instructions. 2.8. De Novo Assembly and Gene Functional Annotation The raw sequencing reads were filtered by applying fastx_toolkit_0.0.14 ([67]http://hannonlab.cshl.edu/fastx_toolkit/, accessed on 30 May 2019) and fastq_0.19.5 ([68]https://github.com/OpenGene/fastp, accessed on 30 May 2019) to obtain clean reads for the subsequent de novo assembly through Trinity Software ([69]https://github.com/trinityrnaseq/trinityrnaseq/wiki, accessed on 30 May 2019). After processing with the TransRate software ([70]http://hibberdlab.com/transrate/, accessed on 30 May 2019), any redundant sequences were removed, and unique transcript fragments (unigenes) were generated. A BLAST search of the unigenes against several databases was conducted, including the non-redundant protein database (NR), the Swiss-Prot database, the protein families database (Pfam), the Clusters of Orthologous Groups database (COG), the Gene Ontology database (GO), and the Kyoto Encyclopedia of Genes and Genomes database (KEGG). 2.9. Analysis of Differentially Expressed Genes The software RSEM v1.3.3 was used to record the mapping results and to quantify gene expression. The results of gene abundance were normalized and correlated with the TPM (Transcript per million reads) value. The DESeq2 1.26.0 software was used for screening differentially expressed genes (DEGs), with fold change (FC) ≥ 2, FC ≤ 2, and p value < 0.05 as screening criteria. These DEGs were displayed with volcanic maps. Then all downregulated DEGs with |Log[2]FC| ≥ 2 were screened and their distribution in KEGG pathways was categorized to explore the function of the most significantly downregulated DEGs. Finally, representative KEGG pathways and relative DEGs were further explored. 2.10. Validation of RNA Sequencing Data Using Quantitative RT-PCR (qPCR) To further validate the gene expression results from the RNA sequencing data, fifteen differently expressed unigenes were randomly selected for q-PCR analysis. The genes included Claudin 10, Myosin 11, Myosin 10, Caspase 3, Cathepsin L, Cathepsin H, alpha-Tubulin 8, etc. Firstly, we conducted reverse transcription of all RNA samples by using the Eastep^® RT Master Mix Kit (Promega, Beijing, China), and then the cDNA were collected for the follow-up qPCR experiment. The Primer premier 5 software (Premier Biosoft International, San Francisco, CA, USA) was applied to design the gene-specific qPCR primers, as shown in [71]Table S1. An amount of 18S mRNA was selected as a housekeeping gene [[72]26]. The qPCR experiment was performed in a 10 μL reaction volume, containing 0.2 μL of forward primer, 0.2 μL of reverse primer, 1 μL of cDNA template, 3 μL of nuclease-free water, and 5 μL of 5 × SYBR^® Premix Ex Taq, according to the manual of the ChamQ Universal SYBR qPCR Master Mix kit (Vazyme, Nanjing, China), on an Eppendorf Mastercycler^® ep Realplex4 instrument (Hamburg, Germany), with the following program: denaturation at 95 °C for 2 min, followed with 40 cycles of 95 °C for 10 s, and annealing at 60 °C for 10 s and at 72 °C for 15 s. The relative gene expression was calculated and normalized with the 2^−ΔΔCt method [[73]28]. All PCR reactions were repeated three times. 2.11. Statistical Analysis The growth indices, serum biochemical parameters, and gene expression were analyzed with the One-way analyses of variance (ANOVA) and the Tukey’s test using SPSS 23.0 (Chigaco, IL, USA). The results are presented as mean ± SE, with p < 0.05 indicating a significant difference. 3. Results 3.1. Growth Performance, Antioxidant Capacity, and Non-Specific Immunity Indices After 8 weeks of PC supplementation, the PC group showed no significant differences in growth when compared with the control group (p > 0.05). The SSI and HSI of groupers from different groups showed no significant differences either (p > 0.05) ([74]Table S2). The antioxidant T-AOC and CAT indices of the PC group were both significantly lower than those in the control group. There was no significant difference in the amounts of SOD. Significantly higher amounts of ALP and LZM were observed in the PC group ([75]Table 1). Table 1. Statistics of the antioxidant capacity and immunity in serum. T-AOC CAT SOD ALP LZM C 7.04 ± 1.42 * 47.42 ± 17.16 * 19.75 ± 4.03 41.19 ± 30.05 112.65 ± 21.52 PC 4.59 ± 1.28 33.763 ± 28.02 21.45 ± 3.59 107.59 ± 24.68 ** 222.65 ± 24.50 ** [76]Open in a new tab The results are presented as mean ± SE (n = 5), and different asterisks above each figure mean significant differences in the Tukey’s test (p < 0.05). Control: control group; and PC: Phellodendri Cortex group. * p < 0.05; and ** p < 0.01. 3.2. De Novo Assembly and Sequence Annotation All raw reads were uploaded and deposited in the National Center for Biotechnology Information (NCBI, PRJNA9008813). After filtration, 73.93 Gb of clean reads, with error rate < 0.025%, Q20 > 98%, Q30 > 94%, and about 48% GC of content ([77]Table S3), were assembled to obtain a total of 67,825 unigenes from 10 samples ([78]Table 2). The length distribution of the unigenes showed that almost 60% of the unigenes were clustered in the range of 0–1500 bp ([79]Figure S1). Table 2. Statistics of the transcriptomic sequences from six libraries. Database Unigene Number Percentage NR 26,128 38.52% Swiss-Prot 20,629 30.42% Pfam 18,636 27.48% COG 11,718 17.28% GO 7745 11.42% KEGG 16,156 23.82% Total annotated unigenes 26,864 39.61% Total unigenes 67,825 100.00% [80]Open in a new tab 3.3. Functional Annotation A total of 67,825 unigenes were annotated in six databases (Nr, Swiss-Prot, Pfam, COG, GO, KEGG), among which 26,864 unigenes were annotated in at least one of the six databases. There were 26,128 (38.52%), 20,629 (30.42%), 18,636 (27.48%), 11,718 (17.28%), 7745 (11.42%), and 16,156 (23.82%) unigenes annotated in the Nr, Swiss-Prot, Pfam, COG, GO, and KEGG databases, respectively ([81]Table 2). The Nr annotation revealed that the top five species showing the highest similarity with the intestine transcriptome of E.fuscoguttatus were Lates calcarifer (4893, 18.97%), Larimichthy crocea (4622, 17.92%), Seriola dumerili (2727, 10.57%), Seriola lalandi (1837, 7.12%), and Stegastes partitus (1269, 4.92%) ([82]Figure S2). The GO annotation revealed that a total of 34,561 unigenes were annotated in three GO categories: molecular function, cellular component, and biological process ([83]Figure S3). In the molecular function category, binding (3790, 42.37%) and catalytic activity (3223, 36.03%) were the two classes with the most annotations; in the cellular component category, most unigenes were annotated in the classes of cell (2364, 18.57%), cell part (2324, 18.25%), membrane (2281, 17.91%), and membrane part (2082, 16.35%); and in the biological process category, most unigenes were annotated in the classes of cellular process (2820, 21.74%), metabolic process (1976, 15.23%), and single-organism process (1925, 14.84%). In addition, the annotation in the KEGG pathways showed that the majority of annotated unigenes belonged to six of KEGG’s first categories (Signal transduction, Cancers: Overview, Immune system, Endocrine system, Transport and catabolism, and Infectious diseases: Viral), which contained twenty of KEGG’s second categories ([84]Figure S4). 3.4. Identification and Annotation of DEGs in the E. fuscoguttatus Intestinal Transcriptome in Response to PC Supplementation When the samples of the PC group were compared to those of the control group, a total of 1356 DEGs, including 298 upregulated DEGs (with an elevated expression in the PC group) and 1058 downregulated DEGs (with a decreased expression in the PC group) were identified (as shown in [85]Figure 1). The top 20 most significant enriched KEGG pathways (sub-categories, as shown in [86]Figure 2) belonged to five KEGG categories: human diseases (7 pathways); metabolism (5 pathways); organismal systems (4 pathways); cellular processes (3 pathways), and environmental information processing (1 pathway). In particular, the highest numbers of DEGs belonged to pathways related to epithelium integrity (Tight junction, 13 DEGs), apoptosis (Apoptosis, 11 DEGs; Parkinson’s disease, 11 DEGs; Alzheimer’s disease, 11 DEGs; and Huntington’s disease, 10 DEGs), and immune regulation (Cytokine–cytokine receptor interaction, 11 DEGs; Antigen processing and presentation, 9 DEGs; Phagosome, 10 DEGs; Influenza A, 10 DEGs; Legionellosis, 8 DEGs; and Rheumatoid arthritis, 6 DEGs). Representative enriched DEGs included caspase-3, cathepsin B, cathepsin L, claudin-11, claudin-1, and 14-3-3 protein beta/alpha-1 ([87]Table 3). Figure 1. [88]Figure 1 [89]Open in a new tab Volcano plot of DEGs. Figure 2. [90]Figure 2 [91]Open in a new tab KEGG pathway enrichment analysis of the annotated DEGs. The Y-axis represents the KEGG pathway; the X-axis represents the rich factor. The dot size represents the number of DEGs of the pathway, and the dot color represents the q-values. Table 3. Selected DEGs involved in the PC stimulation of E. fuscoguttatus. Gene Id Gene Description Gene Symbol KEGG Categories Log2(FC) Regulate Apoptosis/Immune responses TRINITY_DN21229_c6_g3 cathepsin B-like CTSB Immune system; Cell growth and death −2.2 down TRINITY_DN15782_c3_g2 cathepsin L, partial CTSL1 Immune system; Cell growth and death −1.3 down TRINITY_DN16713_c4_g5 cathepsin L, partial CTSL2 Immune system; Cell growth and death −1.2 down TRINITY_DN18426_c2_g3 heat shock protein 90 beta HSP90B Immune system; Cardiovascular disease; Endocrine system −2.6 down TRINITY_DN20083_c0_g1 caspase-3-like CASP3 Cell growth and death; Immune system; Infectious disease: viral; Infectious disease: bacteria 1.0 up Epithelium integrity TRINITY_DN24936_c0_g1 claudin-11 CLAU-11 Cellular community-eukaryotes; Immune system; Signaling molecules and interaction; Infectious disease: viral −2.8 down TRINITY_DN20497_c0_g2 14-3-3 protein beta/alpha-1 14-3-3PB/A1 Aging; Signal transduction; Cell growth and death; Cancer: overview; Infectious disease: viral −2.2 down TRINITY_DN13254_c0_g2 myosin light polypeptide 6 MYOLP6 Cellular community-eukaryotes; endocrine system; circulatory system −1.7 down TRINITY_DN19287_c0_g2 tubulin alpha-1B chain-like TUA-1B Cell growth and death; Cellular community- eukaryotes; Infectious disease: bacterial; Transport and catabolism −1.6 down *TRINITY_DN19417_c2_g1 integrin alpha-6-like IL-6 Signal transduction; Signaling molecules and interaction; Cellular community-eukaryotes; Immune system; Cancer: overview −1.3 down TRINITY_DN19287_c0_g3 tubulin alpha-8 chain-like TUA8 Cell growth and death; Cellular community- eukaryotes; Infectious disease: bacterial; Transport and catabolism −1.2 down TRINITY_DN9279_c0_g2 claudin-10-like isoform X1 CLAU-10 Cellular community-eukaryotes; Immune system; Signaling molecules and interaction; Infectious disease: viral 1.9 up [92]Open in a new tab 3.5. Validation of 15 Putative DEGs using qPCR The gene expression of 15 randomly selected DEGs, including Claudin 10, Myosin 11, Myosin 10, Caspase 3, Cathepsin L, Cathepsin H, Tubulin-alpha 8, etc, were analyzed. The comparative results are shown in [93]Figure 3. It is evident that the RNA sequencing data is valid and the transcriptome analysis is credible. Figure 3. [94]Figure 3 [95]Open in a new tab Validation of the RNA-Seq results and comparison of the gene expression between the control group and the PC supplementation group of E. fuscoguttatus using qPCR. The expressions of genes were examined relative to 18S RNA expressions as the endogenous control gene. The expression values were normalized against those of the control group. The results are the mean ± SEM of intestine tissues derived from five individual fish (n = 3). CASP3: Caspase 3; HYP: Hypothetical protein EH28_09075; CATL: Cathepsin L; TUA8: Tubulin alpha-8; CTSB: Cathepsin B; CLAU-10: Claudin 10; MYO11: Myosin 11; MYO10: Myosin 10; ADGPL3: Adhesion G protein-coupled receptor L3; GJ32P: gap junction Cx32.2 protein; GJD3P: gap junction delta-3 protein; GJ32PL: gap junction Cx32.2 protein-like; rGAP: ras GTPase-activating-like protein; AKTIP: AKT-interacting protein; and AKT1S1: proline-rich AKT1 substrate 1. 3.6. Disease Resistance Due to the presence of strain [96]GDH11385 [[97]27], though mortality was observed for 7 days after the challenge, [98]Figure 4 only shows cumulative mortality records from 0 to 36 h, because fish deaths only occurred within the first 36 h. Figure 4. [99]Figure 4 [100]Open in a new tab Cumulative mortality (%) of E. fuscoguttatus fed with PC for 8 weeks after the Vibrio harveyi challenge. Abbreviations: C, control group; and PC, Phellodendron chinese cortex fed group. Note: Data are expressed as the mean ± SD (n = 3). No significant difference between groups (p > 0.05). The 8-week mortality rate was recorded, and the comparison charts are displayed in [101]Figure 4. The challenge results showed that 8 weeks of supplementation with the PC water extract caused no significant differences in disease resistance against Vibrio harveyi in tiger groupers. 4. Discussion As a traditional Chinese medicine, Phellodendri chinese cortex is commonly used in combination with other traditional Chinese herbs to treat inflammatory disease, diarrhea, dysentery, and fatty liver disease, and to help with weight reduction [[102]29,[103]30]. Studies of phellodendri Chinese cortex used as an immunopotentiator by itself are rare. To our best knowledge, the only study of the phellodendri Chinese cortex in fish was Cai et al. [[104]26], in which the short-term (14-day) supplementation with PC was found to be immunosuppressive and to potentially enhance the disease resistance of E. fuscoguttatus against Vibrio harveyi significantly [[105]26]. Catalase (CAT) and the total antioxidant capacity (T-AOC) are both important antioxidant parameters. Widely distributed in tissues, CAT can exercise an antioxidant function through degradation of hydrogen peroxide [[106]31]. The T-AOC aggregates all antioxidants and is applied mostly for scavenging ROS and iNOS [[107]32]. In this present study, the CAT and the T-AOC of the grouper serum significantly declined after 8 weeks of supplementation with PC. This indicated that antioxidant enzymes and antioxidant non-enzymatic molecules significantly reduced in response to oxidative stress. These changes were different from the results of Cai et al. [[108]26]. Chen et al. [[109]33] found that the T-AOC in muscle was significantly reduced in Rhynchocypris that were fed highly oxidized fish oil compared with a fresh fish oil group after 4 weeks. Lysome (LZM) and Alkaline phosphate (ALP) are widely perceived as immune markers in serum. Lysozyme (LZM) is a crucial enzyme that disturbs the pathogens’ cytoderm, whereas ALP is a classical hydrolase. The increased amount of LZM and ALP in serum normally suggests enhanced immunity in aquatic organisms [[110]34]. The amount of LZM and ALP in the serum of groupers increased significantly after 8 weeks of supplementation with PC. This is different from Cai et al. [[111]26], who found no significant differences in ALP activity after 1 week or 2 weeks of PC supplementation to the diet of groupers. In summary, the long-term exposure to PC seemed to decrease the antioxidant capacity of tiger groupers and to enhance the non-specific immunity of tiger groupers. The intestinal tract is not only involved in food digestion and absorption but it also plays an important role in immune, metabolic, and osmotic regulation [[112]35,[113]36]. The intestine transcriptome DEG analysis of the control vs. the PC group revealed that the top 20 pathways in the KEGG database mainly involved three major categories: neurodegeneration and apoptosis; immune and infectious diseases; and endocrine, digestion, and metabolism. Apoptosis is a form of programmed cell death that is critical to physiologic homeostasis in almost every organ system. The Caspase family plays an important role in apoptosis activation and executions, especially in the mRNA and protein expression of caspase 3, which is an essential regulator of programmed cell death. The Caspase 3 gene encodes a protein that triggers cell death through cleaving specific proteins in the cytoplasm and nucleus. Cell apoptosis is induced by various biotic factors, such as bile acid disorder [[114]37], and various abiotic factors, such as Carnitine [[115]38], alcohol, radiotherapy, severe burns [[116]39], etc. Previous studies of Caspase 3 in aquatic organisms revealed that in grass carp (Ctenopharyngodon idella) [[117]40] and turbot (Scophthalmus maximus (L.) [[118]41] long-term stress and pathogen infection could both bring about organ damage by inducing an apoptosis response through a significant increase in the expression of caspase 3 in the spleen. In our research, the qPCR and RNA-Seq results both indicated that the mRNA expression of caspase 3 was significantly increased after 8 weeks of PC supplementation. This result is also in agreement with Eissa et al. [[119]42], who studied berberine, the major component of PC, in rats. Eissa et al. [[120]42] reported that the Caspae-3 mRNA level was significantly upregulated in rats that received berberine (50 mg/kg/day) via oral gavage for 6 consecutive weeks [[121]42]. Contrary to the pro-apoptosis effects of the long-term administration of PC and berberine, several studies [[122]43,[123]44] have shown that the short-term administration of berberine has significant anti-apoptosis effects. For example, a 7-day treatment of berberine (50 mg/kg or 200 mg/kg) induced significantly downregulated caspase-3 in mice (p < 0.01) [[124]43]; a 10-day administration of berberine (5 mg/mL) resulted in a significantly decreased expression of the pro-apoptosis caspase-3 genes in the epithelial cells of mice [[125]44]. These results indicate that the long-term administration of PC might induce damage to the physical barrier of healthy tiger groupers’ intestines through activation of the apoptosis process. Generally speaking, when the apoptosis response increases, it may be accompanied by changes in other cellular processes of organisms [[126]45]. The DEG analysis of this present study showed that there were a pronounced number of DEGs, which indicated impaired epithelium integrity and increased permeability in the intestines of tiger groupers. For example, our RNA-seq data demonstrated that the expression of Claudin 10 was upregulated and the expression of Claudin 11 was downregulated. These two genes both belong to the claudin family (which are integral membrane proteins and components of tight junction strands) yet have opposite functions. While the expression of Claudin-10 promotes the tight junction cation pore in the intestinal epithelium (to increase permeability), the expression of Claudin 11 enhances barrier integrity, because the Claudin 11 protein is an obligatory protein for tight junction formation [[127]46,[128]47]. In the intestinal tract, the overexpression of Claudin 10 or the insufficient expression of Claudin 11 is usually induced by pathogen infection and considered as a sign of intestinal damage [[129]48]. This damaging effect of long-term PC supplementation on intestine cell structure and permeability was also reflected by the downregulation of many other DEGs that are relevant to structural constituents, such as tubulin alpha-1B (TUBA1B) chain-like, tubulin alpha-8 chain-like, myosin light polypeptide 6, 14-3-3 protein beta/alpha-1, and integrin alpha-6-like. Alpha tubulins are a core protein family that form microtubules. TUBA1B is a structural constituent of the cytoskeleton, which is important in maintaining a stable backbone in both axons and dendrites. Lower levels of TUBA1B were observed in patients with schizophrenia (a serious mental disorder), compared to controls [[130]49]. TUBA1B is also positively correlated with the process of DNA replication and mitotic cell cycle. In mice, the TUBA8 mRNA/protein was significantly higher in the fatty liver and tumor tissues of mice and in hepatocarcinoma (liver cancer) patients [[131]50,[132]51]. The downregulation of the TUBA1B gene in normal groupers reflects both the damaging effect on intestinal integrity and the potential anti-fatty liver and anti-cancer effects of long-term PC supplementation. Myosin light polypeptide 6 is a structural constituent of muscle, regulating the function of myosin. Downregulated myosin gene expressions induce decreased levels of myosin, which may compromise the mitochondrial structure and consequently its function, leading to the loss of cellular energetics, an increase in the amount of free radicals, and eventually cell death [[133]52]. 14-3-3 proteins were expected to play a significant role in cellular osmoregulation. Expression of the 14-3-3 beta/alpha-1 transcripts in fish were most abundant in the gill and to a lesser extent in the intestine. The downregulation of this gene in PC-supplemented groupers indicated the impaired function of osmotic regulation [[134]53]. The anti-inflammatory effect of PC as a traditional herbal medicine used to treat various infectious or inflammatory diseases has long been recognized. In this present study, both the transcriptome and RT-PCR analysis revealed that the PC exhibited strong downregulating effects on the immune- and inflammation-associated genes and pathways of tiger groupers’ intestines. In particular, along with many other immune-associated genes (HSPA1_8 and HSPA90A, for example), the cathepsin B and cathepsin L genes (CTSBs and CTSLs) were significantly downregulated in the PC-supplemented group. CTSBs and CTSLs, with similar structures and functions, are ubiquitously expressed lysosomal cysteine proteinases that play important roles in inflammation, infection, and cancer, where they are highly expressed [[135]54]. In channel catfish infected with Edwardsiella Ictaluri and grouper fish challenged with Singapore grouper iridovirus and Vibrio anguillarum, strongly upregulated CTSBs and CTSLs were observed along with many other pro-inflammatory genes in the intestines [[136]55,[137]56]. Studies have shown that the activation of CTSB and CTSL gene expression after infection is associated with the destruction of elastin-rich tissues during inflammatory responses, and that the inhibition of CTSBs and CTSLs could effectively ameliorate the mucosal damage of inflammatory diseases such as colitis [[138]57]. Therefore, in addition to its anti-inflammatory effects, the supplementation with PC also exhibited a potential cathepsin inhibition effect, which might be helpful in the prevention and therapy of fish enteritis [[139]57]. Like many medicinal herbs, PC has an optimum dietary intervention period to exhibit the best anti-Vibrio effect in grouper fish. Cai et al. [[140]26] and this present study show that, while 7-day and 56-day supplementation with PC did not significantly improve the disease resistance of tiger groupers against V. harveyi, the 14-day supplementation with PC significantly decreased the cumulative mortality of tiger groupers after the V. harveyi challenge. Therefore, the disease resistance-improving effect of PC varied with different supplementation periods. The trial design had some limitations. In future, we need to consider the PC dose treatment, the sampling time, and the challenge dose to verify the effect of the PC. 5. Conclusions This study demonstrated that a long-term (8-weeks) PC diet did not affect the growth performance of tiger groupers, significantly increased the ALP and LZM activity in grouper serum, and had no significant effect on the disease resistance of tiger groupers. The transcriptome analysis revealed that long-term PC supplementation to groupers’ diet may potentially damage the structural and functional integrity of healthy groupers’ intestines. On the other hand, long-term PC supplementation also exhibited potential anti-inflammatory and cathepsin inhibition effects, which might be helpful for the prevention and therapy of fish enteritis. Supplementary Materials The following supporting information can be downloaded at: [141]https://www.mdpi.com/article/10.3390/life13122336/s1, Table S1. Primers used for qPCR; Table S2. Effects of PC diet on growth performance of tiger grouper after 8-week feeding trial; Table S3. Statistics of unigenes assembled; Table S4. Statistics of annotation results for intestine transcriptome unigenes; Figure S1. The length distribution of the unigenes; Figure S2. The species distribution based on Nr annotation; Figure S3. Top 47 annotated GO terms in the intestine transcriptome of E. fuscoguttatus; Figure S4. Annotated KEGG pathways in the intestine transcriptome of E. fuscoguttatus. [142]Click here for additional data file.^ (550.1KB, zip) Author Contributions Y.C. designed the experiments, provided funding, and drafted the article. H.S. conducted the majority of the experiments and data analysis. Y.Z. (Yu Zheng) helped with the fish rearing and sampling experiments. Y.Z. (Yongcan Zhou) provided funding and was in charge of project administration. W.G. helped with the collection of samples and data analysis. J.L. helped with fish rearing system maintenance and the observation of clinical signs. S.W. designed the experiments, provided funding support, and provided critical writing, both review and editing. All authors have read and agreed to the published version of the manuscript. Institutional Review Board Statement The animal study protocol was approved by the Ethics Committee of Hainan University (protocol code: HNUAUCC-2020-00011 and date of approval: 3 January 2019). Informed Consent Statement Not applicable. Data Availability Statement All data generated or analyzed during this study are included in this article. Raw data are available from the authors upon reasonable request. RNA-seq data have been deposited in the Sequence Read Archive (SRA) APP. The accession code is PRJNA908813. Conflicts of Interest The authors declare no conflict of interest. Funding Statement This study was funded by the National Natural Science Foundation of China (Grant No. 31860739 and Grant No. 32173009) and the National Key Research and Development program (2021YFC2600). Footnotes Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. References