Transcriptomic analysis of the head kidney of Topmouth culter (Culter alburnus) infected with Flavobacterium columnare with an emphasis on phagosome pathway

Effects of dietary melatonin on antioxidant and immune function of the Pacific white shrimp (Litopenaeus vannamei), as determined by transcriptomic analysis

Melatonin (MT) is regarded as an antioxidant and immunostimulant that can efficiently scavenge free radicals and activate antioxidant enzymes. The aim of this study was to investigate the effects of dietary MT on the growth performance and immune function of the Pacific white shrimp (Litopenaeus vannamei). Six groups of L. vannamei were supplemented with dietary MT at 0, 22.5, 41.2, 82.7, 165.1, and 329.2 mg/kg levels for 2 months. RNA-Seq analysis was performed to obtain transcriptome data of the control group and the group supplemented with dietary MT at 82.7 mg/kg BW. In total, 1220 DEGs (799 up-regulated and 421 down-regulated) were identified. Pathways and genes related to growth performance and immune function were verified by real-time quantitative polymerase chain reaction. The total hemocyte count, phagocytosis rate, and respiratory burst were significantly increased in the MT (82.7 mg/kg BW) group as compared to the control group. Analysis of antioxidant-related enzymes in the hepatopancreas showed that dietary MT (82.7 mg/kg BW) significantly increased activities of superoxide dismutase (SOD), glutathione peroxidase (GPx), and catalase, while dietary MT at 41.2 mg/kg BW significantly increased activities of glutathione S-transferase, lysozyme (LZM), and phenoloxidase (PO). At the transcriptional level, dietary MT up-regulated expression levels of genes associated with antioxidant immunity and growth, which included PO, SOD, LZM, GPx, chitin synthase, ecdysone receptor, calcium-calmodulin dependent protein kinase I, and retinoid X receptor. In conclusion, dietary MT may improve the growth performance and immune function of L. vannamei to some extent.

Comparative transcriptomes reveal different tolerance mechanisms to Streptococcus agalactiae in hybrid tilapia, nile tilapia, and blue tilapia

Tilapia is one of the most economically important freshwater fish farmed in China. Streptococcosis outbreaks have been extensively documented in farmed tilapia species. Hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂) exhibit greater disease resistance than Nile tilapia (O. niloticus) and blue tilapia (O. aureus). However, the molecular mechanism underlying the enhanced tolerance of hybrid tilapia is still poorly understood. In this study, comparative transcriptome analysis was performed to reveal the different tolerance mechanisms to Streptococcus agalactiae in the three tilapia lines. In total, 1982, 2355, and 2076 differentially expressed genes were identified at 48 h post-infection in hybrid tilapia, Nile tilapia, and blue tilapia, respectively. Functional enrichment analysis indicated that numerous metabolic and immune-related pathways were activated in all three tilapia lines. The differential expression of specific genes associated with phagosome, focal adhesion, cytokine-cytokine receptor interaction, and toll-like receptor signaling pathways contributed to the resistance of hybrid tilapia. Notably, immune response genes in hybrid tilapia, such as P38, TLR5, CXCR3, CXCL12, PSTPIP1, and TFR, were generally suppressed under normal conditions but selectively induced following pathogen challenge. These results expand our knowledge of the molecular mechanisms underlying S. agalactiae tolerance in hybrid tilapia and provide valuable insights for tilapia breeding programs.

17β-Trenbolone activates androgen receptor, upregulates transforming growth factor beta/bone morphogenetic protein and Wnt signaling pathways, and induces masculinization of caudal and anal fins in female guppies (Poecilia reticulata)

Sexually mature female guppies (Poecilia reticulata) were exposed to environmentally relevant concentrations (20, 200, and 2000 ng/L) of 17β-trenbolone for four weeks. As evidenced by the increased caudal fin index and anal fins developing into gonopodium-like structures, exposed females displayed masculinized secondary sexual characteristics. Differential gene expression and subsequent pathway analysis of mRNA sequencing data revealed that the transcription of transforming growth factor beta/bone morphogenetic protein signaling pathway and Wnt signaling pathway were upregulated following 17β-trenbolone exposure. Enzyme-linked immunosorbent assays showed that the bone morphogenetic protein 7 protein content was elevated after 17β-trenbolone exposure. Finally, real-time PCR revealed that 17β-trenbolone treatment significantly increased androgen receptor mRNA levels, and molecular docking showed potent interaction between 17β-trenbolone and guppy androgen receptor. Furthermore, 17β-trenbolone-induced masculinization of caudal and anal fins in female guppies, concomitant to the upregulated expression of differentially expressed genes involved in the above-mentioned two signaling pathways, was significantly inhibited by flutamide (androgen receptor antagonist). These findings demonstrated that 17β-trenbolone masculinized fins of female guppies by activating the androgen receptor. This study revealed that 17β-trenbolone could upregulate signaling pathways related to fin growth and differentiation, and eventually cause caudal and anal fin masculinization in female guppies.

The lysosome-phagosome pathway mediates immune regulatory mechanisms in Mesocentrotus nudus against Vibrio coralliilyticus infection

Sea urchins are a popular model species for studying invertebrate diseases. The immune regulatory mechanisms of the sea urchin Mesocentrotus nudus during pathogenic infection are currently unknown. This study aimed to reveal the potential molecular mechanisms of M. nudus during resistance to Vibrio coralliilyticus infection by integrative transcriptomic and proteomic analyses. Here, we identified a total of 135,868 unigenes and 4,351 proteins in the four infection periods of 0 h, 20 h, 60 h and 100 h in M. nudus. In the I20, I60 and I100 infection comparison groups, 10,861, 15,201 and 8,809 differentially expressed genes (DEGs) and 2,188, 2,386 and 2,516 differentially expressed proteins (DEPs) were identified, respectively. We performed an integrated comparative analysis of the transcriptome and proteome throughout the infection phase and found very a low correlation between transcriptome and proteome changes. KEGG pathway analysis revealed that most upregulated DEGs and DEPs were involved in immune strategies. Notably, "lysosome" and "phagosome" activated throughout the infection process, could be considered the two most important enrichment pathways at the mRNA and protein levels. The significant increase in phagocytosis of infected M. nudus coelomocytes further demonstrated that the lysosome-phagosome pathway played an important immunological role in M. nudus resistance to pathogenic infection. Key gene expression profiles and protein‒protein interaction analysis revealed that cathepsin family and V-ATPase family genes might be key bridges in the lysosome-phagosome pathway. In addition, the expression patterns of key immune genes were verified using qRT‒PCR, and the different expression trends of candidate genes reflected, to some extent, the regulatory mechanism of immune homeostasis mediated by the lysosome-phagosome pathway in M. nudus against pathogenic infection. This work will provide new insights into the immune regulatory mechanisms of sea urchins under pathogenic stress and help identify key potential genes/proteins for sea urchin immune responses.

Integrated Analysis of lncRNA and circRNA Mediated ceRNA Regulatory Networks in Skin Reveals Innate Immunity Differences Between Wild-Type and Yellow Mutant Rainbow Trout (Oncorhynchus mykiss)

Fish skin is a vital immune organ that forms the first protective barrier preventing entry of external pathogens. Rainbow trout is an important aquaculture fish species that is farmed worldwide. However, our knowledge of innate immunity differences between wild-type (WR_S) and yellow mutant rainbow trout (YR_S) remains limited. In this study, we performed whole transcriptome analysis of skin from WR_S and YR_S cultured in a natural flowing water pond. A total of 2448 mRNAs, 1630 lncRNAs, 22 circRNAs and 50 miRNAs were found to be differentially expressed (DE). Among these DEmRNAs, numerous key immune-related genes, including ifih1, dhx58, trim25, atp6v1e1, tap1, tap2, cd209, hsp90a.1, nlrp3, nlrc3, and several other genes associated with metabolism (gstp1, nampt, naprt and cd38) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEmRNAs revealed that many were significantly enriched in innate immune-related GO terms and pathways, including NAD+ADP-ribosyltransferase activity, complement binding, immune response and response to bacterium GO terms, and RIG-I-like receptor signaling, NOD-like receptor signaling and phagosome KEGG pathways. Furthermore, the immune-related competing endogenous RNA networks were constructed, from which we found that lncRNAs MSTRG.11484.2, MSTRG.32014.1 and MSTRG.29012.1 regulated at least three immune-related genes (ifih1, dhx58 and irf3) through PC-5p-43254_34, PC-3p-28352_70 and bta-miR-11987_L-1R-1_1ss8TA, and tap2 was regulated by two circRNAs (circRNA5279 and circRNA5277) by oni-mir-124a-2-p5_1ss13GA. The findings expand our understanding of the innate immune system of rainbow trout, and lay the foundation for further study of immune mechanisms and disease resistance breeding.

The Isolation and Full-Length Transcriptome Sequencing of a Novel Nidovirus and Response of Its Infection in Japanese Flounder (Paralichthys olivaceus)

A novel nidovirus, CSBV Bces-Po19, was isolated from the marine fish, Japanese flounder (Paralichthys olivaceus). The viral genome was 26,597 nucleotides long and shared 98.62% nucleotide identity with CSBV WHQSR4345. PacBio Sequel and Illumina sequencing were used to perform full-length transcriptome sequencing on CSBV Bces-Po19-sensitive (S) and -resistant (R) Japanese flounder. The results of negative staining revealed bacilliform and spherical virions. There were in total 1444 different genes between CSBV Bces-Po19 S and R groups, with 935 being up-regulated and 513 being down-regulated. Metabolism-, immune-, and RNA-related pathways were significantly enriched. Furthermore, CSBV Bces-Po19 infection induced alternative splicing (AS) events in Japanese flounder; the S group had a higher numbers of AS events (12,352) than the R group (11,452). The number of long non-coding RNA (lncRNA) in the S group, on the other hand, was significantly lower than in the R group. In addition to providing valuable information that sheds more light on CSBV Bces-Po19 infection, these research findings provide further clues for CSBV Bces-Po19 prevention and treatment.

Identification of mRNA-miRNA-lncRNA regulatory network associated with the immune response to Aeromonas salmonicides infection in the black rockfish (Sebastes schlegelii)

In aquaculture, Aeromonas salmonicides (A. salmonicida) is a main fish pathogen because of its nearly worldwide distribution, and broad host range. Recently, an increasing number of evidences have uncovered the roles of mRNA-miRNA-lncRNA network in fish diseases. In current study, RNA-seq was conducted in the black rockfish spleen following A. salmonicida infection at 0 h (Sp0 or control) and three different post-infection time-points (2 h: Sp2, 12 h: Sp12 and 24 h: Sp24, respectively) to comprehensively identify differentially expressed (DE) mRNAs, miRNAs and lncRNAs. Enrichment analysis and protein-protein interaction (PPI) analysis of DE mRNAs were performed. Then, expression and correlation analysis for mRNAs and their upstream miRNAs and lncRNAs were conducted. Finally, a total of 1364 mRNAs, 17 miRNAs and 1584 lncRNAs exhibited significantly differential expressions during bacterial infection in the black rockfish spleen. Functional enrichment analysis suggested that they were significantly enriched in several immune-related pathways, including Amino sugar and nucleotide sugar metabolism, Cell adhesion molecules (CAMs), Neuroactive ligand-receptor interaction, Nicotinate and nicotinamide metabolism, Pentose and glucuronate interconversions, Phagosome, Proteasome, etc. Subsequently, 1091 lncRNA-miRNA-mRNA pathways (323 in Sp2, 609 in Sp12 and 207 in Sp24) were constructed including 400 lncRNAs, 69 miRNAs, and 70 mRNAs. Meanwhile, NLRC3/novel-264/LNC_00116154 pathway demonstrated important immune modulating function in the black rockfish against A. salmonicida infection. Finally, the novel mRNA-miRNA-lncRNA sub-networks were established, among which all mRNAs and ncRNAs possessed significant predictive values for further studies for immune responses in the black rockfish.

Transcriptomic analysis of gills in nitrite-tolerant and -sensitive families of Litopenaeus vannamei

Nitrite stress is a major environmental factor that limits aquatic animal growth, reproduction and survival. Even so, some shrimps still can withstand somewhat high concentrations of nitrite environment. However, few studies have been conducted about the tolerance molecular mechanism of Litopenaeus vannamei in the high concentration nitrite. To identify the genes and pathways involved in the regulation of nitrite tolerance, we performed comparative transcriptomic analysis in the L. vannamei nitrite-tolerant (NT) and nitrite-sensitive (NS) families, and untreated shrimps were used as the control group. After 24 h of nitrite exposure (NaNO₂, 112.5 mg/L), a total of 1521 and 868 differentially expressed genes (DEGs) were obtained from NT compared with NS and control group, respectively. Functional enrichment analysis revealed that most of these DEGs were involved in immune defense, energy metabolism processes and endoplasmic reticulum (ER) stress. During nitrite stress, energy metabolism in NT was significantly enhanced by activating the related genes expression of oxidative phosphorylation (OXPHOS) pathway and tricarboxylic acid (TCA) cycle. Meanwhile, some DEGs involved in innate immunity- related genes and pathways, and ER stress responses also were highly expressed in NT. Therefore, we speculate that accelerated energy metabolism, higher expression of immunity and ER related genes might be the important adaptive strategies for NT in relative to NS under nitrite stress. These results will provide new insights on the potential tolerant molecular mechanisms and the breeding of new varieties of nitrite tolerant L. vannamei.

iTRAQ-based quantitative proteomic analysis of Pelteobagrus vachelli liver in response to hypoxia

Dissolved oxygen is one of the determinants in the healthy farming of Pelteobagrus vachelli. This study, we conducted quantitative proteomics on the juvenile P. vachelli livers using iTRAQ. P. vachelli were treated by 3.75 ± 0.25 mg O₂/L (hypoxia group) and 7.25 ± 0.25 mg O₂/L (control group) for 90 days. The results revealed that under hypoxic conditions, P. vachelli grew slower than control group. Proteomic profiling enabled us to identify 2618 proteins, of which 176 were significantly differentially abundant proteins (DAPs). Verification of protein regulation based on qRT-PCR indicated that the proteomics data were reliable. The top 20 significantly DAPs (10 up-regulated, 10 down-regulated) were primarily involved in energy metabolism, apoptosis inhibition, and heavy metal detoxification. KEGG pathway enrichment analysis revealed significant enrichment of 'protein digestion and absorption', 'glycolysis/gluconeogenesis', and 'phagosome'. Combining the proteomics results of short-term hypoxia (treated with 0.70 ± 0.10 mg O₂ /L for 4 h), we screened 36 common DAPs. The analysis of the 36 common DAPs indicated that P. vachelli responded to the hypoxia by regulating energy supply, inhibiting apoptosis, and disturbing defensive system. Our results lay a theoretical foundation for the cultivation of hypoxia-tolerant species and eco-breeding of P. vachelli. SIGNIFICANCE OF THE STUDY: The hypoxia tolerance of Pelteobagrus vachelli is poor, which will seriously lead to its death in high-density culture. This study analysed the liver proteome of P. vachelli under long-term hypoxia stress (treated for 90 days at 3.75 ± 0.25 mg O₂/L), and then combined the proteome results of short-term hypoxia stress (treated for 4 h at 0.70 ± 0.10 mg O₂/L). The results showed P. vachelli responded to the hypoxia by regulating energy supply, inhibiting apoptosis and disturbing defensive system. The study contributes to the breeding of new hypoxia-tolerant species of P. vachelli and lays the theoretical foundation for eco-breeding.

The Dynamic Immune Response of Yellow Catfish (Pelteobagrus fulvidraco) Infected With Edwardsiella ictaluri Presenting the Inflammation Process

Edwardsiella ictaluri is a highly destructive pathogen in cultured yellow catfish, thus it was very necessary to study the immune response of yellow catfish against bacterial infection. In this study, RNA-Seq technology was used to study the immune response in two distinct tissues of yellow catfish at eight different time points (h) after E. ictaluri infection. The number of differentially expressed genes (DEGs) in the spleen and liver was low at 3 h and 6 h post-infection, respectively. Afterwards, the most number of DEGs in the spleen was detected at 72 h, while the number of DEGs in the liver maintained a high level from 24 h to 120 h. The GO and KEGG enrichment analyses of DEGs at different time points uncovered that cytokines were continuously transcribed at 6 h to 120 h; whereas the liver is the main organ that secretes the components of the complement system, and metabolic regulation was activated from 12 h to 120 h. Moreover, an overview of the inflammation response of yellow catfish was exhibited including pattern-recognition receptors, inflammatory cytokines, chemokines, complements, and inflammation-related signal pathways. The similar expression tendency of nine genes by qRT-PCR validated the accuracy of transcriptome analyses. The different transcriptomic profiles obtained from the spleen and liver will help to better understand the dynamic immune response of fish against bacterial infection, and will provide basic information for establishing effective measures to prevent and control diseases in fish.

Identification of reproduction-related genes and pathways in the Culter alburnus H-P-G axis and characterization of their expression differences in malformed and normal gynogenetic ovaries

This study applied RNA-seq technology to discover reproduction-related genes and pathways in female topmouth culter brain (including pituitary) and ovarian tissues. In functional analysis, 2479 and 2605 unigenes in the brain and ovary tissue were assigned to the "reproductive process" subcategory in addition to the 2660 and 2845 unigenes assigned to the "reproduction" subcategory. Twenty-three complete cDNA sequences were identified through the different gene expression (DGE) approach from five reproduction-related pathways (MAPK signaling pathway, neuroactive ligand-receptor interaction pathway, gonadotropin-releasing hormone signaling pathway, oocyte meiosis pathway, and steroid biosynthesis pathway). The expression levels of 16 candidate genes using qPCR in this study were in accordance with the results of transcriptome analysis. In addition, the expression levels of the FSH, 3β-HSD, PGR, and NPYR genes in malformed gynogenetic ovaries were considerably low, which was consistent with the progress of oocytogenesis in the ovaries of topmouth culter. The high expression of these four genes in the ovaries of normal topmouth culter suggested they might involve in the preparation for the shift of oogenesis to ovulation. Hence, our work identified a set of annotated gene products that are candidate factors affecting reproduction in the topmouth culter H-P-G axis. These results could be essential for further research in functional genomics and genetic editing for topmouth culter reproduction.

Transcriptome analysis and immune-related genes expression reveals the immune responses of Macrobrachium rosenbergii infected by Enterobacter cloacae

Macrobrachium rosenbergii is an important cultural species in China and other Southeast Asian countries. However, Enterobacter cloacae infection has caused a great economic loss in M. rosenbergii culture industry. The immune responses of M. rosenbergii to the E. cloacae infection is not fully characterized. To investigate the immune response of M. rosenbergii against E. cloacae, we performed transcriptome analysis of the M. rosenbergii hepatopancreas with and without E. cloacae infection using RNA-seq. After assembly and annotation, 29,731 high quality unigenes were obtained from RNA-seq data. Differential expression analysis revealed the existence of 2498 significantly differently expressed genes (DEGs) at 12 h post infection, with 1365 up-regulated and 1133 down-regulated genes. Among these DEGs, some well-known immune-related genes were up-regulated significantly, including C-type lectin 1, lectin 3, anti-lipopolysaccharide factor 2, Cu/Zn superoxide dismutase and heat shock protein 70. GO analysis demonstrated 24 biological process subcategories, 14 cellular component subcategories, and 12 molecular function subcategories that were enriched among these DEGs, and some DEGs were clustered into immune related subcategories such as immune system process, response to stimulus, biological adhesion, and antioxidant activity. These DEGs were enriched into 216 KEGG pathways including a core set of immune correlated pathways notably in phagosome and lysosome. In addition, 5 up-regulated and 5 down-regulated immune-related DEGs were selected for further validation by quantitative real-time PCR and the results showed consistence with the RNA-seq data. Additionally, the expression level of six selected immune-related genes (ALF2, CLEC1, LEC3, hemocyanin1, HSP70 and SOD) based on the transcriptomic data were monitored at different point of time in hepatopancreas, gill, hemolymph and intestine. Results revealed these immune-related genes were significantly up-regulated in different tissues from 6 to 24 h after E. cloacae infection. Overall, these results provided valuable information for further studying the immune response of M. rosenbergii against E. cloacae infection.

Analysis of the treatment regimen efficacy for columnaris disease in Pterophyllum scalare

The article presents the results of studies on the treatment scheme efficacy for columnaris in Pterophyllum scalare, common under private aquarium husbandry conditions. To establish the diagnosis, the clinical features of the diseased fish, pathological and anatomical changes and the results of microscopic and microbiological studies were taken into account. Separate chemical and microbiological parameters of aquarium water were also studied. It was established that fish disease developed against the background of adverse changes in the chemical composition and microbiocenosis of aquarium water. High alkalinity and excess of phosphates compared to the norm provoked accumulation of opportunistic microbiota, resulting in a balance disorder in the parasite-host system and development of clinical manifestation of the fish disease. During the disease outbreak, bacteriological indices of water indicated a high level of organic contamination and a low intensity of water self-purification processes. Clinically, the disease was manifested in P. scalare by decrease in appetite and motor activity, onset of ulcerative lesions of various shapes and sizes on the surface of the body and on the gill covers. Selected pure cultures of Flavobacterium columnare showed sensitivity to enrofloxacin (growth retardation zone 31.3 ± 1.0 mm); moderate resistance was found to tylosin. The microorganisms were resistant to amoxicillin, doxycycline, benzylpenicillin and tetracycline. Microscopic studies of intestinal specimens of dead P. scalare revealed numerous motile flagellates. It has been shown that an effective treatment regimen that provides recovery for 70% of diseased P. scalare is the use of enroxil 10% solution for five days, metronidazole three times a day, and “API MelaFix” for seven days. It is proved that the following measures are effective to restore the disrupted hydro-balance: periodic water replacement in the amount of 20% of the total volume, providing the aquarium with active aeration systems, planting slow-growing plants and reducing the amount of fish food provided. The measures developed were efficient, they led to elimination of the outbreak of columnaris in the P. scalare and to restoration of biological equilibrium in a closed aquatic ecosystem.

Integrated analysis of mRNA-miRNA expression in Tilapia infected with Tilapia lake virus (TiLV) and identifies primarily immuneresponse genes

We investigated differential gene expression in Tilapia infected with the Tilapia Lake virus (TiLV).We used high-throughput sequencing to identify mRNAs and miRNAs involved in TiLV infection progression We identified 25,359 differentially expressed genes that included 863 new genes. We identified 1770, 4142 and 4947 differently expressed genes comparing non-infected controls with 24 and 120 h infections and between the infected groups, respectively. These genes were enriched to 291 GO terms and 62 KEGG pathways and included immune system progress and virion genes. High-throughput miRNA sequencing identified 316 conserved miRNAs, 525 known miRNAs and 592 novel miRNAs. Furthermore, 138, 198 and 153 differently expressed miRNAs were found between the 3 groups listed above, respectively. Target prediction revealed numerous genes including erythropoietin isoform X2, double-stranded RNA-specific adenosine deaminase isoform X1, bone morphogenetic protein 4 and tapasin-related protein that are involved in immune responsiveness. Moreover, these target genes overlapped with differentially expressed mRNAs obtained from RNA-seq. These target genes were significantly enriched to GO terms and KEGG pathways including immune system progress, virion and Wnt signaling pathways. Expression patterns of differentially expressed mRNA and miRNAs were validated in 20 mRNA and 19 miRNAs by qRT-PCR. We also were able to construct a miRNA-mRNA target network that can further understand the molecular mechanisms on the pathogenesis of TiLV and guide future research in developing effective agents and strategies to combat TiLV infections in Tilapia.

Integrated transcriptomic and proteomic analyses of grass carp intestines after vaccination with a double-targeted DNA vaccine of Vibrio mimicus

Intestinal mucosal immunity plays a vital role against Vibrio mimicus infection because it is an enteric pathogen causing serious vibriosis in fish. In the previous studies, we developed an oral double-targeted DNA vaccine of V. mimicus and demonstrated that the vaccine could elicit significantly higher intestinal mucosal immune response than did naked DNA vaccine. But, little is known underlying regulatory molecular mechanisms of the enhanced intestinal mucosal immunity. Here the transcriptome and proteome in the intestines of the grass carps immunized or not with the double-targeted DNA vaccine were investigated by using RNA-seq and iTRAQ-coupled LC-MS/MS. Compared with the control group, a total of 5339 differentially expressed genes (DEGs) and 1173 differentially expressed proteins (DEPs) were identified in the immunized fish intestines. Subsequently, the integrated analysis between transcriptome and proteome data revealed that 250 DEPs were matched with the corresponding DEGs (named associated DEPs/DEGs) at both transcriptome and proteome levels. Fifty of all the associated DEPs/DEGs were immune-related and mainly enriched in phagosome, antigen-processing and presentation, complement and coagulation cascades, NLRs and MAPK signaling pathways via Gene Ontology and KEGG pathway analyses, which suggested the coordination of the five activated pathways was essential to the enhanced intestinal mucosal immune response in the immunized fish. The protein-protein interaction analysis showed that 60 of the 63 immune-related DEPs to form an integrated network. Additionally, randomly selected DEGs and DEPs were respectively validated by quantitative real-time RT-PCR and multiple reaction monitoring (MRM) assay, indicating that the both RNA-Seq and iTRAQ results in the study were reliable. Overall, our comprehensive transcriptome and proteome data provide some key genes and their protein products for further research on the regulatory molecular mechanisms underlying the enhanced intestinal mucosal immunity.

Transcriptomic analysis of oligochaete immune responses to myxosporeans infection: Branchiura sowerbyi infected with Myxobolus cultus

The Myxozoa are endoparasites characterized by a two-host life cycle that typically involves invertebrates and vertebrates as definitive and intermediate hosts, respectively. However, little is known about invertebrate-myxosporean interactions, particularly about patterns of host immune defense. We used RNA-sequencing to identify genes that are possibly involved in the immune responses of the oligochaete Branchiura sowerbyi naturally infected with Myxobolus cultus. De novo assembly of the B. sowerbyi transcriptome yielded 119,031 unigenes, with an average length of 896 bp and an N50 length of 1754 bp. Comparative transcriptome analysis revealed 4059 differentially expressed genes (DEGs) between M. cultus-infected and uninfected B. sowerbyi groups, including 3802 upregulated genes and 257 downregulated genes. Among the B. sowerbyi immune factors implicated in the responses to M. cultus infection, DEGs related to lectins, ubiquitin-mediated proteolysis, phagocytosis, oxidative-antioxidative responses, proteases, and protease inhibitors were upregulated. The expression of some immune-related molecules such as calmodulin, heat shock proteins, antimicrobial peptides, lysenin, and serum amyoid A protein were also significantly upregulated. The expression patterns of 14 immune-related DEGs identified by RNA-seq were validated by quantitative real-time polymerase chain reaction. This study is the first attempt to characterize the B. sowerbyi transcriptome and identify immune-related molecules possibly associated with M. cultus infection. It is also the first report of invertebrate host-myxosporean interactions at the transcriptomic level. Our results will facilitate the elucidation of adaptive evolution mechanisms of myxosporean parasites in the definitive host and the genetic basis for differences in resistance of invertebrate hosts of different genotypes to a myxosporean species.

Nutrient sensing signaling functions as the sensor and regulator of immunometabolic changes in grass carp during Flavobacteriumcolumnare infection

In order to illustrate the immunometabolic changes of fish during bacterial infection, grass carp (Ctenopharyngodon idellus) was injected with Flavobacteriumcolumnare(F.columnare) and then the immune response, nutrient metabolism and related signaling pathways were assayed from 6 h post injection (hpi) to 7 days post injection (dpi). After F.columnare injection, gill lamellae showed obvious fusion and higher mRNA expression levels of pro-inflammatory cytokines. The mRNA expression levels of TNF-α, IL-1β and IL-8 in the head kidney were also significantly upregulated at 6 hpi and 3 dpi. Moreover, the expression of IgZ in the gill was significantly upregulated at 3 dpi and 7 dpi, while the expression of IgM in the head kidney was significantly upregulated at 1 dpi and 3 dpi after F.columnare injection. During bacterial infection, the systematic nutrient metabolism was also significantly affected. Hepatic glycolysis, indicated by GK mRNA expression and PK activity, was significantly upregulated at 1 dpi, while glucogenesis, indicated by PEPCK mRNA expression and enzyme activity, was significantly increased at later time, which resulted in the decreased hepatic glycogen content at 1dpi but increased glycogen content at 7 dpi in the experimental group. LPL, which catalyzed the lipid catabolism, showed decreased mRNA expression and enzyme activity at 6 hpi, while ACC, which was rate-limiting of FA synthesis, was significantly increased at 6 hpi, 3 dpi and 7 dpi. During this process, the nutrient sensing signaling was also significantly affected. TOR signaling in grass carp was significantly activated while ERK signaling was significantly inhibited after F.columnare infection, both of which might function as the sensor and regulator of fish immunometabolic changes.

The important role of phagocytosis and interleukins for nile tilapia (Oreochromis niloticus) to defense infection of Aeromonas hydrophila based on transcriptome analysis

Tilapia is an important economic fish worldwide. It is vital to understand the mechanism of immune response for the prevention and treatment the infection of Aeromonas hydrophila. Based on high-throughput sequencing of Illumina HiSeq™, we found differentially expressed genes in the immune-related pathway were classified into phagosome, cytokine-cytokine receptor interaction and toll-like receptor signaling pathway. Gene Ontology terms were divided into three categories of transporting function, DNA replication activity and energy supply activity. The first one was related to phagocytosis and the process or transporting of antigen driven by tubulins; the second one was to differentiation and proliferation of lymphocyte activated by cytokines; and the former two both needed energy provided by the third one. According to colchicine assay, cross-immune assay, ELISA of interleukins and classical phagocytosis assay, phagocytosis and interleukins were verified to be most important to defense the infection of A. hydrophila.

Transcriptome analysis reveals novel insights in air-breathing magur catfish (Clarias magur) in response to high environmental ammonia

The facultative air-breathing magur catfish (Clarias magur) frequently face different environmental challenges, such as hyper-ammonia, and desiccation stresses in their natural habitats. All these stresses lead to higher accumulation of body ammonia, thereby causing various harmful effects to the fish due to its toxicity. Nonetheless, the mechanisms underlying ammonia-induced toxicity is yet not clear. In the present study, we used RNA sequencing and utilized a modified method for de novo assembly of the transcriptome to provide an exhaustive study on the transcriptomic alterations of magur catfish in response to high environmental ammonia (HEA; 25 mM NH₄Cl). The final contig assembly produced a total of 311,076 unique transcripts (termed as unigenes) with a GC content of 48.3% and the average length of 599 bp. A considerable number of SSR marker associated with these unigenes were also detected. A total of 279,156 transcripts were successfully annotated by using various databases. Comparative transcriptomic analysis revealed a total of 3453 and 19,455 genes were differentially expressed in the liver and brain tissues, respectively, in ammonia-treated fish compared to the control. Enrichment analysis of the differentially expressed genes (DEGs) showed that several GO and KEGG pathway terms were significantly over-represented. Functional analysis of significantly elevated DEGs demonstrated that ammonia stress tolerance of the magur catfish was associated with quite a few pathways related to immune response, oxidative stress, and apoptosis, as well as few transporter proteins involved with ammonia and urea transport. Both liver and brain tissues showed HEA-mediated oxidative damage with consequent activation of antioxidant machinery. However, elevated ROS levels led to an activation of inflammatory cytokines and thus innate immune response in the liver. Conversely, in the brain ROS-mediated irreversible cell damages activated apoptosis via both p53-Bax-Bcl2 and caspase-mediated pathways. The present study provides a novel understanding of the molecular responses of this air-breathing catfish against the ammonia-induced stressors, which could elucidate the underlying mechanisms of adaptation of this facultative air-breather living under various environmental constraints.

Comparative Study of Immune Reaction Against Bacterial Infection From Transcriptome Analysis

Transcriptome analysis is a powerful tool that enables a deep understanding of complicated physiological pathways, including immune responses. RNA sequencing (RNA-Seq)-based transcriptome analysis and various bioinformatics tools have also been used to study non-model animals, including aquaculture species for which reference genomes are not available. Rapid developments in these techniques have not only accelerated investigations into the process of pathogenic infection and defense strategies in fish, but also used to identify immunity-related genes in fish. These findings will contribute to fish immunotherapy for the prevention and treatment of bacterial infections through the design of more specific and effective immune stimulants, adjuvants, and vaccines. Until now, there has been little information regarding the universality and diversity of immune reactions against pathogenic infection in fish. Therefore, one of the aims of this paper is to introduce the RNA-Seq technique for examination of immune responses in pathogen-infected fish. This review also aims to highlight comparative studies of immune responses against bacteria, based on our previous findings in largemouth bass (Micropterus salmoides) against Nocardia seriolae, gray mullet (Mugil cephalus) against Lactococcus garvieae, orange-spotted grouper (Epinephelus coioides) against Vibrio harveyi, and koi carp (Cyprinus carpio) against Aeromonas sobria, using RNA-seq techniques. We demonstrated that only 39 differentially expressed genes (DEGs) were present in all species. However, the number of specific DEGs in each species was relatively higher than that of common DEGs; 493 DEGs in largemouth bass against N. seriolae, 819 DEGs in mullets against L. garvieae, 909 in groupers against V. harveyi, and 1471 in carps against A. sobria. The DEGs in different fish species were also representative of specific immune-related pathways. The results of this study will enhance our understanding of the immune responses of fish, and will aid in the development of effective vaccines, therapies, and disease-resistant strains.

Bioinformatics Analysis of Key Differentially Expressed Genes in Nonalcoholic Fatty Liver Disease Mice Models

Nonalcoholic fatty liver disease (NAFLD) is a global health problem characterized by excessive accumulation of fat in the liver without effect of other pathological factors including hepatitis infection and alcohol abuse. Current studies indicate that gene factors play important roles in the development of NAFLD. However, the molecular characteristics of differentially expressed genes (DEGs) and associated mechanisms with NAFLD have not been well elucidated. Using two microarray data associated with the gene expression profiling in liver tissues of NAFLD mice models, we identified and selected several common key DEGs that contributed to NAFLD. Based on bioinformatics analysis, we discovered that the DEGs were associated with a variety of biological processes, cellular components, and molecular functions and were also related to several significant pathways. Via pathway crosstalk analysis based on overlapping DEGs, we observed that the identified pathways could form large and complex crosstalk networks. Besides, large and complex protein interaction networks of DEGs were further constructed. In addition, many hub host factors with a high degree of connectivity were identified based on interaction networks. Furthermore, significant modules in interaction networks were found, and the DEGs in the identified modules were found to be enriched with distinct pathways. Taken together, these results suggest that the key DEGs, associated pathways, and modules contribute to the development of NAFLD and might be used as novel molecular targets for the treatment of NAFLD.

Combined transcriptomic/proteomic analysis of crucian carp Carassius auratus gibelio in cyprinid herpesvirus 2 infection

Cyprinid herpesvirus 2 (CyHV-2) is a pathogen of herpesviral hematopoietic necrosis disease of crucian carp. Our study aimed to investigate the molecular mechanisms and immune response at the mRNA and protein levels in head kidney during CyHV-2 infection. Three days after infection with CyHV-2, 7085 differentially expressed genes were identified by transcriptome sequencing, of which 3090 were up-regulated and 3995 were down-regulated. And 338 differentially expressed proteins including 277 up-regulated and 61 down-regulated were identified using tandem mass tag labeling followed by liquid chromatography tandem mass spectrometry. Notably, 128 differentially co-expressed genes at mRNA and protein levels (cDEGs) were reliably quantified, including 86 co-up-regulated and 42 co-down-regulated. In addition, 10 cDGEs in the above pathways were selected for qRT-PCR to confirm the validity of the transcriptome and proteome changes by showing that RIG-I, MDA5, LGP2, FAS, PKR and PKZ up-regulated and Integrin α, Integrin β2, NCF2 and NCF4 down-regulated. This indicated that after CyHV-2 infection, the herpes simplex infection pathway, RIG-I like receptor signaling pathway, necroptosis pathway and p53 signaling pathway were activated and the phagosome pathway was suppressed. Our findings reveal the pathogenesis and the host immune mechanism of CyHV-2 infection of crucian carp.

Transcriptomic analysis of liver from grass carp (Ctenopharyngodon idellus) exposed to high environmental ammonia reveals the activation of antioxidant and apoptosis pathways

High concentration of ammonia in aquatic system leads to detrimental effects on the health of aquatic animals. However, the mechanism underlying ammonia-induced toxicity is still not clear. To better understand the mechanism of ammonia toxicity effects on fish, juvenile grass carp was employed in the present study. RNA high-throughput sequencing technique was applied to analyze the total RNAs extracted from the liver of fish after 8 h post exposure to the water containing 2 mM NH₄HCO₃ which experimentally mimicked the high environmental ammonia (HEA). A total of 49,971,114 and 53,826,986 clean reads were obtained in control and 2 mM HEA group, respectively, in which there were 911 differently expressed genes (DEGs) including 563 up-regulated and 348 down-regulated genes. In addition, 10 DEGs were validated by quantitative PCR. These DEGs were involved in several pathways related with oxidative stress or apoptosis. Further analysis on oxidative stress, histopathology and cellular apoptosis in grass carp liver after HEA exposure revealed interesting findings. Increased reactive oxygen species (ROS) content and superoxide dismutase (SOD) activity together with the decreased catalase (CAT) activity were detected, which may be effected by DEGs and related pathways such as FOXO signaling pathway. The histopathology and TUNEL assays results confirmed that apoptosis was induced in liver when fish had suffered HEA. Combined with the results of transcriptomic experiments, c-Myc-Bax-Caspase9 apoptosis pathway could be involved in grass carp liver apoptosis induced by ammonia stress.