Transport Stress Induces Skin Innate Immunity Response in Hybrid Yellow Catfish (Tachysurus fulvidraco♀ × P. vachellii♂) Through TLR/NLR Signaling Pathways and Regulation of Mucus Secretion

The transport of live fish is a necessary step for commercial production. The skin of teleost fish is the first non-specific immune barrier against exogenous stimuli, and it plays an important protective role under transport stress. Thus, the aim of this study was to explore the skin responses to transport stress in hybrid yellow catfish (Tachysurus fulvidraco♀ × Pseudobagrus vachellii♂) through transcriptome and biochemical analyses. Water samples were collected during a simulated transport treatment. Biochemical indexes and/or gene expression in blood, skin, and mucus in fish in control groups and transport-stress groups (0 h, 2 h, 4 h, 8 h, 16 h) were assayed. The levels of total ammonia–nitrogen and nitrite–nitrogen in the water increased with increasing transport time. Comparison of skin transcriptomes between the control group and the group subjected to 16 h of transport revealed 1547 differentially expressed genes (868 up-regulated and 679 down-regulated). The results of the transcriptome analysis were validated by analyses of the expression levels of selected genes by qRT-PCR. The results indicated that the toll-like receptors and nod-like receptors signaling pathways mediate the skin’s immune response to transport stress: tlr9, mfn2, and ikbke were significantly up-regulated and nfkbia and map3k7cl were significantly down-regulated under transport stress. With increasing transport time, lysozyme activity and the immunoglobulin M content in skin mucus first increased and then decreased. The number of mucous cells peaked at 8 h of transport stress, and then decreased. The mucus cells changed from types II and IV to types I, II, III, and IV. The amounts of red and white blood cells and the levels of hemoglobin and hematocrit first increased and then decreased during 16 h of transport stress. Together, the results showed that the skin responds to transport stress by activating the immune signaling pathway and regulating mucus secretion. These findings have important biological significance for selecting strains that tolerate transport, as well as economic significance for optimizing the transport conditions for scaleless fish.


INTRODUCTION
The transportation of live fish is an essential step to meet the demand for fish in different regions. The extent of harm caused by transport depends on the duration and severity of the stress, as well as the health status of the fish (1). During transportation, fish are in an enclosed environment, and their metabolic wastes affect water quality parameters such as the pH and the concentrations of ammonia, nitrite, and carbon dioxide (2)(3)(4)(5). The accumulation of total ammonia-nitrogen and nitritenitrogen in water can lead to transportation stress (6). With prolonged transit time, transport stress can cause physiological dysfunction in fish, and consequently diminish their resistance to pathogens, even resulting in death (7)(8)(9)(10). Studies have shown that transport stress affects the innate immunity of fish (9).
For teleost fish, the skin is a key component of innate immunity. It ensures that they can survive and adapt to the environment (11,12). Skin with its outermost mucus separates the individual from environment and is the crucial interface for communication and contact with external factors (13). The skin and its mucus often be used to study the immune response of fish under stress in aquaculture. In sea bream (Sparus aurata L.), chronic crowding stress was found to reduce the barrier function of damaged skin and suppress the local immune response to scale removal (14). In Epinephelus coioides, the skin showed a strong immune response to the invasion of Cryptocaryon irritans, and genes encoding Toll-like receptors (tlr1, tlr2, tlr5 and tlr5s) were found to be significantly up-regulated (15). As a skin appendage, mucus is produced by dynamic secretory cells in the epidermis. The change of its composition and rheological properties plays an important role in regulating mucosal immunity in fish (16). The different proteases found in fish surface mucus are involved in innate immunity (17). Among these proteases, immunoglobulin M (IgM) in skin mucus is the first responder to foreign invaders, and mucosal lysozyme (LZM) is a antimicrobial enzyme with antibacterial, antimetastatic, antiviral, and anti-inflammatory properties (18). Mucosal cells and the composition of mucus can be affected by exogenous factors (such as stress, acid, and infections) and endogenous factors (such as sex and developmental stage) (19)(20)(21). In sea bass (Laeolabrax japonicus), vigorous proliferation of skin mucus cells was observed under prolonged hypoxia and high nitrate levels (22). In large yellow croaker (Larimichthys crocea), the cortisol, malondialdehyde and IgM levels in skin mucus increased first and then decreased under high-temperature (23).
Hybrid yellow catfish (Tachysurus fulvidraco♀ × Pseudobagrus vachellii♂) shows good resistance and a fast growth rate, it has a short breeding cycle and a high meat content, and it is widely farmed (24). However, hybrid yellow catfish seem to be highly susceptible to transport stress. In recent years, high-throughput sequencing technology has been used to analyze the responses of aquatic animals to various stresses. Using this technology, researchers have discovered a series of molecular changes that occur under transport stress. For example, a total 1275 differentially expressed genes (DEGs) were identified in the liver of zebrafish (25) and 1000 DEGs were identified in the muscle of rainbow trout under transport stress (26). Those genes were found to be involved in metabolic processes, cellular homeostasis, and the immune response, indicating that these processes rapidly change under transport stress. Thus, transcriptome analysis is a useful tool to explore the responses of fish to transport stress. However, data on skin transcriptome responses of hybrid yellow catfish to transport stress are lacking. Here, we used transcriptome analyses to identify differentially expressed genes (DEGs) in yellow catfish under transport stress. Several randomly selected DEGs were verified by qRT-PCR to confirm the reliability of the RNA-Seq results. We also analyzed skin sections and mucus biochemical factors, including IgM content and LZM activity, in fish under transport stress. Our results revealed significant changes in gene expression levels in the skin and mucus attributes of yellow catfish under transport stress. These results shed light on the responses and regulatory mechanisms in skin of fish under transport stress.

Fish and Experimental Design
Juvenile yellow catfish with an average length (± SD) of 5.06 ± 0.82 cm and an average mass ( ± SD) of 7.06 ± 1.30 g were obtained from the Freshwater Fisheries Research Centre (Yixing, China). The fish were fed with commercial feed (crude protein 38.0%, crude fat 6%) at 2% (w/w) of their body weight twice per day for 7 days in a recirculating aquaculture system (temperature 26 ± 3°C DO > 5.5 mg/L; ammonia nitrogen and nitrite < 0.01 mg/L; pH 7.5 ± 0.3) before experiments. The experimental fish were subjected to a simulated transport treatment after starvation for 24 h. For the simulated transport treatment, 30 fish were placed in each double-layered plastic nylon bag (40 × 82 cm) with about one-third of the volume occupied by freshwater and two-thirds of the volume occupied by pure oxygen (5). The transport bags were sealed with rubber bands and placed on an analog transporter (simulated automotive shaker; Changzhou, China). The vibration frequency was set at 100 RPM (revolutions per minute) according to the method of Wu et al. (5). We established five transport groups (0 h, 2 h, 4 h, 8 h, and 16 h) and one group that was allowed to recover for 96 h after 16 h of transport stress. Each group had three replicates. Control fish of each group were obtained from the recirculating

Sample Collection
At each time point, 12 individual fish were randomly selected from each group and were anesthetized immediately with 200 mg/L MS-222 (27). Water samples were collected at the same time. Blood samples were taken from the tail vein using a 2.5-mL syringe to use for blood cell analysis. Before collecting mucus, water was gently removed from the fish. The fish were then placed in a sterile bag for 1 minute to collect the mucus (28). The skin mucus was immediately homogenized with Tris-buffered saline (TBS, pH 8.0, 50 mM Tris HCl, 150 mM NaCl) on ice and the mixture was centrifuged at 12000 g and 4°C for 15 min (29). The supernatant was collected and stored at −80°C until further analysis. Skin was collected from the left side of each fish behind the dorsal fin, and was separated into two subsamples: one subsample was used for histological analysis, and the second subsample was frozen in liquid nitrogen and then stored at −80°C until transcriptome analysis.

Water Quality Detection
We used a portable DO meter (YSI 556, Yellow Springs Instruments, Yellow Springs, OH, USA) to determine the DO content. According to Qiang et al. (30), we determined total ammonia-nitrogen (TAN) and nitrite-nitrogen (NO 2 -N) by spectrophotometric methods.

Blood Analysis
White blood cells (WBC); hemoglobin (HGB); hematocrit (HCT); and red blood cells (RBC) were quantified using a Blood cell analyzer (Mindray, bc-5300, Shenzhen, China) according to the manufacturer's instructions. All of the reagents were obtained from Shenzhen MINDRAY BioMedical (Shenzhen, China).

Mucus Analysis
The activity of LZM was assessed using a commercially available kit (Jiancheng Institute of Biotechnology). The IgM in skin mucus was detected using enzyme-linked immunoassay (ELISA) (Shanghai Lengton Bioscience Co., Ltd.) as described by Guardiola et al. (29,31).

Histological Analysis
The

Quality Control, RNA-Seq Reads Mapping, and DEGs Analysis
Unqualified raw data were filtered using Cutadapt to obtain clean data. We aligned reads from the control group and 16-h transport group to the UCSC (http://genome.ucsc.edu/) Homo sapiens reference genome using the HISAT package, which initially removed a portion of the reads based on quality information accompanying each read, and then mapped the reads to the reference genome. For each sample, the mapped reads were assembled using StringTie to reconstruct a comprehensive transcriptome using perl scripts. Then, all the transcriptomes were merged to reconstruct a comprehensive transcriptome using perl scripts. After the final transcriptome was generated, StringTie and edgeR were used to estimate gene transcript levels. StringTie was used to calculate the FPKM values, which were indicative of the mRNA levels. The DEGs between control and 16-h transport stressed fish were selected based on log2 (fold change) >1 or log2 (fold change) <-1 and statistical significance of p < 0.05 by the R package. Functional analysis of DEGs was conducted using tools at the GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases. The results of GO and KEGG enrichment analysis are presented as a scatter diagram, which was constructed using GGplot2. Based on the KEGG pathways, we identified the DEGs. We focused on identifying immune-related pathways in the skin of hybrid yellow catfish under transport stress and selected key genes for qRT-PCR verification.

Validation of DEGs in Skin by qRT-PCR
We verified the transcriptional patterns of DEGs under transport stress by qRT-PCR. Total RNA was extracted from skin using the RNAiso Plus kit (TaKaRa, Dalian, China). The quality and concentration of RNA were determined from the absorbance values at 260 nm and 280 nm (A260/A280 ratio, 1.9-2.1). Then, cDNA was synthesized by PrimeScript ™ RT reagent kit with gDNA Eraser (TaKaRa). The relative transcript levels of genes associated with TLRs and NLR signaling pathways were determined by qRT-PCR using the following program: one cycle at 95°C for 30s, then 40 cycles at 95°Cfor 5 s, and 60-63°C for 1 min. The amplification mixture (20 mL) comprised 2 mL cDNA, 10 mL TB Green Premix Ex Taq II (Takara), 0.5 mL each primer, and 7 mL ddH 2 O. b-actin was the reference gene. The relative transcript levels were determined by the 2 −DDCT method. The primers used are shown in Table 1.

Statistical Analysis
Differences within the same treatment group at different times were detected by one-way analysis of variance (ANOVA). Data are presented as mean ± standard error of the mean (SEM). Differences were analyzed by Duncan's multiple range test using SPSS 26.0 (SPSS Inc., Chicago, IL, USA). The independentsamples t test was used to detect significant differences between the control group and transport group at each time point. The level of significance was P <0.05.

Water Quality Parameters During Simulated Transport
As shown in Figure 1, there was no significant change in TAN and NO 2 -N in the water of the control groups. The TAN and NO 2 -N contents differed significantly between the control and the transport groups at all times except for 0 h. In the transport groups, TAN (P< 0.05; Figure 1A) levels clearly increased from 0 h onwards. The NO 2 -N (P < 0.05; Figure 1B) levels increased during transport, but were not significantly different between 2 h and 4 h. The DO levels were higher than 22 mg/L in the transport groups during the simulated transport treatment.

Changes in Blood Parameters in Fish During Simulated Transport
For control groups, there was no significant difference in RBC, WBC, HGB, and HCT among the sampling times. For transport groups, the RBC, WBC and HGB were significantly up-regulated at 4 h and 8 h compared with 0 h, and significantly downregulated at 16 h compared with 4 h and 8 h but still higher than at 0 h (P<0.05; Figures 2B, C). The HCT were significantly increased at 2 h and 4 h of transport compared with 0 h, peaked at 4h and then by 8 h and 16h they decreased to a level similar to that at 2h (P<0.05; Figures 2A, D). In addition, compared with the control groups, the transportation groups showed higher RBC, WBC, HGB, and HCT levels at all sampling times except for 0 h. In the 96h recovery group, the amounts of RBC and the levels of HGB and HCT were not significantly different from those in the control group. However, the amount of WBC was higher than that in the control group, but lower than that in the 16 h transport stress group (P<0.05; Figure 2B).

Mucus Biochemical Analyses and Skin Histological Analyses
The LZM activity and IgM content in mucus did not differ among time points in the control groups. However, the LZM content differed between the control groups and the transportation groups at all sampling times except for 0 h. The IgM content also differed significantly between the control groups and transportation groups at all sampling times except for 0 h. In the transportation groups, LZM activity was higher at 4 h and 8 h than at 0 h and 2 h of transport (P<0.05; Figure 3A) but significantly decreased by 16 h of transport. The IgM content showed a slight upward trend at 2 h and 4 h and was clearly higher at 8 h than at 0 h of transport, but decreased at16h (P<0.05; Figure 3B). The LZM and Igm levels in the 96h recovery group were higher than those in the control group but lower than those in the 16h transport stress group. The mucous cells were stained four different colors by AB-PAS. The results revealed that the skin mucus contained mainly type II and IV mucous cells, and the number of mucus cells was higher in the transportation groups than in the control groups by 4 h. However, after 8 h, the density of mucus cells decreased significantly. The mucus cells were types II and IV at 0 h, 2 h, and 4 h of transport but showed increased proportions of types I and III after 4 h of transport (Figure 4).

Analysis and Verification of DEGs
The skin transcriptome data of yellow catfish have been submitted to NCBI GenBank under the accession number GSE179444.

Identification of DEGs
We compared and analyzed the transcriptome data to identify the genes with significant differences in transcript levels between the con-s and str-s libraries, using the following criteria: | log2foldchange|≥1, P<0.05, and FPKM>10.We identified 1547 DEGs between the con-s and str-s libraries (868 up-regulated and 679 down-regulated) ( Figure 5).

GO Classification and KEGG Enrichment Analysis of DEGs
The DEGs were classified into biological process, cellular component, and molecular function categories by a Gene Ontology enrichment analysis. The GO category enrichment analysis showed that the cellular component subcategories most enriched with DEGs were membrane, integral component of membranes, and nucleus. In the biological process category, the subcategories most enriched with DEGs  were biological process, regulation of transcription, and DNAtemplate. In the molecular function category, the subcategories most enriched with DEGs were metal ion binding and molecular function ( Figure 6). The sequences of the DEGs were annotated using tools at the KEGG database. The KEGG pathway analysis revealed the pathways most enriched with DEGs in the skin of yellow catfish under transport stress (Figure 7), which included the Toll-like receptors (TLRs) signaling pathway, the NOD-like receptors (NLRs) signaling pathway, salmonella infection, cytokine-cytokine receptor interaction, and the AGE-RAGE signaling pathway related to diabetic complications.

Confirmation of Immune-Related DEGs by qRT-PCR
The 20 DEGs in the TLRs/NLRs signal pathway are shown in a heat map analysis (Figure 8). To confirm the reliability of the RNA-Seq data, we selected six immune-related genes for analysis by qRT-PCR. The selected genes were tlr9, encoding toll like receptor 9; mfn2, encoding mitofusin; txnipb, encoding thioredoxin interacting protein b; nfkbia, encoding NFKB inhibitor alpha; map3k7cl, encoding MAP3K7C-terminal like protein; and ikbke, encoding inhibitor of nuclear factor kappa B kinase subunit epsilon. All of these genes are related to the TLRs/NLRs signal pathways. The transcriptome sequencing results of these genes are shown in Table 3. Compared with the control group, the 16-h transportation group showed significantly higher transcript levels of tlr9, mfn2 and ikbke genes in the skin (P<0.05; Figure 9) and significantly lower transcript levels of txnipb, nfkbia and map3k7cl (P<0.05; Figure 9B). The correlation between transcriptome sequencing results and qRT-PCR results is shown in Figure 9A.

DISCUSSION
Previous studies have shown that transportation stress can induce an immune response and impair the health of fish (33). It is very important to explore the mechanism of how transport      (34). Changes in water quality during transportation can cause stress responses in fish. The concentration of TAN, the major excretory product of fish, in water is an important indicator of fish survival during transport (35). Fish were found to excrete high levels of un-ionized ammonia under transport stress, resulting in increased blood ammonia concentrations and reduced blood oxygen-carrying capacity (36,37). Increases in the TAN level promote the production of NO 3 − , which can be transformed into toxic derivatives, alter physiological functions, and elicit an immune response in the body (38). In the present study, the concentrations of TAN and NO 2 -N in the water increased as the transportation time increased. The highest TAN and NO 2 -N levels were 0.394 mg l −1 and 0.123 mg l −1 , respectively, in the 16h transport stress group. These compounds may have been among the factors that induced the immune response in fish under transport stress (6). Similarly, an inhibitory effect of high concentrations of ammonia on the innate immune response was detected in zebrafish during acclimation to transportassociated stress (25). Additionally, transport-induced changes in stress hormones can also affect the immune function of vertebrates. The way that stressors modify immune function is very complex (39), so further studies are required to determine the mechanism by which transport stress affects the immune response.
Innate immunity is the main immune system of teleost fish due to the underdevelopment of adaptive immunity (40,41). In the innate immunity of fish, the TLRs/NLRs signaling pathways are the two most widely studied pattern recognition receptors (42). The TLRs are type I transmembrane proteins that can be classified into two groups: cell surface TLRs and intracellular   TLRs. The NLRs are intracellular cytoplasmic receptors (43,44). Transcription of the genes encoding components of the TLRs/ NLRs pathways is induced by stress, leading to stimulation of the host's immune and inflammation responses (45)(46)(47). Studies have demonstrated that the activation of the TLR9 pathway provides an important link between immune and inflammatory phenotypes (48,49). In juvenile turbot, inflammation and the immune response were found to be related to increased expression of genes encoding TLRs, including tlr2, tlr3, tlr5b, tlr9, tlr21, and tlr22 (50). In rainbow trout, mechanical injury of muscle tissue was found to trigger the expression of immunerelated genes, especially tlr9, which showed significantly increased transcript levels at 4 h and 24 h after injury (2.5-fold and 7.2-fold increases, respectively) (51). In our study, we also detected increased transcript levels of tlr9 in yellow catfish skin at 16 h of transport stress. The transcript level of mfn2 and txnipb, encoding components of the NLRs signaling pathway, were upregulated and down-regulated respectively in the skin of fish under transport stress. mfn2 encodes a multifunctional mitochondrial fusion protein that is involved in the activation of innate immunity. This protein can exacerbate inflammation, while knockdown of mfn2 was shown to ameliorate inflammation (52,53). Similar to our results, another study found that mfn2 transcript levels were elevated in human and rat chondrocytes during metabolic disorders and inflammation (52). The txnipb plays an important role in pathophysiological consequences including elevated inflammatory response, cellular immunity, and tumorigenesis. Therefore, we speculate that transport stress might induce an immune response and aggravate inflammation through tlr9, mfn2 and txnipb mediated TLRs/NLRs pathways. Furthermore, activation of tlr9 could lead to the recruitment of MyD88, resulting in the activation of NF-kB and MAPK kinases (50). In our study, we detected significantly decreased transcript levels of nfkbia and map3k7cl, which also encode components of the TLRs/NLRs pathway, in the skin of fish at 16 h of simulated transport. The protein encoded by nfkbia was identified as a target for glucocorticoid-mediated immunosuppression (54). In mice, nfkbia was poorly expressed and the NF-kB signaling pathway was activated during myocardial infarction, while overexpression of nfkbia reduced the expression of inflammatory cytokines and cardiomyocyte apoptosis, and improved cardiomyocyte viability (55). Because map3k7cl was downregulated in Chinese patients with non-small cell lung cancer, Niu et al. speculated that map3k7cl in the leukocytes may contribute to pathogenesis via reducing inflammatory processes, regulating mitochondrial reactive oxygen species signaling, and mediating the Wnt pathway (56). Studies have shown that hyper-activated or incomplete knock-down of TLRs can lead to tissue injury and inflammatory diseases (57). The up-regulated expression of ikbke is also associated with the inflammatory response. In infected erythrocytes with low pathogenic avian influenza virus, the expression levels of MyD88, CCL5, and IKBKE were found to be increased (58). In summary, we conclude that 16 h of simulated transport stress induced skin innate immune responses and exacerbated inflammation. Fish immunity is associated with stronger mucosal defenses. The immune-related components of skin mucus such as LZM, IgM, and alkaline phosphatase play important roles in nonspecific immunity (59). As a key enzyme involved in innate immunity, LZM plays a role in defense against pathogens during the early developmental stages of bony fish (60). Its activity varies with fish age, size, and sex, and with the season, temperature, water quality, toxins, and different types of stress (61)(62)(63)(64). During the transport of yellow catfish, the activity of LZM in skin mucus might be related to the degree of stress and its duration (65). Previous research showed that in the serum LZM activity was significantly decreased in rainbow trout under transport or water pollution stress (66). IgM is the functional immunoglobulin that is essential for defense mechanisms in the mucosal immunity of fish. It is a common practice to add immunostimulants to enhance innate immune responses, including LZM activity and IgM levels, to relieve stress in fish during transport (67). In the present study, we found that LZM activity and IgM contents peaked at 8 h of transport, and then decreased. This is consistent with the results of Wu et al. (67), who found that LZM activity in the skin mucus and IgM content in the serum tiger grouper first increased and then decreased  during transport. The histological results were consistent with these results, and showed that the number and types of mucous cells changed with increasing transport time. The skin mucus cell can mirror the health status of fish. It is a valuable matrix for monitoring stress, pathogen exposure, and nutritional effects. Skin mucus secretion, and components of the mucus, change when aquatic animals are under stress conditions. For example, prolonged exposure to high temperature resulted in an increase and then a decrease in the secretion of skin mucus in yellow croaker (Larimichthys crocea) (23). Similarly, a significant increase in mucus secretion was detected in gilthead seabream (Sparus aurata L.) under acute crowding and anesthetic stress (68). In sea bass, compared with fish in a control group, those exposed to hypersaline conditions secreted more skin mucus with higher contents of cortisol, glucose, and protein (69). Skin mucus can protect fish against harmful environmental factors and prevent the loss of physiological metabolites and water (70). Increased mucus secretion is a protective mechanism of fish under stress, and is accompanied by enhanced immune function. However, when the degree of the stress exceeds the tolerance range, it can result in immune suppression and decreased mucus secretion. Further research is required to explore the functions of different mucous cell types. At the same time, investigating the difference in mucus production between normal and stressed fish is very important to improve their survival (71). These results may indicate that high-intensity transport stress partially disrupts skin mucosal immune function and causes immunosuppression. After 96 h of recovery, LZM activity and IgM levels were lower than those at 16 h of transport, but still higher than those in the control group, indicating that the immune function of fish skin mucus was still recovering. Notably, some studies have shown that teleost mucosal IgM is generated independently of systemic antibody (72). Therefore, there may be weak relationships between the tissue immune system and the mucosal cell immune system (73). Blood is also an important component of the immune system, and changes in hematological parameters can be indicative of the physiological health of fish (74,75). Under stress conditions, changes in the skin immunity of fish are often accompanied by changes in various hematological indices (76). The WBC, RBC, HGB, and HCT are the hematological variables that are most commonly monitored during stress (77). Ronza et al. (78) confirmed that RBC have important immunological functions through a blood transcriptomics analysis of healthy and diseased turbot (Scophthalmus maximus). WBC are involved in cellular immunity, an increase in the WBC count is often correlated with increased antibody production (79). Thus, on the basis of our results, we speculate that the increase in the number or specific volume of RBC and WBC was due to the transport stressinduced enhanced hematological-immunological response. Consistent with our results, Boaventura et al. (80) reported increased HGB, HCT, and total plasma protein in juveniles of Lophiosilurus alexandri after a 4-h transport treatment. Other studies have reported that exposure to ammonia and nitrite can lead to decreases in WBC, RBC, HGB, and HCT in Takifugu rubripes and yellow catfish (81,82). Therefore, the changes in the number or specific volume of RBC and WBC at 16 h of transport in this study may be due to immunosuppression caused by the increased contents of TAN and NO 2 -N in the water during a long transport period. After 96 h of recovery, the amount of RBC and the levels of HGB and HCT returned to normal, but the amount of WBC remained higher than that in the control group. This might indicate that blood immunity was not fully restored after transport stress. The modulation of blood parameters during stress is related to splenic contraction and hemodilution. The sustained contraction of the spleen increases the amounts of circulating WBC, RBC, and HGB in the blood (83)(84)(85). Therefore, the changes in hematological immunity are probably related to splenic contraction and changes in blood parameters that were secondary to transport stress.
In summary, we explored the effects of transport stress on yellow catfish by conducting physiological, histological, and molecular analyses of the skin. During the transport of hybrid yellow catfish, TAN and NO 2 -N accumulated in the water, and this may be one of the factors that triggered the immune response in the skin. On the one hand, the skin transcriptome results suggested that 16 h of transport stress may induce immune responses and stimulate pro-inflammatory responses by up-regulating tlr9, mfn2, and ikbke and down-regulating txnipb, nfkbia and mapk3cl in the TLRs/NLRs pathways. On the other hand, transport times shorter than 8 h may not seriously damage the mucous immune system, but longer transport times (16 h) may inhibit the immune function of mucus, leading to the decrease of LZM and IgM, and the change of the number and type of mucus cells. In addition, on the basis of the performance of the recovery group, we suggest that the physiological recovery time should be longer than 96 h after 16 h of transport stress.
Minimizing the effects of stress during transport of hybrid yellow catfish is important for successful aquaculture. In addition, a better understanding of the stress response of hybrid yellow catfish to transport is necessary to formulate optimal transport procedures. On the basis of our research, we identified two possible ways to reduce transport stress. Firstly, inadequate or excessive mucus secretion is not conducive to fish survival, so controlling mucus secretion during transport will improve the survival rate (71). Secondly, total ammonianitrogen and nitrite-nitrogen can cause immunosuppression, so controlling their levels in the water during transport will also reduce the adverse effects of transport stress.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

ETHICS STATEMENT
The animal study was reviewed and approved by Freshwater Fisheries Research Center at the Chinese Academy of Fishery Sciences (Wuxi, China). Written informed consent was obtained from the owners for the participation of their animals in this study.