Transcriptome Analysis Based on RNA-Seq in Understanding Pathogenic Mechanisms of Diseases and the Immune System of Fish: A Comprehensive Review

Transcriptome profiling of lumpfish (Cyclopterus lumpus) head kidney to Renibacterium salmoninarum at early and chronic infection stages

Renibacterium salmoninarum causes Bacterial Kidney Disease (BKD) in several fish species. Atlantic lumpfish, a cleaner fish, is susceptible to R. salmoninarum. To profile the transcriptome response of lumpfish to R. salmoninarum at early and chronic infection stages, fish were intraperitoneally injected with either a high dose of R. salmoninarum (1 × 109 cells dose-1) or PBS (control). Head kidney tissue samples were collected at 28- and 98-days post-infection (dpi) for RNA sequencing. Transcriptomic profiling identified 1971 and 139 differentially expressed genes (DEGs) in infected compared with control samples at 28 and 98 dpi, respectively. At 28 dpi, R. salmoninarum-induced genes (n = 434) mainly involved in innate and adaptive immune response-related pathways, whereas R. salmoninarum-suppressed genes (n = 1537) were largely connected to amino acid metabolism and cellular processes. Cell-mediated immunity-related genes showed dysregulation at 98 dpi. Several immune-signalling pathways were dysregulated in response to R. salmoninarum, including apoptosis, alternative complement, JAK-STAT signalling, and MHC-I dependent pathways. In summary, R. salmoninarum causes immune suppression at early infection, whereas lumpfish induce a cell-mediated immune response at chronic infection. This study provides a complete depiction of diverse immune mechanisms dysregulated by R. salmoninarum in lumpfish and opens new avenues to develop immune prophylactic tools to prevent BKD.

Coping with salinity extremes: Gill transcriptome profiling in the black-chinned tilapia (Sarotherodon melanotheron)

Steeper and sometimes extreme salinity gradients increasingly affect aquatic organisms because of climate change. Hypersalinity habitats demand powerful physiological adaptive strategies. Few teleost species have the capacity to spend their whole life cycle in salinities way over seawater levels. Focusing on the multifunctional gill, we unraveled the tilapia S. melanotheron key strategies to cope with different environmental conditions, ranging from freshwater up to hypersaline habitats. De novo transcriptome assembly based on RNAseq allowed for the analysis of 40,967 annotated transcripts among samples collected in three wild populations at 0, 40 and 80 ‰. A trend analysis of the expression patterns revealed responses across the salinity gradient with different gene pathways involved. Genes linked to ion transport, pH regulation and cell surface receptor signaling were mainly upregulated in the high salinity habitat. We identified tight junction proteins that were critical in high salinity habitats and that were different from the well-known tightening junctional proteins identified and expressed in fresh water. Expression profiles also suggest a change in the vascular tone that could be linked to an osmorespiratory compromise not only in fresh water, but also in high salinity environments. A striking downregulation of genes linked to the immune system and to the heat shock response was observed suggesting an energetic trade-off between immunity and acclimation/adaptation in the hypersaline habitat. The high expression of transcripts coding for immune and heat shock response in the freshwater habitat suggests the establishment of powerful mechanisms to protect gills from environmental threats and to maintain protein integrity. Non-directional expression trends were also detected with an upregulation of genes only in the hypersaline habitat (80 ‰) or only in the marine habitat (40 ‰). Unravel physiological strategies in S. melanotheron populations will help to better understand the molecular basis of fish euryhalinity in salinity-contrasted environments.

Effects of phoxim on antibacterial infection of silver carp

The efficacy of phoxim in treating bacterial sepsis in silver carp is significant, yet its underlying mechanism remains elusive. This study aimed to establish a model of Aeromonas veronii infection in silver carp and subsequently treat the infected fish with 10 μg/L phoxim. Kidney and intestine samples from silver carp were collected for transcriptome analysis and assessment of intestinal microbial composition, with the aim of elucidating the mechanism underlying the efficacy of phoxim in treating bacterial sepsis in silver carp. The results of transcriptome and intestinal microbial composition analysis of silver carp kidney indicated that A. veronii infection could up-regulate the expression of il1β, il6, nos2, ctsl, casp3 et al., which means, signifying that the kidney of silver carp would undergo inflammation, induce apoptosis, and alter the composition of intestinal microorganisms. Phoxim immersion might enhance the energy metabolism of silver carp and change its intestinal microbial composition, potentially elevating the antibacterial infection resistance of silver carp. These findings may contribute to an understanding of how phoxim can effectively treat bacterial sepsis in silver carp.

Chromosome level genome assembly and transcriptome analysis of E11 cells infected by tilapia lake virus

The E11 cell line, derived from striped snakehead fish (Channa striata), possesses a distinctive feature: it is persistently infected with a C-type retrovirus. Notably, it exhibits high permissiveness to piscine nodavirus and the emerging tilapia lake virus (TiLV). Despite its popularity in TiLV research, the absence of genome assembly for the E11 cell line and Channa striata has constrained research on host-virus interactions. This study aimed to fill this gap by sequencing, assembling, and annotating the E11 cell line genome. Our efforts yielded a 600.5 Mb genome including 24 chromosomes with a BUSCO score of 98.8%. In addition, the complete proviral DNA sequence of snakehead retrovirus (SnRV) was identified in the E11 cell genome. Comparative genomic analysis between the E11 cell line and another snakehead species Channa argus revealed the loss of many immune-related gene families in the E11 cell genome, indicating a compromised immune response. We also conducted transcriptome analysis of mock- and TiLV-infected E11 cells, unveiling new perspectives on virus-virus and host-virus interactions. The TiLV infection suppressed the high expression of SnRV in E11 cells, and activated some other endogenous retroviruses. The protein-coding gene comparison revealed a pronounced up-regulation of genes involved in immune response, alongside a down-regulation of genes associated with specific metabolic processes. In summary, the genome assembly and annotation of the E11 cell line provide valuable resources to understand the SnRV and facilitate further studies on nodavirus and TiLV. The RNA-seq profiles shed light on the cellular mechanisms employed by fish cells in response to viral challenges, potentially guiding the development of therapeutic strategies against TiLV in aquaculture. This study also provides the first insights into the viral transcriptome profiles of endogenous SnRV and evading TiLV, enhancing our understanding of host-virus interactions in fish.

Unveiling Immunomodulatory Effects of Euglena gracilis in Immunosuppressed Mice: Transcriptome and Pathway Analysis

The immunomodulatory effects of Euglena gracilis (Euglena) and its bioactive component, β-1,3-glucan (paramylon), have been clarified through various studies. However, the detailed mechanisms of the immune regulation remain to be elucidated. This study was designed not only to investigate the immunomodulatory effects but also to determine the genetic mechanisms of Euglena and β-glucan in cyclophosphamide (CCP)-induced immunosuppressed mice. The animals were orally administered saline, Euglena (800 mg/kg B.W.) or β-glucan (400 mg/kg B.W.) for 19 days, and CCP (80 mg/kg B.W.) was subsequently administered to induce immunosuppression in the mice. The mice exhibited significant decreases in body weight, organ weight, and the spleen index. However, there were significant improvements in the spleen weight and the spleen index in CCP-induced mice after the oral administration of Euglena and β-glucan. Transcriptome analysis of the splenocytes revealed immune-related differentially expressed genes (DEGs) regulated in the Euglena- and β-glucantreated groups. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that pathways related with interleukin (IL)-17 and cAMP play significant roles in regulating T cells, B cells, and inflammatory cytokines. Additionally, Ptgs2, a major inflammatory factor, was exclusively expressed in the Euglena-treated group, suggesting that Euglena's beneficial components, such as carotenoids, could regulate these genes by influencing immune lymphocytes and inflammatory cytokines in CCP-induced mice. This study validated the immunomodulatory effects of Euglena and highlighted its underlying mechanisms, suggesting a positive contribution to the determination of phenotypes associated with immune-related diseases and the research and development of immunotherapies.

Mechanism of sturgeon intestinal inflammation induced by Yersinia ruckeri and the effect of florfenicol intervention

The mechanism by which Y. ruckeri infection induces enteritis in Chinese sturgeon remains unclear, and the efficacy of drug prevention and control measures is not only poor but also plagued with numerous issues. We conducted transcriptomic and 16 S rRNA sequencing analyses to examine the differences in the intestinal tract of hybrid sturgeon before and after Y. ruckeri infection and florfenicol intervention. Our findings revealed that Y. ruckeri induced the expression of multiple inflammatory factors, including il1β, il6, and various chemokines, as well as casp3, casp8, and multiple tumor necrosis factor family members, resulting in pathological injury to the body. Additionally, at the phylum level, the relative abundance of Firmicutes and Bacteroidota increased, while the abundance of Plesiomonas and Cetobacterium decreased at the genus level, altering the composition of the intestinal flora. Following florfenicol intervention, the expression of multiple apoptosis and inflammation-related genes was down-regulated, promoting tissue repair. However, the flora became further dysregulated, increasing the risk of infection. In conclusion, our analysis of the transcriptome and intestinal microbial composition demonstrated that Y. ruckeri induces intestinal pathological damage by triggering apoptosis and altering the composition of the intestinal microbiota. Florfenicol intervention can repair pathological damage, but it also exacerbates flora imbalance, leading to a higher risk of infection. These findings help elucidate the molecular mechanism of Y. ruckeri-induced enteritis in sturgeon and evaluate the therapeutic effect of drugs on intestinal inflammation in sturgeon.

Effects of different combinations of koumine and gelsemine on growth performance, intestinal health, and transcriptome of Cyprinus carpio

Koumine (KM) and gelsemine (GS) have shown significant benefits in livestock production, but their potential in aquaculture remains largely unexplored. This study examined the impact of different KM and GS combinations as feed additives on C. carpio (90 fish per group, initial weight 1.95 ± 0.08 g). KM and GS were introduced in ratios of 2:2 (mg/kg), 2:1 (mg/kg), and 2:0.67 (mg/kg) over a 10-week aquaculture experiment. The results demonstrate that the 2:1 (mg/kg) group increases the villus length, muscular layer thickness, crude protein, and crude fat content. Regarding fatty acid content, KM and GS enhance the levels of various fatty acids, including the total saturated fatty acid and total monounsaturated fatty acid. Additionally, KM and GS improve the composition and function of the intestinal microbiota. The 2:1 (mg/kg) group significantly elevates the enzymatic activities of SOD, MDA, CAT and upregulates the expression of immune-related genes such as toll-like receptor 2, transforming growth factor β, and glutathione S-transferase. Transcriptomic analysis suggests that KM and GS may have potential benefits for nutrient utilization and immune regulation in C. carpio. In summary, this study provides valuable insights into the use of KM and GS as feed additives in aquaculture.

Comparative transcriptomic analyses of brain-liver-muscle in channel catfish (Ictalurus punctatus) with differential growth rate

Growth is an important economically trait for aquatic animals. The popularity of farmed channel catfish (Ictalurus punctatus) in China has recently surged, prompting a need for research into the genetic mechanisms that drive growth and development to expedite the selection of fast-growing variants. In this study, the brain, liver and muscle transcriptomes of channel catfish between fast-growing and slow-growing groups were analyzed using RNA-Seq. Totally, 63, 110 and 86 differentially expressed genes (DEGs) were from brain, liver and muscle tissues. DEGs are primarily involved in growth, development, metabolism and immunity, which are related to the growth regulation of channel catfish, such as growth hormone receptor b (ghrb), fibroblast growth factor receptor 4 (fgfr4), bone morphogenetic protein 1a (bmp1a), insulin-like growth factor 2a (igf2a), collagen, type I, alpha 1a (col1a1a), acyl-CoA synthetase long chain family member 2 (acsl2) and caveolin 1 (cav1). This study advances our knowledge of the genetic mechanisms accounting for differences in growth rate and offers crucial gene resources for future growth-related molecular breeding programs in channel catfish.

Comparative transcriptome analysis reveals potential regulatory mechanisms of genes and immune pathways following Vibrio harveyi infection in red drum (Sciaenops ocellatus)

Red drum (Sciaenops ocellatus), as an important economical marine fish, has been affected by various bacterial diseases in recent years. Vibrio harveyi cause fatal vibriosis in S. ocellatus, leading to massive mortality and causing significant setbacks in aquaculture. However, the regulatory mechanisms of S. ocellatus response to V. harveyi infection are poorly understood. In this regard, we performed transcriptomic analysis with head kidney tissues of S. ocellatus after V. harveyi infection from 12 h to 48 h to reveal genes, gene expression profiles, and pathways involved in immune and inflammation responses. Specifically, a total of 9,599, 5,728, and 7144 differentially expressed genes (DEGs) were identified after V. harveyi infection at 12 h, 24 h, and 48 h, respectively, and 1,848 shared DEGs have been identified from the above three comparison groups. Subsequent pathway analysis revealed that the shared DEGs following V. harveyi were involved in complement and coagulation cascades (C1R, C1QC, C3, C4, C5, C7, C8A, C8B, C8G, C9, CFB, CFH, and CFI), MAPK signaling pathway, chemokine signaling pathway (CCL19, CXCL8, CXCL12, CXCL14, CCR4, CCR7, and CXCR2), PPAR signaling pathway (PPAR-α, PPAR-γ and PPAR-β), and TNF signaling pathway. Finally, the expression patterns of DEGs in head kidney tissues and S. ocellatus macrophages were validated by qRT-PCR, suggesting the reliability of RNA sequencing for gene expression analysis. This dynamic transcriptome analyses provided insights into gene expression regulation and immune related pathways involved in S. ocellatus after V. harveyi infection, and provides useful information for further study on the immune defense mechanisms in S. ocellatus as well as other teleost species.

Pivotal Role of GSTO2 in Ferroptotic Neuronal Injury After Intracerebral Hemorrhage

Previous research has found that an adaptive response to ferroptosis involving glutathione peroxidase 4 (GPX4) is triggered after intracerebral hemorrhage. However, little is known about the mechanisms underlying adaptive responses to ferroptosis. To explore the mechanisms underlying adaptive responses to ferroptosis after intracerebral hemorrhage, we used hemin-treated HT22 cells to mimic brain injury after hemorrhagic stroke in vitro to evaluate the antioxidant enzymes and performed bioinformatics analysis based on the mRNA sequencing data. Further, we determined the expression of GSTO2 in hemin-treated hippocampal neurons and in a mouse model of hippocampus-intracerebral hemorrhage (h-ICH) by using Western blot. After hemin treatment, the antioxidant enzymes GPX4, Nrf2, and glutathione (GSH) were upregulated, suggesting that an adaptive response to ferroptosis was triggered. Furthermore, we performed mRNA sequencing to explore the underlying mechanism, and the results showed that 2234 genes were differentially expressed. Among these, ten genes related to ferroptosis (Acsl1, Ftl1, Gclc, Gclm, Hmox1, Map1lc3b, Slc7a11, Slc40a1, Tfrc, and Slc39a14) were altered after hemin treatment. In addition, analysis of the data retrieved from the GO database for the ten targeted genes showed that 20 items on biological processes, 17 items on cellular components, and 19 items on molecular functions were significantly enriched. Based on the GO data, we performed GSEA and found that the glutathione metabolic process was significantly enriched in the hemin phenotype. Notably, the expression of glutathione S-transferase omega (GSTO2), which is involved in glutathione metabolism, was decreased after hemin treatment, and overexpression of Gsto2 decreased lipid reactive oxygen species level in hemin-exposed HT22 cells. In addition, the expression of GSTO2 was also decreased in a mouse model of hippocampus-intracerebral hemorrhage (h-ICH). The decreased expression of GSTO2 in the glutathione metabolic process may be involved in ferroptotic neuronal injury following hemorrhagic stroke.

Effects of sublethal fipronil exposure on cross-generational functional responses and gene expression in Binodoxys communis

The effective systemic insecticide fipronil is widely used on a variety of crops and in public spaces to control insect pests. Binodoxys communis (Gahan) (Hymenoptera: Braconidae) is the dominant natural enemy of Aphis gossypii Glover (Homoptera: Aphididae), an important cotton pest, and has good efficiency in inhibiting aphid populations. The direct effects of environmental residues of sublethal fipronil doses on adult B. communis have not previously been reported. This study therefore aimed to evaluate the side effects and transcriptomic impacts of sublethal fipronil doses on B. communis. The results showed that exposure to the LC10 dose of fipronil significantly reduced the survival rate and parasitism rate of the F0 generation, but did not affect these indicators in the F1 generation. The LC25 dose did not affect the survival or parasitic rates of the F0 generation, but did significantly reduce the survival rate of F1 generation parasitoids. These results indicated that sublethal doses of fipronil affected B. communis population growth. Transcriptome analysis showed that differentially expressed genes (DEGs) in B. communis at 1 h after treatment were primarily enriched in pathways associated with fatty acid elongation, biosynthesis of fatty acids, and fatty acid metabolism. DEGs at 3 days after treatment were mainly enriched in ribosomal functions, glycolysis/gluconeogenesis, and tyrosine metabolism. Six DEGs (PY, ELOVL, VLCOAR, MRJP1, ELOVL AAEL008004-like, and RPL13) were selected for validation with real-time fluorescent quantitative PCR. This is the first report of sublethal, trans-generational, and transcriptomic side effects of fipronil on the dominant parasitoid of A. gossypii. The results of this study show that adaptation of parasitoids to high concentrations of pesticides may be at the expense of their offspring. These findings broaden our overall understanding of the intergenerational adjustments used by insects to respond to pesticide stress and call for risk assessments of the long-term impacts and intergenerational effects of other pesticides.

Pathological observation and transcriptomic analysis of thymus injury in PRRSV-infected piglets

The thymus, the central immune organ in mammals, plays an important role in immune defense. Porcine reproductive and respiratory syndrome virus (PRRSV) infection in piglets can cause thymus injury and immunosuppression. However, the mechanisms of thymus injury remain unknown. This study was aimed at investigating the specific manifestations of thymus injury through the construction of a PRRSV-infected piglet model and histopathological observation. In this study, fourteen 40-day-old PRRSV-free piglets were randomly divided into two groups, eleven of which were intramuscularly injected with 3 mL of PRRSV WUH3 virus suspension (106 PFU /mL) in the infection group, and three of which were sham-inoculated with 3 mL of RPMI-1640 medium in the control group. Clinical necropsy and samples collection were performed on day 8 after artificial infection. With the Illumina platform, the transcriptomes of piglet thymus tissues from infected and control piglets were sequenced to explore the relationships of differentially expressed genes (DEGs) and signaling pathways with thymus injury. The immune organs of PRRSV-infected piglets were severely damaged. The histopathological findings in the thymus indicated that PRRSV infection was associated with a large decrease in lymphocytes, cell necrosis and cell apoptosis; an increase in blood vessels and macrophages; thymic corpuscle hyperplasia; and interstitial widening of the thymic lobules. The transcriptomic analysis results revealed that the Gene Ontology functions of DEGs were enriched primarily in biological processes such as angiogenesis, regulation of angiogenesis and positive regulation of cell migration. Moreover, greater numbers of blood vessels and macrophages were observed in the thymus in PRRSV-infected than control piglets. KEGG pathway enrichment analysis revealed that the DEGs were significantly enriched in the Toll-like receptor signaling pathway, chemokine signaling pathway, IL-17 signaling pathway and TNF signaling pathway. The expression of TLR8, IRF5, the chemokines CCL2, CCL3L1 and CCL5; and their receptors CCR1, CCR2 and CCR5 was significantly up-regulated in PRRSV infection, thus suggesting that these cytokines were associated with the pathological processes of thymus injury.

DNA Methylation Profiling of Ovarian Tissue of Climbing Perch (Anabas testudienus) in Response to Monocrotophos Exposure

Epigenetic modifications like DNA methylation can alter an organism's phenotype without changing its DNA sequence. Exposure to environmental toxicants has the potential to change the resilience of aquatic species. However, little information is available on the dynamics of DNA methylation in fish gonadal tissues in response to organophosphates. In the present work, reduced-representation bisulfite sequencing was performed to identify DNA methylation patterns in the ovarian tissues of Anabas testudienus exposed to organophosphates, specifically monocrotophos (MCP). Through sequencing, an average of 41,087 methylated cytosine sites were identified and distributed in different parts of genes, i.e., in transcription start sites (TSS), promoters, exons, etc. A total of 1058 and 1329 differentially methylated regions (DMRs) were detected as hyper-methylated and hypo-methylated in ovarian tissues, respectively. Utilizing whole-genome data of the climbing perch, the DMRs, and their associated overlapping genes revealed a total of 22 genes within exons, 45 genes at transcription start sites (TSS), and 218 genes in intergenic regions. Through gene ontology analysis, a total of 16 GO terms particularly involved in ovarian follicular development, response to oxidative stress, oocyte maturation, and multicellular organismal response to stress associated with reproductive biology were identified. After functional enrichment analysis, relevant DMGs such as steroid hormone biosynthesis (Cyp19a, 11-beta-HSD, 17-beta-HSD), hormone receptors (ar, esrrga), steroid metabolism (StAR), progesterone-mediated oocyte maturation (igf1ar, pgr), associated with ovarian development in climbing perch showed significant differential methylation patterns. The differentially methylated genes (DMGs) were subjected to analysis using real-time PCR, which demonstrated altered gene expression levels. This study revealed a molecular-level alteration in genes associated with ovarian development in response to chemical exposure. This work provides evidence for understanding the relationship between DNA methylation and gene regulation in response to chemicals that affect the reproductive fitness of aquatic animals.

Subtle transcriptomic response of Eurasian perch (Perca fluviatilis) associated with Triaenophorus nodulosus plerocercoid infection

Determining the physiological effects of parasites and characterizing genes involved in host responses to infections are essential to improving our understanding of host-parasite interactions and their ecological and evolutionary consequences. This task, however, is complicated by high diversity and complex life histories of many parasite species. The use of transcriptomics in the context of wild-caught specimens can help ameliorate this by providing both qualitative and quantitative information on gene expression patterns in response to parasites in specific host organs and tissues. Here, we evaluated the physiological impact of the widespread parasite, the pike tapeworm (Triaenophorus nodulosus), on its second intermediate host, the Eurasian perch (Perca fluviatilis). We used an RNAseq approach to analyse gene expression in the liver, the target organ of T. nodulosus plerocercoids, and spleen which is one of the main immune organs in teleost fishes. We compared perch collected from multiple lakes consisting of individuals with (n = 8) and without (n = 6) T. nodulosus plerocercoids in the liver. Results revealed a small number of differentially expressed genes (DEGs, adjusted p-value ≤0.05) in both spleen (n = 22) and liver (n = 10). DEGs in spleen consisted of mostly upregulated immune related genes (e.g., JUN, SIK1, THSB1), while those in the liver were often linked to metabolic functions (e.g., FABP1, CADM4, CDAB). However, Gene Ontology (GO) analysis showed lack of functional enrichment among DEGs. This study demonstrates that Eurasian perch displays a subtle response at a gene expression level to T. nodulosus plerocercoid infection. Given that plerocercoids are low-metabolic activity transmission stages, our results suggest that moderate T. nodulosus plerocercoid infection most likely does not provoke an extensive host immune response and have relatively low physiological costs for the host. Our findings illustrate that not all conspicuous infections have severe effects on host gene regulation.

Transcriptomic Variation of Amphiprion Percula (Lacepède, 1802) in Response to Infection with Cryptocaryon Irritans Brown, 1951

Cryptocaryon irritans (Brown 1951) frequently infect the Pomacentridae fishes causing severe economic losses. However, the anti-C. irritans' molecular mechanism in these fishes remains largely unknown. To address this issue, we conducted RNA-Seq for C. irrtians-infected gills of the clownfish Amphiprion percula (Lacepède 1802) at the early (day 1) and late (day 3) stages of infection. A total of 1655 differentially expressed genes (DEGs) were identified. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of DEGs showed a vast genetic variation related to the following aspects: ECM-receptor interaction, P13K-Akt signalling, cytokine-cytokine receptor interaction, and endocytosis. During the early phase of infection, key genes involved in ATP production, energy homeostasis, and stress control were abruptly increased. In the late phase, however, acute response molecules of the peripheral nervous system (synaptic transmission and local immunity), metabolic system triggering glycogen synthesis, energy maintenance, and osmoregulation were found to be critical. The highest number of upregulated genes (URGs) recovered during the early phase was included under the 'biological process' category, which primarily functions as response to stimuli, signalling, and biological regulation. In the late phase, most of the URGs were related to gene regulation and immune system processes under 'molecular function' category. The immune-related URGs of early infection include major histocompatibility complex (MHC) class-II molecules apparently triggering CD4+ T-cell-activated Th responses, and that of late infection include MHC class-1 molecules for the possible culmination of CD8+ T-cell triggered cytotoxicity. The high level of genic single nucleotide polymorphisms (SNPs) identified during the late phase of infection is likely to influence their susceptibility to secondary infection. In summary, the identified DEGs and their related metabolic and immune-related pathways and the SNPs may provide new insights into coordinating the immunological events and improving resistance in Pomacentridae fishes against C. irritans.

Transcriptome analysis reveals immune regulation in the spleen of koi carp (Cyprinus carpio Koi) during Aeromonas hydrophila infection

A. hydrophila (Aeromonas hydrophila) is one of the most hazardous pathogenic microorganisms threatening the aquaculture industry and exhibits zoonotic-like characteristics. This study was designed to investigate the differential gene expression and pathway enrichment in the spleen of koi carp (Cyprinus carpio koi) upon A. hydrophila infection. The Illumina NovaSeq 6000 sequencing platform was used to identify 252 DEGs (differentially expressed genes), including 112 upregulated genes and 140 downregulated genes, in the spleens of koi carp challenged with A. hydrophila compared to those in the spleens of koi carp treated with PBS (phosphate-buffered saline). DEGs were shown to be involved in 133 pathways by KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis. Numerous immunological disease-related pathways, such as the immune defense network for IgA production, Staphylococcus aureus infection, and antigen processing and presentation, were enriched in the DEGs. In addition, the expression levels of 10 randomly screened DEGs, including the inflammatory factor nlrp3 (NOD-like receptor family pyrin domain containing 3), cytokine il-8 (interleukin-8), c2 (complement c2), c3 (complement c3), and the lipid mediator cox1 (cyclooxygenase-1), were compared by qPCR. The results showed that six genes, including il-8, cox1, and nlrp3, were upregulated according to both RNA-seq and qPCR validation, while four, including c2 and c3, showed downregulated expression. This result verified a strong correlation between the RNA-seq and qPCR datasets at the expression level. Moreover, this study provided splenic transcriptome data for koi carp during A. hydrophila infection and provided theoretical support for future drug development.

Hemolysin from Aeromonas hydrophila enhances the host's serum enzyme activity and regulates transcriptional responses in the spleen of Cyprinus rubrofuscus

Aeromonas hydrophila is a conditional pathogen impacting public hygiene and safety. Hemolysin is a virulence factor of Aeromonas hydrophila that causes erythrocyte hemolysis, yet its transcriptional response to Cyprinus rubrofuscus remains unknown. Our investigation confirmed the hemolysis of hemolysin from A. hydrophila. Serum enzyme activity was evaluated weekly after C. rubrofuscus were immunized with hemolysin Ahh1. The results showed that the hemolysin enhances the serum superoxide dismutase (SOD), lysozyme (LZM), and catalase (CAT) activity, which reached a maximum on day 14. To elucidate the molecular interaction between hemolysin from A. hydrophila and the host, we performed transcriptome sequencing on the spleen of C. rubrofuscus 14 days post hemolysin infection. The total number of clean reads was 41.37 Gb, resulting in 79,832 unigenes with an N50 length of 1863 bp. There were 1982 significantly differentially expressed genes (DEGs), including 1083 upregulated genes and 899 downregulated genes. Transcript levels of the genes, such as LA6BL, CD2, and NLRC5, were significantly downregulated, while those of IL11, IL1R2, and IL8 were dramatically upregulated. The DEGs were mainly enriched in the immune disease, viral protein interaction with cytokine and cytokine receptor, and toll-like receptor pathways, suggesting that hemolysin stimulation can activate the transcriptional responses. RT-qPCR experiments results of seven genes, IL-8, STAT2, CTSK, PRF1, CXCL9, TLR5, and SACS, showed that their expression was highly concordant with RNA-seq data. We clarified for the first time the key genes and signaling pathways response to hemolysin from A. hydrophila, which offers strategies for treating and preventing diseases.

Transcriptomic Analysis of Liver Tissue of Black Sea Bass (Centropristis striata) Exposed to High Nitrogen Environment

The black sea bass, Centropristis striata, is a potential candidate for commercial aquaculture. Due to inadequate removal of nitrogen in its breeding environment, C. striata exhibits increased nitrate concentration, which can cause acute toxicity, including energy metabolism damage and tissue damage. Therefore, RNA-seq technology was applied to characterize genes associated with toxicity tolerance under nitrate stress. The nitrate treatment caused significant changes in a total of 8920 genes, of which 2949 genes were up-regulated and 5971 genes were down-regulated. It was found that significantly enriched GO terms and KEGG were associated with blood microparticles, inhibitors of enzyme activity, and complement and coagulation cascade pathways. Furthermore, through bioinformatics analysis, it was found that these different pathways obtained in GO and KEGG enrichment analysis were mostly related to the immune and inflammatory response of fish. This study expands our understanding of the mechanism of nitrate stress affecting the liver function of C. striata.

Differentially expressed transcripts of Tetracapsuloides bryosalmonae (Cnidaria) between carrier and dead-end hosts involved in key biological processes: novel insights from a coupled approach of FACS and RNA sequencing

Tetracapsuloides bryosalmonae is a malacosporean endoparasite that infects a wide range of salmonids and causes proliferative kidney disease (PKD). Brown trout serves as a carrier host whereas rainbow trout represents a dead-end host. We thus asked if the parasite adapts to the different hosts by changing molecular mechanisms. We used fluorescent activated cell sorting (FACS) to isolate parasites from the kidney of brown trout and rainbow trout following experimental infection with T. bryosalmonae. The sorted parasite cells were then subjected to RNA sequencing. By this approach, we identified 1120 parasite transcripts that were expressed differentially in parasites derived from brown trout and rainbow trout. We found elevated levels of transcripts related to cytoskeleton organisation, cell polarity, peptidyl-serine phosphorylation in parasites sorted from brown trout. In contrast, transcripts related to translation, ribonucleoprotein complex biogenesis and subunit organisation, non-membrane bounded organelle assembly, regulation of protein catabolic process and protein refolding were upregulated in rainbow trout-derived parasites. These findings show distinct molecular adaptations of parasites, which may underlie their distinct outcomes in the two hosts. Moreover, the identification of these differentially expressed transcripts may enable the identification of novel drug targets that may be exploited as treatment against T. bryosalmonae. We here also describe for the first time how FACS based isolation of T. bryosalmonae cells from infected kidney of fish fosters research and allows to define differentially expressed parasite transcripts in carrier and dead-end fish hosts.

Transcriptome analysis reveals genes connected to temperature adaptation in juvenile antarctic krill Euphausia superba

BACKGROUND

The Antarctic krill, Euphausia superba (E. superba), is a key organism in the Antarctic marine ecosystem and has been widely studied. However, there is a lack of transcriptome data focusing on temperature responses.

METHODS

In this study, we performed transcriptome sequencing of E. superba samples exposed to three different temperatures: -1.19 °C (low temperature, LT), - 0.37 °C (medium temperature, MT), and 3 °C (high temperature, HT).

RESULTS

Illumina sequencing generated 772,109,224 clean reads from the three temperature groups. In total, 1,623, 142, and 842 genes were differentially expressed in MT versus LT, HT versus LT, and HT versus MT, respectively. Moreover, Kyoto Encyclopedia of Genes and Genomes analysis revealed that these differentially expressed genes were mainly involved in the Hippo signaling pathway, MAPK signaling pathway, and Toll-like receptor signaling pathway. Quantitative reverse-transcription PCR revealed that ESG037073 expression was significantly upregulated in the MT group compared with the LT group, and ESG037998 expression was significantly higher in the HT group than in the LT group.

CONCLUSIONS

This is the first transcriptome analysis of E. superba exposed to three different temperatures. Our results provide valuable resources for further studies on the molecular mechanisms underlying temperature adaptation in E. superba.

Transcriptomic analyses reveal the potential antibacterial mechanism of citral against Staphylococcus aureus

Background

The emergence of multi-drug resistant Staphylococcus aureus (S. aureus) has posed a challenging clinical problem for treating its infection. The development of novel or new antibacterial agents becomes one of the useful methods to solve this problem, and has received more attention over the past decade. Citral is reported to have antibacterial activity against S. aureus, but its mechanism is yet entirely clear.

Methods

To reveal the antibacterial mechanism of citral against S. aureus, comparative transcriptomic analysis was carried out to analyze the gene expression differences between the citral-treated and untreated groups. The changes of protein, adenosine triphosphate (ATP) and reactive oxygen species (ROS) content in S. aureus caused by citral were also examined.

Results

Six hundred and fifty-nine differentially expressed genes were obtained according to the comparative transcriptomic analysis, including 287 up-regulated genes and 372 down-regulated genes. The oxidoreductase activity and fatty acid degradation pathway were enriched in up-regulated genes, and ribosome and S. aureus infection pathway were enriched in down-regulated genes. Meanwhile, physiological trials revealed a decline in ATP and protein levels, but an increase in ROS content within the citral-treated group. Thus, it can be inferred that the antibacterial effects of citral against S. aureus were likely due to its ability to decrease ATP content by down-regulating ATP synthase genes (atpD and atpG), reduce protein content, induce cell membrane and cell wall damages, accumulate ROS, and down-regulate virulence factor genes to reduce pathogenicity.

Conclusion

These findings revealed the antibacterial mechanism of citral was likely a type of multi-target mode that affected multiple molecular processes in S. aureus, which lays the groundwork for further exploitation of citral as a therapeutic candidate against S. aureus infections.

MyBrain-Seq: A Pipeline for MiRNA-Seq Data Analysis in Neuropsychiatric Disorders

High-throughput sequencing of small RNA molecules such as microRNAs (miRNAs) has become a widely used approach for studying gene expression and regulation. However, analyzing miRNA-Seq data can be challenging because it requires multiple steps, from quality control and preprocessing to differential expression and pathway-enrichment analyses, with many tools and databases available for each step. Furthermore, reproducibility of the analysis pipeline is crucial to ensure that the results are accurate and reliable. Here, we present myBrain-Seq, a comprehensive and reproducible pipeline for analyzing miRNA-Seq data that incorporates miRNA-specific solutions at each step of the analysis. The pipeline was designed to be flexible and user-friendly, allowing researchers with different levels of expertise to perform the analysis in a standardized and reproducible manner, using the most common and widely used tools for each step. In this work, we describe the implementation of myBrain-Seq and demonstrate its capacity to consistently and reproducibly identify differentially expressed miRNAs and enriched pathways by applying it to a real case study in which we compared schizophrenia patients who responded to medication with treatment-resistant schizophrenia patients to obtain a 16-miRNA treatment-resistant schizophrenia profile.

Transcriptome analysis reveals the host immune response upon LMBV infection in largemouth bass (Micropterus salmoides)

Largemouth bass (Micropterus salmoides) is one of the important economical freshwater aquaculture species in China. However, the outbreak of viral diseases always caused great economic losses in the largemouth bass aquaculture industry. Largemouth bass virus (LMBV), a double-stranded DNA (dsDNA) virus belonging to genus Ranavirus, family Iridoviridae causes high mortality in cultivated largemouth bass. However, host responses, especially the molecular events involved in LMBV infection still remained largely uncertain. Here, we established an in vivo model of LMBV infection, and systematically investigated the mRNA expression profiles of host genes in liver and spleen from infected largemouth bass using RNA sequencing (RNA-seq). Histopathological analysis indicated that necrotic cells and the formed necrotic focus were present in spleen, while numerous basophilic cells, hepatocytes volume shrinkage, nucleus pyknosis, and the disappeared boundary of hepatocytes were observed in the liver of infected largemouth bass. Transcriptomic analysis showed that transcription levels of 5128 genes (2804 up-regulated genes and 2324 down-regulated) in liver and 7008 genes (2603 up-regulated and 4405 down-regulated) in spleen were altered significantly. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that numerous co-regulated differentially expressed genes (DEGs) in liver and spleen were enriched in the pathways related to cell death and immune signaling, such as apoptosis, necroptosis, cytokine-cytokine receptor interaction and JAK-STAT signaling. Moreover, the DEGs specially regulated by LMBV infection in liver were significantly enriched in the KEGG pathways related to metabolism and cell death, while those in spleen were enriched in the immune related pathways. In addition, the expression changes of several randomly selected genes, such as SOCS1, IL-6, CXCL2, CASP8, CYC and TNF from qPCR were consistent with the transcriptomic data. Taken together, our findings will provide new insights into the fundamental patterns of molecular responses induced by LMBV in vivo, but also contribute greatly to understanding the host defense mechanisms against iridoviral pathogens.

Single-Cell RNA Sequencing Reveals Microevolution of the Stickleback Immune System

The risk and severity of pathogen infections in humans, livestock, or wild organisms depend on host immune function, which can vary between closely related host populations or even among individuals. This immune variation can entail between-population differences in immune gene coding sequences, copy number, or expression. In recent years, many studies have focused on population divergence in immunity using whole-tissue transcriptomics. But, whole-tissue transcriptomics cannot distinguish between evolved differences in gene regulation within cells, versus changes in cell composition within the focal tissue. Here, we leverage single-cell transcriptomic approaches to document signatures of microevolution of immune system structure in a natural system, the three-spined stickleback (Gasterosteus aculeatus). We sampled nine adult fish from three populations with variability in resistance to a cestode parasite, Schistocephalus solidus, to create the first comprehensive immune cell atlas for G. aculeatus. Eight broad immune cell types, corresponding to major vertebrate immune cells, were identified. We were also able to document significant variation in both abundance and expression profiles of the individual immune cell types among the three populations of fish. Furthermore, we demonstrate that identified cell type markers can be used to reinterpret traditional transcriptomic data: we reevaluate previously published whole-tissue transcriptome data from a quantitative genetic experimental infection study to gain better resolution relating infection outcomes to inferred cell type variation. Our combined study demonstrates the power of single-cell sequencing to not only document evolutionary phenomena (i.e., microevolution of immune cells) but also increase the power of traditional transcriptomic data sets.

Dual RNA-seq of spleens extracted from channel catfish infected with Aeromonas veronii reveals novel insights into host-pathogen interactions

Interactions between host and pathogen are involving various dynamic changes in transcript expression and critical for understanding host immunity against infections and its associated pathogenesis. Herein, we established a model of channel catfish infected with Aeromonas veronii. The infected fish had prominent body surface bleeding, and the spleen showed hyperemia and swelling. Then, the spleen of channel catfish infected with A. veronii was analyzed by dual RNA sequencing (RNA-seq), and the transcriptome data were compared with uninfected channel catfish spleen or bacteria cultured in vitro. The transcript expression profile of pathogen-host interaction between A. veronii and channel catfish was successfully studied. During infection, the host was enriched for multiple immune-related signaling pathways, such as the Toll-like receptor signaling pathway, Cytokine-cytokine receptor interaction, and T cell receptor signaling pathway; and significantly upregulated for many innate immune-related genes, including IL-8. At the same time, we found that A. veronii mainly harmed the host spleen through hemolysin. Our current findings are of great significance in clarifying the pathogenesis of channel catfish induced by A. veronii and provide gene targets for developing preventive measures.

Transcriptome analyses reveal differentially expressed genes associated with development of the palatal organ in bighead carp (Hypophthalmichthys nobilis)

The palatal organ is a filter-feeding related organ and occupies a considerable proportion of the head of bighead carp (Hypophthalmichthys nobilis), a large cyprinid fish intensive aquaculture in Asia. In this study, we performed RNA-seq of the palatal organ during growth periods of two (M2), six (M6) and 15 (M15) months of age after hatching. The numbers of differentially expressed genes (DEGs) were 1384, 481 and 1837 for M2 VS M6, M6 VS M15 and M2 VS M15 respectively. The following signaling pathways of energy metabolism and cytoskeleton function were enriched, including ECM-receptor interaction, Cardiac muscle contraction, Steroid biosynthesis and PPAR signaling pathway. Several members of collagen family (col1a1, col2a1, col6a2, col6a3, col9a2), Laminin gamma 1 (lamc1), integrin alpha 1 (itga1), Fatty acid binding protein 2 (fads2) and lipoprotein lipase (lpl), and Protein tyrosine kinase 7 (Ptk7) are candidate genes for growth and development of basic tissues of the palatal organ. Furthermore, taste-related genes such as fgfrl1, fgf8a, fsta and notch1a were also identified, which may be involved in the development of taste buds of the palatal organ. The transcriptome data obtained in this study provide insights into the understanding functions and development mechanisms of palatal organ, and potential candidate genes that may be related to the genetic modulation of head size of bighead carp.

Single-cell RNA sequencing reveals homogeneous transcriptome patterns and low variance in a suspension CHO-K1 and an adherent HEK293FT cell line in culture conditions

Recombinant mammalian host cell lines, in particular CHO and HEK293 cells, are used for the industrial production of therapeutic proteins. Despite their well-known genomic instability, the control mechanisms that enable cells to respond to changes in the environmental conditions are not yet fully understood, nor do we have a good understanding of the factors that lead to phenotypic shifts in long-term cultures. A contributing factor could be inherent diversity in transcriptomes within a population. In this study, we used a full-length coverage single-cell RNA sequencing (scRNA-seq) approach to investigate and compare cell-to-cell variability and the impact of standardized and homogenous culture conditions on the diversity of individual cell transcriptomes, comparing suspension CHO-K1 and adherent HEK293FT cells. Our data showed a critical batch effect from the sequencing of four 96-well plates of CHO-K1 single cells stored for different periods of time, which was and may be therefore identified as a technical variable to consider in experimental planning. Besides, in an artificial and controlled culture environment such as used in routine cell culture technology, the gene expression pattern of a given population does not reveal any marker gene capable to disclose relevant cell population substructures, both for CHO-K1 cells and for HEK293FT cells. The variation observed is primarily driven by the cell cycle.

Long non‐coding RNA, a supreme post‐transcriptional immune regulator of bacterial or virus‐driven immune evolution in teleost

The global aquaculture boom, fuelled by a reduction in wild population and detection of novel viruses, has created a demanding market, hence, there is a pressing need to investigate the immune system of fish, further. As the most diverse community of vertebrates and a central contributor to the progressing global aquaculture market, teleost continues to draw vast scientific interest. Recent breakthroughs in multi‐omics technologies have provided a platform to understand the role of long non‐coding RNA (lncRNA) in the host immune system during infection. Emerging evidence shows that teleost lncRNA might have a regulatory role in immune responses, mostly through lncRNA–microRNA (miRNA) sponging. Teleost lncRNA shares a functionally active short sequence complement to target the miRNA which is conserved among the several fish species. Recent report suggests that rhabdovirus exploits a lncRNA in teleost and, to dodge the host immune mechanism and negatively regulate the immune system. This observation reveals the essentiality of lncRNA in pathogen‐driven immunity in teleost. Reports available on the function of teleost lncRNA are still in early stages and experimental verifications are a limiting factor. Unravelling the lncRNA‐mediated immune regulation in fishes could be used against the invading pathogens to strengthen the aquaculture production. This review elaborates on the experimentally identified and functionally characterized lncRNA and its regulatory role in the teleost immune response during infection and pathogen‐driven host immune evolution, which could eventually lead to achieving high standards in aquaculture productivity.

Blood-based gene expression as non-lethal tool for inferring salinity-habitat history of European eel (Anguilla anguilla)

The European eel is a facultative catadromous species, meaning that it can skip the freshwater phase or move between marine and freshwater habitats during its continental life stage. Otolith microchemistry, used to determine the habitat use of eel or its salinity history, requires the sacrifice of animals. In this context, blood-based gene expression may represent a non-lethal alternative. In this work, we tested the ability of blood transcriptional profiling to identify the different salinity-habitat histories of European eel. Eels collected from different locations in Norway were classified through otolith microchemistry as freshwater residents (FWR), seawater residents (SWR) or inter-habitat shifters (IHS). We detected 3451 differentially expressed genes from blood by comparing FWR and SWR groups, and then used that subset of genes in a machine learning approach (i.e., random forest) to the extended FWR, SWR, and IHS group. Random forest correctly classified 100% of FWR and SWR and 83% of the IHS using a minimum of 30 genes. The implementation of this non-lethal approach may replace otolith-based microchemistry analysis for the general assessment of life-history tactics in European eels. Overall, this approach is promising for the replacement or reduction of other lethal analyses in determining certain fish traits.

Gene expression patterns in Atlantic salmon (Salmo salar) with severe nephrocalcinosis

Nephrocalcinosis is a common disorder in farmed Atlantic salmon, but the consequences for the fish physiology are not well understood. We performed a transcriptome study in kidneys of Atlantic salmon (Salmo salar) smolts without and with severe chronic nephrocalcinosis (NC). The study revealed that numerous genes are differentially expressed in fish with NC compared with healthy salmon. The most evident changes in gene expression patterns in the NC group were a massive downregulation of metabolism and energy production, upregulation of signalling pathways important for tissue repair and function maintenance and upregulation of inflammatory responses. Overall, the extensive tissue damage and the gene regulation responses that affect salmon with severe nephrocalcinosis are highly likely to have dramatic consequences on fish survival.

Transcriptome characterization and SSR discovery in the giant spiny frog Quasipaa spinosa

The giant spiny frog Quasipaa spinosa (Amphibia: Ranidae) is a large unique frog species found mainly in southern China with a low amount of fat and high protein, and it has become one of the most important aquaculture animal species in China. To better understand its genetic background and screen potential molecular markers for artificial breeding and species conservation, we constructed an expression profile of Q. spinosa with high-throughput RNA sequencing and acquired potential SSR markers. Approximately 81.7 Gb of data and 93,887 unigenes were generated. The transcriptome contains 2085 (80.7 %) complete BUSCOs, suggesting that our assembly methods were effective and accurate.These unigenes were functionally classified using 7 functional databases, yielding 17,482 Pfam-, 12,752 Sting-, 17,526 KEGG-, 24,341 Swiss-Prot-, 28,604 Nr-, 16,287 GO- and 12,752 COG-annotated unigenes. Among several amphibian species, Q. spinosa unigenes had the highest number of hits to Xenopus tropicalis (35.25 %), followed by Xenopus laevis (12.68 %). 1417 unigenes were assigned to the immune system. In addition, a total of 33,019 candidate SSR markers were identified from the constructed library. Further tests with 20 loci and 118 large-scale breeding specimens gathered from four culture farms in China showed that 15 (75 %) loci were polymorphic, with the number of alleles per locus varying from 3 to 9 (mean of 4.3). The PIC values for the SSR markers ranged from 0.19 to 0.82, with an average value of 0.43, indicating moderate polymorphism in Q. spinosa. The transcriptomic profile and SSR repertoire obtained in the present study will facilitate population genetic studies and the selective breeding of amphibian species.

Transcriptomic Analysis in Marine Medaka Gill Reveals That the Hypo-Osmotic Stress Could Alter the Immune Response via the IL17 Signaling Pathway

Fish gills are the major osmoregulatory tissue that contact the external water environment and have developed an effective osmoregulatory mechanism to maintain cellular function. Marine medaka (Oryzias melastigma) has the ability to live in both seawater and fresh water environments. The present study performed a seawater (SW) to 50% seawater (SFW) transfer, and the gill samples were used for comparative transcriptomic analysis to study the alteration of hypo-osmotic stress on immune responsive genes in this model organism. The result identified 518 differentiated expressed genes (DEGs) after the SW to SFW transfer. Various pathways such as p53 signaling, forkhead box O signaling, and the cell cycle were enriched. Moreover, the immune system was highlighted as one of the top altered biological processes in the enrichment analysis. Various cytokines, chemokines, and inflammatory genes that participate in the IL-17 signaling pathway were suppressed after the SW to SFW transfer. On the other hand, some immunoglobulin-related genes were up-regulated. The results were further validated by real-time qPCR. Taken together, our study provides additional gill transcriptome information in marine medaka; it also supports the notion that osmotic stress could influence the immune responses in fish gills.

Comparative transcriptomic analysis reveals different host cell responses to Singapore grouper iridovirus and red-spotted grouper nervous necrosis virus

Singapore grouper iridovirus (SGIV) and red-spotted grouper nervous necrosis virus (RGNNV) are important pathogens that cause high mortality and heavy economic losses in grouper aquaculture. Interestingly, SGIV infection in grouper cells induces paraptosis-like cell death, while RGNNV infection induces autophagy and necrosis characterized morphologically by vacuolation of lysosome. Here, a comparative transcriptomic analysis was carried out to identify the different molecular events during SGIV and RGNNV infection in grouper spleen (EAGS) cells. The functional enrichment analysis of DEGs suggested that several signaling pathways were involved in CPE progression and host immune response against SGIV or RGNNV. Most of DEGs featured in the KEGG "lysosome pathway" were up-regulated in RGNNV-infected cells, indicating that RGNNV induced lysosomal vacuolization and autophagy might be due to the disturbance of lysosomal function. More than 100 DEGs in cytoskeleton pathway and mitogen-activated protein kinase (MAPK) signal pathway were identified during SGIV infection, providing additional evidence for the roles of cytoskeleton remodeling in cell rounding during CPE progression and MAPK signaling in SGIV induced cell death. Of note, consistent with changes at the transcriptional levels, the post-translational modifications of MAPK signaling-related proteins were also detected during RGNNV infection, and the inhibitors of extracellular signal-regulated kinase (ERK) and p38 MAPK significantly suppressed viral replication and virus induced vacuoles formation. Moreover, the majority of DEGs in interferon and inflammation signaling were obviously up-regulated during RGNNV infection, but down-regulated during SGIV infection, suggesting that SGIV and RGNNV differently manipulated host immune response in vitro. In addition, purine and pyrimidine metabolism pathways were also differently regulated in SGIV and RGNNV-infection cells. Taken together, our data will provide new insights into understanding the potential mechanisms underlying different host cell responses against fish DNA and RNA virus.

Potential biomarkers of endocrine and habitat disruption identified via RNA-Seq in Salvelinus fontinalis with proximity to fracking operations in Pennsylvania headwater stream ecosystems

Unconventional natural gas development (fracking) has been a rapidly expanding technique used for the extraction of natural gas from the Marcellus Shale formation in Pennsylvania. There remains a knowledge gap regarding the ecological impacts of fracking, especially regarding the long-term health of native Brook trout (Salvelinus fontinalis) populations. During the summer of 2015, Brook trout were sampled from twelve streams located in forested, northwestern Pennsylvania in order to evaluate the impacts of fracking on Brook trout. Four stream sites were undisturbed (no fracking activity), three had a developed well pad without fracking activity, and five had active fracking with natural gas production. Liver tissue was isolated from two to five fish per stream and underwent RNA-Seq analysis to identify differentially expressed genes between ecosystems with differing fracking status. Data were analyzed individually and with samples pooled within-stream to account for hierarchical data structure and variation in sample coverage within streams. Differentially expressed and differentially alternatively spliced genes had functions related to lipid and steroid metabolism, mRNA processing, RNA polymerase and protein regulation. Unique to our study, genes related to xenobiotic and stress responses were found as well as potential markers for endocrine disruption and saline adaptation that were identified in watersheds with active fracking activity. These results support the utility of RNA-Seq to assess trout health and suggest detrimental impacts of fracking on sensitive trout populations.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Genome-wide alternative splicing profile in the posterior kidney of brown trout (Salmo trutta) during proliferative kidney disease

BACKGROUND

The cnidarian myxozoan parasite Tetracapsuloides bryosalmonae causes chronic proliferative kidney disease (PKD) in salmonids. This parasite is a serious threat to wild and cultured salmonids. T. bryosalmonae undergoes intra-luminal sporogonic development in the kidney of brown trout (Salmo trutta) and the viable spores are released via urine. We investigated the alternative splicing pattern in the posterior kidney of brown trout during PKD.

RESULTS

RNA-seq data were generated from the posterior kidney of brown trout collected at 12 weeks post-exposure to T. bryosalmonae. Subsequently, this data was mapped to the brown trout genome. About 153 significant differently expressed alternatively spliced (DEAS) genes, (delta PSI = 5%, FDR P-value < 0.05) were identified from 19,722 alternatively spliced events. Among the DEAS genes, the least and most abundant alternative splicing types were alternative 5' splice site (5.23%) and exon skipping (70.59%), respectively. The DEAS genes were significantly enriched for sodium-potassium transporter activity and ion homeostasis (ahcyl1, atp1a3a, atp1a1a.1, and atp1a1a.5). The protein-protein interaction network analysis enriched two local network clusters namely cation transporting ATPase C-terminus and Sodium/potassium ATPase beta chain cluster, and mixed inclusion of Ion homeostasis and EF-hand domain cluster. Furthermore, the human disease-related salmonella infection pathway was significantly enriched in the protein-protein interaction network.

CONCLUSION

This study provides the first baseline information about alternative splicing in brown trout during PKD. The generated data lay a foundation for further functional molecular studies in PKD - brown trout infection model. The information generated from the present study can help to develop therapeutic strategies for PKD in the future.

Transcriptome analysis provides insight into deltamethrin-induced fat accumulation in 3T3-L1 adipocytes

Previously, deltamethrin (a Type-II pyrethroid) has been reported to increase triglyceride (fat) accumulation in adipocytes, while its underlying molecular mechanism is not fully determined. The aim of this study was to further investigate the molecular mechanisms of deltamethrin induced fat accumulation in murine 3T3-L1 adipocytes. Consistent to previous reports, deltamethrin (10 μM) significantly promoted adipogenesis in 3T3-L1 adipocytes. RNA sequencing (RNA-seq) results showed that 721 differentially expressed genes (DEGs) were identified after deltamethrin treatment, involving in 58 Functional groups of Gene Ontology (GO) and 255 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Several key functional groups regulating adipogenesis, such as fat cell differentiation (Igf1, Snai2, Fgf10, and Enpp1) and cytosolic calcium ion concentration (Nos1, Cxcl1, and Ngf) were significantly enriched. Collectively, these results suggest that the promotion of adipogenesis by deltamethrin was attributed to an obesogenic global transcriptomic response, which provides further understanding of the underlying mechanisms of deltamethrin-induced fat accumulation.

Molecular, biochemical and behavioral responses of Daphnia magna under long-term exposure to polystyrene nanoplastics

The presence and potential toxicity of nanoplastics (NPs) in aquatic ecosystems is an issue of growing concern. Although many studies have investigated the adverse effects of short-term exposure to high concentrations of NPs to aquatic organisms, the information on the consequences caused by the administration of low NPs concentrations over long-term exposure is limited. The present study aimed at investigating the effects induced by a long-term exposure (21-days) to two sub-lethal concentrations of polystyrene nanoplastics (PS-NPs; 0.05 and 0.5 µg/mL) on Daphnia magna. A multi-level approach was performed to assess potential sub-individual (i.e., molecular and biochemical) and individual (i.e., behavioural) adverse effects. At molecular level, the modulation of the expression of genes involved in antioxidant defence, response to stressful conditions and specific physiological pathways was investigated. Oxidative stress (i.e., the amount of pro-oxidants, the activity of antioxidant and detoxifying enzymes and lipid peroxidation) and energetic (i.e., protein, carbohydrate, lipid and total caloric content) biomarkers were applied to assess effects at the biochemical level, while swimming activity was measured to monitor changes in individual behavior. Although the 21-days exposure to PS-NPs induced a slight modulation of gene involved in oxidative stress response, biochemical analyses showed that D. magna individuals did not experience an oxidative stress condition. Significant changes in energy reserves of individuals exposed for 21 days to both the PS-NPs concentrations were observed, but no alterations of swimming activity occurred. Our results highlighted that the exposure to low concentrations of PS-NPs could pose a limited risk to D. magna individuals and suggested the importance of a multi-level approach to assess the risks of NPs on aquatic organisms.

Transcriptome analysis in the spleen of Northern Snakehead (Channa argus) challenged with Nocardia seriolae

Northern snakehead (Channa argus) is an indigenous fish species and is one of popularly cultured snakeheads in China and other Asian countries. Unfortunately, Nocardia seriolae infections have caused considerable losses in the snakehead aquaculture industry. However, the infectivity and the immune response induced by N. seriolae in snakehead are unclear. In order to better understand the immune response of Northern snakehead in a series of time points after N. seriolae challenge, we conducted the transcriptomic comparison in snakehead spleen at 48, 96, and 144 h after the challenge of N. seriola against their control counterparts. Gene annotation and pathway analysis of differentially expressed genes (DEGs) were carried out to understand the functions of the DEGs. Additionally, protein-protein interaction networks were conducted to obtain the interaction relationships of immune-related DEGs. These results revealed the expression changes of multiple DEGs and signaling pathways involved in immunity during N. seriolae infection, which will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection in the northern snakehead.

Time-Course Transcriptome Analysis of Lungs From Mice Infected With Hypervirulent Klebsiella pneumoniae via Aerosolized Intratracheal Inoculation

Hypervirulent Klebsiella pneumoniae (hvKp) can cause life-threatening community-acquired infections among healthy young individuals and is thus of concern for global dissemination. In this study, a mouse model of acute primary hvKp pneumonia was established via aerosolized intratracheal (i.t.) inoculation, laying the foundation for conducting extensive studies related to hvKp. Subsequently, a time-course transcriptional profile was created of the lungs from the mouse model at 0, 12, 24, 48 and 60 hours post-infection (hpi) using RNA Sequencing (RNA-Seq). RNA-Seq data were analyzed with the use of Mfuzz time clustering, weighted gene co-expression network analysis (WGCNA) and Immune Cell Abundance Identifier for mouse (ImmuCellAI-mouse). A gradual change in the transcriptional profile of the lungs was observed that reflected expected disease progression. At 12 hpi, genes related to acute phase inflammatory response increased in expression and lipid metabolism appeared to have a pro-inflammatory effect. At 24 hpi, exacerbation of inflammation was observed and active IFN-γ suggested that signaling promoted activation and recruitment of macrophages occurred. Genes related to maintaining the structural integrity of lung tissues showed a sustained decrease in expression after infection and the decrease was especially marked at 48 hpi. TNF, IL-17, MAPK and NF-kB signaling pathways may play key roles in the immunopathogenesis mechanism at all stages of infection. Natural killer (NK) cells consistently decreased in abundance after infection, which has rarely been reported in hvKp infection and could provide a new target for treatment. Genes Saa1 and Slpi were significantly upregulated during infection. Both Saa1, which is associated with lipopolysaccharide (LPS) that elicits host inflammatory response, and Slpi, which encodes an antimicrobial protein, have not previously been reported in hvKp infections and could be important targets for subsequent studies. To t our knowledge, this paper represents the first study to investigate the pulmonary transcriptional response to hvKp infection. The results provide new insights into the molecular mechanisms underlying the pathogenesis of hvKp pulmonary infection that can contribute to the development of therapies to reduce hvKp pneumonia.

Global gene expression responses of Atlantic salmon skin to Moritella viscosa

Moritella viscosa is a Gram-negative pathogen that causes large, chronic ulcers, known as winter-ulcer disease, in the skin of several fish species including Atlantic salmon. We used a bath challenge approach to profile the transcriptome responses of M. viscosa-infected Atlantic salmon skin at the lesion (Mv-At) and away from the lesion (Mv-Aw) sites. M. viscosa infection was confirmed through RNA-based qPCR assays. RNA-Seq identified 5212 and 2911 transcripts differentially expressed in the Mv-At compared to no-infection control and Mv-Aw groups, respectively. Also, there were 563 differentially expressed transcripts when comparing the Mv-Aw to control samples. Our results suggest that M. viscosa caused massive and strong, but largely infection site-focused, transcriptome dysregulations in Atlantic salmon skin, and its effects beyond the skin lesion site were comparably subtle. The M. viscosa-induced transcripts of Atlantic salmon were mainly involved in innate and adaptive immune response-related pathways, whereas the suppressed transcripts by this pathogen were largely connected to developmental and cellular processes. As validated by qPCR, M. viscosa dysregulated transcripts encoding receptors, signal transducers, transcription factors and immune effectors playing roles in TLR- and IFN-dependent pathways as well as immunoregulation, antigen presentation and T-cell development. This study broadened the current understanding of molecular pathways underlying M. viscosa-triggered responses of Atlantic salmon, and identified biomarkers that may assist to diagnose and combat this pathogen.

Transcriptomic Analysis of Fish Hosts Responses to Nervous Necrosis Virus

Nervous necrosis virus (NNV) has been responsible for mass mortalities in the aquaculture industry worldwide, with great economic and environmental impact. The present review aims to summarize the current knowledge of gene expression responses to nervous necrosis virus infection in different fish species based on transcriptomic analysis data. Four electronic databases, including PubMed, Web of Science, and SCOPUS were searched, and more than 500 publications on the subject were identified. Following the application of the appropriate testing, a total of 24 articles proved eligible for this review. NNV infection of different host species, in different developmental stages and tissues, presented in the eligible publications, are described in detail, revealing and highlighting genes and pathways that are most affected by the viral infection. Those transcriptome studies of NNV infected fish are oriented in elucidating the roles of genes/biomarkers for functions of special interest, depending on each study’s specific emphasis. This review presents a first attempt to provide an overview of universal host reaction mechanisms to viral infections, which will provide us with new perspectives to overcome NNV infection to build healthier and sustainable aquaculture systems.

Genome-Wide Expression and Physiological Profiling of Pearl Millet Genotype Reveal the Biological Pathways and Various Gene Clusters Underlying Salt Resistance

Pearl millet (Pennisetum glaucum L.) is a vital staple food and an important cereal crop used as food, feed, and forage. It can withstand heat and drought due to the presence of some unique genes; however, the mechanism of salt stress has been missing in pearl millet until now. Therefore, we conducted a comparative transcriptome profiling to reveal the differentially expressed transcripts (DETs) associated with salt stress in pearl millet at different time points, such as 1, 3, and 7 h, of salt treatment. The physiological results suggested that salt stress significantly increased proline, malondialdehyde (MDA) content, and hydrogen peroxide (H₂O₂) in pearl millet at 1, 3, and 7 h of salt treatment. In addition, pearl millet plants regulated the activities of superoxide dismutase, catalase, and peroxidase to lessen the impact of salinity. The transcriptomic results depicted that salt stress upregulated and downregulated the expression of various transcripts involved in different metabolic functions. At 1 and 7 h of salt treatment, most of the transcripts were highly upregulated as compared to the 3 h treatment. Moreover, among commonly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, the mitogen-activated protein kinase (MAPK) signaling pathway and peroxisome pathway were significantly enriched. The DETs related to hormone signaling (auxins, ethylene, gibberellin, and abscisic acid), kinases, protein modifications, and degradation were also identified, depicting the possible role of hormones and kinases to enhance plant tolerance against salt stress. Furthermore, the transcription factors, such as ethylene-responsive element binding factors (ERF), basic helix-loop-helix (bHLH), HMG box-containing protein (HBP), MADS, myeloblastosis (MYB), and WRKY, were predicted to significantly regulate different transcripts involved in salt stress responses at three different time points. Overall, this study will provide new insights to better understand the salt stress regulation mechanisms in pearl millet to improve its resistance against salinity and to identify new transcripts that control these mechanisms in other cereals.

Meta-Analysis of Esophageal Cancer Transcriptomes Using Independent Component Analysis

Independent Component Analysis is a matrix factorization method for data dimension reduction. ICA has been widely applied for the analysis of transcriptomic data for blind separation of biological, environmental, and technical factors affecting gene expression. The study aimed to analyze the publicly available esophageal cancer data using the ICA for identification and comprehensive analysis of reproducible signaling pathways and molecular signatures involved in this cancer type. In this study, four independent esophageal cancer transcriptomic datasets from GEO databases were used. A bioinformatics tool « BiODICA-Independent Component Analysis of Big Omics Data» was applied to compute independent components (ICs). Gene Set Enrichment Analysis (GSEA) and ToppGene uncovered the most significantly enriched pathways. Construction and visualization of gene networks and graphs were performed using the Cytoscape, and HPRD database. The correlation graph between decompositions into 30 ICs was built with absolute correlation values exceeding 0.3. Clusters of components-pseudocliques were observed in the structure of the correlation graph. The top 1,000 most contributing genes of each ICs in the pseudocliques were mapped to the PPI network to construct associated signaling pathways. Some cliques were composed of densely interconnected nodes and included components common to most cancer types (such as cell cycle and extracellular matrix signals), while others were specific to EC. The results of this investigation may reveal potential biomarkers of esophageal carcinogenesis, functional subsystems dysregulated in the tumor cells, and be helpful in predicting the early development of a tumor.

Full-Length Transcriptome: A Reliable Alternative for Single-Cell RNA-Seq Analysis in the Spleen of Teleost Without Reference Genome

Fish is considered as a supreme model for clarifying the evolution and regulatory mechanism of vertebrate immunity. However, the knowledge of distinct immune cell populations in fish is still limited, and further development of techniques advancing the identification of fish immune cell populations and their functions are required. Single cell RNA-seq (scRNA-seq) has provided a new approach for effective in-depth identification and characterization of cell subpopulations. Current approaches for scRNA-seq data analysis usually rely on comparison with a reference genome and hence are not suited for samples without any reference genome, which is currently very common in fish research. Here, we present an alternative, i.e. scRNA-seq data analysis with a full-length transcriptome as a reference, and evaluate this approach on samples from Epinephelus coioides-a teleost without any published genome. We show that it reconstructs well most of the present transcripts in the scRNA-seq data achieving a sensitivity equivalent to approaches relying on genome alignments of related species. Based on cell heterogeneity and known markers, we characterized four cell types: T cells, B cells, monocytes/macrophages (Mo/MΦ) and NCC (non-specific cytotoxic cells). Further analysis indicated the presence of two subsets of Mo/MΦ including M1 and M2 type, as well as four subsets in B cells, i.e. mature B cells, immature B cells, pre B cells and early-pre B cells. Our research will provide new clues for understanding biological characteristics, development and function of immune cell populations of teleost. Furthermore, our approach provides a reliable alternative for scRNA-seq data analysis in teleost for which no reference genome is currently available.

Comparative Transcriptome Analysis of the Expression of Antioxidant and Immunity Genes in the Spleen of a Cyanidin 3-O-Glucoside-Treated Alzheimer's Mouse Model

Cyanidin 3-O-glucoside (C3G) is a well-known antioxidant found as a dietary anthocyanin in different fruits and vegetables. It has protective and therapeutic effects on various diseases. It can reduce neuronal death from amyloid-beta (Aβ)-induced toxicity and promote the inhibition of Aβ fibrillization. Antioxidant and immune modulation might play a critical role in the properties of C3G against Alzheimer's disease (AD) and other diseases. However, limited studies have been performed on the mechanism involved in the effect of C3G through transcriptome analysis. Thus, the objective of this study was to perform comparative transcriptome analysis of the spleen to determine gene expression profiles of wild-type mice (C57BL/6J Jms), an Alzheimer's mouse model (APPswe/PS1dE9 mice), and a C3G-treated Alzheimer's mouse model. Differentially expressed antioxidant, immune-related, and AD pathways genes were identified in the treated group. The validation of gene expression data via RT-PCR studies further supported the current findings. Six important antioxidant genes (S100a8, S100a9, Prdx2, Hp, Mpst, and Prxl2a) and a high number of immune-related genes were found to be upregulated in the treatment groups, suggesting the possible antioxidant and immunomodulatory mechanisms of C3G, respectively. Further studies are strongly recommended to elucidate the precise role of these essential genes and optimize the therapeutic function of C3G in AD and other disease conditions.

RNA-Seq-Based Analysis of Cortisol-Induced Differential Gene Expression Associated with Piscirickettsia salmonis Infection in Rainbow Trout (Oncorhynchus mykiss) Myotubes

Salmonid rickettsial septicemia (SRS) is the major infectious disease of the Chilean salmonid aquaculture industry caused by Piscirickettsia salmonis. Intensive farming conditions generate stress and increased susceptibility to diseases, being skeletal muscle mainly affected. However, the interplay between pathogen infection and stress in muscle is poorly understood. In this study, we perform an RNA-seq analysis on rainbow trout myotubes that are pretreated for 3 h with cortisol (100 ng/mL) and then infected with P. salmonis strain LF-89 for 8 h (MOI 50). Twelve libraries are constructed from RNA samples (n = 3 per group) and sequenced on Illumina HiSeq 4000. A total of 704,979,454 high-quality reads are obtained, with 70.25% mapped against the reference genome. In silico DETs include 175 total genes-124 are upregulated and 51 are downregulated. GO enrichment analysis reveals highly impacted biological processes related to apoptosis, negative regulation of cell proliferation, and innate immune response. These results are validated by RT-qPCR of nine candidate transcripts. Furthermore, cortisol pretreatment significantly stimulated bacterial gene expression of ahpC and 23s compared to infection. In conclusion, for the first time, we describe a transcriptomic response of trout myotubes infected with P. salmonis by inducing apoptosis, downregulating cell proliferation, and intrinsic immune-like response that is differentially regulated by cortisol.

Genome Analysis of Phytophthora nicotianae JM01 Provides Insights into Its Pathogenicity Mechanisms

Phytophthora nicotianae is a widely distributed plant pathogen that can cause serious disease and cause significant economic losses to various crops, including tomatoes, tobacco, onions, and strawberries. To understand its pathogenic mechanisms and explore strategies for controlling diseases caused by this pathogen, we sequenced and analyzed the whole genome of Ph. nicotianae JM01. The Ph. nicotianae JM01 genome was assembled using a combination of approaches including shotgun sequencing, single-molecule sequencing, and the Hi-C technique. The assembled Ph. nicotianae JM01 genome is about 95.32 Mb, with contig and scaffold N50 54.23 kb and 113.15 kb, respectively. The average GC content of the whole-genome is about 49.02%, encoding 23,275 genes. In addition, we identified 19.15% of interspersed elements and 0.95% of tandem elements in the whole genome. A genome-wide phylogenetic tree indicated that Phytophthora diverged from Pythium approximately 156.32 Ma. Meanwhile, we found that 252 and 285 gene families showed expansion and contraction in Phytophthora when compared to gene families in Pythium. To determine the pathogenic mechanisms Ph. nicotianae JM01, we analyzed a suite of proteins involved in plant-pathogen interactions. The results revealed that gene duplication contributed to the expansion of Cell Wall Degrading Enzymes (CWDEs) such as glycoside hydrolases, and effectors such as Arg-Xaa-Leu-Arg (RXLR) effectors. In addition, transient expression was performed on Nicotiana benthamiana by infiltrating with Agrobacterium tumefaciens cells containing a cysteine-rich (SCR) protein. The results indicated that SCR can cause symptoms of hypersensitive response. Moreover, we also conducted comparative genome analysis among four Ph. nicotianae genomes. The completion of the Ph. nicotianae JM01 genome can not only help us understand its genomic characteristics, but also help us discover genes involved in infection and then help us understand its pathogenic mechanisms.

Dual-Organ Transcriptomic Analysis of Rainbow Trout Infected With Ichthyophthirius multifiliis Through Co-Expression and Machine Learning

Ichthyophthirius multifiliis is a major pathogen that causes a high mortality rate in trout farms. However, systemic responses to the pathogen and its interactions with multiple organs during the course of infection have not been well described. In this study, dual-organ transcriptomic responses in the liver and head kidney and hemato-serological indexes were profiled under I. multifiliis infection and recovery to investigate systemic immuno-physiological characteristics. Several strategies for massive transcriptomic interpretation, such as differentially expressed genes (DEGs), Poisson linear discriminant (PLDA), and weighted gene co-expression network analysis (WGCNA) models were used to investigate the featured genes/pathways while minimizing the disadvantages of individual methods. During the course of infection, 6,097 and 2,931 DEGs were identified in the head kidney and liver, respectively. Markers of protein processing in the endoplasmic reticulum, oxidative phosphorylation, and the proteasome were highly expressed. Likewise, simultaneous ferroptosis and cellular reconstruction was observed, which is strongly linked to multiple organ dysfunction. In contrast, pathways relevant to cellular replication were up-regulated in only the head kidney, while endocytosis- and phagosome-related pathways were notably expressed in the liver. Moreover, interestingly, most immune-relevant pathways (e.g., leukocyte trans-endothelial migration, Fc gamma R-mediated phagocytosis) were highly activated in the liver, but the same pathways in the head kidney were down-regulated. These conflicting results from different organs suggest that interpretation of co-expression among organs is crucial for profiling of systemic responses during infection. The dual-organ transcriptomics approaches presented in this study will greatly contribute to our understanding of multi-organ interactions under I. multifiliis infection from a broader perspective.