Identification and characterization of caspases genes in rainbow trout (Oncorhynchus mykiss) and their expression profiles after Aeromonas salmonicida and Vibrio anguillarum infection

Characterization of caspase gene family in Sebastes schlegelii and their expression profiles under Aeromonas salmonicida and Vibrio anguillarum infection

The caspase, functioning as a proteinase, plays a crucial role in eukaryotic cell apoptosis, regulation of apoptosis, cellular growth, differentiation, and immunity. The identification of caspase gene family in Sebastes schlegelii is of great help to understand its antimicrobial research. In S. schlegelii, we totally identified nine caspase genes, including four apoptosis initiator caspases (caspase 2, caspase 8, caspase 9 and caspase 10), four apoptosis executioners (caspase 3a, caspase 3b, caspase 6, and caspase 7) and one inflammatory executioner (caspase 1). The duplication of caspase 3 genes on chr3 and chr8 may have been facilitated by whole genome duplication (WGD) events or other complex evolutionary processes. In general, the number of caspase genes relatively conserved in high vertebrates, while exhibiting variation in teleosts. Furthermore, syntenic analysis and phylogenetic relationships analysis supported the correct classification of these caspase gene family in S. schlegelii, especially for genes with duplicated copies. Additionally, the expression patterns of these caspase genes in different tissues of S. schlegelii under healthy conditions were assessed. The results revealed that the expression levels of most caspase genes were significantly elevated in the intestine, spleen, and liver. To further investigate the potential immune functions of these caspase genes in S. schlegelii, we challenged individuals with A. salmonicida and V. anguillarum, respectively. After infection with A. salmonicida, the expression levels of caspase 1 in the liver and spleen of S. schlegelii remained consistently elevated throughout the infection time points. The expression levels of most caspase family members in the intestine exhibited significant divergence following V. anguillarum infection. This study provides a comprehensive understanding of the caspase gene families in S. schlegelii, thereby establishing a solid foundation for further investigations into the functional roles of these caspase genes.

Impact of different levels of handling on Solea senegalensis culture: effects on growth and molecular markers of stress

Aquaculture routine practices may cause stress induction on the fish and compromise their welfare affecting the production. This experiment aimed to evaluate the potential links between handling during culture with stress responses and growth on Senegalese sole (Solea senegalensis). We worked with two fish cohorts in terms of initial body weight and culture stage: Trial 1 included specimens in the fattening stage (226 ± 4.96 g) and Trial 2 animals in the pre-fattening stage (27.20 ± 0.44 g). The tested culture protocol, which lasted 6 and 4 months for Trial 1 and 2, respectively, mainly reduced handling-derived stressors in the experimental tanks via lowering routine samplings to a minimum. This decrease of the handling-derived stress was reflected in both trials with lower concentration of circulating cortisol in blood plasma from the experimental fish when compared to controls. Moreover, the proposed protocol promoted higher growth in the fish cultured in the less disturbing protocol in Trial 2. Higher specific growth rates and mean body weight and length were reported. In order to further explore the potential beneficial effects of our protocol, we studied the musculoskeletal from Trial 2 gene expression of key genes regulating glucocorticoid signaling pathway and apoptosis: glucocorticoid receptors 1 and 2 (gr1, gr2), heat shock protein 90 AA (hsp90aa), and caspase 6 (casp6). In line with the cortisol reduced level in this trial, gr1, hsp90aa, and casp6 genes showed lower expression in the samples coming from the experimental group. The findings of this study provide valuable information to the aquaculture industry for the management of Solea senegalensis stress and welfare.

Transcriptome analysis of response mechanism to Microcystin-LR and microplastics stress in Asian clam (Corbicula fluminea)

In this study, we analyzed the hepatopancreas tissues of Asian Clam (Corbicula fluminea) exposed to three different adverse environmental conditions from the same batch using RNA-seq. The four treatment groups included the Asian Clam group treated with Microcystin-LR (MC), the Microplastics-treated group (MP), the Microcystin-LR and Microplastics-treated group (MP-MC), and the Control group. Our Gene Ontology analysis revealed 19,173 enriched genes, and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis identified 345 related pathways. The KEGG pathway analysis demonstrated that the MC vs control group and the MP vs control group were significantly enriched in immune and catabolic pathways such as Antigen processing and presentation, Rheumatoid arthritis, Lysosome pathway, Phagosome pathway, and Autophagy pathway. We also evaluated the effects of Microplastics and Microcystin-LR on the activities of eight antioxidant enzymes and immune enzymes in Asian clams. Our study enriched the genetic resources of Asian clams and provided valuable information for understanding the response mechanism of Asian clams to microplastics and microcystin in the environment, through the identification of differentially expressed genes and related pathway analyses from the large number of transcriptome sequences obtained.

Comparative Transcriptome Analysis of Head Kidney of Aeromonas hydrophila-infected Hypoxia-tolerant and Normal Large Yellow Croaker

The large yellow croaker (Larimichthys crocea) is one of the most economically important marine fish on the southeast coast of China and much of its yield is usually lost by hypoxia. To address this problem and lay a foundation for culturing a new strain of large yellow croaker with hypoxia tolerance, our research group screened a hypoxia-tolerant population of L. crocea. Surprisingly, we also found that hypoxia-tolerant population exhibited higher survival when infected with pathogens compared to the normal population during the farming operation. In order to understand the mechanism underlying the higher survival rate of the hypoxia-tolerant population and enrich the head kidney immune mechanism of L. crocea infected with pathogens, we compared and analyzed the head kidney transcriptome of the hypoxia-tolerant and normal individuals under Aeromonas hydrophila infection. We obtained 159.68 GB high-quality reads, of which more than 87.61% were successfully localized to the reference genome of L. crocea. KEGG analysis revealed differentially expressed genes in the signaling pathways involving immunity, cell growth and death, transport and catabolism, and metabolism. Among these, the toll-like receptor signaling pathway, Nod-like receptor signaling pathway, cytokine-cytokine receptor interaction, phagosome, apoptosis, and OXPHOS pathways were enriched in both groups after infection compared to before, and were enriched in infected tolerant individuals compared to normal individuals. In addition, we found that the expression of hif1α and its downstream genes were higher in the hypoxia-sensitive group of fish than in the normal group. In conclusion, our results showed some signaling pathways and hub genes, which may participate in A. hydrophila defense in the head kidney of two populations, and may contribute to the higher survival rate in the hypoxia-tolerant population. Overall, these findings increase our understanding of the defense mechanism within the head kidney of L. crocea under A. hydrophila infection, and suggest a preliminary hypothesis for why hypoxia-tolerant individuals may exhibit a higher survival rates after infection. Our study provides scientific evidence for the breeding of a new hypoxia-tolerant strain of L. crocea for aquaculture.

Identification of three inflammatory Caspases in common carp (Cyprinus carpio L.) and its role in immune response against bacterial infection

Inflammatory Caspases are key effectors of the inflammasomes and play an important role in innate immune response. However, there are few studies on the homologs of inflammatory Caspases in bony fish. In the present study, three inflammatory Caspase genes were cloned from common carp and named CcCaspase-A1, CcCaspase-A2 and CcCaspase-B. Nucleotide sequences alignment revealed that the three Caspases were very similar in structure, which contained a PYD domain in the N-terminal, and a CASc domain in the C-terminal. In the phylogenetic tree, CcCaspase-A1 and CcCaspase-A2 were close to the Caspase-A of grass carp, and CcCaspase-B was close to the DrCaspase-B of zebrafish. In healthy common carp, the expression levels of CcCaspase-A1 and CcCaspase-A2 were the highest in the gills, and CcCaspase-B was the highest in the spleen. After immune stimulation with Edwardsiella tarda or Aeromonas hydrophila, the expression levels of all CcCaspases increased significantly. The fluorescence localization assays showed that all these CcCaspases were expressed in the cytoplasm, and were involved in the assembly of CcNLRP1 inflammasome. These results suggest that the inflammatory CcCaspases play a key role in immune response of common carp against bacterial infection, which may enrich the knowledge of inflammasome in fish, and provide basic data for the prevention and treatment of fish infectious diseases.

Time-course transcriptome analyses of spleen in rainbow trout (Oncorhynchus mykiss) post-Flavobacterium psychrophilum infection

Flavobacterium psychrophilum, the etiological agent of bacterial coldwater disease and rainbow trout fry syndrome, causes considerable losses in salmonid aquaculture globally. Systemic F. psychrophilum infections in rainbow trout (Oncorhynchus mykiss) lead to a range of clinical signs, including ulcerative lesions in the skin and muscle and splenitis. Previous studies offered an integrative analysis of the skeletal muscle response to F. psychrophilum infection in rainbow trout. However, little is known about the molecular mechanism of immune response in the spleen, which is an important immune organ of rainbow trout. Here, we investigated the time-course splenic transcriptome profiles in uninfected rainbow trout (CK) and F. psychrophilum–infected rainbow trout at day 3 and day 7 (D3, D7) by RNA-seq analyses. Among the 7,170 differentially expressed genes (DEGs) in the three comparisons (D3 vs. CK, D7 vs. CK, D3 vs. D7), 1,286 DEGs showed consistent upregulation or downregulation at D3 and D7 and were associated with pattern recognition, acute-phase response, complement cascade, chemokine and cytokine signaling, and apoptosis. The Real time quantitative PCR (RT-qPCR) analysis of eight DEGs confirmed the accuracy of the RNA-Sequencing (RNA-seq) data. Our results reflected a general process from pathogen recognition to inflammatory cytokine generation and delineated a putative Toll-like receptor signaling pathway in rainbow trout spleen, following F. psychrophilum infection. Taken together, these results provide new insights into the molecular mechanism of the immune response to F. psychrophilum infection and are a valuable resource for future research on the prevention and control of bacterial coldwater disease during salmon culture.

Crosstalk between Growth and Osmoregulation of GHRH-SST-GH-IGF Axis in Triploid Rainbow Trout (Oncorhynchus mykiss)

Smolting is an important development stage of salmonid, and an energy trade-off occurs between osmotic regulation and growth during smolting in rainbow trout (Oncorhynchus mykiss). Growth hormone releasing hormone, somatostatin, growth hormone and insulin-like growth factor (GHRH-SST-GH-IGF) axis exhibit pleiotropic effects in regulating growth and osmotic adaptation. Due to salmonid specific genome duplication, increased paralogs are identified in the ghrh-sst-gh-igf axis, however, their physiology in modulating osmoregulation has yet to be investigated. In this study, seven sst genes (sst1a, sst1b, sst2, sst3a, sst3b, sst5, sst6) were identified in trout. We further investigated the ghrh-sst-gh-igf axis of diploid and triploid trout in response to seawater challenge. Kidney sst (sst1b, sst2, sst5) and sstr (sstr1b1, sstr5a, sstr5b) expressions were changed (more than 2-fold increase (except for sstr5a with 1.99-fold increase) or less than 0.5-fold decrease) due to osmoregulation, suggesting a pleiotropic physiology of SSTs in modulating growth and smoltification. Triploid trout showed significantly down-regulated brain sstr1b1 and igfbp2a1 (p < 0.05), while diploid trout showed up-regulated brain igfbp1a1 (~2.61-fold, p = 0.057) and igfbp2a subtypes (~1.38-fold, p < 0.05), suggesting triploid trout exhibited a better acclimation to the seawater environment. The triploid trout showed up-regulated kidney igfbp5a subtypes (~6.62 and 7.25-fold, p = 0.099 and 0.078) and significantly down-regulated igfbp5b2 (~0.37-fold, p < 0.05), showing a conserved physiology of teleost IGFBP5a in regulating osmoregulation. The IGFBP6 subtypes are involved in energy and nutritional regulation. Distinctive igfbp6 subtypes patterns (p < 0.05) potentially indicated trout triggered energy redistribution in brain and kidney during osmoregulatory regulation. In conclusion, we showed that the GHRH-SST-GH-IGF axis exhibited pleiotropic effects in regulating growth and osmoregulatory regulation during trout smolting, which might provide new insights into seawater aquaculture of salmonid species.

A Comparative Review of Pyroptosis in Mammals and Fish

Pyroptosis is a form of programmed cell death, which is executed by gasdermin family proteins. Under the stimulation of pathogen- and/or damage-associated molecular patterns, pattern recognition receptors (PRRs) such as Nod like receptors could recruit apoptosis-associated speck-like protein containing a CARD (ASC) and pro-caspases to form inflammasomes and then activate caspases through various pathways. The activated caspases then cleave gasdermin family proteins, and N-terminal (NT) domains of gasdermins were released to form oligomeric pores, resulting in the increased membrane permeability, cell swelling, and final pyroptosis. During this process, caspases also promote the maturation and release of inflammatory cytokines such as IL-1β and IL-18, thus pyroptosis is also named inflammatory cell death. Unlike numerous gasdermin family proteins in mammals, only gasdermin E (GSDME) has been identified in fish. GSDME in fish can be cleaved by caspase-a/-b to release its NT domain and induce pyroptosis. Studies indicated that pyroptosis in fish mainly depends on NLR family pyrin domain-containing 3 (NLRP3) inflammasome. ASC and different caspase proteins also were identified in different fish species. The influences of pathogenic microorganism infection and environmental pollutants on fish pyroptosis were studied in recent years. Considering that fish living environment is affected by multiple factors such as water salinity, temperature, oxygen supply, and highly fluctuating food supply, the in-depth research about fish pyroptosis will contribute to revealing the mechanism of pyroptosis during evolution.

Identification and characterization of mkk genes and their expression profiles in rainbow trout (Oncorhynchus mykiss) symptomatically or asymptomatically infected with Vibrio anguillarum

Mitogen-activated protein kinase kinases (MKKs) are intermediate kinases of mitogen-activated protein kinases (MAPKs) signaling pathways. MKKs are activated by mitogen-activated protein kinase kinase kinase (MKKK) and then the activated MKKs trigger the activation of downstream MAPKs. MAPK signaling pathways play an important role in regulating immune functions including apoptosis and inflammation. However, studies on identification and characterization of mkk repertoire in rainbow trout (Oncorhynchus mykiss) are still limited. Trout experienced 4 rounds (4R) of whole genome duplication (WGD), thus exhibiting increased paralogs of mkks with potentially functional diversity. In this study, we identified 17 mkk genes in trout and the following bacterial challenge (Vibrio anguillarum) studies showed functional diversity of different mkk subtypes. Vibrio anguillarum infection resulted in significantly up-regulated mkk2 subtypes in spleen and liver, and mkk4b3 in spleen, suggesting immunomodulation was regulated by activation of ERK, p38 and JNK pathways. Compared to other mkk subtypes, mkk6s were down-regulated in symptomatic group, rather than asymptomatic group. The organisms present negative feedback on MAPK activation, thus reducing extra damage to cells. We observed down-regulated mkk6s with up-regulated genes (dusp1 & dusp2) involved in negative feedback of MAPK activation. Based on these results, we might propose the distinct expression patterns of genes associated with MAPK pathways resulted in different phenotypes and symptoms of trout in response to bacterial challenge.

Transcriptional Signatures of Immune, Neural, and Endocrine Functions in the Brain and Kidney of Rainbow Trout (Oncorhynchus mykiss) in Response to Aeromonas salmonicida Infection

Rainbow trout (Oncorhynchus mykiss) serves as one of the most important commercial fish with an annual production of around 800,000 tonnes. However, infectious diseases, such as furunculosis caused by Aeromonas salmonicida infection, results in great economic loss in trout culture. The brain and kidney are two important organs associated with “sickness behaviors” and immunomodulation in response to disease. Therefore, we worked with 60 trout and investigated transcriptional responses and enrichment pathways between healthy and infected trout. We observed that furunculosis resulted in the activation of toll-like receptors with neuroinflammation and neural dysfunction in the brain, which might cause the “sickness behaviors” of infected trout including anorexia and lethargy. We also showed the salmonid-specific whole genome duplication contributed to duplicated colony stimulating factor 1 (csf-1) paralogs, which play an important role in modulating brain immunomodulation. Enrichment analyses of kidneys showed up-regulated immunomodulation and down-regulated neural functions, suggesting an immune-neural interaction between the brain and kidney. Moreover, the kidney endocrine network was activated in response to A. salmonicida infection, further convincing the communications between endocrine and immune systems in regulating internal homeostasis. Our study provided a foundation for pathophysiological responses of the brain and kidney in response to furunculosis and potentially offered a reference for generating disease-resistant trout strains.

Immune correlates of NF-κB and TNFα promoter DNA methylation in Japanese flounder (Paralichthys olivaceus) muscle and immune parameters change response to vibrio anguillarum infection

Vibrio anguillarum infection can activate NF-κB/TNFα pathway in the immune organs of fish. Fish muscle is also an important immune organ, but the research on its immune function is few. Our aim was to study regulating mechanism of NF-κB and TNFα gene expressions in the muscle of Japanese flounder (Paralichthys olivaceus) which was under Vibrio anguillarum infection (0, 24, 48, 72 and 96 h). The results showed that the expressions of NF-κB and TNFα increased significantly at 48 h, and there was a significant positive correlation between them. In situ hybridization confirmed the co-existence of NF-κB and TNFα genes in Japanese flounder muscle. Interestingly, the expression of the TNFα gene was regulated by the DNA methylation and its methylation level was negatively correlated with the expression. The lowest methylation level of TNFα occurred at 48 h under Vibrio anguillarum infection (P < 0.05). And more, when the fragment (-2122 ∼ -730) was deleted on TNFα gene promoter, double luciferase activity was the highest, indicating that fragment (-730-0) was the transcription factor binding region. The site (-78 ~ -69) on the fragment (-730-0) binding NF-κB was mutated, and double luciferase activity decreased significantly. The results confirmed that the site (-78 ~ -69) was indeed an important binding site for NF-κB. In addition, the activity of TNFα in the serum of Japanese flounder changed with the prolongation of vibrio anguillarum infection, and the concentration of other immune factors such as ALP, ALT, AST and LDH also changed in the muscle under vibrio anguillarum infection. They all showed a trend of first increasing and then decreasing. Above studies implied that Japanese flounder responded to Vibrio anguillarum infection at the immune level with the change of its methylation status and the activation of transcription factor. By studying the mechanism of immune pathways, understanding the response to immune stress is great significant to the research of fish breeding for disease resistance.

Transcriptome Analysis in the Head Kidney of Rainbow Trout (Oncorhynchus mykiss) Immunized with a Combined Vaccine of Formalin-Inactivated Aeromonas salmonicida and Vibrio anguillarum

This study aimed to identify the molecular mechanisms regulated by a combined vaccine against Aeromonas salmonicida and Vibrio anguillarum (O1 serotype). These bacteria cause furunculosis and vibriosis, respectively, and are associated with a high mortality in rainbow trout in Korea. The vaccine upregulated gene expression of TCRα, T-bet, sIgM, and mIgM, markers of an activated adaptive immune response. On days 1, 3, and 5, transcriptome analysis revealed 862 (430 up- and 432 downregulated), 492 (204 up- and 288 downregulated), and 741 (270 up- and 471 downregulated) differentially expressed genes (DEGs), respectively. Gene ontology (GO) enrichment analysis identified 377 (108 MF, 132 CC, 137 BP), 302 (60 MF, 180 CC, 62 BP), and 314 (115 MF, 129 CC, 70 BP) GOs at days 1, 3, and 5, respectively. Kyoto Encyclopedia of Genetic and Genomic enrichment analysis identified eight immune system-related pathways like cytokine-cytokine receptor interaction, NF-kappaB signaling pathway, TNF signaling pathway, NOD-like receptor signaling pathway, cytosolic DNA sensing pathway, cell adhesion molecule, complement and coagulation cascade, and antigen processing and presentation. In the analysis of the protein–protein interaction of immune-related DEGs, a total of 59, 21, and 21 interactional relationships were identified at days 1, 3, and 5, respectively, with TNF having the highest centrality at all three time points.

Transcriptional Profiles of Genes Related to Stress and Immune Response in Rainbow Trout (Oncorhynchus mykiss) Symptomatically or Asymptomatically Infected With Vibrio anguillarum

Rainbow trout (Oncorhynchus mykiss) is one of the most common aquaculture fish species worldwide. Vibriosis disease outbreaks cause significant setbacks to aquaculture. The stress and immune responses are bidirectionally modulated in response to the health challenges. Therefore, an investigation into the regulatory mechanisms of the stress and immune responses in trout is invaluable for identifying potential vibriosis treatments. We investigated the transcriptional profiles of genes associated with stress and trout immune functions after Vibrio anguillarum infection. We compared the control trout (CT, 0.9% saline injection), asymptomatic trout (AT, surviving trout with minor or no symptoms after bacteria injection), and symptomatic trout (ST, moribund trout with severe symptoms after bacteria injection). Our results showed activated immunomodulatory genes in the cytokine network and downregulated glucocorticoid and mineralocorticoid receptors in both AT and ST, indicating activation of the proinflammatory cytokine cascade as a common response in AT and ST. Moreover, the AT specifically activated the complement- and TNF-associated immune defenses in response to V. anguillarum infection. However, the complement and coagulation cascades, as well as steroid hormone homeostasis in ST, were disturbed by V. anguillarum. Our studies provide new insights toward understanding regulatory mechanisms in stress and immune functions in response to diseases.

Molecular characterization of a novel p38 MAPK cDNA from Cyclina sinensis and its potential immune-related function under the threat of Vibrio anguillarum

The p38 mitogen-activated protein kinase (MAPK) is one important member of MAPK family and reported to serve a predominant function in regulating innate immunity after the occurrence of certain infection. In the present study, one novel p38 MAPK gene was acquired from Cyclina sinensis based on the RNA-seq analysis and designated as Csp38 MAPK. This novel gene contained a full length of 1781 bp, 1104 bp of which was deemed as open reading frames and gave rise to a peptide of 367 amino acids with a predicted molecular weight of 42.31 KDa. A conserved serine/threonine protein kinase (S_Tkc) region along with a Thr-Gly-Tyr motif was discovered in the deduced sequence. According to the phylogenetic analysis, there was a close relationship between this kinase and Meretrix petechialis p38 MAPK. As for the expression pattern, this newly-identified Csp38 MAPK was ubiquitously distributed in several tissues throughout the body but with varied abundance. After the challenge of Vibrio anguillarum, both the transcription and phosphorylation level of Csp38 MAPK in hemolymph were coordinately altered with a time-dependent manner. Besides, with the application of double strand RNA homologous to myeloid differentiation factor 88 (MyD88) of C. sinensis, the activation of Csp38 MAPK was found to obviously decrease in hemolymph after the pathogen stimulation. Hence, our experimental data presented evidence for the potential involvement of p38 MAPK in response to bacterial invaders in C. sinensis, possibly facilitating the understanding for pathogen-induced innate immunity in clams.