- 1Yangtze River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Wuhan, China
- 2Fishery Resource and Environment Research Center, Chinese Academy of Fishery Sciences, Beijing, China
- 3College of Animal Science and Technology, Hunan Agricultural University, Changsha, China
Deltamethrin (Del) is a widely used pyrethroid insecticide and a dangerous material that has brought serious problems to the healthy breeding of aquatic animals. However, the toxicological mechanisms of Del on channel catfish remain unclear. In the present study, we exposed channel catfish to 0, 0.5, and 5 μg/L Del for 6 h, and analyzed the changes in histopathology, trunk kidney transcriptome, and intestinal microbiota composition. The pathological analyses showed that a high concentration of Del damaged the intestine and trunk kidney of channel catfish in the early stage. The transcriptome analysis detected 32 and 1837 differentially expressed genes (DEGs) in channel catfish trunk kidneys after exposure to 0.5 and 5 μg/L Del, respectively. Moreover, the KEGG pathway and GO enrichment analyses showed that the apoptosis signaling pathway was significantly enriched, and apoptosis-related DEGs, including cathepsin L, p53, Bax, and caspase-3, were also detected. These results suggested that apoptosis occurs in the trunk kidney of channel catfish in the early stage of acute exposure to Del. We also detected some DEGs and signaling pathways related to immunity and drug metabolism, indicating that early exposure to Del can lead to immunotoxicity and metabolic disorder of channel catfish, which increases the risk of pathogenic infections and energy metabolism disorders. Additionally, 16S rRNA gene sequencing showed that the composition of the intestinal microbiome significantly changed in channel catfish treated with Del. At the phylum level, the abundance of Firmicutes, Fusobacteria, and Actinobacteria significantly decreased in the early stage of Del exposure. At the genus level, the abundance of Romboutsia, Lactobacillus, and Cetobacterium decreased after Del exposure. Overall, early exposure to Del can lead to tissue damage, metabolic disorder, immunotoxicity, and apoptosis in channel catfish, and affect the composition of its intestinal microbiota. Herein, we clarified the toxic effects of Del on channel catfish in the early stage of exposure and explored why fish under Del stress are more vulnerable to microbial infections and slow growth.
Introduction
Pyrethroids are synthetic pesticides that were first isolated from the flower extract of Chrysanthemum. Due to its selectivity and low toxicity to non-target organisms, such as mammals and birds, pyrethroids are widely used to control agricultural and residential pests (1–4). However, pyrethroids are highly toxic to aquatic animals since they lack the corresponding degrading enzymes (5). Additionally, due to the lipophilicity of pyrethroids (6, 7), in the aquatic environment, they can easily enter into fish’s bodies through gills and are transported to each tissue through blood circulation and accumulate. Previous studies have shown that the LD50 of pyrethroids to aquatic animals is 10 to 1000 times lower compared to mammals and birds (8–11). Therefore, the pollution of pyrethroids can harm aquatic animals. Deltamethrin (Del) is an important pyrethroid with a broad-spectrum insecticidal activity that is widely used in aquatic animal pest control (12). Del can flow into rivers and lakes through surface runoff and domestic wastewater, threatening the health of aquatic organisms, including fish (13, 14). For example, Feoet al. reported that the concentration of Del in water samples from the Ebro River Delta in Spain was between 0.73 - 58.80 ng/L, and Amin et al. indicated a concentration in water samples of 0.20-2.00 μg/mL (15, 16). There have also been reports of a large number of deaths of aquatic animals caused by Del pollution (17). Hence, considering the high toxicity and environmental residues of Del, its harmful effects on aquatic organisms have attracted increasing attention.
Furthermore, Del has been proved to be highly toxic to different aquatic animals, including Chinese mitten crab (Eriocheir sinensis), Chinese rare minnow (Gobiocypris rarus), black tiger shrimp (Penaeus monodon), zebrafish (Danio rerio), and Tilapia mossambica (Oreochromis niloticus) (12, 18–24). Meanwhile, Del is an important stress source for freshwater ecosystems (20, 25). The potential health risks and toxicity of Del are mainly related to neurotoxicity and hepatorenal toxicity (26, 27). According to previous studies, Del can inhibit the activities of acid phosphatase, lysozyme, and phenoloxidase in crayfish (26). Long-term exposure to Del might also induce inflammation, oxidative stress, and apoptosis in different organs of carp in a dose-dependent manner (28). Del might also lead to developmental toxicity, neurotoxicity, and cardiovascular toxicity in zebrafish (29, 30). Additionally, subacute exposure to environmentally relevant concentrations of Del can inhibit the host immune response and disease resistance of snakehead fish (Channa argus) (24). Besides, Del exposure leads to histopathological changes and immunotoxicity of Chinese rare minnows, making them vulnerable to pathogen infections (31), and can cause gill tissue damage and immunotoxicity in crucian carp (12). These results demonstrated that different degrees of exposure to Del can have toxic effects on aquatic animals, especially on the immune system. Therefore, Del is a dangerous material, but its harmful molecular mechanisms remain unclear.
In the past decade, RNA-seq has been widely used in aquatic toxicology research to explore potential functional gene networks and gene mining (32). RNA-seq has also been used as a visualization tool to reveal the whole gene expression of immune response and signaling pathways of fish in specific physiological processes. Thus, we can explore the transcriptional landscape at the overall level and better describe the potential molecular mechanisms of Del-induced toxicity (33). The trunk kidney is an important organ for fish to transform and excrete exogenous substances. Meanwhile, it also plays a key role in immune response and maintaining homeostasis (34, 35). Hence, acute exposure to Del will inevitably cause a strong trunk kidney response. Transcriptome changes in the trunk kidney are of great significance to reveal the toxic mechanisms of Del in fish. The intestinal microbiota is also crucial to fish health and plays a key role in maintaining host physiology, nutrient supply, and intestinal immune response (36). Previous studies have shown that various environmental pollutants can cause changes in the composition of the intestinal microbial community, thereby affecting the health of aquatic animals (37, 38). Therefore, under acute exposure to Del, the composition of fish intestinal microbiota will change to adapt to the new environment and have an important impact on the host life activity processes. Thus, it is necessary to clarify the changes in fish intestinal microbiota under Del exposure to better understand its toxic effects.
Channel catfish (Ictalurus punctatus) was introduced in China from the United States in 1984. Due to their adaptation to the breeding environment in China, the breeding output and area have explosive growth, with an annual output of 308000 tons, becoming one of the important varieties of freshwater fish cultured with Chinese characteristics. However, diseases have always been the bottleneck restricting the further development of the channel catfish industry. Pathogenic microorganisms that infect channel catfish, including Aeromonas hydrophila, Yersinia ruckeri, Edwardsiella ictalurid, and Streptococcus iniae, have caused great economic losses to the channel catfish breeding industry (39, 40). Additionally, the inhibition of the immune system of channel catfish is an important cause of pathogenic microbial infections. Therefore, in the present study, we aimed to explore the toxic mechanisms of Del on channel catfish by analyzing the changes in trunk kidney transcriptome and the intestinal microbiome. We focused on the toxic effects of Del on channel catfish immune system to provide the theoretical basis for disease prevention and control.
Materials and Methods
Reagents and Fish
Deltamethrin (purity > 98%) was obtained from Sigma Aldrich (Chemical Co., USA), then dissolved in dimethyl sulfoxide (DMSO) as a stock solution. All healthy juvenile channel catfish (100 ± 20 g) were purchased in the Wuhan Baishazhou aquatic market. Before the experiment, these fish were domesticated for one week under laboratory conditions free of specific pathogens at 25 ± 1°C.
Experimental Design and Samples Collection
Ninety healthy channel catfish were randomly selected and divided into three groups, with three replicates in each group and 10 fish in each replicate. According to the pre-test results about lethal concentration of deltamethrin on channel catfish and the residual amount of Del in the nature water environment, 0.5 (L group) and (H group) 5 μg/L of Del were selected to study its toxic effect on channel catfish. The control group (C group) was soaked with DMSO. The final concentration of DMSO in experimental and control groups was under 0.001% (v/v). Experimental fish were not fed throughout the exposure period. After 6 h of exposure, 15 channel catfish were selected from each group to collect intestinal contents. Five intestinal contents were mixed to form a sample. The intestinal contents were squeezed into sterile centrifuge tubes, immediately frozen in liquid nitrogen, then stored at - 80°C for DNA extraction. Meanwhile, trunk kidneys were taken from 15 randomly selected channel catfish in each group, and five trunk kidneys were mixed to form a sample. Then, samples were immediately frozen in liquid nitrogen and stored at - 80°C for further transcriptomics. Additionally, the trunk kidneys and intestines of channel catfish in each group were cut and put into 10% neutral formaldehyde to make paraffin sections and stained with Hematoxylin and Eosin (HE) to study the pathological damage of Del to the intestines and trunk kidneys of channel catfish.
RNA Extraction, Library Construction, and Sequencing
Nine trunk kidney samples (three samples per treatment) were selected for RNA sequencing (RNA-seq). After total RNA was extracted, enriched mRNA by Oligo(dT) beads were fragmented using fragmentation buffer and reverse transcripted into cDNA with random primers using the NEBNext® Ultra™RNA library prep kit for Illumina following the manufacturer’s instructions. For the first-strand cDNA synthesis, the fragmented and primed mRNA was reversed to cDNA in a 20 μL reaction using protoscript II reverse transcriptase at 25°C for 10 min, then at 42°C for 15 min, and at 70°C for 15 min. For the synthesis of the second strand cDNA, the second strand synthase mixture was added to the first strand synthesis reaction at 16°C for 1 h with 80 μL. Then, the cDNA fragment was purified using the QiaQuick PCR Extraction Kit (Qiagen, Germany), and the end repair and poly (a) were performed. Finally, the Illumina sequencing adapter was connected. The ligation products were selected by size using agarose gel electrophoresis, PCR-amplified, and sequenced using Illumina HiSeq2500 by Genedenovo Biotechnology Co, Ltd (Guangzhou, China). For RNA-seq of fish samples under sterile conditions, libraries were sequenced on the Illumina HiSeq Xten platform.
Bioinformatic Analyses
Adaptor sequences were trimmed by cutadapt with at least 30 nt of remaining length. Clean reads were mapped to the channel catfish reference genome using HISAT2 (41). Gene expression levels were estimated by the number of fragments per million fragments per thousand base transcripts (FPKM). DESeq2 was used for differential expression analysis. The threshold of differentially expressed genes (DEGs) was set as FDR < 0.05 and |log2 [fold change (FC)]| > 1 (42). The gene ontology (GO) enrichment analysis of DEGs was performed with the “GOseq” R package, while the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed with KOBAS software (43).
Quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR)
We used quantitative reverse transcription (qRT)-PCR to verify the expression level of DEGs identified by RNA-Seq. The primers used are shown in Table 1 (39, 44–50). Channel catfish EF-1 α was used as the internal parameter to homogenize the expression levels. The qRT-PCR was performed using quantum Studio 6 flex (life technologies, USA), and melting curves were generated at the end of the run to confirm the specificity. The 2−ΔΔCt method was used to calculate the relative levels of DEGs (51).
DNA Extraction, PCR Amplification, and Sequencing
Bacterial DNA was purified using the TIANamp Bacteria DNA kit (Tiangen Biotech Inc., Beijing, China), following the manufacturer’s protocols. A region encompassing the V3–V4 hypervariable regions of the 16S rRNA gene was amplified using the primers 341F (CCTACGGGNGGCWGCAG) and 806R (GGACTACHVGGGTATCTAAT) (52). The PCR amplification system comprehended 20 μL, including 10 μL premixed Taq (Takara), 1 μL of each primer (10 μM), 2 μL gDNA, and 6 μL ddH2O. The PCR reaction procedure comprehended a pre-incubation at 95°C for 3 min, 30 cycles of 95°C for 30 s, 50°C for 30 s, 72°C for 1 min, and finally 72°C for 7 min. The amplified PCR products were separated by agarose gel electrophoresis and recovered by gel cutting and quantified by QuantiFluorTM fluorimetry. Purified amplification products were mixed in equal quantities, the sequenced connector was joined, the sequencing library was constructed, and Illumina PE250 was used for sequencing.
Biodiversity Analysis
For microbiome data analysis, FASTP (https://github.com/OpenGene/fastp) was used to filter the original data containing adapters or low-quality readings, and paired-end clean reads were merged as raw tags using FLSAH (53), with a minimum overlap of 10 bp and mismatch error rates of 2%. Clean labels were searched according to the reference database, and the reference-based chimerism check was performed using the uchime algorithm (http://www.drive5.com/usearch/manual/uchime_algo.html). All chimeric labels were removed to finally obtain valid labels for further analysis. Effective tags were clustered into operational taxonomic units (OTUs) of ≥ 97% similarity using UPARSE (54). The tag sequence with the highest abundance was selected as the representative sequence within each cluster. A between-groups Venn analysis was performed using R (version 3.4.1) to identify unique and common OTUs. Representative sequences were classified into organisms through a naive Bayesian model using RDP based on the Silva database (https://www.arb-silva.de/). Chao1, Simpson, and other alpha diversity indices were calculated in QIIME. The OTU sparse curve and rank abundance curve were also drawn in QIIME. The microbial communities of different samples were evaluated by principal coordinate analysis (PCoA) and arithmetic mean group method (UPGMA) β Diversity analysis. Calculation and plot multivariate statistics were performed in R. The KEGG pathway analysis of OTU was inferred by picrust (55) and tax4fun (56).The microbiome phenotypes were classified using BugBase.
Statistical Analyses
Data are expressed as means ± standard deviations. SPSS 22 (American SPSS company) was used for data analyses. Independent sample t-test, one-way analysis of variance (ANOVA), and Dunnett test were used to detect differences between groups. A p < 0.05 was considered statistically significant.
Results
Pathological Changes of Channel Catfish in the Early Stage of Del Exposure
In the present study, we found that the trunk kidney and intestine of channel catfish in the H group presented different degrees of pathological changes, including a disordered arrangement of intestinal structure, disordered shedding of some intestinal villous epithelial cells, infiltration of inflammatory cells, and atrophy and degeneration of trunk kidney glomerulus. However, the trunk kidney glomerulus in controls (C group) and the L group was not significantly denatured, and the tissue was normal. Meanwhile, the intestine of L group fish showed edema and degeneration of some villous epithelial cells and infiltration of tissue inflammatory cells. On the other hand, the intestine of the C group was normal (Figure 1). The clinical symptoms of channel catfish in the C and L groups had no significant change in the activity state. Meanwhile, fish in the H group were belly up, unable to swim, floating, and spot bleeding on the body surface after Del exposure for 6 h.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/916100/fimmu-13-916100-HTML-r1/image_m/fimmu-13-916100-g001.jpg)
Figure 1 Pathological changes of the intestine and trunk kidney in different groups. Black arrows represent inflammatory cell infiltration; Red arrows represent edema and degeneration of intestinal chorionic epithelial cells; Green arrows represent the shedding of intestinal chorionic epithelial cells; Yellow arrows represent tubular atrophy and denaturation. C: 0 μg/L Del; L: 0.5 μg/L Del; H: 5 μg/L Del.
Transcriptomes
Nine cDNA libraries were constructed from channel catfish trunk kidneys and sequenced. The characteristics of these libraries are summarized in Table 2. After quality control of sequencing data, the number of clean reads of each library ranged from 40330662 to 56410298, and the Q30 value of each sequencing library was higher than 93.41%. Then, clean reads were mapped to the reference genome (http://ftp.ensembl.org/pub/release-105/fasta/ictalurus_punctatus/) and presented mapping rates ranging from 88.70 to 89.92% in different libraries.
Identification of DEGs
Based on the differential expression analysis of DESeq2, DEGs between controls (C group) and experimental groups (L and H) were identified and visualized using volcano maps (Figures 2A, B). A total of 32 DEGs were detected in the L group and 1837 were detected in the H group (Figure 2C; Tables S1, S2). Among these DEGs, 30 were upregulated and 2 were downregulated in the L group. Meanwhile, 693 DEGs were upregulated and 1144 were downregulated in the H group (Figure 2C). The differential gene expression map is shown in Figure 2D. Moreover, apoptosis-, immune-, and drug metabolism-related DEGs were also detected in the trunk kidney of channel catfish after exposure to Del (Tables S3–S5), which were mainly concentrated in the H group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/916100/fimmu-13-916100-HTML-r1/image_m/fimmu-13-916100-g002.jpg)
Figure 2 Summary of differentially expressed genes (DEGs) in channel catfish after exposure to 0.5 and 5 μg/L Del. Volcano plots of DEGs in the trunk kidney treated with (A) 0.5 and (B) 5 μg/L Del. Red and blue dots indicate upregulated and downregulated genes, respectively. (C) Statistical results of DEGs in the C (control) vs L (0.5 μg/L), and C (control) vs H (5 μg/L) group comparisons; (D) Heat map of gene expression differences.
Next, to better characterize the biological function of DEGs, GO and KEGG pathway enrichment analyses were performed between control and experimental groups. The DEGs in L and H groups were successfully classified into three major groups, including biological processes (BP), cellular components (CC), and molecular functions (MF) (Figure 3). The most abundant GO terms in the BP category in the L group were oxidation-reduction process (GO: 0055114), followed by organic hydroxyl compound transport (GO: 0015850) and cholesteric metallic process (GO: 0008203); in the CC category, only four GO terms were detected: plasma lipoprotein particle (GO: 0034358), lipoprotein particle (GO: 1990777), protein-lipid complex (GO: 0032994), and high-density lipoprotein particle (GO: 0034364). In the MF category, the most abundant GO terms were oxidation activity (GO: 0016491) and coenzyme binding (GO: 0050662) (Table S6). Only one CC-related term was significantly enriched in the H group, but many terms were significantly enriched in both BP and MF. The GO term significantly enriched in the CC category was MCM complex (GO: 0042555), and the most abundant GO terms in the BP category were negative regulation of biological process (GO: 0048519), followed by negative regulation of cellular process (GO: 0048523), response to organic substances (GO: 0010033), cellular response to organic substances (GO: 0071310) and DNA metallic process (GO: 0006259). In the MF category, the most abundant GO terms were DNA binding (GO: 0003677) and DNA binding transcription factor activity (GO: 0003700) (Table S6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/916100/fimmu-13-916100-HTML-r1/image_m/fimmu-13-916100-g003.jpg)
Figure 3 Gene ontology (GO) annotation of differentially expressed genes (DEGs) in the three main GO categories: biological process (BP), cellular component (CC), and molecular function (MF). C: 0 μg/L, L: 0.5 μg/L and H: 5 μg/L.
The DEGs were also mapped to the KEGG database to analyze their biological functions and important pathways based on the whole transcriptome background. The KEGG pathway database is mainly divided into five categories: metabolism, genetic information processing, environmental information processing, cell process, and biological system. The KEGG pathway enrichment analysis showed that 15 and 175 pathways were enriched in L and H groups, respectively (Table S7, 8). The most significantly enriched pathway in the L group L was Autophagy-animal (KO04140, 3 DEGs) (Figure 4A), while in the H group was DNA replication (KO03030, 16 DEGs) (Figure 4B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/916100/fimmu-13-916100-HTML-r1/image_m/fimmu-13-916100-g004.jpg)
Figure 4 Top 20 KEGG enriched DEGs after exposure to (A) 0.5 and (B) 5 μg/L Del compared to the control group (0 μg/L Del). The vertical and horizontal axes represent different pathways and rich factors, respectively.
Verification With qPCR
To verify the RNA-seq results, the same RNA samples used for sequencing database construction were used to detect the expression of 18 DEGs by qRT-PCR. The primer sequences of all detected DEGs are shown in Table 1. The qRT-PCR results of the 18 DEGs were consistent with the RNA-seq results. These results indicated that the results of RNA-seq were generally reliable (Figure 5).
Sequencing Analysis and Taxonomic Annotation
A total of 1114159 effective tags were obtained from all intestinal samples of channel catfish. The effective ratio was 92.93-97.38% and the average sequence length was 303 to 478 (Tables S9, 10). The sequence sparse curve shows that the sample was sufficient to reflect the community richness and could be used for further data analysis (Figure 6A). The PCoA showed that the consistency of the sample was relatively good (Figure 6B). The α diversity index of intestinal microbiota after exposure to different concentrations of Del is shown in Table S11 and did not differ (p > 0.05). Additionally, a total of 1488 OTUs were detected using the 97% sequence similarity threshold. There were 336 common OTUs between C, L, and H groups. On the other hand, compared with the C group, 314 unique OTUs were detected in the L group and 493 in the H group (Figure 6C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/916100/fimmu-13-916100-HTML-r1/image_m/fimmu-13-916100-g006.jpg)
Figure 6 Rarefaction curves of channel catfish intestinal microbial samples (A); Principal component analysis of samples from two groups (B), n = 3; Venn diagram of the OTU distribution (C). C: 0 μg/L Del; L: 0.5 μg/L Del; H: 5 μg/L Del.
Intestinal Microbial Composition
Among all intestinal samples, Firmicutes was the main phylum in the top 10, followed by Cyanobacteria, Bacteroidetes, Proteobacteria, Verrucomicrobia, and Fusobacteria (Figure 7A). Compared with controls, the abundance of Firmicutes, Cyanobacteria, and Fusobacteria decreased in the L group, while the abundance of Bacteroides, Proteobacteria, and Verrucomicrobia increased. In the H group, the abundances of Firmicutes and Fusobacteria seriously decreased, the abundances of Proteobacteria and Verrucomicrobia slightly decreased, while the abundances of Cyanobacteria and Bacteroidetes significantly increased (Table S10). Moreover, the abundances of Firmicutes, Fusobacteria, Actinobacteria, and Epsilonbacteraeota decreased with increasing Del concentrations. At the genus level, the dominant bacteria in the gut of channel catfish included 10 species, such as Romboutsia, Sediminibacterium, Ralstonia, Akkermansia, and Lactobacillus (Figure 7B). After Del exposure, the abundances of Romboutsia, Lactobacillus, and Cetobacterium decreased (Table S11). The functional prediction of the bacterial structure showed that the intestinal tract of channel catfish in L and H groups increased the relative abundance of Gram-negative bacteria and decreased the relative abundance of Gram-positive bacteria (Figure 8A). Additionally, Del treatment reduced the abundance of anaerobic bacteria in L and H groups. The abundance of aerobic bacteria increased in the L group, while decreased in the H group. The phenotype prediction of the microbiome indicated that these alterations in the bacterial diversity might be involved in different metabolisms as the dominant function (Figure 8B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/916100/fimmu-13-916100-HTML-r1/image_m/fimmu-13-916100-g007.jpg)
Figure 7 Relative abundances of dominant bacterial phyla (A) and genera (B) in the intestine of channel catfish after exposure to different concentrations of Del. C: 0 μg/L Del; L: 0.5 μg/L Del; H: 5 μg/L Del.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/916100/fimmu-13-916100-HTML-r1/image_m/fimmu-13-916100-g008.jpg)
Figure 8 Organism level microbiome phenotypes were predicted by Bugbase (A). Heat map of functional prediction by PICRUSt (B).
Discussion
Large-scale diseases can occur from time to time, resulting in significant economic losses in channel catfish culture (39, 40). Del is not only an important pollutant in water but also an insecticide used in aquaculture, becoming an important factor affecting their healthy cultivation. Therefore, clarifying the toxic mechanisms of Del in channel catfish can provide a theoretical basis for guiding their cultivation and avoiding water pollution. Hence, our current research has important theoretical and practical significance. In the present study, exposure to different concentrations of Del for 6 h led to pathological damage to the intestinal tract and trunk kidney of channel catfish. In the H group, the intestinal villous epithelial cells were disordered and shed, and inflammatory cell infiltration and glomerular atrophy and degeneration were the most serious characteristics, consistent with previous studies (57). The pathological results not only confirmed the toxicity of Del with high concentration to channel catfish but also demonstrated that the immersion experiment was successful, which laid a foundation for subsequent transcriptome and intestinal microbial composition analyses.
Then, we sequenced the transcriptome of the trunk kidney of channel catfish after Del exposure for 6 h. The results showed that the early gene expression of channel catfish was significantly altered after exposure to Del at different concentrations. Only 32 DEGs were enriched in the L group, while 1837 DEGs were enriched in the H group. The enriched DEGs could be divided into three categories: apoptosis, immunity, and metabolism. The L group was mainly enriched in metabolism-related DEGs, such as adobe and p4ha1b, and only two DEGs, such as lgi2a, were significantly downregulated. In the H group, 693 DEGs were upregulated and 1144 DEGs were downregulated, covering the three categories (apoptosis, immunity, and metabolism). The apoptosis-related DEGs enriched in the H group included Ctsl, Gadd45b, Gadd45aa, Caspase-3, Bax, Smad1, Smad6, Smad7, Bcl-2, and TNF superfamily receptors. Moreover, the enriched metabolism-related DEGs included cytochrome P450 and the immune-related ones included the toll-like receptor, chemokine family genes, and IL-1β. These results showed that the low concentration of Del (0.5μg/L) had little effect on channel catfish, but 5μg/L induced a great change in the gene transcription level of the trunk kidney. Hence, the exposure to a high concentration of Del had a strong stimulation on channel catfish in the early stage.
Apoptosis or programmed cell death is considered a key component of homeostasis maintenance and can be induced by various environmental stressors in aquatic organisms, including chemicals (58). Here, we found that channel catfish after 5μg/L Del exposure presented enriched apoptosis-related GO terms, such as tumor necrosis factor receptor binding, tumor necrosis factor receptor superfamily binding, and enriched apoptosis and p53 signaling pathways. Meanwhile, after 0.5μg/L Del exposure the related GO terms and signaling pathways were rarely enriched. Previous studies have shown that Del can induce several types of apoptosis by upregulating the expression of p53 and Bax and inhibiting the expression of Bcl-2 (59, 60). In the current study, p53 and Bax were significantly upregulated, and Bcl-2 was significantly downregulated in the trunk kidney of channel catfish in the H group. Additionally, Del exposure can induce Caspase-3-dependent apoptosis in carp (28). Caspase-3 was also significantly upregulated in the H group. Some studies have also pointed out that cathepsin L (Ctsl) is an important participant in various apoptosis-related pathways by promoting the expression of cytochrome c and bid, and inhibiting the expression of the anti-apoptotic factor Bcl-2 from activating the apoptosis pathway (61, 62). We also found that the Ctsl gene was significantly upregulated but failed to enrich cytochrome c. Moreover, the expression difference of bid was not significant, which might be the problem with the depth of transcriptome sequencing (63, 64). The TGF-β signaling pathway was also enriched in the present study, and the expressions of Smad1, Smad6, and Smad7 were significantly upregulated, which might indicate channel catfish renal fibrosis (65). These results suggested that in the early stage of acute exposure to high concentrations of Del, apoptosis occurred in the trunk kidney of channel catfish, and fibrosis might occur further.
The metabolism comprehends a series of orderly chemical reactions that occur in organisms to maintain life. These reactions enable organisms to grow and reproduce, maintain their structure and respond to the environment. In our research, many metabolism-related DEGs were enriched in channel catfish after Del exposure. Cytochrome P450 (CYP) is involved in the metabolism of exogenous substances (66), especially chemical drugs, and is considered to be a good biomarker to evaluate pollutants in water environments (67). Previous studies have shown that Del exposure can lead to the upregulation of CYP1A1 expression in tuna and rainbow trout (68, 69), and it is time-dependent in rainbow trout (69). Similarly, in the present study, after exposure to 5 μg/L Del, CYP1A1 was also significantly upregulated in channel catfish. CYP1A1 was mainly involved in metabolic pathways and metabolism of xenobiology by cytochrome P450, indicating that CYP1A1 was involved in the metabolism of Del in the trunk kidney of channel catfish. To resist external stimuli, such as pathogenic microorganisms and pollutants, and maintain tissue homeostasis, organisms need to invest a lot of energy, which can disorder energy metabolism (70–72). In the present study, some pathways related to energy metabolism were enriched in channel catfish at the early stage of Del exposure, including lipid metabolism, carbohydrate metabolism, amino acid metabolism, startch and sucrose metabolism, glycolysis/gluconeogenesis, pentose and gluconate interconversion and fructose and mannose metabolism. The results showed that channel catfish provided the required energy through carbohydrate, fatty acid, and amino acid metabolisms to resist the requirements to compensate other functional damage caused by Del toxicity. It also explains that the growth rate of fish is very slow under long-term stress such as pollutants, consistent with previous studies (47, 73). We also found that the expression of myosin, which regulates muscle strength (74), was significantly downregulated, that is, the muscle strength and exercise of channel catfish decreased, consistent with their clinical symptoms. However, we also found that the expression of glucose-6-phosphatase (G6P) was significantly downregulated. G6P can dephosphorylate glucose to promote circulation (47), which means that the energy metabolism of channel catfish was disordered in the early stage of Del acute exposure. However, Hsp70, as a stress protein involved in cell protection and repair (75), was significantly up-regulated. At the same time, it is also enriched in autophagy - animal and mitophagy - animal signal pathways (76, 77), indicating that channel catfish produces an acute stress response under high concentration Del exposure, so as to strengthen cell repair and protection, maintain body homeostasis and improve survival ability.
The fish immune system is very important to the host’s defense, comprehending an important barrier for the body to resist external stimuli such as pathogenic microbial infections (78). Meanwhile, it is also affected by various environmental pollutants, including drugs (79, 80). We showed that in the early stage of acute exposure to a high concentration of Del, the transcriptome of channel catfish trunk kidney was enriched for many immune-related DEGs, including toll-like receptor, chemokine family, and IL-1β. It was also enriched in many immune-related signaling pathways, such as the Toll-like receptor, RIG-I-like receptor, and C-type lectin receptor signaling pathways. The TLR family senses the molecular signatures of microbial pathogens and plays a fundamental role in innate immune responses (81). The expressions of TLR9, TLR21, TLR18, TLR1, TLR5, and TLR19 in channel catfish were significantly downregulated in the early stage of Del acute exposure, which led to the decline of the body’s ability to deal with pathogenic microorganisms. Additionally, alkaline phosphatase (ALP), chemokines, nucleotide-binding oligomeric domain (NOD), and integrin play an important role in innate immunity by promoting inflammation and antibacterial responses and recruiting leukocyte phagocytosis (48, 82–86). However, we found that the expressions of alkaline phosphatase (ALP), CCL34, CXCR5, NOD1, NLRC3 and integrin alpha in were significantly downregulated channel catfish after Del exposure. Along with the enriched immune-related signaling pathways, these results indicated that the immune system was seriously affected by Del exposure. Further, we showed that Del could inhibit the immune defense system of channel catfish, especially the innate immunity. Similar results were detected for snakehead, Chinese rare minnow, and goldfish (12, 24, 31, 87). Therefore, Del can reduce fish immunity and make it more vulnerable to pathogen infections, consistent with the result that copper sulfate stimulation can make channel catfish more vulnerable to Edwardsiella ictalurid (88).
Surprisingly, in the early stage of high concentration Del exposure, C3, iNOS, and IL-1β were significantly upregulated. IL-1β is a powerful proinflammatory factor (89) and iNOS is an important index to measure inflammation (28, 90). Their upregulation expression can represent the body’s inflammation to protect itself in response to Del toxicity. However, the pathological results showed that there was no clear inflammation in the trunk kidney, which might be related to the fact that fish were still in the early stage of exposure and the inflammation did not result in pathological changes yet. Additionally, the complement system can mediate host immune and inflammatory responses and play a key role in preventing pathogen infections (91), indicating that under Del toxicity, the complement system can enhance the innate immunity of channel catfish through compensation, similar to gilthead seabream (Sparus aurata) (68). However, after subacute exposure to Del, the levels of C3, C4, and other immune parameters significantly decreased in snakeheads (24). A potential explanation for this difference might be that the effects of Del on aquatic animals are species and dose-dependent (92), that is, Del can present considerable toxicity differences among different species (24, 93, 94). Overall, in the present study, the conflicting results regarding the transcription level of immune-related genes are likely to reflect the complexity of the immune defense response of channel catfish after acute exposure to high concentrations of Del.
Under normal circumstances, beneficial and harmful bacteria are in a dynamic balance in the intestine (95). Here, in the early stage of Del exposure, significant changes were detected in the composition of the intestinal microbiota of channel catfish, with a significant decrease in Gram-positive bacteria and a significant increase in Gram-negative bacteria. The intestinal microbiome composition analysis of channel catfish showed that the dominant bacteria at the phylum level were Firmicutes, Cyanobacteria, Bacteroides, Proteobacteria, Verrucomicrobia, Fusobacteria, and Actinobacteria. Previous studies have also reported similar results, showing that the intestinal flora of fish is mostly composed of Proteobacteria, Firmicutes, Fusobacteria, and Actinobacteria (96, 97).
Furthermore, Del exposure significantly reduced the abundance of Firmicutes, Fusobacteria, and Actinobacteria, and the Firmicutes/Bacteroidetes ratio also decreased. Bacteroides and Firmicutes play an important role in energy metabolism in the fish intestine. Bacteroides can encode carbohydrate-related enzymes involved in the hydrolysis of glycoconjugates. Firmicutes promote the intestinal absorption and metabolism of fatty acids and play an active role in the growth performance, immunity, digestion, and disease resistance of aquatic animals (98–101). The microbiome phenotype prediction showed that changing bacterial species might participate in many metabolisms as the dominant function. The changes in the number of Bacteroides and Firmicutes caused by Del exposure might also lead to energy metabolism disorders in channel catfish. The transcriptome data showed that channel catfish provided the required energy to resist Del toxicity through carbohydrate, fatty acid, and amino acid metabolisms. The changes in intestinal microbial composition showed that they mainly rely on nutrients in vivo to make up for the massive consumption of energy, rather than exogenous intake. We also found that the abundance of Fusobacteria and Actinobacteria significantly decreased, which can increase the risk of infections (102, 103). Hence, Del exposure in the early stage can make the body more vulnerable to pathogen infections and cause metabolic disorders.
At the genus level, there were 10 dominant bacteria, including Romboutsia, Sediminibacterium, Ralstonia, Akkermansia, Lactobacillus, Olleya, and Cetobacterium. Among them, the abundance of Romboutsia, Lactobacillus, and Cetobacterium significantly decreased with increasing Del concentrations. However, Lactobacillus and Cetobacterium are important beneficial bacteria in the intestine (104). For example, Lactobacillus can reduce the intestinal pH, which might contribute to overcoming pathogen challenges (105). Cetobacterium might be involved in the metabolism of the host and provide vitamin B12 (106). Therefore, the intestinal function of channel catfish might also be affected in the early stage of Del exposure. These results suggested that Del might adversely affect the metabolism and immune function of channel catfish by changing the intestinal microbial community in the early stage of acute exposure.
In summary, the pathological results showed that a high concentration of Del can damage the intestine and trunk kidney of channel catfish in the early stage. Through transcriptome analysis, we detected several DEGs related to immunity, metabolism, and apoptosis in the trunk kidney of channel catfish, indicating that Del exposure in the early stage can also lead to immunotoxicity, metabolic disorders, and tissue damage. Additionally, the 16S rRNA gene sequencing showed that exposure to Del can significantly change the composition of the intestinal microbial community, interfering with the normal metabolism and immune function of the intestine of channel catfish. Overall, we analyzed the adverse effects of Del on the early exposure of channel catfish and showed why it can trigger infectious diseases and slow growth in aquaculture.
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 below: https://www.ncbi.nlm.nih.gov/genbank/, PRJNA820141 https://www.ncbi.nlm.nih.gov/genbank/, PRJNA820136.
Ethics Statement
The animal study was reviewed and approved by the Animal Welfare and Research Ethics Committee of Yangtze River Fisheries Research Institute (YFI2021YYB002). Written informed consent was obtained from the owners for the participation of their animals in this study.
Author Contributions
YY and XA conceived and designed the study. YY and XZ performed most of the experiments. XA, XZ, YH, HZ, NX, YL, and GF provided bioinformatics assistance and support. HZ and YY wrote the manuscript. All authors proofed approved the manuscript.
Funding
This work was supported by Central Public-Interest Scientific Institution Basal Research Fund CAFS(2020TD11), Central Public-interest Scientific Institution Basal Research Fund (NO.YFI202209) and Yancheng fishery high quality development project (No.YCSCYJ2021026).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We sincerely thank anonymous reviewers for their constructive comments which lead to an improved manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.916100/full#supplementary-material
References
1. Velíšek J, Jurčikova J, Dobšikova R, Svobodova Z, Piačkova V, Máchova J, et al. Effects of Deltamethrin on Rainbow Trout (Oncorhynchus Mykiss). Environ Toxicol Pharmacol (2007) 23(3):297–301. doi: 10.1016/j.etap.2006.11.006
2. Milam CD, Farris JL, Wilhide JD. Evaluating Mosquito Control Pesticides for Effect on Target and Nontarget Organisms. Arch Environ Contam Toxicol (2000) 39(3):324–8. doi: 10.1007/s002440010111
3. Parvez S, Raisuddin S. Copper Modulates non-Enzymatic Antioxidants in the Freshwater Fish Channa Punctata (Bloch) Exposed to Deltamethrin. Chemosphere (2006) 62(8):1324–32. doi: 10.1016/j.chemosphere.2005.07.025
4. Prasamthi K, Muralidhara and PS. Rajini: Fenvalerate-Induced Oxidative Damage in Rat Tissues and Its Attenuation by Dietary Sesame Oil. Food Chem Toxicol (2005) 43(2):299–306. doi: 10.1016/j.fct.2004.10.005
5. Glickman AH, Weitman SD, Lech JJ. Differential Toxicity of Trans-Permethrin in Rainbow Trout and Mice. I. Role of Biotransformation. Toxicol Appl Pharmacol (1982) 66(2):153–61. doi: 10.1016/0041-008x(82)90280-0
6. Zhao M, Liu W. Enantioselectivity in the Immunotoxicity of the Insecticide Acetofenate in an In Vitro Model. Environ Toxicol Chem (2009) 28(3):578–85. doi: 10.1897/08-246.1
7. Clasen B, Loro VL, Murussi CR, Tiecher TL, Moraes B, Zanella R. Bioaccumulation and Oxidative Stress Caused by Pesticides in Cyprinus Carpio Reared in a Rice-Fish System. Sci Total Environ (2018) 626:737–43. doi: 10.1016/j.scitotenv.2018.01.154
8. Bradbury SP, Coats JR. Comparative Toxicology of the Pyrethroid Insecticides. Rev Environ Contam Toxicol (1989) 108:133–77. doi: 10.1007/978-1-4613-8850-0_4
9. Casida JE, Gammon DW, Glickman AH, Lawrence LJ. Mechanisms of Selective Action of Pyrethroid Insecticides. Annu Rev Pharmacol Toxicol (1983) 23:413–38. doi: 10.1146/annurev.pa.23.040183.002213
10. Edwards R, Millburn P, Hutson DH. Comparative Toxicity of Cis-Cypermethrin in Rainbow Trout, Frog, Mouse, and Quail. Toxicol Appl Pharmacol (1986) 84(3):512–22. doi: 10.1016/0041-008x(86)90256-5
11. Paul EA, Simonin HA. Toxicity of Three Mosquito Insecticides to Crayfish. Bull Environ Contam Toxicol (2006) 76(4):614–21. doi: 10.1007/s00128-006-0964-4
12. Wu H, Gao J, Xie M, Xiang J, Zuo Z, Tian X, et al. Histopathology and Transcriptome Analysis Reveals the Gills Injury and Immunotoxicity in Gibel Carp Following Acute Deltamethrin Exposure. Ecotoxicol Environ Saf (2022) 234:113421. doi: 10.1016/j.ecoenv.2022.113421
13. Ensminger M, Bergin R, Spurlock F, Goh KS. Pesticide Concentrations in Water and Sediment and Associated Invertebrate Toxicity in Del Puerto and Orestimba Creeks, California, 2007-2008. Environ Monit Assess (2011) 175(1-4):573–87. doi: 10.1007/s10661-010-1552-y
14. Thomas CR, Hose GC, Warne MS, Lim RP. Effects of River Water and Salinity on the Toxicity of Deltamethrin to Freshwater Shrimp, Cladoceran, and Fish. Arch Environ Contam Toxicol (2008) 55(4):610–8. doi: 10.1007/s00244-008-9147-0
15. Feo ML, Ginebreda A, Eljarrat E, Hydrology D. Presence of Pyrethroid Pesticides in Water and Sediments of Ebro River Delta. J Hydrol (2010) 393(3-4):156–62. doi: 10.1016/j.jhydrol.2010.08.012
16. Amin AS, Moalla S, Ali A, Gouda AA. Sensitive Spectrophotometric Determination of Deltamethrin Insecticide in its Formulation and Environmental Samples. IJARCS (2015) 10(2):72–9. doi: 10.20431/2349-0403.0210009
17. Bille L, Binato G, Gabrieli C, Manfrin A, Pascoli F, Pretto T, et al. First Report of a Fish Kill Episode Caused by Pyrethroids in Italian Freshwater. Forensic Sci Int (2017) 281:176–82. doi: 10.1016/j.forsciint.2017.10.040
18. Zhang C, Zhang Q, Pang Y, Song X, Zhou N, Wang J, et al. The Protective Effects of Melatonin on Oxidative Damage and the Immune System of the Chinese Mitten Crab (Eriocheir Sinensis) Exposed to Deltamethrin. Sci Total Environ (2019) 653:1426–34. doi: 10.1016/j.scitotenv.2018.11.063
19. Zhang L, Zhao X, Yan S, Zha J, Ma X. The Immune Responses of the Chinese Rare Minnow (Gobiocypris Rarus) Exposed to Environmentally Relevant Concentrations of Cypermethrin and Subsequently Infected by the Bacteria Pseudomonas Fluorescens. Environ Pollut (2019) 250:990–7. doi: 10.1016/j.envpol.2019.03.126
20. Hong Y, Huang Y, Yan G, Huang Z. Effects of Deltamethrin on the Antioxidant Defense and Heat Shock Protein Expression in Chinese Mitten Crab, Eriocheir Sinensis. Environ Toxicol Pharmacol (2019) 66:1–6. doi: 10.1016/j.etap.2018.12.012
21. Tu HT, Silvestre F, Meulder BD, Thome JP, Phuong NT, Kestemont P. Combined Effects of Deltamethrin, Temperature and Salinity on Oxidative Stress Biomarkers and Acetylcholinesterase Activity in the Black Tiger Shrimp (Penaeus Monodon). Chemosphere (2012) 86(1):83–91. doi: 10.1016/j.chemosphere.2011.09.022
22. Awoyemi OM, Kumar N, Schmitt C, Subbiah S, Crago J. Behavioral, Molecular and Physiological Responses of Embryo-Larval Zebrafish Exposed to Types I and II Pyrethroids. Chemosphere (2019) 219:526–37. doi: 10.1016/j.chemosphere.2018.12.026
23. Abdel-Daim MM, Abdelkhalek NK, Hassan AM. Antagonistic Activity of Dietary Allicin Against Deltamethrin-Induced Oxidative Damage in Freshwater Nile Tilapia; Oreochromis Niloticus. Ecotoxicol Environ Saf (2015) 111:146–52. doi: 10.1016/j.ecoenv.2014.10.019
24. Kong Y, Li M, Shan X, Wang G, Han G. Effects of Deltamethrin Subacute Exposure in Snakehead Fish, Channa Argus: Biochemicals, Antioxidants and Immune Responses. Ecotoxicol Environ Saf (2021) 209:111821. doi: 10.1016/j.ecoenv.2020.111821
25. Antwi FB, Reddy GV. Toxicological Effects of Pyrethroids on non-Target Aquatic Insects. Environ Toxicol Pharmacol (2015) 40(3):915–23. doi: 10.1016/j.etap.2015.09.023
26. Hong Y, Huang Y, Yan G, Yin H, Huang Z. DNA Damage, Immunotoxicity, and Neurotoxicity Induced by Deltamethrin on the Freshwater Crayfish, Procambarus clarkii. Environ Toxicol (2020) 36:16–23. doi: 10.1002/tox.23006
27. Liu X, Gao Q, Feng Z, Tang Y, Zhao X, Chen D, et al. Protective Effects of Spermidine and Melatonin on Deltamethrin-Induced Cardiotoxicity and Neurotoxicity in Zebrafish. Cardiovasc Toxicol (2021) 21(1):29–41. doi: 10.1007/s12012-020-09591-5
28. Arslan H, Altun S, Özdemir S. Acute Toxication of Deltamethrin Results in Activation of iNOS, 8-OHdG and Up-Regulation of Caspase 3, iNOS Gene Expression in Common Carp (Cyprinus Carpio L.). Aquat Toxicol (2017) 187:90–9. doi: 10.1016/j.aquatox.2017.03.014
29. Li M, Liu X, Feng X. Cardiovascular Toxicity and Anxiety-Like Behavior Induced by Deltamethrin in Zebrafish (Danio Rerio) Larvae. Chemosphere (2019) 219:155–64. doi: 10.1016/j.chemosphere.2018.12.011
30. Liu X, Zhang Q, Li S, Mi P, Chen D, Zhao X, et al. Developmental Toxicity and Neurotoxicity of Synthetic Organic Insecticides in Zebrafish (Danio Rerio): A Comparative Study of Deltamethrin, Acephate, and Thiamethoxam. Chemosphere (2018) 199:16–25. doi: 10.1016/j.chemosphere.2018.01.176
31. Zhang L, Hong X, Zhao X, Yan S, Ma X, Zha J. Exposure to Environmentally Relevant Concentrations of Deltamethrin Renders the Chinese Rare Minnow (Gobiocypris Rarus) Vulnerable to Pseudomonas Fluorescens Infection. Sci Total Environ (2020) 715:136943. doi: 10.1016/j.scitotenv.2020.136943
32. Xiong G, Deng Y, Li J, Cao Z, Liao X, Liu Y, et al. Immunotoxicity and Transcriptome Analysis of Zebrafish Embryos in Response to Glufosinate-Ammonium Exposure. Chemosphere (2019) 236:124423. doi: 10.1016/j.chemosphere.2019.124423
33. Hrdlickova R, Toloue M, Tian B. RNA-Seq Methods for Transcriptome Analysis. Wiley Interdiscip Rev RNA (2017) 8(1):e1364. doi: 10.1002/wrna.1364
34. Atamaniuk TM, Kubrak OI, Husak VV, Storey KB, Lushchak VI. The Mancozeb-Containing Carbamate Fungicide Tattoo Induces Mild Oxidative Stress in Goldfish Brain, Liver, and Kidney. Environ Toxicol (2014) 29(11):1227–35. doi: 10.1002/tox.21853
35. Stosik MP, Tokarz-Deptuła B, Deptuła W. Specific Humoral Immunity in Osteichthyes. Cent Eur J Immunol (2018) 43(3):335–40. doi: 10.5114/ceji.2018.80054
36. Pérez T, Balcázar JL, Ruiz-Zarzuela I, Halaihel N, Vendrell D, de Blas I, et al. Host-Microbiota Interactions Within the Fish Intestinal Ecosystem. Mucosal Immunol (2010) 3(4):355–60. doi: 10.1038/mi.2010.12
37. Buerger AN, Dillon DT, Schmidt J, Yang T, Zubcevic J, Martyniuk CJ, et al. Gastrointestinal Dysbiosis Following Diethylhexyl Phthalate Exposure in Zebrafish (Danio Rerio): Altered Microbial Diversity, Functionality, and Network Connectivity. Environ Pollut (2020) 265(Pt B):114496. doi: 10.1016/j.envpol.2020.114496
38. Zhang Y, Sun K, Li Z, Chai X, Fu X, Kholodkevich S, et al. Effescts of Acute Diclofenac Exposure on Intestinal Histology, Antioxidant Defense, and Microbiota in Freshwater Crayfish (Procambarus clarkii). Chemosphere (2021) 263:128130. doi: 10.1016/j.chemosphere.2020.128130
39. Yang Y, Zhu X, Zhang H, Chen Y, Song Y, Ai X. Dual RNA-Seq of Trunk Kidneys Extracted From Channel Catfish Infected With Yersinia Ruckeri Reveals Novel Insights Into Host-Pathogen Interactions. Front Immunol (2021) 12:775708. doi: 10.3389/fimmu.2021.775708
40. Liu T, Wang KY, Wang J, Chen DF, Huang XL, Ouyang P, et al. Genome Sequence of the Fish Pathogen Yersinia Ruckeri SC09 Provides Insights Into Niche Adaptation and Pathogenic Mechanism. Int J Mol Sci (2016) 17(4):557. doi: 10.3390/ijms17040557
41. Kim D, Langmead B, Salzberg SL. HISAT: A Fast Spliced Aligner With Low Memory Requirements. Nat Methods (2015) 12(4):357–60. doi: 10.1038/nmeth.3317
42. Love MI, Huber W, Anders S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol (2014) 15(12):550. doi: 10.1186/s13059-014-0550-8
43. Mao X, Cai T, Olyarchuk JG, Wei L. Automated Genome Annotation and Pathway Identification Using the KEGG Orthology (KO) as a Controlled Vocabulary. Bioinformatics (2005) 21(19):3787–93. doi: 10.1093/bioinformatics/bti430
44. Xu DH, Zhang D, Shoemaker C, Beck B. Immune Response of Channel Catfish (Ictalurus Punctatus) Against Ichthyophthirius Multifiliis Post Vaccination Using DNA Vaccines Encoding Immobilization Antigens. Fish Shellfish Immunol (2019) 94:308–17. doi: 10.1016/j.fsi.2019.08.071
45. Ren Y, Zhao H, Su B, Peatman E, Li C. Expression Profiling Analysis of Immune-Related Genes in Channel Catfish (Ictalurus Punctatus) Skin Mucus Following Flavobacterium Columnare Challenge. Fish Shellfish Immunol (2015) 46(2):537–42. doi: 10.1016/j.fsi.2015.07.021
46. Li C, Beck BH, Peatman E. Nutritional Impacts on Gene Expression in the Surface Mucosa of Blue Catfish (Ictalurus Furcatus). Dev Comp Immunol (2014) 44(1):226–34. doi: 10.1016/j.dci.2013.12.014
47. Paul JS, Small BC. Chronic Exposure to Environmental Cadmium Affects Growth and Survival, Cellular Stress, and Glucose Metabolism in Juvenile Channel Catfish (Ictalurus Punctatus). Aquat Toxicol (2021) 230:105705. doi: 10.1016/j.aquatox.2020.105705
48. Li M, Wang QL, Lu Y, Chen SL, Li Q, Sha ZX. Expression Profiles of NODs in Channel Catfish (Ictalurus Punctatus) After Infection With EdwardsiellaTarda, Aeromonas Hydrophila, Streptococcus Iniae and Channel Catfish Hemorrhage Reovirus. Fish Shellfish Immunol (2012) 33(4):1033–41. doi: 10.1016/j.fsi.2012.06.033
49. Weber TE, Bosworth BG. Effects of 28 Day Exposure to Cold Temperature or Feed Restriction on Growth, Body Composition, and Expression of Genes Related to Muscle Growth and Metabolism in Channel Catfish. Aquaculture (2005) 246(1-4):483–92. doi: 10.1016/j.aquaculture.2005.02.032
50. Pridgeon JW, Russo R, Shoemaker CA, Klesius PH. Expression Profiles of Toll-Like Receptors in Anterior Kidney of Channel Catfish, Ictalurus Punctatus (Rafinesque), Acutely Infected by Edwardsiella Ictaluri. J Fish Dis (2010) 33(6):497–505. doi: 10.1111/j.1365-2761.2010.01159.x
51. Livak KJ, Schmittgen TD. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods (2001) 25(4):402–8. doi: 10.1006/meth.2001.1262
52. Zhou F, Liu B, Liu X, Li Y, Wang L, Huang J, et al. The Impact of Microbiome and Microbiota-Derived Sodium Butyrate on Drosophila Transcriptome and Metabolome Revealed by Multi-Omics Analysis. Metabolites (2021) 11(5):402–8. doi: 10.3390/metabo11050298
53. Magoč T, Salzberg SL. FLASH: Fast Length Adjustment of Short Reads to Improve Genome Assemblies. Bioinformatics (2011) 27(21):2957–63. doi: 10.1093/bioinformatics/btr507
54. Edgar RC. UPARSE: Highly Accurate OTU Sequences From Microbial Amplicon Reads. Nat Methods (2013) 10(10):996–8. doi: 10.1038/nmeth.2604
55. Langille MG, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, et al. Predictive Functional Profiling of Microbial Communities Using 16S rRNA Marker Gene Sequences. Nat Biotechnol (2013) 31(9):814–21. doi: 10.1038/nbt.2676
56. Aßhauer KP, Wemheuer B, Daniel R, Meinicke P. Tax4Fun: Predicting Functional Profiles From Metagenomic 16S rRNA Data. Bioinformatics (2015) 31(17):2882–4. doi: 10.1093/bioinformatics/btv287
57. Cengiz EI. Gill and Kidney Histopathology in the Freshwater Fish Cyprinus Carpio After Acute Exposure to Deltamethrin. Environ Toxicol Pharmacol (2006) 22(2):200–4. doi: 10.1016/j.etap.2006.03.006
58. AnvariFar H, Amirkolaie AK, Miandare HK, Ouraji H, Jalali MA, Üçüncü S. Apoptosis in Fish: Environmental Factors and Programmed Cell Death. Cell Tissue Res (2017) 368(3):425–39. doi: 10.1007/s00441-016-2548-x
59. Wu A, Liu Y. Apoptotic Cell Death in Rat Brain Following Deltamethrin Treatment. Neurosci Lett (2000) 279(2):85–8. doi: 10.1016/s0304-3940(99)00973-8
60. Wu A, Li L, Liu Y. Deltamethrin Induces Apoptotic Cell Death in Cultured Cerebral Cortical Neurons. Toxicol Appl Pharmacol (2003) 187(1):50–7. doi: 10.1016/s0041-008x(02)00032-7
61. Hsu KF, Wu CL, Huang SC, Wu CM, Hsiao JR, Yo YT, et al. Cathepsin L Mediates Resveratrol-Induced Autophagy and Apoptotic Cell Death in Cervical Cancer Cells. Autophagy (2009) 5(4):451–60. doi: 10.4161/auto.5.4.7666
62. Droga-Mazovec G, Bojic L, Petelin A, Ivanova S, Romih R, Repnik U, et al. Cysteine Cathepsins Trigger Caspase-Dependent Cell Death Through Cleavage of Bid and Antiapoptotic Bcl-2 Homologues. J Biol Chem (2008) 283(27):19140–50. doi: 10.1074/jbc.M802513200
63. Sims D, Sudbery I, Ilott NE, Heger A, Ponting CP. Sequencing Depth and Coverage: Key Considerations in Genomic Analyses. Nat Rev Genet (2014) 15(2):121–32. doi: 10.1038/nrg3642
64. Mirsafian H, Ripen AM, Leong WM, Manaharan T, Mohamad SB, Merican AF. Transcriptome Landscape of Human Primary Monocytes at Different Sequencing Depth. Genomics (2017) 109:463–70. doi: 10.1016/j.ygeno.2017.07.003
65. Li X, Chu G, Zhu F, Zheng Z, Wang X, Zhang G, et al. Epoxyeicosatrienoic Acid Prevents Maladaptive Remodeling in Pressure Overload by Targeting Calcineurin/NFAT and Smad-7. Exp Cell Res (2020) 386(1):111716. doi: 10.1016/j.yexcr.2019.111716
66. Peter Guengerich F, Avadhani NG. Roles of Cytochrome P450 in Metabolism of Ethanol and Carcinogens. Adv Exp Med Biol (2018) 1032:15–35. doi: 10.1007/978-3-319-98788-0_2
67. Uno T, Ishizuka M, Itakura T. Cytochrome P450 (CYP) in Fish. Environ Toxicol Pharmacol (2012) 34(1):1–13. doi: 10.1016/j.etap.2012.02.004
68. Guardiola FA, Gónzalez-Párraga P, Meseguer J, Cuesta A, Esteban MA. Modulatory Effects of Deltamethrin-Exposure on the Immune Status, Metabolism and Oxidative Stress in Gilthead Seabream (Sparus Aurata L. ). Fish Shellfish Immunol (2014) 36(1):120–9. doi: 10.1016/j.fsi.2013.10.020
69. Erdoğan O, Ceyhun SB, Ekinci D, Aksakal E. Impact of Deltamethrin Exposure on mRNA Expression Levels of Metallothionein A, B and Cytochrome P450 1A in Rainbow Trout Muscles. Gene (2011) 484(1-2):13–7. doi: 10.1016/j.gene.2011.05.026
70. He X, Qi Z, Hou H, Gao J, Zhang XX. Effects of Chronic Cadmium Exposure at Food Limitation-Relevant Levels on Energy Metabolism in Mice. J Hazard Mater (2020) 388:121791. doi: 10.1016/j.jhazmat.2019.121791
71. Nhan JD, Turner CD, Anderson SM, Yen CA, Dalton HM, Cheesman HK, et al. Redirection of SKN-1 Abates the Negative Metabolic Outcomes of a Perceived Pathogen Infection. Proc Natl Acad Sci U.S.A. (2019) 116(44):22322–30. doi: 10.1073/pnas.1909666116
72. Ganeshan K, Chawla A. Metabolic Regulation of Immune Responses. Annu Rev Immunol (2014) 32:609–34. doi: 10.1146/annurev-immunol-032713-120236
73. McGeer JC, Szebedinszky C, McDonald DG, Wood CM. Effects of Chronic Sublethal Exposure to Waterborne Cu, Cd or Zn in Rainbow Trout. 1: Iono-Regulatory Disturbance and Metabolic Costs. Aquat Toxicol (2000) 50(3):231–43. doi: 10.1016/s0166-445x(99)00105-8
74. Manstein DJ, Preller M. Small Molecule Effectors of Myosin Function. Adv Exp Med Biol (2020) 1239:61–84. doi: 10.1007/978-3-030-38062-5_5
75. Roberts RJ, Agius C, Saliba C, Bossier P, Sung YY. Heat Shock Proteins (Chaperones) in Fish and Shellfish and Their Potential Role in Relation to Fish Health: A Review. J Fish Dis (2010) 33(10):789–801. doi: 10.1111/j.1365-2761.2010.01183.x
76. Levine B, Kroemer G. Biological Functions of Autophagy Genes: A Disease Perspective. Cell (2019) 176(1-2):11–42. doi: 10.1016/j.cell.2018.09.048
77. Pradeepkiran JA, Reddy PH. Defective Mitophagy in Alzheimer’s Disease. Ageing Res Rev (2020) 64:101191. doi: 10.1016/j.arr.2020.101191
78. Magnadottir B. Immunological Control of Fish Diseases. Mar Biotechnol (NY) (2010) 12(4):361–79. doi: 10.1007/s10126-010-9279-x
79. Lee JW, Choi H, Hwang UK, Kang JC, Kang YJ, Kim KI, et al. Toxic Effects of Lead Exposure on Bioaccumulation, Oxidative Stress, Neurotoxicity, and Immune Responses in Fish: A Review. Environ Toxicol Pharmacol (2019) 68:101–8. doi: 10.1016/j.etap.2019.03.010
80. Zhang H, Shen L, Fang W, Zhang X, Zhong Y. Perfluorooctanoic Acid-Induced Immunotoxicity via NF-Kappa B Pathway in Zebrafish (Danio Rerio) Kidney. Fish Shellfish Immunol (2021) 113:9–19. doi: 10.1016/j.fsi.2021.03.004
81. Kawai T, Akira S. TLR Signaling. Semin Immunol (2007) 19(1):24–32. doi: 10.1016/j.smim.2006.12.004
82. Lallès JP. Biology, Environmental and Nutritional Modulation of Skin Mucus Alkaline Phosphatase in Fish: A Review. Fish Shellfish Immunol (2019) 89:179–86. doi: 10.1016/j.fsi.2019.03.053
83. Nibbs RJ, Graham GJ. Immune Regulation by Atypical Chemokine Receptors. Nat Rev Immunol (2013) 13(11):815–29. doi: 10.1038/nri3544
84. Borroni EM, Savino B, Bonecchi R, Locati M. Chemokines Sound the Alarmin: The Role of Atypical Chemokine in Inflammation and Cancer. Semin Immunol (2018) 38:63–71. doi: 10.1016/j.smim.2018.10.005
85. Karsi A, Wolters WR, Waldbieser GC. Assignment of Immune-Related Genes to the Channel Catfish, Ictalurus Punctatus, Genetic Map. Anim Genet (2005) 36(6):502–6. doi: 10.1111/j.1365-2052.2005.01345.x
86. Huhtala M, Heino J, Casciari D, de Luise A, Johnson MS. Integrin Evolution: Insights From Ascidian and Teleost Fish Genomes. Matrix Biol (2005) 24(2):83–95. doi: 10.1016/j.matbio.2005.01.003
87. Zhou S, Dong J, Liu Y, Yang Q, Xu N, Yang Y, et al. Effects of Acute Deltamethrin Exposure on Kidney Transcriptome and Intestinal Microbiota in Goldfish (Carassius Auratus). Ecotoxicol Environ Saf (2021) 225:112716. doi: 10.1016/j.ecoenv.2021.112716
88. Griffin BR, Mitchell AJ. Susceptibility of Channel Catfish, Ictalurus Punctatus (Rafinesque), to Edwardsiella Ictaluri Challenge Following Copper Sulphate Exposure. J Fish Dis (2007) 30(10):581–5. doi: 10.1111/j.1365-2761.2007.00838.x
89. Migliorini P, Italiani P, Pratesi F, Puxeddu I, Boraschi D. The IL-1 Family Cytokines and Receptors in Autoimmune Diseases. Autoimmun Rev (2020) 19(9):102617. doi: 10.1016/j.autrev.2020.102617
90. Wu XX, Huang XL, Chen RR, Li T, Ye HJ, Xie W, et al. Paeoniflorin Prevents Intestinal Barrier Disruption and Inhibits Lipopolysaccharide (LPS)-Induced Inflammation in Caco-2 Cell Monolayers. Inflammation (2019) 42(6):2215–25. doi: 10.1007/s10753-019-01085-z
91. Li Q, Cao M, Fu Q, Yang N, Yan X, Song L, et al. Complement Genes in Black Rockfish (Sebastods Schlegelii): Genome-Wide Identification, Evolution and Their Potential Functions in Response to Vibrio Anguillarum Infection. Fish Shellfish Immunol (2021) 114:119–31. doi: 10.1016/j.fsi.2021.04.020
92. Lu Q, Sun Y, Ares I, Anadón A, Martínez M, Martínez-Larrañaga MR, et al. Deltamethrin Toxicity: A Review of Oxidative Stress and Metabolism. Environ Res (2019) 170:260–81. doi: 10.1016/j.envres.2018.12.045
93. Jiang Q, Jiang Z, Ao S, Gao X, Zhu X, Zhang Z, et al. Multi-Biomarker Assessment in the Giant Freshwater Prawn Macrobrachium Rosenbergii After Deltamethrin Exposure. Ecotoxicol Environ Saf (2021) 214:112067. doi: 10.1016/j.ecoenv.2021.112067
94. Singh S, Tiwari RK, Pandey RS. Evaluation of Acute Toxicity of Triazophos and Deltamethrin and Their Inhibitory Effect on AChE Activity in Channa Punctatus. Toxicol Rep (2018) 5:85–9. doi: 10.1016/j.toxrep.2017.12.006
95. Wang S, Guo C, Zhou L, Zhong Z, Zhu W, Huang Y, et al. Effects of Dietary Supplementation With Epidermal Growth Factor-Expressing Saccharomyces Cerevisiae on Duodenal Development in Weaned Piglets. Br J Nutr (2016) 115(9):1509–20. doi: 10.1017/s0007114516000738
96. Liu C, Zhao LP, Shen YQ. A Systematic Review of Advances in Intestinal Microflora of Fish. Fish Physiol Biochem (2021) 47(6):2041–53. doi: 10.1007/s10695-021-01027-3
97. Bledsoe JW, Waldbieser GC, Swanson KS, Peterson BC, Small BC. Comparison of Channel Catfish and Blue Catfish Gut Microbiota Assemblages Shows Minimal Effects of Host Genetics on Microbial Structure and Inferred Function. Front Microbiol (2018) 9:1073. doi: 10.3389/fmicb.2018.01073
98. Semova I, Carten JD, Stombaugh J, Mackey LC, Knight R, Farber SA, et al. Microbiota Regulate Intestinal Absorption and Metabolism of Fatty Acids in the Zebrafish. Cell Host Microbe (2012) 12(3):277–88. doi: 10.1016/j.chom.2012.08.003
99. Shi X, Ng KK, Li XR, Ng HY. Investigation of Intertidal Wetland Sediment as a Novel Inoculation Source for Anaerobic Saline Wastewater Treatment. Environ Sci Technol (2015) 49(10):6231–9. doi: 10.1021/acs.est.5b00546
100. Miao S, Zhao C, Zhu J, Hu J, Dong X, Sun L. Dietary Soybean Meal Affects Intestinal Homoeostasis by Altering the Microbiota, Morphology and Inflammatory Cytokine Gene Expression in Northern Snakehead. Sci Rep (2018) 8(1):113. doi: 10.1038/s41598-017-18430-7
101. Dong J, Li X, Zhang R, Zhao Y, Wu G, Liu J, et al. Comparative Analysis of the Intestinal Bacterial Community and Expression of Gut Immunity Genes in the Chinese Mitten Crab (Eriocheir Sinensis). AMB Express (2018) 8(1):192. doi: 10.1186/s13568-018-0722-0
102. Knutie SA, Wilkinson CL, Kohl KD, Rohr JR. Early-Life Disruption of Amphibian Microbiota Decreases Later-Life Resistance to Parasites. Nat Commun (2017) 8(1):86. doi: 10.1038/s41467-017-00119-0
103. Knutie SA, Gabor CR, Kohl KD, Rohr JR. Do Host-Associated Gut Microbiota Mediate the Effect of an Herbicide on Disease Risk in Frogs? J Anim Ecol (2018) 87(2):489–99. doi: 10.1111/1365-2656.12769
104. Xu Y, Li Y, Xue M, Yang T, Luo X, Fan Y, et al. Effects of Dietary Saccharomyces Cerevisiae YFI-SC2 on the Growth Performance, Intestinal Morphology, Immune Parameters, Intestinal Microbiota, and Disease Resistance of Crayfish (Procambarus Clarkia). Anim (Basel) (2021) 11(7):1963. doi: 10.3390/ani11071963
105. Fang H, Wang B, Jiang K, Liu M, Wang L. Effects of Lactobacillus Pentosus HC-2 on the Growth Performance, Intestinal Morphology, Immune-Related Genes and Intestinal Microbiota of Penaeus Vannamei Affected by Aflatoxin B1. Aquaculture (2020) 525:735289. doi: 10.1016/j.aquaculture.2020.735289
Keywords: channel catfish, deltamethrin, immunotoxicity, metabolic disorder, apoptosis
Citation: Yang Y, Zhu X, Huang Y, Zhang H, Liu Y, Xu N, Fu G and Ai X (2022) RNA-Seq and 16S rRNA Analysis Revealed the Effect of Deltamethrin on Channel Catfish in the Early Stage of Acute Exposure. Front. Immunol. 13:916100. doi: 10.3389/fimmu.2022.916100
Received: 08 April 2022; Accepted: 09 May 2022;
Published: 03 June 2022.
Edited by:
Brian Peterson, National Cold Water Marine Aquaculture Center (USDA-ARS), United StatesReviewed by:
Jaime Eugenio Figueroa, Austral University of Chile, ChileSudipta Maitra, Visva-Bharati University, India
Copyright © 2022 Yang, Zhu, Huang, Zhang, Liu, Xu, Fu and Ai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hongyu Zhang, zhanghy@cafs.ac.cn; Xiaohui Ai, aixh@yfi.ac.cn
†These authors have contributed equally to this work