Transcriptome analysis of spleen reveals the signal transduction of toll-like receptors after Aeromonas hydrophila infection in Schizothorax prenanti

Role of the Pseudomonas plecoglossicida fliL gene in immune response of infected hybrid groupers (Epinephelus fuscoguttatus ♀ × Epinephelus lanceolatus ♂)

Pseudomonas plecoglossicida, a gram-negative bacterium, is the main pathogen of visceral white-point disease in marine fish, responsible for substantial economic losses in the aquaculture industry. The FliL protein, involved in torque production of the bacterial flagella motor, is essential for the pathogenicity of a variety of bacteria. In the current study, the fliL gene deletion strain (ΔfliL), fliL gene complement strain (C-ΔfliL), and wild-type strain (NZBD9) were compared to explore the influence of the fliL gene on P. plecoglossicida pathogenicity and its role in host immune response. Results showed that fliL gene deletion increased the survival rate (50%) and reduced white spot disease progression in the hybrid groupers. Moreover, compared to the NZBD9 strain, the ΔfliL strain was consistently associated with lower bacterial loads in the grouper spleen, head kidney, liver, and intestine, coupled with reduced tissue damage. Transcriptomic analysis identified 2 238 differentially expressed genes (DEGs) in the spleens of fish infected with the ΔfliL strain compared to the NZBD9 strain. Based on Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, the DEGs were significantly enriched in seven immune system-associated pathways and three signaling molecule and interaction pathways. Upon infection with the ΔfliL strain, the toll-like receptor (TLR) signaling pathway was activated in the hybrid groupers, leading to the activation of transcription factors (NF-κB and AP1) and cytokines. The expression levels of proinflammatory cytokine-related genes IL-1β, IL-12B, and IL-6 and chemokine-related genes CXCL9, CXCL10, and CCL4 were significantly up-regulated. In conclusion, the fliL gene markedly influenced the pathogenicity of P. plecoglossicida infection in the hybrid groupers. Notably, deletion of fliL gene in P. plecoglossicida induced a robust immune response in the groupers, promoting defense against and elimination of pathogens via an inflammatory response involving multiple cytokines.

Roles of eyestalk in salinity acclimatization of mud crab (Scylla paramamosain) by transcriptomic analysis

Salinity acclimatization refers to the physiological and behavioral adjustments made by crustaceans to adapt to varying salinity environments. The eyestalk, a neuroendocrine organ in crustaceans, plays a crucial role in salinity acclimatization. To elucidate the molecular mechanisms underlying eyestalk involvement in mud crab (Scylla paramamosain) acclimatization, we employed RNA-seq technology to analyze transcriptomic changes in the eyestalk under low (5 ppt) and standard (23 ppt) salinity conditions. This analysis revealed 5431 differentially expressed genes (DEGs), with 2372 upregulated and 3059 downregulated. Notably, these DEGs were enriched in crucial biological pathways like metabolism, osmoregulation, and signal transduction. To validate the RNA-seq data, we further analyzed 15 DEGs of interest using qRT-PCR. Our results suggest a multifaceted role for the eyestalk: maintaining energy homeostasis, regulating hormone synthesis and release, PKA activity, and downstream signaling, and ensuring proper ion and osmotic balance. Furthermore, our findings indicate that the crustacean hyperglycemic hormone (CHH) may function as a key regulator, modulating carbonic anhydrase expression through the activation of the PKA signaling pathway, thereby influencing cellular osmoregulation, and associated metabolic processes. Overall, our study provides valuable insights into unraveling the molecular mechanisms of mud crab acclimatization to low salinity environments.

Identifying low salinity adaptation gene expression in the anterior and posterior gills of the mud crab (Scylla paramamosain) by transcriptomic analysis

In the present study, BGISEQ-500 RNA-Seq technology was adopted to investigate how Scylla paramamosain adapts to salinity tolerance at the molecular level and explores changes in gene expression linked to salinity adaptation following exposure to both low salinity (5 ‰) and standard salinity (23 ‰) conditions. A total of 1100 and 520 differentially expressed genes (DEGs) were identified in the anterior and posterior gills, respectively, and their corresponding expression patterns were visualized in volcano plots and a heatmap. Further analysis highlighted significant enrichment of well-established gene functional categories and signaling pathways, including those what associated with cellular stress response, ion transport, energy metabolism, amino acid metabolism, H2O transport, and physiological stress compensation. We also selected key DEGs within the anterior and posterior gills that encode pivotal stress adaptation and tolerance modulators, including AQP, ABCA1, HSP 10, A35, CAg, NKA, VPA, CAc, and SPS. Interestingly, A35 in the gills might regulate osmolality by binding CHH in response to low salinity stress or serve as a mechanism for energy compensation. Taken together, our findings elucidated the intricate molecular mechanism employed by S. paramamosain for salinity adaptation, which involved distinct gene expression patterns in the anterior and posterior gills. These findings provide the foothold for subsequent investigations into salinity-responsive candidate genes and contribute to a deeper understanding of S. paramamosain's adaptation mechanisms in low-salinity surroundings, which is crucial for the development of low-salinity species cultivation and the establishment of a robust culture model.

Stress-associated β -glucan administration stimulates the TLR - MYD88 - NFKB1 signaling pathway in Nile tilapia (Oreochromis niloticus)

There is evidence that the administration of β-glucan can effectively activate several defense mechanisms, such as the Tlr-Myd88-Nfkb1 pathway that induces the expression of immune cytokines. Thus, the objective of this work was to evaluate whether β-glucan acts on the mechanisms of gene transcription via the Tlr-Myd88-Nfkb1 pathway in Nile tilapia under stress after challenge with Streptococcus agalactiae. Therefore, we evaluated the expression of immune system genes such as toll-like receptors 1 (tlr1), toll-like receptors 2 (tlr2), primary myeloid differentiation response gene (myd88) and nuclear factor kappa B1 (nfkb1). A total of 408 fish were distributed in 24 polyethylene boxes and randomly divided into eight groups with 3 replications each: C15: Tilapias received a control diet (free of β-glucan) for 15 days and were sampled after the 15th day of the experiment; C15D: Tilapias received a control diet (free of β-glucan) for 15 days, were challenged on the 14th day and were sampled at the 15th day of the experiment; β15: Tilapias received experimental diet (1g kg-1 of β-glucan) for 15 days and were sampled after 15 days; β15D: Tilapias received an experimental diet (1g kg-1 of β-glucan) for 15 days, were challenged on the 14th day and were sampled at the 15th day of the experiment; C30: Tilapias received a control diet (free of β-glucan) for 30 days and were sampled on the 30th day of the experiment; C30D: Tilapias received a control diet (free of β-glucan) for 30 days, were challenged on the 29th day and were sampled at the 30th day of the experiment; β30: Tilapias received experimental diet (1g kg-1 of β-glucan) for 30 days and were sampled after 30 days and β30D: Tilapias received experimental diet (1g kg-1 of β-glucan) for 30 days, were challenged on the 29th day and were sampled at 30 of the experiment. In the fish sampled at 15 and 30 days of the experiment, after being anesthetized and killed by brain section, cranial kidney and spleen were collected for gene expression analysis. The analyzes showed that the association of β-glucan and stressful management modulated the immune system, using the Tlr-Myd88-Nfkb1 signaling pathway, indicating that this compound can be used to promote early defense and protect fish against diseases.

Effects of TLR2/4 signalling pathway in western mosquitofish (Gambusia affinis) after Edwardsiella tarda infection

Gambusia affinis is regarded as an important animal model. Edwardsiella tarda is one of the most serious pathogens affecting aquaculture. The study explores the effects of TLR2/4 partial signalling pathway in G. affinis of E. tarda infection. The study collected the brain, liver, and intestine after E. tarda LD₅₀ and 0.85% NaCl solution challenge at different times (0 h, 3 h, 9 h, 18 h, 24 h, and 48 h). In these three tissues, the mRNA levels of PI3K, AKT3, IRAK4, TAK1, IKKβ, and IL-1β were substantially enhanced (p < .05) then returned to normal levels. Additionally, Rac1 and MyD88 in liver had different trend with other genes in brain and intestine, which displayed significantly indifference. The overexpression of IKKβ, and IL-1β indicated that E. tarda also caused immune reaction in intestine and liver, which would be consistent with delayed edwardsiellosis, which causes intestinal lesions and liver and kidney necrosis. Additionally, MyD88 plays a smaller role than IRAK4 and TAK1 in this signalling pathways. This study could enrich the understanding of the immune mechanism of the TLR2/4 signalling pathway in fish and might help to prescribe preventive measures against E. tarda to prevent infectious diseases in fish.

UNC93B1 facilitates TLR18-mediated NF-κB signal activation in Schizothorax prenanti

Toll-like receptor 18 (TLR18), a non-mammalian TLR, has been believed to play an important role in anti-bacterial immunity of teleost fishes. UNC93B1 is a classical molecular chaperone that mediates TLRs transport from endoplasmic reticulum to the located membrane. However, TLR18-mediated signal transduction mechanism and the regulatory effect of UNC93B1 to TLR18 are still unclear in teleost fishes. In this study, the coding sequences of TLR18 and UNC93B1 were cloned from Schizothorax prenanti, named spTLR18 and spUNC93B1, respectively. The spTLR18 and spUNC93B1 are 2583 bp and 1878 bp in length, encode 860 and 625 amino acids, respectively. The spTLR18 widely expressed in various tissues with the highest expression level in liver. After stimulation of Aeromonas hydrophila, lipopolysaccharide (LPS) and Poly(I:C), the expression levels of spTLR18 were significantly increased in spleen and head kidney. The spTLR18 located in the cell membrane, while spUNC93B1 located in the cytoplasm. Luciferase and overexpression analysis showed that spTLR18 activated NF-κB and type I IFN signal pathways, and spTLR18-mediated NF-κB activation might depend on the adaptor molecule MyD88. Besides, spUNC93B1 positively regulates spTLR18-mediated NF-κB signal. Our study first uncovers TLR18-UNC93B1-mediated signal transduction mechanism, which contributes to the understanding of TLR signaling pathway in teleost fishes.

Cloning of Toll-like Receptor 3 Gene from Schizothorax prenanti (SpTLR3), and Expressions of Seven SpTLRs and SpMyD88 after Lipopolysaccharide Induction

Toll-like receptor 3 (SpTLR3) from Schizothorax prenanti (S. prenanti) was cloned and identified, and the tissue distribution of the SpTLR3 gene was examined in this study. Moreover, the relative mRNA expression levels of myeloid differentiation factor 88 gene (SpMyD88) and seven TLR genes (SpTLR2, SpTLR3, SpTLR4, SpTLR18, SpTLR22-1, SpTLR22-2 and SpTLR22-3) from S. prenanti after lipopolysaccharide (LPS) challenge were analyzed through quantitative real-time polymerase chain reaction (qRT-PCR). The full length of SpTLR3 gene is 3097 bp, and complete coding sequence (CDS) is 2715 bp, which encodes 904 amino acids. The SpTLR3 amino acid sequence shared 43.94−100% identity with TLR3 sequences from other vertebrates; SpTLR3 was expressed in all eight tissues examined; and the highest level appeared in the liver, which was significantly higher than in all other tissues (p < 0.05), followed by the levels in the heart and muscles. LPS significantly up-regulated all eight genes in the S. prenanti tissues at 12 or 24 h (p < 0.05). Compared with the PBS control group, no significant transcripts changes were found in SpTLR2 or SpTLR3 at 12 h after LPS induction, but they were significantly up-regulated at 24 h (p < 0.001). The most abundant transcripts were found in the head kidney SpTLR22 genes after 24 h LPS induction, with high to low levels, which were SpTLR22-1 (564-fold), SpTLR22-3 (508-fold) and SpTLR22-2 (351-fold). Among these eight genes, the expression level of SpTLR4 was the least up-regulated. Overall, SpTLR4 in the head kidney was involved in the antibacterial immune response earlier, and the level was increased at 12 h with extreme significance after LPS stimulation (p < 0.001), while the other seven genes were the most significantly up-regulated at 24 h post injection. Taken together, the results suggest that SpMyD88, SpTLR2, SpTLR3, SpTLR4, SpTLR18, SpTLR22-1, SpTLR22-2 and SpTLR22-3 participate in an innate immune response stimulated by LPS, and the response intensity of the genes was organ-specific, with differing kinetics. Our findings will contribute to a more complete understanding of the roles of these TLR genes in antibacterial immunity.

Transcriptomic analysis reveals the gene expression profiles in the spleen of spotted knifejaw (Oplegnathus punctatus) infected by Vibrio harveyi

As one of the most valuable maricultured species, spotted knifejaw (Oplegnathus punctatus) has high popularity in eastern Asia. In recent years, diseases caused by Vibrio harveyi have brought huge economic losses in spotted knifejaw industry. To better understand the molecular mechanisms of immune response about V. harveyi resistance in spotted knifejaw, a comparative transcriptome analysis was performed on spleen tissues at five different time points post-infection (0, 12, 24, 48 and 72 hpi). A total of 4279 differentially expressed genes (DEGs) were identified. KEGG pathways analysis showed that multiple immune-related pathways were significant regulated, including Toll-like receptor signaling pathway, ECM-receptor interaction pathway, cytokine-cytokine receptor interaction pathway and hematopoietic cell lineage pathway. Weighted gene co-expression network analysis showed that several immune-related pathways of the highest correlation with 12 hpi (cor = 0.89, P = 7e-06) were significantly enriched. In addition, 12 hpi was a turning point for 7 gene clusters out of 9 that were divided according to gene expression patterns. Therefore, we speculated that 12 hpi might be a very critical time point for spotted knifejaw against V. harveyi infection. Additionally, qRT-PCR was carried out to validate the expressions of 12 DEGs. This study provided the first systematical transcriptome analysis of spotted knifejaw against V. harveyi. The results could help us better understand the dynamic immune responses of spotted knifejaw against bacterial infection, and provide useful information for antibacterial defense in spotted knifejaw industry as well.

Transcriptome analysis of turbot (Scophthalmus maximus) head kidney and liver reveals immune mechanism in response to Vibrio anguillarum infection

The diseases triggered by Vibrio anguillarum infection have created huge economic losses to the turbot (Scophthalmus maximus) farming industry. However, the immune mechanism of turbot to V. anguillarum infection has not been deeply investigated. To better understand the immune response of turbot to V. anguillarum infection, transcriptome analysis of the head kidney and liver of turbot was performed. A total of 15,948 and 11,494 differentially expressed genes (DEGs) were obtained from the turbot head kidney and liver, respectively. Transcriptome analysis revealed that the head kidney and liver of turbot have some differences in the gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis of the DEGs for the different functions of these two organs. Although there are many uncertain factors in this immune process, such as the occurrence of alternative splicing (AS) events and the differences in the protein structure of the DEGs, the NFκB signalling pathway, MKK-dependent AP-1 activation, JAK-STAT signalling pathway, the signal transmission of MHC Ⅰ and a series of DEGs including HSP90 driving NLRP3 to produce inflammatory factors (IL-1β, IL-8, TNFα, etc.) were possible important immune response pathways for turbot to V. anguillarum infection. Overall, our research has conducted a preliminary exploration of the immune mechanism of turbot in response to V. anguillarum infection.

Comprehensive analysis of full-length transcriptomes of Schizothorax prenanti by single-molecule long-read sequencing

Schizothorax prenanti (S. prenanti) is one of the most important aquaculture species in the southwest of China. However, information of the full-length transcripts in S. prenanti remains unknown. In this study, single-molecule real-time (SMRT) sequencing was performed to generate full-length transcriptomes of S.prenanti. In total, 23.26 Gb of clean reads were generated. A total of 312,587 circular consensus sequences (CCS) were obtained with average lengths of 2634 bp and 84.16% (270,662) of CCS were full-length non-chimeric reads. After being corrected with Illumina library sequencing, 18,005 contigs were obtained, with 17,797 (98.81%) successfully annotated in eight public databases, including 15,839 complete open reading frames (ORFs) with an average length of 1330 bp. Furthermore, a total of 4152 alternative splicing (AS) events and 250 long non-coding RNA (lncRNA) transcripts were detected. Additionally, a total of 1129 putative transcription factors (TFs) members from 56 TF families and 11,660 simple sequence repeats (SSRs) were identified. This study provided a valuable resource of full-length transcripts for further research on S. prenanti.

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.

Immune-related genes expression analysis of Western mosquitofish (Gambusia affinis) challenged with Aeromonas hydrophila

The great Gambusia affinis (G. affinis) is considered as an important animal model to study the endocrine disruption, ecological behavior, and environmental pollutant. The present study aims to build a new promising infection model with Aeromonas hydrophila (A. hydrophila) in aquaculture. The mRNA expression of Rac1, MyD88, IRAK4, TAK1, IKKβ, and IL-1β in G. affinis were significance higher (P < 0.05) in the liver of G. affinis than that of brain and intestine. And the PI3K mRNA expression level was significant lower (P < 0.05) in the intestine than that of brain and liver. The mRNA levels of AKT3 were significant higher (P < 0.05) in the brain than that of liver and intestine. And then the brain, liver, and intestine were collected at different time points (0 h, 3 h, 9 h, 18 h, 24 h, 48 h) after post injection of LD₅₀ of A. hydrophila. The 0.85% NaCl was used as a negative control for the LD₅₀ of A. hydrophila. The RT-PCR results showed that mRNA expressions of TLR2/4 pathway downstream genes MyD88, IRAK4, TAK1, Rac1, IKKβ, and IL-1β were firstly significantly up-regulated (P < 0.05) and were then backed to the 0 h group levels in three tissues. In contrast, mRNA expressions of TLR2/4 pathway downstream genes PI3K and AKT3 were firstly significantly decreased (P < 0.05) and were then increased to the 0 h group levels in brain and intestine. In summary, the results indicated that A. hydrophila could cause inflammatory reaction in intestinal and brain. In addition, the liver showed a provocative reaction when infected with A. hydrophila.

Transcriptome-based evaluation and validation of suitable housekeeping gene for quantification real-time PCR under specific experiment condition in teleost fishes

Quantification real-time PCR (qRT-PCR) is a common method in analysis of gene expression, but the stable reference genes for the normalization analysis have not been appreciated before identifying expression pattern of genes in teleost fishes. In this study, we selected eight candidate reference genes (18S, Actin, EF-1α, 40S, B2M, TUBA, UBCE and GAPDH) basing on transcriptome analysis and the traditional housekeeping genes, and analyzed the stability of the reference genes in spleen, head kidney and head kidney leukocytes (HKL) after pathogen challenge in Schizothorax prenanti (S. prenanti). Three common programs (geNorm, NormFinder and Bestkeeper) were used to evaluate the stability of the candidate reference genes. Two reference genes, Actin and EF-1α presented higher stability, while 18S and GAPDH were the lower stable genes, both in in vitro and in vivo. An important immune gene, toll-like receptor 22a (TLR22a), was selected to validate the stability of the proposed reference genes (Actin and EF-1α) across different experiment treatments. The results reveal that Actin and EF-1α are quite suitable reference genes for the normalization analysis. Otherwise, using the most stable gene Actin to validate the reliable of transcriptome data showed the high correlation between the fold change of transcriptome data and qRT-PCR data. In conclusion, our study not only acquired the suitable reference gene for the qRT-PCR assay under specific experiment condition, but also provided a comprehensive method to evaluate and validate the reference gene based on transcriptome analysis in teleost fishes.

The potential sensing molecules and signal cascades for protecting teleost fishes against lipopolysaccharide

Lipopolysaccharide (LPS) is a classical pathogen-associated molecular pattern that can trigger strong inflammatory response mainly by TLR4-mediated signaling pathway in mammals, but the molecular mechanism of anti-LPS immunity is unclear in teleost fishes. In this study, we analyzed the gene expression features based on transcriptome analysis in Schizothorax prenanti (S. prenanti), after stimulation with two sources of LPS from Aeromonas hydrophila and Escherichia coli (Ah. LPS and Ecoli. LPS). 921 different expression genes (DEGs) after Ah. LPS stimulation and 975 DEGs after Ecoli.LPS stimulation were acquired, but only 706 and 750 DEGs were successfully annotated into the databases, respectively. Both of two groups of DGEs were significantly enriched into immune-related pathways by KEGG enrichment analysis, such as "Toll-like receptor signaling pathway", "Cytokine-cytokine receptor interaction" and "JAK-STAT signaling pathway". The annotated DEGs from Ah. LPS and Ecoli. LPS stimulation shared 470 DEGs, including 88 immune-related DEGs (IRGs) identified mainly by KEGG enrichment to immune-related signaling pathways. Among the shared IRGs, four pattern-recognition genes (TLR5, TLR25, PTX3 and C1q) were induced with high expression foldchange, and IFN-γ and relative genes also showed higher expression levels than control. Meanwhile, inflammatory signals were highlighted by upregulating the expression of inflammatory cytokines (IL-1β, IL-10 and IL-8). Moreover, some non-shared IRGs (including TLR2 and TLR4) were identified, suggesting that different sources of LPS own different potentials for the induction of immune gene expression. In conclusion, TLR5, TLR25, PTX3 and C1q may function as the sensing molecules to catch the invasion signal of LPS. The anti-LPS immune response may be involved into TLR25/TLR5-mediated inflammatory signals that regulate subsequently the activation of PTX3/C1q-modulated complement pathway upon the induction of PTX3 expression, and the crosstalk between IFN-γ and TLR signaling pathways in teleost fishes. This study will contribute to further explore the molecular mechanism of LPS-induced immunity in teleost fishes.

Clone, identification and functional character of two toll-like receptor 5 molecules in Schizothorax prenanti

Mammal Toll-like receptor 5 (TLR5) can directly recognize bacterial flagellin, initiate the inflammatory signaling cascades and trigger body immune system to clear the "non-self" substances. In teleosts, TLR5 has presented more complexes not only in increasing the molecular types, but also in elevating the functional diversity. In this study, we identified two TLR5 family members in Schizothorax prenanti, named as spTLR5-1 and spTLR5-2. The complete coding sequence (CDS) of spTLR5-1 is 2622 bp, encoding 873 amino acids, while the complete CDS of spTLR5-2 is 2640 bp, encoding 879 amino acids. Phylogenetic analysis showed that spTLR5-1 and spTLR5-2 were clustered to the TLR5 of schizothorax richardsonii and Cyprinus carpio respectively. The 3D structure analysis exhibited that the α-helix, β-sheet, and the ligand binding site of spTLR5-1, spTLR5-2 and human TLR5 have large differences. The spTLR5-1 and spTLR5-2 had extensively expressed in various tissues, including the higher expression in liver, spleen and head kidney. Both the expression levels of spTLR5-1 and spTLR5-2 were significantly up-regulated after Aeromonas hydrophila (A. hydrophila) challenge. And, the downstream genes, such as AP-1, IKK-α, NF-kB, IL-1β, IL-8 and TNF-α, were also significantly up-regulated after A. hydrophila challenge. Apart from that, the luciferase reporter assay demonstrated that the co-transfection of spTLR5-1 or spTLR5-2 into HEK293T cells showed the significantly increased NF-kB luciferase activity after flagellin stimulation. In conclusion, our results reveal that both two molecular types of fish TLR5 may commonly mediate the recognition of flagellin and the activation of the downstream inflammatory signaling molecules.

Comparative transcriptomics and histopathological analysis of crucian carp infection by atypical Aeromonas salmonicida

Aeromonas salmonicida is a ubiquitous fish pathogen known to cause furunculosis. With the emergence of new subtypes and the expansion of the host range, it has threatened the health of a variety of marine and freshwater fish, particularly the non-salmonids, manifesting differently from the classical furunculosis. Although there have been reports of infection by atypical strains on the crucian carp, the pathogenesis and tissue pathology remain unclear. In this study, transcriptomics and histopathology were used to analyze the immune response and lesions of crucian carp infected with A. salmonicida. Comparative analysis showed 6579 differentially expressed genes (DEGs) (3428 down-regulated and 3151 up-regulated) were identified on day 5 post-infection (5 dpi). Further annotation and analysis revealed that the DEGs were enriched in enzyme regulator activity, response to oxidative stress, iron ion homeostasis and other functions, and mitogen-activated protein kinase (MAPK), nuclear factor-κB (NF-κB), toll-like receptor (TLR), and nucleotide-binding oligomerization domain (NOD)-like receptor (NLR) etc., and immune-related signaling pathways. Meanwhile, the four C-type lysozyme genes found in all DEGs were significantly up-regulated after infection. In addition, there was severe bleeding on the body of the infected fish. Also, the intestine, liver, spleen, and kidney showed varying degrees of inflammatory damage, especially the goblet cell hyperplasia of intestinal mucosa epithelium and degeneration and necrosis of renal tubular epithelium cells. Additionally, with the increase in pathogen concentration, the cumulative mortality increased, the severity of lesions in the hindgut and head-kidney tissues increased. The relative expression levels of four immune-related genes (TNF-α, IL-1β, IL-11, C-lysozyme) were also significantly upregulated, compared with the control (P < 0.05). In conclusion, this study provides a scientific basis for further study on the immune response, pathological diagnosis, and prevention of crucian carp infection caused by atypical A. salmonicida.

Multiple subtypes of TLR22 molecule from Schizothorax prenanti present the functional diversity in ligand recognition and signal activation

Evolutionary development has increased the diversity of genotypes and the complexity of gene functions in fish. TLR22 has been identified as a teleost-specific gene, but its functions are tremendously different among different fish species. Whether the functional diversity relates to the difference of genotypes remains poorly understand. In this study, we cloned and identified three TLR22 molecules from Schizothorax prenanti (S. prenanti), named as spTLR22-1, spTLR22-2 and spTLR22-3. The full-length coding regions of spTLR22s are 2841 bp, 2805 bp and 2868 bp and coding 946 aa, 934 aa and 955 aa, respectively. All spTLR22s are composed of multiple leucine-rich repeat (LRR) domains, a transmembrane structure and a Toll/IL-1 receptor (TIR) region. The phylogenetic analysis showed that three spTLR22s were close to Cyprinus carpio TLR22-1, TLR22-2 and TLR22-3, respectively. Among the spTLR22s, they presented not close relationship but remained to belong to TLR22 subfamily. All spTLR22s were ubiquitously expressed in all tested tissues, but the expression levels of spTLR22s were dominant in immune-related tissues, such as gill and spleen. The expression levels of spTLR22-1 and spTLR22-3 were significantly increased after treatment with bacteria, LPS and Poly(I:C). However, spTLR22-2 seems like no response to these treatments. The luciferase reporter assay demonstrated that all spTLR22s could activate NF-κB signaling pathway, but only spTLR22-1 and spTLR22-2 could activate IFN-β signaling pathway. Interestingly, in the ligand recognition analysis, spTLR22-1 and spTLR22-3 but not spTLR22-2 had the recognized potential to Poly(I:C), and all spTLR22s could not recognize LPS. Both spTLR22-1 and spTLR22-3 significantly up-regulated the expression of anti-viral-related genes (Mx, IFN and ISG15) and down-regulated the expression of anti-inflammatory factor IL-10 after the overexpression in carp EPC cell line, but spTLR22-2 failed to impact the expression of these genes. Moreover, we found that all spTLR22s localized to the intracellular region. Taken together, our results reveal that spTLR22-1 and spTLR22-3 but not spTLR22-2 may be involved into the anti-viral immune response via IFN-β signaling pathway, and all spTLR22s can activate NF-κB signaling pathway but only spTLR22-1 and spTLR22-3 response to the stimulation of bacteria and LPS.

The key differentially expressed genes and proteins related to immune response in the spleen of pufferfish (Takifugu obscurus) infected by Aeromonas hydrophila

The immune mechanism elicited in pufferfish (Takifugu obscurus) against the invasion of Aeromonas hydrophila is still poorly understood. We examined the spleen of pufferfish at the transcriptome and proteome levels by using Illumina-seq and TMT coupled mass spectrometry after 12 h infection by A. hydrophila, respectively. A total of 2,339 genes (1,512 up-regulated and 827 down-regulated) and 537 (237 up-regulated and 300 down-regulated) proteins were identified. GO and KEGG analyses revealed that the responses to stimulus were the main biological processes, intestinal immune network for IgT production and calcium signaling pathway. Fourteen genes (8 up-regulated and 6 down-regulated) and proteins (5 up-regulated and 9 down-regulated) involved immune responses or signal transduction were validated by qRT-PCR and parallel reaction monitoring to confirm the reliability of the transcriptomic and proteomic analyses, respectively. Moreover, qRT-PCR and flow cytometry were used to detect dynamics of the genes in calcium signaling pathway and changes of concentration of cytoplasm Ca²⁺ in spleen cells within a 72 h challenge. This study provides the findings regarding immune response, especially intestinal immune network for IgT production pathway and calcium signaling pathway at the molecular, protein and cellular in pufferfish after infection by A. hydrophila. These results would provide a new insight and molecular targets into the response to pathogenic infection in pufferfish.