RNA-seq analysis of mucosal immune responses reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus

Early divergent responses to virulent and attenuated vaccine isolates of Flavobacterium covae sp. nov. In channel catfish, Ictalurus punctatus

Columnaris disease continues to inflict substantial losses among freshwater cultured species since its first description one hundred years ago. The experimental and anecdotal evidence suggests an expanded range and rising virulence of columnaris worldwide due to the warming global climate. The channel catfish (Ictalurus punctatus) are particularly vulnerable to columnaris. A recently developed live attenuated vaccine (17-23) for Flavobacterium columnare (now Flavobacterium covae sp. nov.) demonstrated superior protection for vaccinated catfish against genetically diverse columnaris isolates. In this study, we aimed to elucidate the molecular mechanisms and patterns of immune evasion and host manipulation linked to virulence by comparing gene expression changes in the host after the challenge with a virulent (BGSF-27) or live attenuated F. covae sp. nov. vaccine (17-23). Thirty-day-old fry were accordingly challenged with either virulent or vaccine isolates. Gill tissues were collected at 0 h (control), 1 h, and 2 h post-infection, which are two critical time points in early host-pathogen interactions. Transcriptome profiling of the gill tissues revealed a larger number (518) of differentially expressed genes (DEGs) in vaccine-exposed fish than those exposed to the virulent pathogen (321). Pathway analyses suggested potent suppression of early host immune responses by the virulent isolate through a higher expression of nuclear receptor corepressors (NCoR) responsible for antagonizing macrophage and T-cell signaling. Conversely, in vaccinated fry, we observed induction of Ca2+/calmodulin-dependent protein kinase II (CAMKII), responsible for clearing NCoR, and commensurate up-regulation of transcription factor AP-1 subunits, c-Fos, and c-Jun. As in mammalian systems, AP-1 expression was connected with a broad immune activation in vaccinated fry, including induction of CC chemokines, proteinases, iNOS, and IL-12b. Relatedly, divergent expression patterns of Src tyrosine kinase Lck, CD44, and CD28 indicated a delay or suppression of T-cell adhesion and activation in fry exposed to the virulent isolate. Broader implications of these findings will be discussed. The transcriptomic differences between virulent and attenuated bacteria may offer insights into how the host responds to the vaccination or infection and provide valuable knowledge to understand the early immune mechanisms of columnaris disease in aquaculture.

First report of L1 cell adhesion molecule (L1CAM) in flat fish (Starry flounder): Molecular features and expression analysis

In this study, the starry flounder L1 cell adhesion molecule (L1CAM) sequence was obtained using next-generation sequencing, and the integrity of the sequence was verified by cloning and sequencing. First, the amino acid sequence was predicted using the cDNA sequence, and the gene was then identified through multiple sequence alignment analysis with related sequences and phylogenetic analysis. Thus, homogeneity was confirmed. The expression level of PsL1CAM (Platichthys stellatus L1CAM) mRNA in healthy starry flounder was detected in all tissues used in the experiment, and tissue- and gene-specific expression levels were confirmed. In addition, as a result of mRNA expression analysis after artificial infection with viral hemorrhagic septicemia virus (VHSV) and Streptococcus parauberis PH0710, significant expression changes and characteristics were confirmed following infection with VHSV and S. parauberis PH0710. After artificial infection with VHSV, the expression level of PsL1CAM mRNA was significantly upregulated in almost all major tissues of the starry flounder, whereas it was significantly downregulated in mucosal-associated lymphoid tissues, such as the gills and intestine. Infection with S. parauberis PH0710 significantly upregulated the expression of PsL1CAM mRNA in almost all major tissues of the starry flounder, whereas it was significantly downregulated in the heart after infection. Our results indicate that PsL1CAM may be involved in the host immune response to starry flounders.

Catfishes: A global review of the literature

This study aims to elucidate the evolution of catfish research publications over recent decades, identify emerging research clusters, examine keyword patterns, determine major contributors (including authors, organizations, and funding agencies), and analyze their collaborative networks and citation bursts on a global scale. The USA, Brazil, China, and India collectively contribute approximately 67% of the total catfish research publications, with a marked increase in prevalence since 2016. The most frequently occurring and dominant keywords are "channel catfish" and "responses," respectively. Intriguingly, our findings reveal 28 distinct article clusters, with prominent clusters including "yellow catfish," "channel catfish", "pectoral girdle," "African catfish", "Rio Sao Francisco basin," "Edwardsiella ictaluri," and "temperature mediated". Concurrently, keyword clustering generates seven main clusters: "new species", "growth performance", "heavy metal", "gonadotropin-releasing", "essential oil", and "olfactory receptor". This study further anticipates future research directions, offering fresh perspectives on the catfish literature landscape. To the best of our knowledge, this is the first article to conduct a comprehensive mapping review of catfish research publications worldwide.

Transcriptome analysis of Macrobrachium rosenbergii hemocytes in response to Staphylococcus aureus infection

The aquaculture industry has suffered significant financial losses as a result of disease outbreaks. In particular, disease outbreaks have become a major problem that can seriously affect the sustainable development of the Macrobrachium rosenbergii aquaculture industry. It is crucial to determine the defense mechanism of the host after pathogenic invasion in order to provide effective defense measures after disease outbreaks. Shrimp, like other invertebrates, primarily depend on their innate immune systems to defend against pathogens, and recognize and resist pathogens through humoral and cellular immune responses. In this investigation, we used RNA-seq technology to investigate the transcriptome of hemocytes from M. rosenbergii induced by Staphylococcus aureus. Our main targets were immune pathways and genes related to innate immunity. RNA-seq identified 209,069 and 204,775 unigenes in the control and experimental groups, respectively. In addition, we identified 547 and 1734 differentially expressed genes (DEGs) following S. aureus challenge after 6 and 12 h (h), respectively. GO and KEGG enrichment analysis revealed that the DEGs were significantly enriched in several biological signalling pathways, including NOD-like receptor, PI3K-Akt, Toll and Imd, IL-17, TGF-beta, RIG-I-like receptor, cAMP, apoptosis, and C-type lectin receptor. Sixteen DEGs were chosen at random for qPCR verification; these results concurred with those from sequencing. Our findings revealed that immune-related genes play an important role in antibacterial activities and have specific functions for gram-positive bacteria. These results provide more data for the prevention of M. rosenbergii diseases and offer a basis for the better prevention of diseases.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Parasites and the neuroendocrine control of fish intestinal function: an ancient struggle between pathogens and host

Most individual fish in wild and farmed populations can be infected with parasites. Fish intestines can harbour protozoans, myxozoans and helminths, which include several species of digeneans, cestodes, nematodes and acanthocephalans. Enteric parasites often induce inflammation of the intestine; the pathogen provokes changes in the host physiology, which will be genetically selected for if they benefit the parasite. The host response to intestinal parasites involves neural, endocrine and immune systems and interaction among these systems is coordinated by hormones, chemokines, cytokines and neurotransmitters including peptides. Intestinal fish parasites have effects on the components of the enteric nervous and endocrine systems; mechanical/chemical changes impair the activity of these systems, including gut motility and digestion. Investigations on the role of the neuroendocrine system in response to fish intestinal parasites are very few. This paper provides immunohistochemical and ultrastructural data on effects of parasites on the enteric nervous system and the enteric endocrine system in several fish–parasite systems. Emphasis is on the occurrence of 21 molecules including cholecystokinin-8, neuropeptide Y, enkephalins, galanin, vasoactive intestinal peptide and serotonin in infected tissues.

The early function of cortisol in liver during Aeromonas hydrophila infection: Dynamics of the transcriptome and accessible chromatin landscapes

In China, channel catfish (Ictalurus punctatus) is an important aquaculture species; however, haemorrhagic disease (Aeromonas hydrophila induced disease) in these fish has caused tremendous economic loss due to high morbidity and mass mortality in the breeding industry. The role of cortisol in bacterial diseases, particularly in the acute phase, remains unclear. In this study, liver transcriptome (RNA-seq) and chromatin accessibility (ATAC-seq) analyses were employed to investigate the early functional role of cortisol in Aeromonas hydrophila-stimulated responses. Our experiments confirmed that A. hydrophila infection can initially significantly increase serum cortisol levels at 1 h after infection. At this time point, the increased serum cortisol levels can significantly regulate A. hydrophila-regulated genes by affecting both transcriptome and chromatin accessibility. Cross-analysis of RNA-seq and ATAC-seq revealed that a certain gene group (92 target_DEGs) was regulated at an early time point by cortisol. KEGG enrichment analysis revealed that the top three pathways according to target_DEGs were cancer, glutathione metabolism, and the Notch signalling pathway. The protein-protein interaction analysis of target_DEGs revealed that they may be primarily involved in cell proliferation, CD8⁺ T cell function, glutathione synthesis, and activation of the NF-κB signalling pathway. This suggests that after the emergence of immune stress, the early regulation of cortisol is positive against the immune response. It is possible that in this situation, the animal is attempting to avoid dangerous situations and risks and then cope with the imbalance produced by the stressor to ultimately restore homeostasis. Our results will contribute to future research on fish and provide valuable insight regarding the mechanism of immune regulation by cortisol and the study of bacterial haemorrhagic disease in channel catfish.

Molecular characterization and expression analysis of septin gene family and phagocytic function of recombinant septin 2, 3 and 8 of starry flounder (Platichthys stellatus)

Septin is an evolutionarily conserved family of GTP-binding proteins. Septins are known to be involved in a variety of cellular processes, including cell division, chromosome separation, cell polarity, motility, membrane dynamics, exocytosis, apoptosis, phagocytosis, DNA damage responses, and other immune responses. In this study, the sequences of the septin gene family of starry flounder were obtained using NGS sequencing, and the integrity of the sequences was verified through cloning and sequencing. At first, the amino acid sequence was annotated using the cDNA sequence, and then, the gene sequence was verified through multiple sequence alignment and phylogenetic analyses using the related conserved sequences. The septin gene family was classified into three subgroups based on the phylogenetic analysis. High conservation within the domain and homology between the genes reported in different species were confirmed. The expression level of septin gene family mRNA in each tissue of healthy starry flounder was evaluated to confirm the tissue- and gene-specific expression levels. Additionally, as a result of the analysis of mRNA expression after simulated pathogen infection, significant expression changes and characteristics were confirmed upon infection with bacteria (Streptococcus parauberis PH0710) and virus (VHSV). Based on the current results and that of previous studies, to confirm the immunological function, Septin 2, 3, and 8 were produced as recombinant proteins based on the amino acid sequences, and their role in phagocytosis was further investigated. The results of this study indicate that septin gene family plays a complex and crucial role in the host immune response to pathogens of starry flounder.

Integrative omics analysis highlights the immunomodulatory effects of the parasitic dinoflagellate hhematodinium on crustacean hemocytes

Parasitic dinoflagellates in genus Hematodinium have caused substantial economic losses to multiple commercially valuable marine crustaceans around the world. Recent efforts to better understand the life cycle and biology of the parasite have improved our understanding of the disease ecology. However, studies on the host-parasite interaction, especially how Hematodinium parasites evade the host immune response are lacking. To address this shortfall, we used the comprehensive omics approaches (miRNA transcriptomics, iTRAQ-based proteomics) to get insights into the host-parasite interaction between hemocytes from Portunus trituberculatus and Hematodinium perezi in the present study. The parasitic dinoflagellate H. perezi remodeled the miRNome and proteome of hemocytes from challenged hosts, modulated the host immune response at both post-transcriptional and translational levels and caused post-transcriptional regulation to the host immune response. Multiple important cellular and humoral immune-related pathways (ex. Apoptosis, Endocytosis, ECM-receptor interaction, proPO activation pathway, Toll-like signaling pathway, Jak-STAT signaling pathway) were significantly affected by Hematodinium parasites. Through modulation of the host miRNome, the host immune responses of nodulation, proPO activation and antimicrobial peptides were significantly suppressed. Cellular homeostasis was imbalanced via post-transcriptional dysregulation of the phagosome and peroxisome pathways. Cellular structure and communication was seriously impacted by post-transcriptional downregulation of ECM-receptor interaction and focal adhesion pathways. In conclusion, H. perezi parasites could trigger striking changes in the miRNome and proteome of crustacean hemocytes, and this parasite exhibited multifaceted immunomodulatory effects and potential immune-suppressive mechanisms in crustacean hosts.

Allelically and Differentially Expressed Genes After Infection of Edwardsiella ictaluri in Channel Catfish as Determined by Bulk Segregant RNA-Seq

Identification of genetic markers associated with resistance against enteric septicemia of catfish (ESC) is of great interest for genetic enhancement programs of catfish. In the present study, bulk segregant RNA-Seq analysis was applied to determine differentially expressed genes and alleles after ESC infection. Here we report three genomic regions on LG1, LG12, and LG26, containing significant single-nucleotide polymorphisms (SNPs). These genomic regions aligned well with quantitative trait loci (QTL) previously identified. Within the QTL regions, eleven genes were found to be differentially regulated between phenotypic bulks. Importantly, the QTL on linkage group 1 (LG1) were found to be expressed in the liver, whereas the QTL on LG12 and LG26 were expressed in the intestine, suggesting multiple mechanisms of ESC resistance. It is apparent that apolipoproteins may be important for ESC resistance as the QTL on LG1 included the 14-kDa apolipoprotein genes that are both allelically expressed and differentially expressed between the resistant and susceptible bulks. Traf2 and NCK-interacting protein kinase (TNIK) were found in the QTL on LG12, and it was downregulated in resistant fish, suggesting the importance of NCK downregulation in ESC resistance, as previously reported. In addition, we observed divergent gene expression patterns between the liver and intestine after infection. Immune/inflammatory-related processes were overrepresented from liver DEGs, while those DEGs identified from intestine were enriched for proteolysis and wounding processes. Taken together, the BSR-Seq analysis presented here advanced the knowledge of ESC resistance, providing information of not only positions of QTL but also genes and their differential expression between resistant and susceptible fish, making it one step closer to the identification of the causal genes for ESC resistance.

RNA-Sequencing Analysis of the Spleen and Gill of Takifugu rubripes in Response to Vibrio harveyi Infection

Takifugu rubripes is commonly subjected to the disease-causing bacterium, Vibrio harveyi. However, the mechanism involved in the immune response of T. rubripes to V. harveyi infection is unclear. We conducted a transcriptomic analysis of the spleen and gill from T. rubripes infected with V. harveyi. We obtained 60,981,357 and 60,760,550 clean reads from the control and infected spleens, and 57,407,586 and 57,536,651 clean reads from the control and infected gills, respectively. We also identified 1,560 and 1,213 differentially expressed genes in the spleen and gill, respectively. Gene ontology analysis revealed that the most enriched biological process in both the spleen and gill was "immune response". The most enriched Kyoto Encyclopedia of Genes and Genomes immune response-related pathways were the NOD-like receptor signaling pathway in the spleen and cytokine-cytokine receptor interaction in the gill. We found 10 candidate immune-related genes in the spleen and gill. These putative immune pathways and candidate genes will provide insight into the immune response mechanisms of T. rubripes against V. harveyi.

The Programming of Antioxidant Capacity, Immunity, and Lipid Metabolism in Dojo Loach (Misgurnus anguillicaudatus) Larvae Linked to Sodium Chloride and Hydrogen Peroxide Pre-treatment During Egg Hatching

Non-nutritional stress during early life period has been reported to promote the metabolic programming in fish induced by nutritional stimulus. Sodium chloride (NaCl) and hydrogen peroxide (H₂O₂) have been widely applied during fish egg hatching, but the influences on health and metabolism of fish in their later life remain unknown. In the present study, H₂O₂ treatment at 400mg/L but not 200mg/L significantly increased the loach hatchability and decreased the egg mortality, while NaCl treatment at 1,000 and 3,000mg/L showed no significant influences on the loach hatchability nor egg mortality. Further studies indicated that 400mg/L H₂O₂ pre-treatment significantly enhanced the antioxidant capacity and the mRNA expression of genes involved in immune response of loach larvae, accompanied by the increased expression of genes involved in fish early development. However, the expression of most genes involved in lipid metabolism, including catabolism and anabolism of loach larvae, was significantly upregulated after 200mg/L H₂O₂ pre-treatment. NaCl pre-treatment also increased the expression of antioxidant enzymes; however, only the expression of C1q within the detected immune-related genes was upregulated in loach larvae. One thousand milligram per liter NaCl pre-treatment significantly increased the expression of LPL and genes involved in fish early development. Thus, our results suggested the programming roles of 400mg/L H₂O₂ pre-treatment during egg hatching in enhancing antioxidant capacity and immune response of fish larvae via promoting fish early development.

Immune Response and Apoptosis-Related Pathways Induced by Aeromonas schubertii Infection of Hybrid Snakehead (Channa maculata♀ × Channa argus♂)

Aeromonas schubertii is the etiological pathogen of internal organ nodules in snakehead fish. Infections with A. schubertii produce a significant economic loss in aquaculture. Therefore, it is important to examine the immune mechanisms by which snakeheads defend against A. schubertii infection. In this study, we established a hybrid snakehead infection model by intraperitoneal injection of A. schubertii that produced internal organ nodules. The splenic immune response of infected fish was examined at the transcriptome level by Illumina-seq analysis. Results showed 14,796 differentially expressed genes (DEGs) following A. schubertii infection, including 4441 up-regulated unigenes and 10,355 down-regulated unigenes. KEGG analysis showed 2084 DEGs to be involved in 192 pathways, 14 of which were immune-related. Twelve DEGs were used to validate quantitative real-time PCR results with RNA-seq data. Time-course expression analysis of six genes demonstrated modulation of the snakehead immune response by A. schubertii. Furthermore, transcriptome analysis identified a substantial number of DEGs that were involved in the apoptosis signaling pathway. TUNEL analysis of infected spleens confirmed the presence of apoptotic cells. This study provided new information for a further understanding of the pathogenesis of A. schubertii in snakeheads, which can be used to prevent and possibly treat A. schubertii infections.

Effects of dietary silkrose of Antheraea yamamai on gene expression profiling and disease resistance to Edwardsiella tarda in Japanese medaka (Oryzias latipes)

We previously identified a novel acidic polysaccharide, silkrose-AY, from the Japanese oak silkmoth (Antheraea yamamai), which can activate an innate immune response in mouse macrophage cells. However, innate immune responses stimulated by silkrose-AY in teleosts remain unclear. Here, we show the influence of dietary silkrose-AY in medaka (Oryzias latipes), a teleost model, in response to Edwardsiella tarda infection. Dietary silkrose-AY significantly improved the survival of fish and decreased the number of bacteria in their kidneys after the fish were artificially infected with E. tarda by immersion. We also performed a microarray analysis of the intestine, which serves as a primary barrier against microbial infection, to understand the profiles of differentially expressed genes (DEGs) evoked by silkrose-AY. The dietary silkrose-AY group showed differential expression of 2930 genes when compared with the control group prior to E. tarda infection. Gene ontology and pathway analysis of the DEGs highlighted several putative genes involved in pathogen attachment/recognition, the complement and coagulation cascade, antimicrobial peptides/enzymes, opsonization/phagocytosis, and epithelial junctional modification. Our findings thus provide fundamental information to help understand the molecular mechanism of bacterial protection offered by insect-derived immunostimulatory polysaccharides in teleosts.

Dietary Phytogenics and Galactomannan Oligosaccharides in Low Fish Meal and Fish Oil-Based Diets for European Sea Bass (Dicentrarchus labrax) Juveniles: Effects on Gill Structure and Health and Implications on Oxidative Stress Status

An effective replacement for fish meal (FM) and fish oil (FO) based on plant-based raw materials in the feed of marine fish species is necessary for the sustainability of the aquaculture sector. However, the use of plant-based raw materials to replace FM and FO has been associated with several negative health effects, some of which are related to oxidative stress processes that can induce functional and morphological alterations in mucosal tissues. This study aimed to evaluate the effects of dietary oligosaccharides of plant origin (5,000 ppm; galactomannan oligosaccharides, GMOS) and a phytogenic feed additive (200 ppm; garlic oil and labiatae plant extract mixture, PHYTO) on the oxidative stress status and mucosal health of the gills of juvenile European sea bass (Dicentrarchus labrax). The experimental diets, low FM and FO diets (10%FM/6%FO) were supplemented with GMOS from plant origin and PHYTO for 63 days. GMOS and PHYTO did not significantly affect feed utilization, fish growth, and survival. GMOS and PHYTO downregulated the expression of β-act, sod, gpx, cat, and gr in the gills of the fish compared with that in fish fed the control diet. The expression of hsp70 and ocln was upregulated and downregulated, respectively, in the GMOS group compared with that in the control group, whereas the expression of zo-1 was downregulated in the PHYTO group compared with that in the GMOS group. The morphological, histopathological, immunohistochemical, and biochemical parameters of the fish gills were mostly unaffected by GMOS and PHYTO. However, the PHYTO group had lower incidence of lamellar fusion than did the control group after 63 days. Although the tissular distribution of goblet cells was unaffected by GMOS and PHYTO, goblet cell size showed a decreasing trend (−11%) in the GMOS group. GMOS and PHYTO significantly reduced the concentration of PCNA+ in the epithelium of the gills. The above findings indicated that GMOS and PHYTO in low FM/FO-based diets protected the gill epithelia of D. labrax from oxidative stress by modulating the expression of oxidative enzyme-related genes and reducing the density of PCNA+ cells in the gills of the fish.

Identification and Characterization of Long Non-coding RNAs in the Intestine of Olive Flounder (Paralichthys olivaceus) During Edwardsiella tarda Infection

Long non-coding RNAs (lncRNAs) play widespread roles in fundamental biological processes, including immune responses. The olive flounder (Paralichthys olivaceus), an important economical flatfish widely cultured in Japan, Korea, and China, is threatened by infectious pathogens, including bacteria, viruses, and parasites. However, the role of lncRNAs in the immune responses of this species against pathogen infections is not well-understood. Therefore, in this study, we aimed to identify lncRNAs in the intestine of olive flounder and evaluate their differential expression profiles during Edwardsiella tarda infection, which is an important zoonotic and intestinal pathogen. A total of 4,445 putative lncRNAs were identified, including 3,975 novel lncRNAs and 470 annotated lncRNAs. These lncRNAs had shorter lengths and fewer exons compared with mRNAs. In total, 115 differentially expressed lncRNAs (DE-lncRNAs) were identified during E. tarda infection. To validate the expression pattern of lncRNAs, six DE-lncRNAs were randomly selected for quantitative real-time PCR. The co-located and co-expressed mRNAs of DE-lncRNAs were predicted, which were used to conduct the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The target genes of DE-lncRNAs enriched numerous immune-related processes and exhibited a strong correlation with immune-related signaling pathways. To better understand the extensive regulatory functions of lncRNAs, the lncRNA–miRNA–mRNA regulatory networks were constructed, and two potential competing endogenous RNA (ceRNA) networks, LNC_001979-novel_171-Potusc2 and LNC_001979-novel_171-Podad1, were preliminarily identified from the intestine of olive flounders for the first time. In conclusion, this study provides an invaluable annotation and expression profile of lncRNAs in the intestine of olive flounder infected with E. tarda; this forms a basis for further studies on the regulatory function of lncRNAs in the intestinal mucosal immune responses of olive flounder.

RGD-binding integrins mediated phagocytosis involved in the entry of Edwardsiella tarda into mudskipper MO/MФ

The versatile fish pathogen Edwardsiella tarda is an intracellular pathogen with the ability to invade and replicate in host phagocytes. However, the mechanism mediating the uptake of E. tarda in fish monocytes/macrophages (MO/MΦ) is not yet understood. Generating mudskipper kidney-derived MO/MФ transcriptomic resources from mudskipper challenged by E. tarda is crucial for understanding the molecular mechanisms underlying the mudskipper invasion process. In the present study, a total of 1185 up-regulated and 885 down-regulated differentially expressed genes (DEGs) were identified using RNA-seq. Enrichment and pathway analysis of DEGs revealed the centrality of the phagosome and regulation of actin cytoskeleton pathways in pathogen entry. The progress of phagosome formation was observed by transmission electron microscopy. Eight conserved integrin (ITG) subunit genes, belonging to the phagocytic receptors, were found in the transcriptomic sequence data. Additionally, quantitative real-time PCR showed that the mRNA expressions of most ITG subunit genes were related to the different infection times of E. tarda and the different bacterial pathogens. Further assays demonstrated that phagocytosis of FITC-labeled E. tarda by mudskipper MO/MФ was significantly reduced by the tetrapeptide Asp-Gly-Arg-Ser (RGDS). In summary, phagocytosis is one of the entry pathways into mudskipper MO/MΦ, and RGD-binding ITGs are involved in the phagosome formation process.

A tale of two fish: Comparative transcriptomics of resistant and susceptible steelhead following exposure to Ceratonova shasta highlights differences in parasite recognition

Diseases caused by myxozoan parasites represent a significant threat to the health of salmonids in both the wild and aquaculture setting, and there are no effective therapeutants for their control. The myxozoan Ceratonova shasta is an intestinal parasite of salmonids that causes severe enteronecrosis and mortality. Most fish populations appear genetically fixed as resistant or susceptible to the parasite, offering an attractive model system for studying the immune response to myxozoans. We hypothesized that early recognition of the parasite is a critical factor driving resistance and that susceptible fish would have a delayed immune response. RNA-seq was used to identify genes that were differentially expressed in the gills and intestine during the early stages of C. shasta infection in both resistant and susceptible steelhead (Oncorhynchus mykiss). This revealed a downregulation of genes involved in the IFN-γ signaling pathway in the gills of both phenotypes. Despite this, resistant fish quickly contained the infection and several immune genes, including two innate immune receptors were upregulated. Susceptible fish, on the other hand, failed to control parasite proliferation and had no discernible immune response to the parasite, including a near-complete lack of differential gene expression in the intestine. Further sequencing of intestinal samples from susceptible fish during the middle and late stages of infection showed a vigorous yet ineffective immune response driven by IFN-γ, and massive differential expression of genes involved in cell adhesion and the extracellular matrix, which coincided with the breakdown of the intestinal structure. Our results suggest that the parasite may be suppressing the host’s immune system during the initial invasion, and that susceptible fish are unable to recognize the parasite invading the intestine or mount an effective immune response. These findings improve our understanding of myxozoan-host interactions while providing a set of putative resistance markers for future studies.

Integrated Analysis of circRNA-miRNA-mRNA Regulatory Networks in the Intestine of Sebastes schlegelii Following Edwardsiella tarda Challenge

Sebastes schlegelii, an important aquaculture species, has been widely cultured in East Asian countries. With the increase in the cultivation scale, various diseases have become major threats to the industry. Evidence has shown that non-coding RNAs (ncRNAs) have remarkable functions in the interactions between pathogens and their hosts. However, little is known about the mechanisms of circular RNAs (circRNAs) and coding RNAs in the process of preventing pathogen infection in the intestine in teleosts. In this study, we aimed to uncover the global landscape of mRNAs, circRNAs, and microRNAs (miRNAs) in response to Edwardsiella tarda infection at different time points (0, 2, 6, 12, and 24 h) and to construct regulatory networks for exploring the immune regulatory mechanism in the intestine of S. schlegelii. In total, 1,794 mRNAs, 87 circRNAs, and 79 miRNAs were differentially expressed. The differentially expressed RNAs were quantitatively validated using qRT-PCR. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that most of the differentially expressed mRNA genes and the target genes of ncRNAs were related to immune signaling pathways, such as the NF-κB signal pathway, pathogen recognition receptors related to signaling pathways (Toll-like receptors and Nod-like receptors), and the chemokine signaling pathway. Based on these differentially expressed genes, 624 circRNA-miRNA pairs and 2,694 miRNA-mRNA pairs were predicted using the miRanda software. Integrated analyses generated 25 circRNA-miRNA-mRNA interaction networks. In a novel_circ_0004195/novel-530/IκB interaction network, novel_530 was upregulated, while its two targets, novel_circ_0004195 and IκB, were downregulated after E. tarda infection. In addition, two circRNA-miRNA-mRNA networks related to apoptosis (novel_circ_0003210/novel_152/apoptosis-stimulating of p53 protein 1) and interleukin (novel_circ_0001907/novel_127/interleukin-1 receptor type 2) were also identified in our study. We thus speculated that the downstream NF-κB signaling pathway, p53 signaling pathway, and apoptosis pathway might play vital roles in the immune response in the intestine of S. schlegelii. This study revealed a landscape of RNAs in the intestine of S. schlegelii during E. tarda infection and provided clues for further study on the immune mechanisms and signaling networks based on the circRNA-miRNA-mRNA axis in S. schlegelii.

Time-resolved RNA-seq provided a new understanding of intestinal immune response of European eel (Anguilla anguilla) following infection with Aeromonas hydrophila

No studies systematically examined the intestinal immune response for yellow stage of European eel (Anguilla anguilla) with Aeromonas hydrophila infection by time-resolved RNA-seq. Here, we examined transcriptional profiles of the intestines at three-time points following infection with A. hydrophila. Intraperitoneal injections caused mortalities within 48 h post-injection (hpi), with the survival rate 87.5% at 24 hpi and 83.9% at 48 hpi. The result from KEGG pathway enrichment analysis showed that the immune related "cytosolic DNA-sensing pathway" was significantly enriched at the first and second time points (6 hpi and 18 hpi), with the up-regulated expression of irf3, il1b, tnfaip3, cxcl8a, ap1-2, c-fos, polr3d, polr3g and polr3k both at 6 hpi and 18 hpi, but not at the third time point (36 hpi). According to the KEGG annotation, 326 immune and inflammation-related DEGs were found. The co-expression network of those 326 DEGs revealed the existence of three modules, and tlr1 was found to be in the center of the biggest module which contained massive DEGs from "signal transduction" and "transport and catabolism". The c3 isoforms showed different expression pattern among the three time points, indicating a unique activation of complement systems at 18 hpi. Furthermore, two cathelicidins (aaCATH_1 and aaCATH_2) were highly up-regulated at the first two time points, and the bacterial growth inhibition assay revealed their antibacterial properties against A. hydrophila. Our data indicated the important roles of cytosolic DNA-sensing pathway, as well as transcripts including tlr1, c3, polr and cathelicidins in the intestine of A. anguilla in response to A. hydrophila infection. The present study will provide leads for functional studies of host-pathogen interactions.

RNA-Seq analysis of the guppy immune response against Gyrodactylus bullatarudis infection

Gyrodactylids are ubiquitous ectoparasites of teleost fish, but our understanding of the host immune response against them is fragmentary. Here, we used RNA-Seq to investigate genes involved in the primary response to infection with Gyrodactylus bullatarudis on the skin of guppies, Poecilia reticulata, an important evolutionary model, but also one of the most common fish in the global ornamental trade. Analysis of differentially expressed genes identified several immune-related categories, including IL-17 signalling pathway and Th17 cell differentiation, cytokine-cytokine receptor interaction, chemokine signalling pathway, NOD-like receptor signalling pathway, natural killer cell-mediated cytotoxicity and pathways involved in antigen recognition, processing and presentation. Components of both the innate and the adaptive immune responses play a role in response to gyrodactylid infection. Genes involved in IL-17/Th17 response were particularly enriched among differentially expressed genes, suggesting a significant role for this pathway in fish responses to ectoparasites. Our results revealed a sizable list of genes potentially involved in the teleost-gyrodactylid immune response.

Differentially expressed proteins in the intestine of Cynoglossus semilaevis Günther following a Shewanella algae challenge

Fish intestine is an important constituent of the mucosal immune system. The gut and gut-associated lymphoid tissue construct a local immune environment. A Shewanella algae strain was previously reported to be a pathogen causing ascitic disease accompanied with intestinal inflammation in Cynoglossus semilaevis. This study aimed to investigate the intestine immune response in C. semilaevis to S. algae infection at the protein level. Two-dimensional electrophoresis coupled with mass spectrometry proteomics was utilized to compare protein expression in the intestines from normal and S. algae-infected C. semilaevis. A total of 70 differentially expressed proteins (DEPs), consisting of 16 upregulated and 54 downregulated proteins, were identified in the intestine tissue of C. Semilaevis. These protein expression changes were further validated using western blot analysis and quantitative real-time PCR. Gene ontology enrichment analysis showed that these 70 DEPs could be assigned across three categories: "cellular components", "molecular function", and "biological process". Forty-one DEPs (six up-regulated and 35 down-regulated proteins) related to metabolic processes were identified. In addition, 20 DEPs (eight up-regulated and 12 down-regulated proteins) related to stress and immune responses were identified. A protein-protein interaction network generated by the STRING (Search Tool for the Retrieval of Interacting Genes/protein) revealed that 30 DEPs interacted with one another to form an integrated network. Among them, 29 DEPs were related to stress, immune, and metabolism processes. In the network, some of the immune related proteins (C9, FGB, KNG1, apolipoprotein A-IV-like, and PDIA3) were up-regulated and most DEPs involved in metabolism processes were down-regulated. These results indicate that the immune defense response of the intestine was activated and the intestinal function associated with metabolism processes was disturbed. This study provides valuable information for further research into the functions of these DEPs in fish.

Advancing the fathead minnow (Pimephales promelas) as a model for immunotoxicity testing: Characterization of the renal transcriptome following Yersinia ruckeri infection

Recent studies have utilized the fathead minnow (Pimephales promelas) to explore the immunotoxic effects associated with a variety of environmental contaminants in the absence of immunological stimuli. Though this approach allows for alterations in the resting immune system to be detected, previous evidence suggests that many immunotoxic effects may only manifest in the activated immune system. However, basic immune responses to pathogens have not been well described in this species. To expand the utility of the fathead minnow as a model for immunotoxicity testing, a more comprehensive understanding of the activated immune system is required. As such, the main goal of this study was to characterize the transcriptomic response to pathogen infection in the fathead minnow using RNA sequencing. To achieve this goal, female fathead minnows were intraperitoneally injected with either Hank's Balanced Salt Solution (sham-injected) or Yersinia ruckeri (pathogen-injected). Eight hours following injection, fish were sacrificed for the assessment of general morphological (i.e., mass, length, condition factor, hepatic index) and immunological (i.e., leukocyte counts, spleen index) endpoints. To assess the molecular immune response to Y. ruckeri, kidney tissue was collected for transcriptomic analysis. A comparison of sham- and pathogen-injected fish revealed that >1800 genes and >500 gene networks were differentially expressed.Gene networks associated with inflammation, innate immunity, complement, hemorrhaging and iron absorption are highlighted and their utility within the context of immunotoxicity is discussed. These data reveal pathogen-related molecular endpoints to improve data interpretation of future studies utilizing the fathead minnow as a model for immunotoxicity.

Effects of inflammation and/or infection on the neuroendocrine control of fish intestinal motility: A review

Food is the largest expense in fish farms. On the other hand, the fish health and wellbeing are determining factors in aquaculture production where nutrition is a vital process for growing animals. In fact, it is important to remember that digestion and nutrition are crucial for animals' physiology. However, digestion is a very complex process in which food is processed to obtain necessary nutrients and central mechanisms of this process require both endocrine and neuronal regulation. In this context, intestinal motility is essential for the absorption of the nutrients (digestive process determining nutrition). An imbalance in the intestinal motility due to an inadequate diet or an infectious process could result in a lower use of the food and inefficiency in obtaining nutrients from food. Very frequently, farmed fish are infected with different pathogenic microorganism and this situation could alter gastrointestinal physiology and, indirectly reduce fish growth. For these reasons, the present review focuses on analysing how different inflammatory molecules or infections can alter conventional modulators of fish intestinal motility.

Liver transcriptome analysis and cortisol immune-response modulation in lipopolysaccharide-stimulated in channel catfish (Ictalurus punctatus)

Channel catfish (Ictalurus punctatus) is an important aquaculture species in China. In channel catfish, diseases such as haemorrhagic, sepsis and tail-rot disease are all caused by bacteria as general in China. Most of the pathogenic bacteria are Gram-negative bacteria. Liver transcriptome analysis of the co-injection of cortisol and lipopolysaccharide (LPS) was performed in this study. Preliminary evidence from the results suggest that after the emergence of immune stress, cortisol will up-regulate the complement cascade pathway, down-regulate the coagulation cascade pathway, down-regulate the platelet activation pathway, down-regulate antigen presentation pathway, and show complex regulation relationship to inflammatory factors. At 12 h, the number of differential genes regulated by cortisol was about half less than the number of differential genes regulated by LPS. At 24 h, there was no significant difference between the number of differential genes regulated by cortisol and LPS, but the types of differential genes vary widely. KEGG enrichment analysis found that cortisol regulated LPS-stimulated immune responses mainly focus on cytokines, complement and coagulation cascades pathways, antigen presentation pathways, haematopoiesis, and inflammation. It is suggested that there may be some strategic choice in the regulation of immune response by cortisol. These results will help understand the pathogenesis and host defence system in bacterial disease caused by Gram-negative bacteria.

Identification of Multigene Biomarker for Shrimp White Feces Syndrome by Full-Length Transcriptome Sequencing

The pacific white shrimp, Litopenaeus vannamei, with the largest shrimp industry production in the world, is currently threatened by a severe disease, white feces syndrome (WFS), which cause devastating losses globally, while its causal agents remain largely unknown. Herein, compared to the Control shrimp by metagenomic analysis, we firstly investigated that the altered functions of intestinal microbial community in WFS shrimp were the enrichment of bacterial chemotaxis and flagellar assembly pathways, hinting at a potential role of pathogenic bacteria for growth and development, which might be related to WFS occurrence. Single-molecule real-time (SMRT) sequencing was to further identify the gene structure and gene regulation for more clues in WFS aetiology. Totally 50,049 high quality transcripts were obtained, capturing 39,995 previously mapped and 10,054 newly detected transcripts, which were annotated to 30,554 genes. A total of 158 differentially expressed genes (DEGs) were characterized in WFS shrimp. These DEGs were strongly associated with various immune related genes that regulated the expression of multiple antimicrobial peptides (e.g., antilipopolysaccharide factors, penaeidins, and crustin), which were further experimentally validated using quantitative PCR on transcript level. Collectively, multigene biomarkers were identified to be closely associated with WFS, especially those functional alterations in microbial community and the upregulated immune related gene with antibacterial activities. Our finding not only inspired our cogitation on WFS aetiology from both microbial and host immune response perspectives with combined metagenomic and full-length transcriptome sequencing, but also provided valuable information for enhancing shrimp aquaculture.

Transcriptome analysis of human brain microvascular endothelial cells response to Neisseria meningitidis and its antigen MafA using RNA-seq

Interaction of Neisseria meningitidis (NM) with human brain microvascular endothelial cells (hBMECs) initiates of multiple cellular processes, which allow bacterial translocation across the blood-brain barrier (BBB). NM is equipped with several antigens, which interacts with the host cell receptors. Recently we have shown that adhesin MafA (UniProtKB-X5EG71), relatively less studied protein, is one of those surface exposed antigens that adhere to hBMECs. The present study was designed to comprehensively map the undergoing biological processes in hBMECs challenged with NM or MafA using RNA sequencing. 708 and 726 differentially expressed genes (DEGs) were identified in hBMECs exposed to NM and MafA, respectively. Gene ontology analysis of the DEGs revealed that several biological processes, which may alter the permeability of BBB, were activated. Comparative analysis of DEGs revealed that MafA, alike NM, might provoke TLR-dependent pathway and augment cytokine response. Moreover, both MafA and NM were able to induce genes involved in cell surface modifications, endocytosis, extracellular matrix remodulation and anoikis/apoptosis. In conclusion, this study for the first time describes effect of NM on the global gene expression in hBMECs using high-throughput RNA-seq. It also presents ability of MafA to induce gene expression, which might aid NM in breaching the BBB.

Genome-Wide Association Mapping and Gene Expression Analyses Reveal Genetic Mechanisms of Disease Resistance Variations in Cynoglossus semilaevis

The sustainable development of aquaculture has been impeded by infectious diseases worldwide. However, the genomic architecture and the genetic basis underlying the disease resistance remain poorly understood, which severely hampers both the understanding of the evolution of fish disease resistance traits and the prevention of these diseases in the aquaculture community. Cynoglossus semilaevis is a representative and commercially-important flatfish species. Here we combined genome-wide association study and Fst and nucleotide diversity filtration to identify loci important for the disease resistance. Based on 1,016,774 single-nucleotide polymorphisms (SNPs) identified from 650 Gb genome resequencing data of 505 individuals, we detected 33 SNPs significantly associated with disease resistance and 79 candidate regions after filtration steps. Both the allele frequencies and genotype frequencies of the associated loci were significantly different between the resistant and susceptible fish, suggesting a role in the genetic basis of disease resistance. The SNP with strongest association with disease resistance was located in Chr 17, at 145 bp upstream of fblx19 gene, and overlapped with the major quantitative trait locus previously identified. Several genes, such as plekha7, nucb2, and fgfr2, were also identified to potentially play roles in the disease resistance. Furthermore, the expression of some associating genes were likely under epigenetic regulations between the bacterial resistant and susceptible families. These results provide insights into the mechanism that enable variation of disease resistance to bacterial pathogen infection. The identified polymorphisms and genes are valuable targets and molecular resources for disease resistance and other traits, and for advanced breeding practice for superior germplasm in fish aquaculture.

Transcriptomic characterization of adult zebrafish infected with Streptococcus agalactiae

Streptococcus agalactiae is a major aquaculture pathogen infecting various saltwater and freshwater fish. To better understand the mechanism of the immune responses to S. agalactiae in wildtype zebrafish, the transcriptomic profiles of two organs containing mucosal-associated lymphoid tissues from S. agalactiae-infected and non-infected groups were obtained using RNA-seq techniques. In the intestines, 6735 and 12908 differently expressed genes (DEGs) were identified at 24 hpi and 48 hpi, respectively. Among 66 and 116 significantly enriched pathways, 15 and 21 pathways were involved in immune system or signal transduction at 24 hpi and 48 hpi, respectively. A number of genes involved in Toll-like receptor signaling pathway, RIG-I-like receptor signaling pathway, NOD-like receptor signaling pathway, T cell receptor signaling pathway, B cell receptor signaling pathway, Antigen processing and presentation, NF-kappa B signaling pathway and PI3K-Akt signaling pathway were significantly downregulated. In the skins, 3113 and 4467 DEGs were identified at 24 hpi and 48 hpi, respectively. Among 24 and 56 significantly enriched pathways, 4 and 13 pathways were involved in immune system or signal transduction at 24 hpi and 48 hpi, respectively. More immune-related signaling pathways including Leukocyte transendothelial migration, Cytokine-cytokine receptor interaction, PI3K-Akt signaling pathway, IL-17 signaling pathway, MAPK signaling pathway, TNF signaling pathway, Complement and coagulation cascades, Hematopoietic cell lineage and Jak-STAT signaling pathway were differently enriched for upregulated DEGs at 48 hpi, which were completely different from that in the intestines. Furthmore, comparative transcriptome analysis revealed that the downregulated 1618 genes and upregulated 1622 genes existed both at 24 hpi and 48 hpi for the intestine samples. In the skins, the downregulated 672 genes and upregulated 428 genes existed both at 24 hpi and 48 hpi. Three pathways related to immune processes were significantly enriched for downregulated DEGs both in the intestines and skins collected at 24 hpi and 48 hpi, which included Antigen processing and presentation, Intestinal immune network for IgA production and Hematopoietic cell lineage. Interaction network analysis of DEGs identified the main DEGs in the sub-network of complement and coagulation cascades both in the intestines and skins. Twenty of DEGs involved in complement and coagulation cascades were further validated by Real-time quantitative PCR. Altogether, the results obtained in this study will provide insight into the immune response of zebrafish against S. agalactiae XQ-1 infection in fatal conditions, and reveal the discrepant expression pattern of complement and coagulation cascades in the intestines and skins.

Identification of Potential Immune-Related circRNA-miRNA-mRNA Regulatory Network in Intestine of Paralichthys olivaceus During Edwardsiella tarda Infection

Olive flounder (Paralichthys olivaceus) is an important economical flatfish in Japan, Korea, and China, but its production has been greatly threatened by disease outbreaks. In this research, we aimed to explore the immune responsive mechanism of P. olivaceus against Edwardsiella tarda infection by profiling the expression of circRNA, miRNA, and mRNA by RNA-seq and constructing a regulatory circular circRNA–miRNA–mRNA network. Illumina sequencing of samples from normal control (H0), 2 h (H2), 8 h (H8), and 12 h (H12) post-challenge was conducted. Differentially expressed (DE) circRNA (DE–circRNAs), miRNAs (DE–miRNAs), and mRNAs [differential expression genes (DEGs)] between challenge and control groups were identified, resulting in a total of 62 DE–circRNAs, 39 DE–miRNAs, and 3,011 DEGs. Based on the differentially expressed gene results, miRNA target interactions (circRNA–miRNA pairs and miRNA–mRNA pairs) were predicted by MiRanda software. Once these paired were combined, a preliminary circRNA–miRNA–mRNA network was generated with 198 circRNA–miRNA edges and 3,873 miRNA–mRNA edges, including 44 DE–circRNAs, 32 DE–miRNAs, and 1,774 DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed to evaluate the function of the DEGs in this network, and we focused and identified two important intestinal immune pathways (herpes simplex infection and intestinal immune network for IgA production) that showed statistical significance between the challenge and control groups. Furthermore, three critical DEGs (nectin2, MHC II α-chain, and MHC II β-chain) were identified, mapped into the preliminary circRNA–miRNA–mRNA network, and new circRNA–miRNA–mRNA regulatory networks were constructed. In conclusion, we, for the first time, identified circRNA–miRNA–mRNA network from P. olivaceus in the pathogenesis of E. tarda and provided valuable resources for further analyses of the molecular mechanisms and signaling networks.

FOXO genes in channel catfish and their response after bacterial infection

FOXO proteins are a subgroup of the forkhead family of transcription factors that play crucial roles in lifespan regulation. In addition, FOXO proteins are also involved in immune responses. After a systematic study of FOXO genes in channel catfish, Ictalurus punctatus, seven FOXO genes were identified and characterized, including FOXO1a, FOXO1b, FOXO3a, FOXO3b, FOXO4, FOXO6a and FOXO6b. Through phylogenetic and syntenic analyses, it was found that FOXO1, FOXO3 and FOXO6 were duplicated in the catfish genome, as in the zebrafish genome. Analysis of the relative rates of nonsynonymous (dN) and synonymous (dS) substitutions revealed that the FOXO genes were globally strongly constrained by negative selection. Differential expression patterns were observed in the majority of FOXO genes after Edwardsiella ictaluri and Flavobacterium columnare infections. After E. ictaluri infection, four FOXO genes with orthologs in mammal species were significantly upregulated, where FOXO6b was the most dramatically upregulated. However, after F. columnare infection, the expression levels of almost all FOXO genes were not significantly affected. These results suggested that either a pathogenesis-specific pattern or tissue-specific pattern existed in catfish after these two bacterial infections. Taken together, these findings indicated that FOXO genes may play important roles in immune responses to bacterial infections in catfish.

Comparative transcriptome analysis reveals the different roles between hepatopancreas and intestine of Litopenaeus vannamei in immune response to aflatoxin B1 (AFB1) challenge

Aflatoxin B1 (AFB1) is a mycotoxin mainly produced by Aspergillus flavus and Aspergillus parasiticus contaminating food, feed ingredients and products of animal origin. In mammals, this toxin causes widespread organ-specific damage; it is immunotoxicity and could promote hepatotoxicity, alter intestinal functions and so on. In this study, we conducted transcriptome and histomorphology analyses of hepatopancreas and intestinal in Litopenaeus vannamei (L. vannamei) challenged with AFB1. Totally 12,014 and 1387 differentially expression genes (DEGs) were identified in the hepatopancreas and intestine, respectively. In hepatopancreas, a total of 1995 DEGs were mainly annotated and grouped into 18 processes or pathways related to animal immune system. With respect to intestine, a total of 152 DEGs were mainly annotated to 7 processes or pathways related to animal immune system. Meanwhile, we determined the relative mRNA expression of several crucial representative immune genes including Toll, immune deficiency (IMD), prophenoloxidase (proPO), Rab and glutathione S-transferase (GST) in the hepatopancreas and intestines of shrimp at 3-, 6-, 12-, 18-, 24- and 30-d after challenged by AFB1. Exposure to AFB1 increased mortality, decrease weight gain rate, severely destroyed the histomorphology of hepatopancreas and intestine, and resulted in the damaged of immune system of shrimp. The present data reveals the different roles between hepatopancreas and intestine of L. vannamei in immune response to AFB1 challenge, and provides insight into the molecular basis of the relationship between hepatopancreas and intestinal immunity during either homeostasis or inflammation.

Modulation of posterior intestinal mucosal proteome in rainbow trout (Oncorhynchus mykiss) after Yersinia ruckeri infection

Yersinia ruckeri is the causative agent of enteric redmouth disease in salmonids. In fish, the intestine represents an important site of nutrient uptake, host–pathogen interactions, and defense. The posterior intestine can be inflamed, reddened, and filled with an opaque, yellowish fluid during Y. ruckeri infection. Herein, we report an investigation on the proteome alteration in the posterior intestinal mucosa of rainbow trout (Oncorhynchus mykiss) after exposure to Y. ruckeri. The intestinal mucosal proteins were identified and quantified by a shotgun proteomic approach by applying data-independent quantification with sequential windowed acquisition of all theoretical mass spectra (SWATH). A total of 437 proteins were found to be differentially up- or downregulated in the posterior intestine. Gene ontology of upregulated proteins pointed to their involvement into exopeptidase, endopeptidase, and hydrolase activities, while the downregulated proteins were involved in lipid metabolism, actin binding, and translation processes. Additionally, upregulated proteins were predicted to be involved in lysosome, oxidative phosphorylation, and metabolic pathways, while downregulated proteins were implicated in focal adhesion, regulation of actin cytoskeleton, protein digestion and absorption pathways. This study showed that Y. ruckeri infection can alter protein abundance involved in serine-type carboxypeptidase, cysteine and aspartic-type endopeptidases, metallopeptidases, antioxidant defense, calcium ion binding, glycolytic and carbohydrate metabolic processes in the proteome of the intestinal mucosa of rainbow trout.

Evaluation of a Recombinant Flavobacterium columnare DnaK Protein Vaccine as a Means of Protection Against Columnaris Disease in Channel Catfish (Ictalurus punctatus)

Flavobacterium columnare causes substantial losses among cultured finfish species. The Gram-negative bacterium is an opportunistic pathogen that manifests as biofilms on the host's mucosal surfaces as the disease progresses. We previously demonstrated that the dominant mucosal IgM antibody response to F. columnare is to the chaperone protein DnaK that is found in the extracellular fraction. To establish the efficacy of using recombinant protein technology to develop a new vaccine against columnaris disease, we are reporting on two consecutive years of vaccine trials using a recombinant F. columnare DnaK protein (rDnaK). In year one, three groups of channel catfish (n = 300) were immunized by bath immersion with a live attenuated F. columnare isolate, rDnaK or sham immunized. After 6 weeks, an F. columnare laboratory challenge showed a significant increase in survival (>30%) in both the live attenuated and rDnaK vaccines when compared to the non-immunized control. A rDnaK-specific ELISA revealed significant levels of mucosal IgM antibodies in the skin of catfish immunized with rDnaK at 4- and 6-weeks post immunization. In the second year, three groups of channel catfish (n = 300) were bath immunized with rDnaK alone or with rDnaK after a brief osmotic shock or sham immunized. After 6 weeks a laboratory challenge with F. columnare was conducted and showed a significant increase in survival in the rDnaK (> 25%) and in rDnaK with osmotic shock (>35%) groups when compared to the non-immunized control. The rDnaK-specific ELISA demonstrated significant levels of mucosal IgM antibodies in the skin of catfish groups immunized with rDnaK at 4- and 6-weeks post immunization. To further understand the processes which have conferred immune protection in the rDnaK group, we conducted RNA sequencing of skin samples from the non-immunized (n = 6) and rDnaK treated channel catfish at 1-week (n = 6) and 6 weeks (n = 6) post immunization. Significantly altered gene expression was identified and results will be discussed. Work to further enhance the catfish immune response to F. columnare rDnaK is underway as this protein remains a promising candidate for additional optimization and experimental trials in a production setting.

Transcriptome analysis of the spleen provides insight into the immunoregulation of Mastacembelus armatus under Aeromonas veronii infection

Mastacembelus armatus, also known as the zigzag eel, is an economically important species of freshwater fish that is very popular with consumers as a high-grade table fish in China. Recently, the wild population of this fish has declined gradually due to overfishing and various types of ecological imbalance. Meanwhile, the aquaculture of this spiny eel has flourished in southern China. To understand the immune response of zigzag eel to Aeromonas veronii, we carried out transcriptome sequencing of zigzag eel spleens after artificial bacterial infection. After assembly, 110,328 unigenes were obtained with 44.42% GC content. A total of 27,098 unigenes were successfully annotated by four public protein databases, namely, Nr, UniProt, KEGG and KOG. Differential expression analysis revealed the existence of 1278 significantly differentially expressed unigenes at 24 h post infection, with 767 unigenes upregulated and 511 unigenes downregulated. After GO and KEGG enrichment analyses, many immune-related GO categories and pathways were significantly enriched. The typical significantly enriched pathways included toll-like receptor signaling pathway, cytokine-cytokine receptor interaction and TNF signaling pathway. In addition, 40,027 microsatellites (SSRs) and 52,716 candidate single nucleotide polymorphisms (SNPs) were identified from the infection and control transcriptome libraries. Overall, this transcriptomic analysis provided valuable information for studying the immune response of zigzag eels against bacterial infection.

Polyadenylation sites and their characteristics in the genome of channel catfish (Ictalurus punctatus) as revealed by using RNA-Seq data

Polyadenylation plays important roles in gene expression regulation in eukaryotes, which typically involves cleavage and poly(A) tail addition at the polyadenylation site (PAS) of the pre-mature mRNA. Many eukaryotic genes contain more than one PASs, termed as alternative polyadenylation (APA). As a crucial post-transcriptional regulation, polyadenylation affects various aspects of RNA metabolism such as mRNA stability, translocation, and translation. However, polyadenylation has been rarely studied in teleosts. Here we conducted polyadenylation analysis in channel catfish, a commercially important aquaculture species around the world. Using RNA-Seq data, we identified 20,320 PASs which were classified into 14,500 clusters by merging adjacent PASs. Most of the PASs were found in 3' UTRs, followed by intron regions based on the annotation of channel catfish reference genome. No apparent difference in PAS distribution was observed between the sense and antisense strand of the channel catfish genome. The sequence analysis of nucleotide composition and motif around PASs yielded a highly similar profile among various organisms, suggesting the conservation and importance of polyadenylation in evolution. Using APA genes with more than two PASs, gene ontology enrichment revealed genes particularly involved in RNA binding. Reactome pathway analysis showed the enrichment of the innate immune system, especially neutrophil degranulation.

Novel insights into the immune regulatory effects of ferritins from blunt snout bream, Megalobrama amblycephala

Ferritins play vital roles in maintenance of iron homeostasis as iron storage proteins. Recently, the immune function of ferritins have attracted increasing attention, especially their roles in defense against pathogenic infections. However, the immune regulatory mechanism of fish ferritins are not well known. In the present study, comparative digital gene expression (DGE) profiling was performed to explore the regulatory effects of the Megalobrama amblycephala ferritins (MamFers) using MamFers overexpressed and control L8824 cells (Ctenopharyngodon idella hepatic cell line). Clean reads were aligned to the C. idella genome and differential expression analysis was conducted with representative differentially expressed genes pointed out. On that basis, further studies were performed to verify two pivotal regulated pathways in L8824 and EPC (Epithelioma Papulosum Cyprini cell line) cells, respectively. The results showed that NLRC5 (NOD-like Receptor Family CARD Domain Containing 5) mediated the regulation of MamFers on expression of MHC I (Major Histocompatibility Complex Class I) and its chaperone β2M (Beta-2-Microglobulin) in L8824 cells. Then, β2M further mediated the regulation of MamFers on hepcidin expression, indicating that MamFers regulated the expression of hepcidin via NLRC5/MHC I/β2M axis. In addition, MamFers regulated the adhesion of Aeromonas hydrophila to EPC cells by regulating the expression of two extracellular matrix proteins Intgβ1 (integrin β1) and FN (fibronectin). In a word, the present study provided novel insights into the immune regulatory functions of fish ferritins.

Transcriptome profiles in the spleen of African catfish (Clarias gariepinus) challenged with Aeromonas veronii

The African catfish, Clarias gariepinus, an important cultured freshwater species in many countries, possess the characteristic of high disease resistance. However, little genomic information for this character of the fish is available up to now. To address the shortfall and to better understand C. gariepinus immune response to pathogen infection at molecular level, C. gariepinus were challenged with potent A. veronii and the high-throughput RNA sequencing (RNA-seq) technology were employed to produce transcriptomes from spleen. In total, an average of 46,073,372 clean reads obtained were de novo assembled into 156,955 unigenes with an average length of 1082 bp. All of unigenes were annotated to seven public databases. Three comparisons were separately conducted between the infected groups at 3 h, 24 h, 48 h post-challenge and control group. A total of 2482 differentially expressed unigenes (DEGs) were identified. Among these, 114 immune-related DEGs were captured, including 88, 42, and 31 genes at 3 h, 24 h and 48 h after infection respectively, for analysis of expression pattern and enrichment. The 114 DEGs displayed four expression patterns by cluster analysis and they were significantly enriched in 38 pathways (q < 0.01) related to the immune or disease, five of which were NF-kappa B, TNF, NLR, TLR and RLR pathways. Finally, the expression levels of twelve selected immune-related DEGs involved in above five pathways were scrutinized. Seven of which were up-regulated at 3 h after infection, afterward, their expression dropped to control level. In summary, this study provides valuable transcriptome resource for understanding the defense mechanisms of C. gariepinus in resistance to pathogens from the gene expression viewpoint, which also open up the possibility to study the immune complexity and to better comprehend the interrelationships between some immune pathways in C. gariepinus.

Detection and quantification of Aeromonas schubertii in Channa maculata by TaqMan MGB probe fluorescence real-time quantitative PCR

Aeromonas schubertii is a major epidemiological agent that threatens cultured snakeheads (Channidae) and has caused great economic losses in fish-farming industries in China in recent years. In present study, a specific TaqMan minor groove binder (MGB) probe fluorescence real-time quantitative PCR (qPCR) assay was developed to rapidly detect and quantify A. schubertii. A pair of qPCR primers and a TaqMan MGB probe were selected from the rpoD gene, which were shown to be specific for A. schubertii. A high correlation coefficient (R²  = 0.9998) in a standard curve with a 103% efficiency was obtained. Moreover, the qPCR method's detection limit was as low as 18 copies/μl, which was 100 times more sensitive than that of conventional PCR. The detection results for the A. schubertii in pond water and fish tissue were consistent with those of the viable counts. Bacterial load changes detected by qPCR in different tissues of snakeheads infected with A. schubertii showed that the gills and intestines may be the entry for A. schubertii, and the spleen and kidney are major sites for A. schubertii replication. The established method in present study should be a useful tool for the early surveillance and quantitation of A. schubertii.

Increased Alternative Splicing as a Host Response to Edwardsiella ictaluri Infection in Catfish

Alternative splicing is the process of generating multiple transcripts from a single pre-mRNA used by eukaryotes to regulate gene expression and increase proteomic complexity. Although alternative splicing profiles have been well studied in mammalian species, they have not been well studied in aquatic species, especially after biotic stresses. In the present study, genomic information and RNA-Seq datasets were utilized to characterize alternative splicing profiles and their induced changes after bacterial infection with Edwardsiella ictaluri in channel catfish (Ictalurus punctatus). A total of 27,476 alternative splicing events, derived from 9694 genes, were identified in channel catfish. Exon skipping was the most abundant while mutually exclusive exon was the least abundant type of alternative splicing. Alternative splicing was greatly induced by E. ictaluri infection with 21.9% increase in alternative splicing events. Interestingly, genes involved in RNA binding and RNA splicing themselves were significantly enriched in differentially alternatively spliced genes after infection. Sequence analyses of splice variants of a representative alternatively spliced gene, splicing factor srsf2, revealed that certain spliced transcripts may undergo nonsense-mediated decay (NMD), suggesting functional significance of the induced alternative splicing. Although statistical analysis was not possible with such large datasets, results from quantitative real-time PCR from representative differential alternative splicing events provided general validation of the bacterial infection-induced alternative splicing. This is the first comprehensive study of alternative splicing and its changes in response to bacterial infection in fish species, providing insights into the molecular mechanisms of host responses to biotic stresses.

A Review of Molecular Responses of Catfish to Bacterial Diseases and Abiotic Stresses

Catfish is one of the major aquaculture species in the United States. However, the catfish industry is threatened by several bacterial diseases such as enteric septicemia of catfish (ESC), columnaris disease and Aeromonas disease, as well as by abiotic stresses such as high temperature and low oxygen. Research has been conducted for several decades to understand the host responses to these diseases and abiotic stresses. With the development of sequencing technologies, and the application of genome-wide association studies in aquaculture species, significant progress has been made. This review article summarizes recent progress in understanding the molecular responses of catfish after bacterial infection and stress challenges, and in understanding of genomic and genetic basis for disease resistance and stress tolerance.

Enterotoxic effects of Aeromonas hydrophila infection in the catfish, Clarias gariepinus: Biochemical, histological and proteome analyses

Aeromonas hydrophila is considered as a potential risk to fish populations in the aquaculture industry and could also pose a serious threat to humans. In this study, the impact of A. hydrophila infection in the air-breathing catfish, Clarias gariepinus was analyzed using a multidimensional approach. Aeromonas hydrophila (1 × 10⁷ cells) was injected into C. gariepinus intraperitoneally and maintained at an ambient temperature and photoperiod with periodical monitoring for morphological changes. After 7 days post-infection, tissue samples of the gills, liver, intestine, and kidney were subjected to biochemical, histological, transmission electron microscope (TEM) and proteomic analyses. Observed results indicated distinct morphological changes with the significant increase of ROS and oxidative stress enzymes (CAT and SOD) in tissues of the infected group when compared to the control. Histological analysis in infected fish revealed the presence of pyknotic nuclei, early stages of necrosis in the liver, degradation of renal tubules and widened sinusoidal space in kidneys along with enlargement of the epithelial region in the intestine. TEM analysis of the infected intestine showed degeneration of villi and the presence of multinucleated erythrocytes. Two-dimensional proteomic and mass spectrometry analysis of intestine and liver displayed up-regulation of several immune regulatory proteins such as proteasome subunit 3 protein, prolactin and intermediated filament protein; and down-regulation of proteins including actin, serine/arginine-rich splicing factor and carbonic anhydrase. Taken together, these results suggest that the identified proteins may have a role in immune regulation against A. hydrophila infection in C. gariepinus and support further investigations of host-pathogen interactions.

Phytogenic feed-additive effects on channel catfish rhamnose-binding lectin levels, and susceptibility to Edwardsiella ictaluri

We investigated the effects of a phytogenic feed additive on disease susceptibility to Edwardsiella ictaluri in channel catfish Ictalurus punctatus and regulation of 6 rhamnose-binding lectin (RBL) genes. Juvenile catfish (n = 250, 13.4 ± 0.1 g) were allotted to the following treatments: control (floating diet) or EO (floating diet supplemented with essential oils; Digestarom® P.E.P. MGE). The fish were fed their respective diets for 6 wk. Following subjection to different feed treatments, all fish were exposed to pathogenic E. ictaluri by bath immersion. Another group of fish were not challenged (non-challenged controls, fed control feed). Mucosal tissue samples were taken to quantify gene expression levels of RBL on Days 1 and 2 post-challenge. After challenge, survival was higher (64.4 vs. 48.0%) in fish fed EO compared to controls (p < 0.05). Relative to non-challenged controls, gill RBL1a mRNA was higher in fish fed EO (p < 0.05) on Day 1 while gill RBL3b was higher in fish fed EO (p < 0.01) on Days 1 and 2, respectively. RBL5a in the skin and proximate small intestine did not change significantly relative to non-challenged fish on Days 1 and 2 of the disease challenge. Results demonstrate that Digestarom® P.E.P. MGE improved survival of channel catfish challenged with E. ictaluri. One of the mechanisms through which essential oils may improve survival is through upregulation of RBL1a and RBL3b in the gill.

Genome-wide association analysis of intra-specific QTL associated with the resistance for enteric septicemia of catfish

Disease resistance is one of the most important traits for aquaculture industry. For catfish industry, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the most severe disease, causing enormous economic losses every year. In this study, we used three channel catfish families with 900 individuals (300 fish per family) and the 690K catfish SNP array, and conducted a genome-wide association study to detect the quantitative trait loci (QTL) associated with ESC resistance. Three significant QTL, with two of located on LG1 and one on LG26, and three suggestive QTL located on LG1, LG3, and LG21, respectively, were identified to be associated with ESC resistance. With a well-assembled- and -annotated reference genome sequence, genes around the involved QTL regions were identified. Among these genes, 37 genes had known functions in immunity, which may be involved in ESC resistance. Notably, nlrc3 and nlrp12 identified here were also found in QTL regions of ESC resistance in the channel catfish × blue catfish interspecific hybrid system, suggesting this QTL was operating within both intra-specific channel catfish populations and interspecific hybrid backcross populations. Many of the genes of the Class I MHC pathway, for mediated antigen processing and presentation, were found in the QTL regions. The positional correlation found in this study and the expressional correlation found in previous studies indicated that Class I MHC pathway was significantly associated with ESC resistance. This study validated one QTL previously identified using the second and fourth generation of the interspecific hybrid backcross progenies, and identified five additional QTL among channel catfish families. Taken together, it appears that there are only a few major QTL for ESC disease resistance, making marker-assisted selection an effective approach for genetic improvements of ESC resistance.

Mucus-Pathogen Interactions in the Gastrointestinal Tract of Farmed Animals

Gastrointestinal infections cause significant challenges and economic losses in animal husbandry. As pathogens becoming resistant to antibiotics are a growing concern worldwide, alternative strategies to treat infections in farmed animals are necessary in order to decrease the risk to human health and increase animal health and productivity. Mucosal surfaces are the most common route used by pathogens to enter the body. The mucosal surface that lines the gastrointestinal tract is covered by a continuously secreted mucus layer that protects the epithelial surface. The mucus layer is the first barrier the pathogen must overcome for successful colonization, and is mainly composed of densely glycosylated proteins called mucins. The vast array of carbohydrate structures present on the mucins provide an important setting for host-pathogen interactions. This review summarizes the current knowledge on gastrointestinal mucins and their role during infections in farmed animals. We examine the interactions between mucins and animal pathogens, with a focus on how pathogenic bacteria can modify the mucin environment in the gut, and how this in turn affects pathogen adhesion and growth. Finally, we discuss analytical challenges and complexities of the mucus-based defense, as well as its potential to control infections in farmed animals.

Applications of transcriptomics and proteomics in understanding fish immunity

With the development of intensive aquaculture, economic losses increasingly result from fish mortality due to pathogen infection. In recent years, a growing number of researchers have used transcriptomic and proteomic analyses to study fish immune responses to exogenous pathogen infection. Integrating transcriptomic and proteomic analyses provides a better understanding of the fish immune system including gene expression, regulation, and the intricate biological processes underlying immune responses against infection. This review focuses on the recent advances in the fields of transcriptomics and proteomics, which have contributed to our understanding of fish immunity to exogenous pathogens.

Comparison of Channel Catfish and Blue Catfish Gut Microbiota Assemblages Shows Minimal Effects of Host Genetics on Microbial Structure and Inferred Function

The microbiota of teleost fish has gained a great deal of research attention within the past decade, with experiments suggesting that both host-genetics and environment are strong ecological forces shaping the bacterial assemblages of fish microbiomes. Despite representing great commercial and scientific importance, the catfish within the family Ictaluridae, specifically the blue and channel catfish, have received very little research attention directed toward their gut-associated microbiota using 16S rRNA gene sequencing. Within this study we utilize multiple genetically distinct strains of blue and channel catfish, verified via microsatellite genotyping, to further quantify the role of host-genetics in shaping the bacterial communities in the fish gut, while maintaining environmental and husbandry parameters constant. Comparisons of the gut microbiota among the two catfish species showed no differences in bacterial species richness (observed and Chao1) or overall composition (weighted and unweighted UniFrac) and UniFrac distances showed no correlation with host genetic distances (Rst) according to Mantel tests. The microbiota of environmental samples (diet and water) were found to be significantly more diverse than that of the catfish gut associated samples, suggesting that factors within the host were further regulating the bacterial communities, despite the lack of a clear connection between microbiota composition and host genotype. The catfish gut communities were dominated by the phyla Fusobacteria, Proteobacteria, and Firmicutes; however, differential abundance analysis between the two catfish species using analysis of composition of microbiomes detected two differential genera, Cetobacterium and Clostridium XI. The metagenomic pathway features inferred from our dataset suggests the catfish gut bacterial communities possess pathways beneficial to their host such as those involved in nutrient metabolism and antimicrobial biosynthesis, while also containing pathways involved in virulence factors of pathogens. Testing of the inferred KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways by DESeq2 revealed minor difference in microbiota function, with only two metagenomic pathways detected as differentially abundant between the two catfish species. As the first study to characterize the gut microbiota of blue catfish, our study results have direct implications on future ictalurid catfish research. Additionally, our insight into the intrinsic factors driving microbiota structure has basic implications for the future study of fish gut microbiota.

GWAS analysis using interspecific backcross progenies reveals superior blue catfish alleles responsible for strong resistance against enteric septicemia of catfish

Infectious diseases pose significant threats to the catfish industry. Enteric septicemia of catfish (ESC) caused by Edwardsiella ictaluri is the most devastating disease for catfish aquaculture, causing huge economic losses annually. Channel catfish and blue catfish exhibit great contrast in resistance against ESC, with channel catfish being highly susceptible and blue catfish being highly resistant. As such, the interspecific backcross progenies provide an ideal system for the identification of quantitative trait locus (QTL). We previously reported one significant QTL on linkage group (LG) 1 using the third-generation backcrosses, but the number of founders used to make the second- and third-generation backcross progenies was very small. Although the third-generation backcross progenies provided a greater power for fine mapping than the first-generation backcrosses, some major QTL for disease resistance may have been missing due to the small numbers of founders used to produce the higher generation backcrosses. In this study, we performed a genome-wide association study using first-generation backcrosses with the catfish 690 K SNP arrays to identify additional ESC disease resistance QTL, especially those at the species level. Two genomic regions on LG1 and LG23 were determined to be significantly associated with ESC resistance as revealed by a mixed linear model and family-based association test. Examination of the resistance alleles indicated their origin from blue catfish, indicating that at least two major disease resistance loci exist among blue catfish populations. Upon further validation, markers linked with major ESC disease resistance QTL should be useful for marker-assisted introgression, allowing development of highly ESC resistant breeds of catfish.

JAK and STAT members in channel catfish: Identification, phylogenetic analysis and expression profiling after Edwardsiella ictaluri infection

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway is one of the main pleiotropic cascades used to transmit information from extracellular receptors to the nucleus, which results in DNA transcription and expression of genes involved in immunity, proliferation, differentiation, migration, apoptosis, and cell survival. Members of JAK family and STAT family have been extensively studied in different mammalian species because of their important roles in innate and adaptive immune responses. However, they have not been systematically studied among teleost fish species. In this study, five JAK family members and eight STAT family members were identified and characterized from channel catfish. Phylogenetic analysis was conducted to properly annotate these genes. Syntenic analysis was also conducted to establish orthology, and confirm the results from phylogenetic analysis. Compared to mammals, more members of the JAK and STAT family were identified in channel catfish genome. Expression of JAK and STAT family members was detected in healthy catfish tissues, but was induced in gill, liver, and intestine after bacterial challenge. Notably, the significant upregulation of STAT1b gene in catfish liver, gill and intestine after Edwardsiella ictaluri infection supported the notion that high STAT1 expression are involved in defense against pathogens. Collectively, the increased expression of JAK and STAT members in tested tissues suggested their crucial function in defending the host against pathogen invasion.