Co-Expression Network Analysis of Spleen Transcriptome in Rock Bream (Oplegnathus fasciatus) Naturally Infected with Rock Bream Iridovirus (RBIV)

Red sea bream iridovirus infection downregulates inflammation-related genes in the spleen of rock bream (Oplegnathus fasciatus)

This study investigated the kinetics of red sea bream iridovirus and host gene expression during infection in rock bream (Oplegnathus fasciatus), a species highly sensitive to the virus. After intraperitoneal injection of the viral solution at 104 TCID50/fish, the viral genome copy number in the spleen was 104.7 ± 0.2 and 105.9 ± 0.4 copies/μg DNA at 3 and 5 days post-injection (dpi), respectively. Using transcriptomic analyses via MiSeq, viral gene transcripts were detected at 3 and 5 dpi. Six genes including RING-finger domain-containing protein and laminin-type epidermal growth factor-like domain genes were significantly expressed at 5 dpi. Further, 334 host genes were differentially expressed compared with those before infection. Genes were clustered into four groups based on their expression profiles. Interferon-stimulated genes were more prevalent in groups showing upregulation at 5 dpi and 3 and 5 dpi. In contrast, the group showing downregulation at 3 dpi included inflammation-related genes, such as granzyme and eosinophil peroxidase genes. Downregulation of certain inflammation-related genes may contribute to the susceptibility of this fish to the virus.

Inoculation of Freund's adjuvant in European eel (Anguilla anguilla) revealed key KEGG pathways and DEGs of host anti-Edwardsiella anguillarum infection

Freund's complete (FCA) and incomplete adjuvants (FIA), generally applied in subunit fishery vaccine, have not been explored on the molecular mechanism of the nonspecific immune enhancement. In this study, we examined the RNA-seq in the spleen of European eel (Anguilla anguilla) inoculated with FCA and FIA (FCIA group) to elucidate the key KEGG pathways and differential expressed genes (DEGs) in the process of Edwardsiella anguillarum infection and A. anguilla anti-E. anguillarum infection using genome-wide transcriptome. After eels were challenged by E. anguillarum at 28 d post the first inoculation (dpi), compared to the control uninfected eels (Con group), the control infected eels (Con_inf group) showed severe pathological changes in the liver, kidney and spleen, although infected eels post the inoculation of FCIA (FCIA_inf group) also formed slight bleeding. Compared to the FCIA_inf group, there was more than 10 times colony forming unit (cfu) in the Con_inf group per 100 μg spleen, kidney or blood, and the relative percent survival (RPS) of eels was 44.4% in FCIA_inf vs Con_inf. Compared to the Con group, the SOD activity in the FCIA group increased significantly in the liver and spleen. Using high-throughput transcriptomics, DEGs were identified and 29 genes were verified using fluorescence real-time polymerase chain reaction (qRT-PCR). The result of DEGs clustering showed 9 samples in 3 groups of Con, FCIA and FCIA_inf were similar, contrast to distinct differences of 3 samples in the Con_inf group. We found 3795 up and 3548 down regulated DEGs in the compare of FCIA_inf vs Con_inf, of which 5 enriched KEGG pathways of "Lysosome", "Autophagy", "Apoptosis", "C-type lectin receptor signaling" and "Insulin signaling" were ascertained, and 26 of 30 top GO terms in the compare were significantly enriched. Finally, protein-protein interactions between the DEGs of the 5 KEGG pathways and other DEGs were explored using Cytoscape 3.9.1. The compare of FCIA_inf vs Con_inf showed 110 DEGs from the 5 pathways and 718 DEGs from other pathways formed total of 9747° in a network, of which 9 hub DEGs play vital roles in anti-infection or apoptosis. Together, the interaction networks revealed that 9 DEGs involved in the 5 pathways underlies the key process of A. anguilla anti-E. anguillarum infection or host cell apoptosis.

Integrative time series of cellular, humoral and molecular response revealed immunotoxicity of bifenthrin to Chinese rare minnow (Gobiocypris rarus) following Pseudomonas fluorescens challenge

Bifenthrin is a common pesticide that is widespread in aquatic environments. Although it has been shown to be toxic to aquatic organisms, its immunotoxicity and mechanism are unclear. Herein, we reported the immunotoxicity of bifenthrin on adult Chinese rare minnow (Gobiocypris rarus) after 28 days of exposure to different concentrations of bifenthrin (0.1, 0.3, and 1.0 μg/L) and 36-h Pseudomonas fluorescens challenge. Bifenthrin inhibited the fish humoral immune response to bacteria by altering the lymphocyte and neutrophil ratios and decreasing the production of lysozyme, complement component 3, immunoglobulin M, and C-reactive protein, particularly were 1.0 μg/L. Bifenthrin caused intestinal damage and significantly reduced the volume of intestinal mucus at 12 and 36 hours postinjection (hpi) (p < 0.05). Moreover, 1.0 μg/L bifenthrin significantly increased the fish mortality and bacterial loads at 12 and 36 hpi (p < 0.05). RNA-seq analysis revealed several enriched genes involved in pathogen attachment and recognition, inflammatory responses, and complement system at the early-to-mid stage of infection (4-12 hpi). Overall, our results corroborated that bifenthrin induced immunotoxicity in Gobiocypris rarus, resulting in immune dysfunction of fish and increasing their sensitivity to bacterial infection and accelerating mortality. Moreover, 4-12 hpi was better than 36 hpi for analyzing immune responses against pathogen infection in fish exposed to bifenthrin.

Improved biomarker discovery through a plot twist in transcriptomic data analysis

Background

Transcriptomic analysis is crucial for understanding the functional elements of the genome, with the classic method consisting of screening transcriptomics datasets for differentially expressed genes (DEGs). Additionally, since 2005, weighted gene co-expression network analysis (WGCNA) has emerged as a powerful method to explore relationships between genes. However, an approach combining both methods, i.e., filtering the transcriptome dataset by DEGs or other criteria, followed by WGCNA (DEGs + WGCNA), has become common. This is of concern because such approach can affect the resulting underlying architecture of the network under analysis and lead to wrong conclusions. Here, we explore a plot twist to transcriptome data analysis: applying WGCNA to exploit entire datasets without affecting the topology of the network, followed with the strength and relative simplicity of DEG analysis (WGCNA + DEGs). We tested WGCNA + DEGs against DEGs + WGCNA to publicly available transcriptomics data in one of the most transcriptomically complex tissues and delicate processes: vertebrate gonads undergoing sex differentiation. We further validate the general applicability of our approach through analysis of datasets from three distinct model systems: European sea bass, mouse, and human.

Results

In all cases, WGCNA + DEGs clearly outperformed DEGs + WGCNA. First, the network model fit and node connectivity measures and other network statistics improved. The gene lists filtered by each method were different, the number of modules associated with the trait of interest and key genes retained increased, and GO terms of biological processes provided a more nuanced representation of the biological question under consideration. Lastly, WGCNA + DEGs facilitated biomarker discovery.

Conclusions

We propose that building a co-expression network from an entire dataset, and only thereafter filtering by DEGs, should be the method to use in transcriptomic studies, regardless of biological system, species, or question being considered.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01398-w.

Comparative Transcriptome Analysis of Spleen Reveals Potential Regulation of Genes and Immune Pathways Following Administration of Aeromonas salmonicida subsp. masoucida Vaccine in Atlantic Salmon (Salmo salar)

Aeromonas salmonicida is a global fish pathogen. Aeromonas salmonicida subsp. masoucida (ASM) is classified as atypical A. salmonicida and caused huge losses to salmonid industry in China. Hence, it is of great significance to develop ASM vaccine and explore its protection mechanism in salmonids. In this regard, we conducted RNA-seq analysis with spleen tissue of Atlantic salmon after ASM vaccination to reveal genes, their expression patterns, and pathways involved in immune protections. In our results, a total of 441.63 million clean reads were obtained, and 389.37 million reads were mapped onto the Atlantic salmon reference genome. In addition, 1125, 2126, 1098, 820, and 1351 genes were significantly up-regulated, and 747, 2626, 818, 254, and 908 genes were significantly down-regulated post-ASM vaccination at 12 h, 24 h, 1 month, 2 months, and 3 months, respectively. Subsequent pathway analysis revealed that many differentially expressed genes (DEGs) following ASM vaccination were involved in cytokine-cytokine receptor interaction (TNFRSF11b, IL-17RA, CCR9, and CXCL11), HTLV-I infection (MR1 and HTLV-1), MAPK signaling pathway (MAPK, IL8, and TNF-α-1), PI3K-Akt signaling pathway (PIK3R3, THBS4, and COL2A1), and TNF signaling pathway (PTGS2, TNFRSF21-l, and CXCL10). Finally, the results of qRT-PCR showed a significant correlation with RNA-seq results, suggesting the reliability of RNA-seq for gene expression analysis. This study provided insights into regulation of gene expression and their involved pathways in Atlantic salmon spleen in responses to vaccine, and set the foundation for further study on the vaccine protective mechanism in Atlantic salmon as well as other teleost species.

Genetics Responses to Hypoxia and Reoxygenation Stress in Larimichthys crocea Revealed via Transcriptome Analysis and Weighted Gene Co-Expression Network

Simple Summary

Hypoxia, which occurs frequently in aquaculture, can cause serious harm to all aspects of the growth, reproduction and metabolism of cultured fish. Due to the intolerance of Larimichthys crocea to hypoxia, Larimichthys crocea often floats head or even dies under hypoxic environment. However, the molecular mechanism of hypoxia tolerance in Larimichthys crocea has not been fully described. Therefore, the aim of this study was to explore the hub regulatory genes under hypoxic stress environment by transcriptome analysis of three key tissues (liver, blood and gill) in Larimichthys crocea. We identified a number of important genes that exercise different regulatory functions. Overall, this study will provide important clues to the molecular mechanisms of hypoxia tolerance in Larimichthys crocea.

Abstract

The large yellow croaker (Larimichthys crocea) is an important marine economic fish in China; however, its intolerance to hypoxia causes widespread mortality. To understand the molecular mechanisms underlying hypoxia tolerance in L. crocea, the transcriptome gene expression profiling of three different tissues (blood, gills, and liver) of L. crocea exposed to hypoxia and reoxygenation stress were performed. In parallel, the gene relationships were investigated based on weighted gene co-expression network analysis (WGCNA). Accordingly, the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that several pathways (e.g., energy metabolism, signal transduction, oxygen transport, and osmotic regulation) may be involved in the response of L. crocea to hypoxia and reoxygenation stress. In addition, also, four key modules (darkorange, magenta, saddlebrown, and darkolivegreen) that were highly relevant to the samples were identified by WGCNA. Furthermore, some hub genes within the association module, including RPS16, EDRF1, KCNK5, SNAT2, PFKL, GSK-3β, and PIK3CD, were found. This is the first study to report the co-expression patterns of a gene network after hypoxia stress in marine fish. The results provide new clues for further research on the molecular mechanisms underlying hypoxia tolerance in L. crocea.

Production of a Bacteria-like Particle Vaccine Targeting Rock Bream (Oplegnathus fasciatus) Iridovirus Using Nicotiana benthamiana

Viral diseases are extremely widespread infections that change constantly through mutations. To produce vaccines against viral diseases, transient expression systems are employed, and Nicotiana benthamiana (tobacco) plants are a rapidly expanding platform. In this study, we developed a recombinant protein vaccine targeting the major capsid protein (MCP) of iridovirus fused with the lysine motif (LysM) and coiled-coil domain of coronin 1 (ccCor1) for surface display using Lactococcus lactis. The protein was abundantly produced in N. benthamiana in its N-glycosylated form. Total soluble proteins isolated from infiltrated N. benthamiana leaves were treated sequentially with increasing ammonium sulfate solution, and recombinant MCP mainly precipitated at 40-60%. Additionally, affinity chromatography using Ni-NTA resin was applied for further purification. Native structure analysis using size exclusion chromatography showed that recombinant MCP existed in a large oligomeric form. A minimum OD₆₀₀ value of 0.4 trichloroacetic acid (TCA)-treated L. lactis was required for efficient recombinant MCP display. Immunogenicity of recombinant MCP was assessed in a mouse model through enzyme-linked immunosorbent assay (ELISA) with serum-injected recombinant MCP-displaying L. lactis. In summary, we developed a plant-based recombinant vaccine production system combined with surface display on L. lactis.

Transcriptome and co-expression network analyses reveal the regulatory pathways and key genes associated with temperature adaptability in the yellow drum (Nibea albiflora)

The yellow drum (Nibea albiflora) is an important marine economy fish, that is widely distributed in the coastal waters of the Northwest Pacific. To understand the molecular regulatory mechanism of the yellow drum under temperature stress, transcriptome analysis was performed under five temperature conditions (10 °C, 15 °C, 20 °C, 24 °C, 28 °C) in the present study. Compared with 20 °C, 163, 401, 276, and 372 differentially expressed genes (DEGs) were obtained at 10 °C, 15 °C, 24 °C and 28 °C, respectively. Gene Ontology (GO) analysis indicated that the DEGs were mainly involved in cellular processes, metabolic processes, catalytic activity, membrane and binding. Meanwhile, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the temperature adaptive regulation of the yellow drum was mainly involved in signal transduction, metabolism, genetic information and protein processing. Weighted gene co-expression network analysis (WGCNA) showed that HMGB1, STAT4, Noct, C1q and CRT may be the key hub genes in the response of the yellow drum to temperature stress. In addition, 20 genes that may be associated with temperature stress were identified based on comparative analysis between the KEGG enrichment and the WGCNA results. Ten DEGs were selected for further validation using quantitative real-time PCR (qRT-PCR), and the results were consistent with the RNA-seq data. This study explored the transcriptional patterns of the yellow drum under temperature stress and provided fundamental information on the temperature adaptability of this species.

Proteome profile of spleen in rock bream (Oplegnathus fasciatus) naturally infected with rock bream iridovirus (RBIV)

BACKGROUND

Rock bream iridovirus (RBIV) is one of the most dangerous pathogens that causes the highest mortality in the aquaculture of rock bream (Oplegnathus fasciatus). Even though RBIV infection leads to huge economic loss, proteome studies on RBIV-infected rock bream have not been conducted to provide information about the differential protein expression pattern by the host protection system.

OBJECTIVE

The purpose of this study was to investigate the protein expression patterns in spleens of rock bream olive after infection by RBIV or mixed infection by RBIV and bacteria.

METHODS

Depending on the infection intensity and sampling time point, fish were divided into five groups: uninfected healthy fish at week 0 as the control (0C), heavily infected fish at week 0 (0H), heavily mixed RBIV and bacterial infected fish at week 0 (0MH), uninfected healthy fish at week 3 (3C), and lightly infected fish at week 3 (3L). Proteins were extracted from the spleens of infected rock bream. We used 2-DE analysis with LC-MS/MS to investigate proteome changes in infected rock bream.

RESULTS

The results of the LC-MS/MS analyses showed different protein expression profiles after infection. Proteins related to oxygen transport and energy generation, such as hemoglobin, beta-globin, and ATP synthase, were mostly expressed in the infected spleen. Whereas proteins involved in structure and cell movement, such as tubulin, myosin, actin binding proteins, and intermediate filament proteins, were down-regulated in the infected spleens. The protein expression profiles between infection by RBIV and mixed infection by RBIV and bacteria showed similar patterns.

CONCLUSIONS

Our results indicated that infection by RBIV or mixed infection by RBIV and bacteria triggered energy generation and oxygen-transport, but cell migration and constructional changes in the spleen were extremely decreased.

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.

Identifying Potential Candidate Hub Genes and Functionally Enriched Pathways in the Immune Responses to Quadrivalent Inactivated Influenza Vaccines in the Elderly Through Co-Expression Network Analysis

Insights into the potential candidate hub genes may facilitate the generation of safe and effective immunity against seasonal influenza as well as the development of personalized influenza vaccines for the elderly at high risk of influenza virus infection. This study aimed to identify the potential hub genes related to the immune induction process of the 2018/19 seasonal quadrivalent inactivated influenza vaccines (QIVs) in the elderly ≥60 years by using weighted gene co-expression network analysis (WGCNA). From 63 whole blood samples from16 elderly individuals, a total of 13,345 genes were obtained and divided into eight co-expression modules, with two modules being significantly correlated with vaccine-induced immune responses. After functional enrichment analysis, genes under GO terms of vaccine-associated immunity were used to construct the sub-network for identification and functional validation of hub genes. MCEMP1 and SPARC were confirmed as the hub genes with an obvious effect on QIVs-induced immunity. The MCEMP1 expression was shown to be negatively correlated with the QIVs-associated reactogenicity within 7 days after vaccination, which could be suppressed by the CXCL 8/IL-8 and exacerbated by the Granzyme-B cytotoxic mediator. Meanwhile, the SPARC expression was found to increase the immune responses to the QIVs and contribute to the persistence of protective humoral antibody titers. These two genes can be used to predict QIVs-induced adverse reaction, the intensity of immune responses, and the persistence of humoral antibody against influenza. This work has shed light on further research on the development of personalized QIVs with appropriate immune responses and long-lasting immunity against the forthcoming seasonal influenza.