Transcriptome profiling based on protein-protein interaction networks provides a core set of genes for understanding blood immune response mechanisms against Edwardsiella tarda infection in Japanese flounder (Paralichthys olivaceus)

Network-based integration of omics, physiological and environmental data in real-world Elbe estuarine Zander

Coastal and estuarine environments are under endogenic and exogenic pressures jeopardizing survival and diversity of inhabiting biota. Information of possible synergistic effects of multiple (a)biotic stressors and holobiont interaction are largely missing in estuaries like the Elbe but are of importance to estimate unforeseen effects on animals' physiology. Here, we seek to leverage host-transcriptional RNA-seq and gill mucus microbial 16S rRNA metabarcoding data coupled with physiological and abiotic measurements in a network analysis approach to decipher the impact of multiple stressors on the health of juvenile Sander lucioperca along one of the largest European estuaries. We find mesohaline areas characterized by gill tissue specific transcriptional responses matching osmosensing and tissue remodeling. Liver transcriptomes instead emphasized that zander from highly turbid areas were undergoing starvation which was supported by compromised body condition. Potential pathogenic bacteria, including Shewanella, Acinetobacter, Aeromonas and Chryseobacterium, dominated the gill microbiome along the freshwater transition and oxygen minimum zone. Their occurrence coincided with a strong adaptive and innate transcriptional immune response in host gill and enhanced energy demand in liver tissue supporting their potential pathogenicity. Taken together, we show physiological responses of a fish species and its microbiome to abiotic factors whose impact is expected to increase with consequences of climate change. We further present a method for the close-meshed detection of the main stressors and bacterial species with disease potential in a highly productive ecosystem.

Whole genome identification of olive flounder (Paralichthys olivaceus) cathepsin genes: Provides insights into its regulation on biotic and abiotic stresses response

Cathepsins are major lysosomal enzymes involved in essential physiological processes, including protein degradation, tissue differentiation, and innate or adaptive responses. Several kinds of cathepsins have been reported in teleost fishes, but no characterization have been performed for the inflammatory response of cathepsin family in olive flounder until now. In our current study, a total of 17 cathepsins in olive flounder were systematically identified and characterized. Phylogenetic analysis clearly indicated that the cathepsin genes was highly conserved. Analysis of structure and motifs exhibited high sequence similarity of cathepsin genes in olive flounder. Expression profiles of cathepsin genes in different tissues and developmental stages showed that cathepsins were temporally and spatially specific. RNA-seq analysis of bacteria and temperature stresses revealed that members of cathepsin were involved in inflammatory responses. Collectively, our findings would provide a further reference for understanding the molecular mechanisms of cathepsins in olive flounder.

Gene network analyses of Sepia esculenta larvae exposed to copper and cadmium: A comprehensive investigation of oxidative stress, immune response, and toxicological mechanisms

Copper (Cu) and Cadmium (Cd), prevalent heavy metals in marine environments, have known implications in oxidative stress, immune response, and toxicity in marine organisms. Sepia esculenta, a cephalopod of significant economic value along China's eastern coastline, experiences alterations in growth, mobility, and reproduction when subjected to these heavy metals. However, the specific mechanisms resulting from heavy metal exposure in S. esculenta remain largely uncharted. In this study, we utilized transcriptome and four oxidative, immunity, and toxicity indicators to assess the toxicological mechanism in S. esculenta larvae exposed to Cu and Cd. The measurements of Superoxide Dismutase (SOD), Malondialdehyde (MDA), Glutathione S-Transferase (GST), and Metallothioneins (MTs) revealed that Cu and Cd trigger substantial oxidative stress, immune response, and metal toxicity. Further, we performed an analysis on the transcriptome data through Weighted Gene Co-expression Network Analysis (WGCNA) and Protein-Protein Interaction (PPI) network analysis. Our findings indicate that exposure methods and duration influence the type and the extent of toxicity and oxidative stress within the S. esculenta larvae. We took an innovative approach in this research by integrating WGCNA and PPI network analysis with four significant physiological indicators to closely examine the toxicity and oxidative stress profiles of S. esculenta upon exposure to Cu and Cd. This investigation is vital in decoding the toxicological, immunological, and oxidative stress mechanisms within S. esculenta when subjected to heavy metals. It provides foundational insights capable of advancing invertebrate environmental toxicology and informs S. esculenta artificial breeding practices.

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

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

Preliminary study on toxicological mechanism of golden cuttlefish (Sepia esculenta) larvae exposed to cd

BACKGROUND

Cadmium (Cd) flows into the ocean with industrial and agricultural pollution and significantly affects the growth and development of economic cephalopods such as Sepia esculenta, Amphioctopus fangsiao, and Loligo japonica. As of now, the reasons why Cd affects the growth and development of S. esculenta are not yet clear.

RESULTS

In this study, transcriptome and four oxidation and toxicity indicators are used to analyze the toxicological mechanism of Cd-exposed S. esculenta larvae. Indicator results indicate that Cd induces oxidative stress and metal toxicity. Functional enrichment analysis results suggest that larval ion transport, cell adhesion, and some digestion and absorption processes are inhibited, and the cell function is damaged. Comprehensive analysis of protein-protein interaction network and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used to explore S. esculenta larval toxicological mechanisms, and we find that among the 20 identified key genes, 14 genes are associated with neurotoxicity. Most of them are down-regulated and enriched to the neuroactive ligand-receptor interaction signaling pathway, suggesting that larval nervous system might be destroyed, and the growth, development, and movement process are significantly affected after Cd exposure.

CONCLUSIONS

S. esculenta larvae suffered severe oxidative damage after Cd exposure, which may inhibit digestion and absorption functions, and disrupt the stability of the nervous system. Our results lay a function for understanding larval toxicological mechanisms exposed to heavy metals, promoting the development of invertebrate environmental toxicology, and providing theoretical support for S. esculenta artificial culture.

Transcriptome Profiling Based at Different Time Points after Hatching Deepened Our Understanding on Larval Growth and Development of Amphioctopus fangsiao

As the quality of life improves, there is an increasing demand for nutrition-rich marine organisms like fish, shellfish, and cephalopods. To address this, artificial cultivation of these organisms is being explored along with ongoing research on their growth and development. A case in point is Amphioctopus fangsiao, a highly valued cephalopod known for its tasty meat, nutrient richness, and rapid growth rate. Despite its significance, there is a dearth of studies on the A. fangsiao growth mechanism, particularly of its larvae. In this study, we collected A. fangsiao larvae at 0, 4, 12, and 24 h post-hatching and conducted transcriptome profiling. Our analysis identified 4467, 5099, and 4181 differentially expressed genes (DEGs) at respective intervals, compared to the 0 h sample. We further analyzed the expression trends of these DEGs, noting a predominant trend of continuous upregulation. Functional exploration of this trend entailed GO and KEGG functional enrichment along with protein-protein interaction network analyses. We identified GLDC, DUSP14, DPF2, GNAI1, and ZNF271 as core genes, based on their high upregulation rate, implicated in larval growth and development. Similarly, CLTC, MEF2A, PPP1CB, PPP1R12A, and TJP1, marked by high protein interaction numbers, were identified as hub genes and the gene expression levels identified via RNA-seq analysis were validated through qRT-PCR. By analyzing the functions of key and core genes, we found that the ability of A. fangsiao larvae to metabolize carbohydrates, lipids, and other energy substances during early growth may significantly improve with the growth of the larvae. At the same time, muscle related cells in A. fangsiao larvae may develop rapidly, promoting the growth and development of larvae. Our findings provide preliminary insights into the growth and developmental mechanism of A. fangsiao, setting the stage for more comprehensive understanding and broader research into cephalopod growth and development mechanisms.

Investigation of Gene Networks in Three Components of Immune System Provides Novel Insights into Immune Response Mechanisms against Edwardsiella tarda Infection in Paralichthys olivaceus

As a quintessential marine teleost, Paralichthys olivaceus demonstrates vulnerability to a range of pathogens. Long-term infection with Edwardsiella tarda significantly inhibits fish growth and even induces death. Gills, blood, and kidneys, pivotal components of the immune system in teleosts, elicit vital regulatory roles in immune response processes including immune cell differentiation, diseased cell clearance, and other immunity-related mechanisms. This study entailed infecting P. olivaceus with E. tarda for 48 h and examining transcriptome data from the three components at 0, 8, and 48 h post-infection employing weighted gene co-expression network analysis (WGCNA) and protein-protein interaction (PPI) network analysis. Network analyses revealed a series of immune response processes after infection and identified multiple key modules and key, core, and hub genes including xpo1, src, tlr13, stat1, and mefv. By innovatively amalgamating WGCNA and PPI network methodologies, our investigation facilitated an in-depth examination of immune response mechanisms within three significant P. olivaceus components post-E. tarda infection. Our results provided valuable genetic resources for understanding immunity in P. olivaceus immune-related components and assisted us in further exploring the molecular mechanisms of E. tarda infection in teleosts.

Identification of the matrix metalloproteinase (MMP) gene family in Japanese flounder (Paralichthys olivaceus): Involved in immune response regulation to temperature stress and Edwardsiella tarda infection

The Matrix metalloproteinase (MMP) gene family is responsible for regulating the degradation of Extra Cellular Matrix (ECM) proteins, which are important for physiological processes such as wound healing, tissue remodeling, and stress response. Although MMPs have been studied in many species, their role in immune response in Japanese flounder (Paralichthys olivaceus) is still not fully understood. This study conducted a comprehensive analysis of MMPs in flounder, including gene structures, evolutionary relationships, conserved domains, molecular evolution, and expression patterns. Analysis revealed that MMP genes could be grouped into 17 subfamilies and were evolutionarily conserved and functionally-constrained. Meanwhile, MMP genes were found to express in different embryonic and larval stages and might play the role of sentinel in healthy tissues. Furthermore, expression profiling showed that MMPs had diverse functions in environmental stress, with 60% (9/15) and 73% (11/15) of MMPs showing differential expression patterns under temperature stress and Edwardsiella tarda (E. tarda) infection, respectively. These findings provide a useful resource for understanding the immune functions of MMP genes in Japanese flounder.

Genome-wide identification of TNF receptor-associated factors in Japanese flounder (Paralichthys olivaceus) and functional analysis of resistance to temperature and Edwardsiella tarda stress

Tumor necrosis factor receptor-associated factors (TRAFs), as the signaling mediators of the tumor necrosis factor (TNFR) superfamily, toll-like receptors (TLR) and interleukin-1 receptor (IL-1R) superfamily, can activate downstream signal transduction pathways and play an important role in the body's immune process. In this study, six TRAF genes, namely PoTRAF2a, PoTRAF2b, PoTRAF3, PoTRAF4, PoTRAF6 and PoTRAF7, were identified and annotated in Japanese flounder by using bioinformatics methods. Phylogenetic analysis confirmed that TRAF genes can be divided into seven groups. Analysis of motif composition and gene structure demonstrated that all PoTRAF members were evolutionarily conserved. The expression patterns of PoTRAF genes were then further investigated in six different developmental stages and eleven tissues of healthy fish, and it was found that there were spatial and tissue specificities among the members. To investigate the immune response of Japanese flounder to abiotic and biotic stresses, we further analyzed the expression profile of PoTRAFs after temperature stress and pathogen challenge. The result showed that PoTRAF3 and PoTRAF4 were observably differentially expressed under temperature stress, indicating that they were involved in the immune response after temperature stress. The expression of PoTRAF2a, PoTRAF2b and PoTRAF4 was significantly different after E. tarda infection, suggesting that they might have antibacterial effects. These results would help to clarify the molecular roles of PoTRAF genes in the regulation of immune and inflammatory responses in Japanese flounder.

Characterization of Japanese flounder (Paralichthys olivaceus) STAT members: An immune-related gene family involved in Edwardsiella tarda and temperature stress

The signal transducer and activator of transcription (STAT) family members are not only the transcriptional activators, but also play important roles in regulating inflammatory response. Some members have been reported to be involved in innate bacterial and antiviral immunity in aquatic organisms. However, no systematic research on STATs has been found in teleost. In this present study, we characterized six STAT genes in Japanese flounder based on bioinformatics methods, namely PoSTAT1, PoSTAT2, PoSTAT3, PoSTAT4, PoSTAT5 and PoSTAT6. The phylogenetic analysis of STATs in fish indicated that STATs were highly conserved and revealed an absence of STAT5 in a few species. Further analysis of gene structures and motifs showed STAT proteins shared a similar structure and probably had similar functionality in Japanese flounder. The expression profiles of different development stages and tissues demonstrated that PoSTATs exhibited specificity in temporality and spatiality as well as PoSTAT4 was highly expressed in gill. The transcriptome data analysis of E. tarda and temperature stress showed that PoSTAT1 and PoSTAT2 were more respective to these two kinds of stress. In addition, the results also demonstrated that these PoSTATs might regulate immune response in different ways, manifested by up-regulation in E. tarda infection and down-regulation in temperature stress. In a word, this systematic analysis of PoSTATs would provide valuable information about the phylogenetic relationship of STATs in fish species and help understand the role of STAT genes in the immune response of Japanese flounder.

Gene network analyses of larvae under different egg-protecting behaviors provide novel insights into immune response mechanisms of Amphioctopus fangsiao

Amphioctopus fangsiao was a representative economic species in cephalopods, which was vulnerable to marine bacteria. Vibrio anguillarum was a highly infectious pathogen that have recently been found to infect A. fangsiao and inhibit its growth and development. There were significant differences in the immune response mechanisms between egg-protected and egg-unprotected larvae. To explore larval immunity under different egg-protecting behaviors, we infected A. fangsiao larvae with V. anguillarum for 24 h and analyzed the transcriptome data about egg-protected and egg-unprotected larvae infected with 0, 4, 12, and 24 h using weighted gene co-expression networks (WGCNA) and protein-protein interaction (PPI) networks. Network analyses revealed a series of immune response processes after infection, and identified six key modules and multiple immune-related hub genes. Meanwhile, we found that ZNF family, such as ZNF32, ZNF160, ZNF271, ZNF479, and ZNF493 might play significant roles in A. fangsiao immune response processes. We first creatively combined WGCNA and PPI network analysis to deeply explore the immune response mechanisms of A. fangsiao larvae with different egg-protecting behaviors. Our results provided further insights into the immunity of V. anguillarum infected invertebrates, and laid the foundation for exploring the immune differences among cephalopods with different egg protecting behaviors.

Genome-wide identification of olive flounder (Paralichthys olivaceus) SOCS genes: Involvement in immune response regulation to temperature stress and Edwardsiella tarda infection

The suppressors of cytokine signaling (SOCS) gene family participates in development and immunity through negative regulation of cytokine signaling pathways. Although the immune response of SOCS gene family members has been extensively characterized in teleost, no similar study has been reported in olive flounder yet. In our present study, a total of 13 SOCSs in olive flounder were identified and characterized systematically. By querying the SOCS sequences of ten teleost fish species, we found there were exactly more members of SOCSs in fish than mammals, which indicated that there were more duplication events occurred in fish than in higher vertebrates. Phylogenetic analysis clearly illuminated that SOCS genes were highly conserved. The analysis of gene structure and motif showed SOCS proteins of olive flounder shared a high level of sequence similarity strikingly. The expression profiles of tissues and developmental stages indicated that SOCS members had a kind of specificity in temporality and spatiality. RNA-Seq analysis of temperature stress and E. Tarda infection demonstrated SOCS members were involved in inflammatory response. In a word, our results would provide a further reference for understanding the mechanism of SOCS genes in olive flounder.

Transcriptome analysis preliminary reveals the immune response mechanism of golden cuttlefish (Sepia esculenta) larvae exposed to Cd

As a well-known marine metal element, Cd can significantly affect bivalve mollusk life processes such as growth and development. However, the effects of Cd on the molecular mechanisms of the economically important cephalopod species Sepia esculenta remain unclear. In this study, S. esculenta larval immunity exposed to Cd is explored based on RNA-Seq. The analyses of GO, KEGG, and protein-protein interaction (PPI) network of 1,471 differentially expressed genes (DEGs) reveal that multiple immune processes are affected by exposure such as inflammatory reaction and cell adhesion. Comprehensive analyses of KEGG signaling pathways and the PPI network are first used to explore Cd-exposed S. esculenta larval immunity, revealing the presence of 16 immune-related key and hub genes involved in exposure response. Results of gene and pathway functional analyses increase our understanding of Cd-exposed S. esculenta larval immunity and improve our overall understanding of mollusk immune functions.

High-quality Japanese flounder genome aids in identifying stress-related genes using gene coexpression network

The Japanese flounder is one of the most economically important marine flatfish. However, due to the increased frequency of extreme weather events and high-density industrial farming, an increasing number of environmental stresses have become severe threats to the healthy development of the Japanese flounder culture industry. Herein, we produced a high-quality chromosome-scale Japanese flounder genome using PacBio Circular Consensus Sequencing technologies. The assembled Japanese flounder genome spanned 588.22 Mb with a contig N50 size of 24.35 Mb. In total, 105.89 Mb of repetitive sequences and 22,565 protein-coding genes were identified by genome annotation. In addition, 67 candidate genes responding to distinct stresses were identified by gene coexpression network analysis based on 16 published stress-related RNA-seq datasets encompassing 198 samples. A high-quality chromosome-scale Japanese flounder genome and candidate stress-related gene set will not only serve as key resources for genomics studies and further research on the underlying stress responsive molecular mechanisms in Japanese flounder but will also advance the progress of genetic improvement and comprehensive stress-resistant molecular breeding of Japanese flounder.

Measurement(s)	genome assembly
Technology Type(s)	PacBio RS II

Splenic protection network revealed by transcriptome analysis in inactivated vaccine-immunized flounder (Paralichthys olivaceus) against Edwardsiella tarda infection

The protective immune response produced by fish after vaccination is crucial for vaccine effectiveness. Our previous studies have shown inactivated vaccine against Edwardsiella tarda can induce immune response in flounder (Paralichthys olivaceus). To elucidate the protective immune response at the genetic level, in this study, flounder was immunized with inactivated E. tarda for 5 weeks, and then they were challenged with E. tarda. The spleen was dissected at 7th day post immunization, 1st and 7th day post challenge, respectively. Transcriptome analysis showed that average of 46 million clean reads were obtained per library, while percentage of clean reads being mapped to reference genome was more than 89% in all cases, which suggested good quality of samples. As for differentially expressed genes (DEGs) identification in inactivated E. tarda groups, at 7th day post immunization, 1422 DEGs were identified and significantly enriched in innate immune-related pathways, such as Phagosome, Cell adhesion molecules and NF-kappa B signaling pathway; At 1st post challenge day, 1210 DEGs were identified and enriched to Antigen processing and presentation and Cell adhesion molecules, indicating that the pathogen was rapidly recognized and delivered; At 7th post challenge day, 1929 DEGs were identified, belonged to Toll-like receptor signaling pathway, Antigen processing and presentation, Th1 and Th2 cell differentiation and Th17 cell differentiation. Compared to 7th post immunization day, 73 immune-associated DEGs were identified at 1st post challenge day. Protein-protein interaction networks analysis revealed 11 hub genes (TLR7, TLR3, CXCR4, IFIH1, TLR8 etc), associated with recognition of pathogens and activation of innate immunity; while for 7th post challenge day, 141 immune-associated DEGs were identified. 30 hub genes (IL6, STAT1, HSP90A.1, TLR7, IL12β etc) were associated with stimulation of lymphocyte differentiation and activation of cellular immunity. Ten immune-related genes were randomly selected for RT-qPCR validation at each time point. In conclusion, data revealed protection of flounder against E. tarda infection by inactivated vaccine is mediated via immediate recognition of pathogen and subsequently activation of cellular immunity. Results give new aspect for vaccine protection cascades, is good references for vaccine evaluation.

Identification and functional analysis of dual-specificity phosphatases (DUSP) genes in Japanese flounder (Paralichthys olivaceus) against temperature and Edwardsiella tarda stress

Dual-specificity Phosphatases (DUSPs) are not only the key regulators of dephosphorylating and inactivating mitogen-activated protein kinases (MAPKs), but play a crucial role in the immune response. However, the role of DUSP genes in Japanese flounder (PoDUSPs) is still unclear. In this study, 28 DUSP genes in Japanese flounder were identified and classified based on the whole genome database. Phylogenetic analysis and protein structure analysis revealed that DUSPs had highly conserved domains in teleosts. Molecular evolution analysis indicated that the PoDUSP genes were conservative during evolution and were functional-constrained. Meanwhile, PoDUSP genes were found to express in different embryonic and larval stages which might play the role of sentinel in healthy organisms. Furthermore, PoDUSP genes' expression profiles after temperature stress and Edwardsiella tarda (E. tarda) infection were determined in Japanese flounder without precedent, and the results demonstrated that Podusp1, Podusp2 and Podusp16 were more respective to temperature variation whereas Podusp1 and Podusp6 were more respective to E. tarda infection. In summary, our results provide useful resources for understanding the immune responsibilities of DUSP genes in flatfish.

Immune responses to Vibrio vulnificus formalin-killed vaccine and ghost vaccine in Scophthalmus maximus

In this research, Vibrio vulnificus formalin-killed (FKCs) vaccine and ghost (VVGs) vaccine were successfully developed, and shown to prevent vibriosis of Scophthalmus maximus resulting from V. vulnificus. The antibody titre of FKCs and VVGs vaccine was 1: 2⁸ and 1: 2¹¹ . The RPS of FKCs and VVGs vaccine was 60% and 80%. In order to improve the understanding of vaccine protection mechanism, transcriptome data was used to analyse the immune response of S. maximus infected with V. vulnificus after vaccination with FKCs and VVGs vaccine. In the SmCon and SmIV groups, a series of innate immune-related genes were upregulated (such as, TLR5, Tp12, AP-1 and IL-1β) or downregulated (such as, CASP6 and CASP8), which suggested that the immune protection mechanism induced by inactivated vaccine was similar to that of autoimmune response. In the SmIV and SmGho group, a number of innate and adaptive immune-related genes (such as, STAT1, IFN-γ and MHC Ia) were activated, in which the expression of these genes was higher in SmGho, and VVGs vaccine induced stronger innate and acquired immune responses. In conclusion, the results lay a foundation for further study on the molecular mechanisms of immune protection induced by VVGs vaccine and FKCs vaccine.

Genome-wide identification and analysis of scavenger receptors and their expression profiling in response to Edwardsiella tarda infection in Japanese flounder (Paralichthys olivaceus)

The scavenger receptors (SRs) gene family, as one of pattern recognition receptors, participates in the innate immune response in diverse lineages. However, the systematic identification, characteristics and functions of SRs family are lacking in teleost. Here, we identified all 19 SRs family members in Japanese flounder (Paralichthys olivaceus) based on the genome and transcriptome data. Phylogenetic and Ka/Ks analysis demonstrated that these SRs genes were divided into five classes and all exhibited pronounced purified selection pressures. Whole genome duplication event was found in colec12, scarb2, and lamp1. Gene structure, functional domain and motif distribution analyses indicated that SRs within the different subfamilies are severely conservative. SRs genes showed diverse expression patterns in the embryogenesis and unchanged tissues. The regulations of 14 SRs genes in blood, gill and kidney after E. tarda infection suggested their roles in innate immune response. Meanwhile, ten SRs genes were differentially expressed after E. tarda stimulation in macrophages in vitro. Then we proved that PoSCARA3 could suppress the activity of NF-κB and AP-1 in HEK 293T cells by dual-luciferase assays. In summary, this study provided valuable basis for further functional characterization and immune functions of SRs genes in P. olivaceus.

Transcriptome profiling based on larvae at different time points after hatching provides a core set of gene resource for understanding the immune response mechanisms of the egg-protecting behavior against Vibrio anguillarum infection in Amphioctopus fangsiao

Mollusks have recently received increasing attention because of their unique immune systems. Mollusks such as Amphioctopus fangsiao are economically important cephalopods, and the effects of their egg-protecting behavior on the larval immune response are unclear. Meanwhile, little research has been done on the resistance response of cephalopod larvae infected with pathogenic bacteria such as Vibrio anguillarum. In this study, V. anguillarum was used to infect the primary hatching A. fangsiao larvae under different egg-protecting behaviors for 24 h, and a total of 7156 differentially expressed genes (DEGs) were identified at four time points after hatching based on transcriptome analysis. GO and KEGG enrichment analyses showed that multiple immune-related GO terms and KEGG signaling pathways were enriched. Protein-protein interaction networks (PPI networks) were used to search functional relationships between immune-related DEGs. Finally, 20 hub genes related to multiple gene functions or involved in multiple signaling pathways were identified, and their accuracy was verified using quantitative RT-PCR. PPI networks were first used to study the effects A. fangsiao larvae after infection with V. anguillarum under different egg-protecting behaviors. The results provide significant genetic resources for exploring invertebrate larval immune processes. The data lays a foundation for further study the immune response mechanisms for invertebrates after infection.

Transcriptome analysis and protein-protein interaction in resistant and susceptible families of Japanese flounder (Paralichthys olivaceus) to understand the mechanism against Edwardsiella tarda

Edwardsiella tarda is one of the most harmful bacterial pathogens for aquaculture flatfish. After artificial infection of 47 Japanese flounder (Paralichthys olivaceus) families, resistant and susceptible families were identified in this study. High-throughput sequencing was performed on the liver transcriptome of uninfected groups (PoRU and PoSU) and infected groups (PoRC and PoSC). Through assembly and annotation, a total of 3012 and 1386 differentially expressed genes (DEGs) were identified in PoRU vs. PoSU and PoRC vs. PoSC. The significant enrichment pathways between PoRU and PoSU were mainly in metabolic and biosynthesis pathways. A total of thirty dominant enrichment pathways between PoRC and PoSC mainly focused on some immune-related pathways, including the hematopoietic cell lineage, cytokine-cytokine receptor interaction, complement and coagulation cascades, antigen processing and presentation, the intestinal immune network for immunoglobulin A (IgA) production and T/B cell receptor signaling pathway. Under the protein-protein interaction (PPI) analysis, hub genes, including CD molecules, complement component factors and chemokines, were identified in the network, and 16 core genes were differentially expressed in resistant and sustainable families in quantitative polymerase chain reaction (qPCR) validation. This study represents the first transcriptome analysis based on resistant and susceptible families and provides resistant genes to understand the potential molecular mechanisms of antibacterial function in marine fish. The results obtained in this study provide crucial information on gene markers for resistant breeding of Japanese flounder.

Transcriptome Profiling Based on Larvae at Different Time Points After Hatching Provides a Core Set of Gene Resource for Understanding the Metabolic Mechanisms of the Brood-Care Behavior in Octopus ocellatus

The metabolic processes of organisms are very complex. Each process is crucial and affects the growth, development, and reproduction of organisms. Metabolism-related mechanisms in Octopus ocellatus behaviors have not been widely studied. Brood-care is a common behavior in most organisms, which can improve the survival rate and constitution of larvae. Octopus ocellatus carried out this behavior, but it was rarely noticed by researchers before. In our study, 3,486 differentially expressed genes (DEGs) were identified based on transcriptome analysis of O. ocellatus. We identify metabolism-related DEGs using GO and KEGG enrichment analyses. Then, we construct protein-protein interaction networks to search the functional relationships between metabolism-related DEGs. Finally, we identified 10 hub genes related to multiple gene functions or involved in multiple signal pathways and verified them using quantitative real-time polymerase chain reaction (qRT-PCR). Protein-protein interaction networks were first used to study the effects of brood-care behavior on metabolism in the process of growing of O. ocellatus larvae, and the results provide us valuable genetic resources for understanding the metabolic processes of invertebrate larvae. The data lay a foundation for further study the brood-care behavior and metabolic mechanisms of invertebrates.

Identification, evolution and expression analyses of mapk gene family in Japanese flounder (Paralichthys olivaceus) provide insight into its divergent functions on biotic and abiotic stresses response

Mitogen-activated protein kinase (MAPK) are a series of serine/threonine protein kinases showing evolutionary conservation, which can be activated by many stimulus signals and then transfer them from cell membrane to nucleus. MAPKs regulate a variety of biological processes, such as apoptosis, hormone signaling and immune response. In this study, 14 putative mapk genes in Japanese flounder were identified, and their basic physical and chemical properties were characterized. Phylogenetic analysis showed that mapk genes were divided into three main subfamilies, including ERK, JNK and the p38 MAPK. Selection pressure analysis revealed they were evolutionarily-constrained and undergone strong purifying selection. Gene structure and conserved protein motif comparison suggested high levels of conservation in members of mapk gene family. The expression patterns were further investigated in each embryonic and larval development stages and different tissues. In addition, RNA-seq analyses after bacteria and temperature stresses suggested mapk genes had different expression patterns. Three mapk genes showed significant differences in response to E. tarda challenge and five were induced significantly after temperature stress, indicating their potential functions. This systematic analysis provided valuable information for further understanding of the regulation mechanism of mapk gene family under different stresses in Japanese flounder.

Transcriptome profiling based on protein-protein interaction networks provides a set of core genes for understanding the immune response mechanisms of the egg-protecting behavior in Octopus ocellatus

Protection via of the immune system is indispensable to the life of organisms. Within an immune network, problems with a given link will affect the normal life activities of the organism. Octopus ocellatus is cephalopod widely distributed throughout the world's oceans. Because of its unique nervous system and locomotive organs, research on this species has gradually increased in recent years. Many immune response mechanisms associated with behaviors of O. ocellatus are still unclear. Moreover, as a factor affecting the normal growth of O. ocellatus, egg protection has rarely been considered in previous behavioral studies. In this study, we analyzed the transcriptome profile of gene expression in O. ocellatus larvae, and identified 5936 differentially expressed genes (DEGs). GO and KEGG enrichment analyses were used to search for immune-related DEGs. Protein-protein interaction networks were constructed to examine the interactions between immune-related genes. Fifteen hub genes involved in multiple KEGG signaling pathways or with multiple protein-protein interaction relationships were obtained and verified by quantitative RT-PCR. We first studied the effects of egg protection on the immunity of O. ocellatus larvae by means of protein-protein interaction networks, and the results provide valuable genetic resources for understanding the immunity of invertebrate larvae. The data serve as a foundation for further research on the egg-protecting behavior of invertebrates.

Blood Transcriptomics of Turbot Scophthalmus maximus: A Tool for Health Monitoring and Disease Studies

Simple Summary

The analysis of blood gene expression is emerging as a relevant source of information about the health status of an organism. While these investigations are increasingly performed in human and terrestrial animals, their potential is still underexplored in fish pathology. The aim of this work was to analyze the blood transcriptional profile of a commercially important flatfish species, turbot (Scophthalmus maximus), in healthy and diseased specimens. The analysis of the most expressed genes in healthy fish indicated that turbot red blood cells have important immunological functions. In diseased fish, parasitized by a myxozoan, the blood analysis reflected a broad inhibition of the immune response followed by intense inflammatory activation in heavy infections. The results showed that turbot response appears delayed, dysregulated and ineffective in stopping the infection. Particularly, a proper development of the adaptive immune response was lacking. This study points out that blood gene expression profiling is a reliable tool for health monitoring, as well as to advance in the knowledge of fish immunity and diseases.

Abstract

Blood transcriptomics is emerging as a relevant tool to monitor the status of the immune system and assist in diagnosis, prognosis, treatment and pathogenesis studies of diseases. In fish pathology, the potential of transcriptome profiling of blood is still poorly explored. Here, RNA sequencing was applied to analyze the blood transcriptional profile of turbot (Scophthalmus maximus), the most important farmed flatfish. The study was conducted in healthy specimens and specimens parasitized by the myxozoan Enteromyxum scophthalmi, which causes one of the most devastating diseases in turbot aquaculture. The blood of healthy turbot showed a transcriptomic profile mainly related to erythrocyte gas transportation function, but also to antigen processing and presentation. In moderately infected turbot, the blood reflected a broad inhibition of the immune response. Particularly, down-regulation of the B cell receptor signaling pathway was shared with heavily parasitized fish, which showed larger transcriptomic changes, including the activation of the inflammatory response. Turbot response to enteromyxosis proved to be delayed, dysregulated and ineffective in stopping the infection. The study evinces that blood transcriptomics can contribute to a better understanding of the teleost immune system and serve as a reliable tool to investigate the physiopathological status of fish.

Blood Will Tell: What Hematological Analyses Can Reveal About Fish Welfare

Blood analyses provide substantial information about the physiological aspects of animal welfare assessment, including the activation status of the neuroendocrine and immune system, acute and long-term impacts due to adverse husbandry conditions, potential diseases, and genetic predispositions. However, fish blood is still not routinely analyzed in research or aquaculture for the assessment of health and/or welfare. Over the years, the investigative techniques have evolved from antibody-based or PCR-based single-parameter analyses to now include transcriptomic, metabolomic, and proteomic approaches and from hematological observations to fluorescence-activated blood cell sorting in high-throughput modes. The range of testing techniques established for blood is now broader than for any other biogenic test material. Evaluation of the particular characteristics of fish blood, such as its cell composition, the nucleation of distinct blood cells, or the multiple isoforms of certain immune factors, requires adapted protocols and careful attention to the experimental designs and interpretation of the data. Analyses of fish blood can provide an integrated picture of the endocrine, immunological, reproductive, and genetic functions under defined environmental conditions and treatments. Therefore, the scarcity of high-throughput approaches using fish blood as a test material for fish physiology studies is surprising. This review summarizes the wide range of techniques that allow monitoring of informative fish blood parameters that are modulated by different stressors, conditions, and/or treatments. We provide a compact overview of several simple plasma tests and of multiparametric analyses of fish blood, and we discuss their potential use in the assessment of fish welfare and pathologies.

RGD-binding integrins mediated phagocytosis involved in the entry of Edwardsiella tarda into mudskipper MO/MФ

The versatile fish pathogen Edwardsiella tarda is an intracellular pathogen with the ability to invade and replicate in host phagocytes. However, the mechanism mediating the uptake of E. tarda in fish monocytes/macrophages (MO/MΦ) is not yet understood. Generating mudskipper kidney-derived MO/MФ transcriptomic resources from mudskipper challenged by E. tarda is crucial for understanding the molecular mechanisms underlying the mudskipper invasion process. In the present study, a total of 1185 up-regulated and 885 down-regulated differentially expressed genes (DEGs) were identified using RNA-seq. Enrichment and pathway analysis of DEGs revealed the centrality of the phagosome and regulation of actin cytoskeleton pathways in pathogen entry. The progress of phagosome formation was observed by transmission electron microscopy. Eight conserved integrin (ITG) subunit genes, belonging to the phagocytic receptors, were found in the transcriptomic sequence data. Additionally, quantitative real-time PCR showed that the mRNA expressions of most ITG subunit genes were related to the different infection times of E. tarda and the different bacterial pathogens. Further assays demonstrated that phagocytosis of FITC-labeled E. tarda by mudskipper MO/MФ was significantly reduced by the tetrapeptide Asp-Gly-Arg-Ser (RGDS). In summary, phagocytosis is one of the entry pathways into mudskipper MO/MΦ, and RGD-binding ITGs are involved in the phagosome formation process.

Impact of Ringer's Solution Challenge Stress to Immunostimulatory Experiment, Insights From Japanese Flounder

Ringer's or phosphate buffer saline (PBS) solution buffer usually was used as dilution butter in intraperitoneal injection. Stress could activate immune response, inflammatory response and glycogen metabolic process. The impact of solution buffer as a stressor to immune system was ignored in immunostimulatory experiment. In this report, we tested the hypothesis that the innate immune response and glycogen metabolic process were altered when it were challenged with Ringer's in Japanese flounder (Paralichthys olivaceus). RNA-seq was performed after challenge with Ringer's at 8 h and 48 h. The data revealed that the expression profiles of blood, gill, and kidney were significantly changed. Differentially expressed genes (DEGs) were identified, and energy metabolic and immune-related genes were up-regulated or down-regulated obviously. GO and KEGG analyses showed that DEGs were mainly enriched in innate immune terms and pathways. Weighted gene co-expression networks analysis (WGCNA) also indicated the highest association module with stress. A total of 16 genes were detected in the gray module, which were immune-related and metabolic-related genes. These results provided fundamental information on intraperitoneal injection with solution buffer. It offered useful clues to further explore the functional mechanism of stress and immunity.

Transcriptome Analysis of Paralichthys olivaceus Erythrocytes Reveals Profound Immune Responses Induced by Edwardsiella tarda Infection

Unlike mammalian red blood cells (RBCs), fish RBCs are nucleated and thus capable of gene expression. Japanese flounder (Paralichthys olivaceus) is a species of marine fish with important economic values. Flounder are susceptible to Edwardsiella tarda, a severe bacterial pathogen that is able to infect and survive in flounder phagocytes. However, the infectivity of and the immune response induced by E. tarda in flounder RBCs are unclear. In the present research, we found that E. tarda was able to invade and replicate inside flounder RBCs in both in vitro and in vivo infections. To investigate the immune response induced by E. tarda in RBCs, transcriptome analysis of the spleen RBCs of flounder challenged with E. tarda was performed. Six sequencing libraries were constructed, and an average of 43 million clean reads per library were obtained, with 85% of the reads being successfully mapped to the genome of flounder. A total of 1720 differentially expressed genes (DEGs) were identified in E. tarda-infected fish. The DEGs were significantly enriched in diverse Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, especially those associated with immunity, disease, and infection. Ninety-one key DEGs involved in 12 immune-related pathways were found to form extensive interaction networks. Twenty-one genes that constituted the hub of the networks were further identified, which were highly regulated by E. tarda and involved in a number of immune processes, notably pathogen recognition and signal transduction, antigen processing and presentation, inflammation, and splicing. These results provide new insights into the immune role of flounder RBCs during bacterial infection.

Gene network analysis reveals a core set of genes involved in the immune response of Japanese flounder (Paralichthys olivaceus) against Vibrio anguillarum infection

Japanese flounder (Paralichthys olivaceus) is one of the most economically important marine fish cultured in north Asia. Vibrio anguillarum is a severe bacterial pathogen to Japanese flounder and many other aquaculture species. In order to understand the immune response of flounder during bacterial infection, we systematically examined the transcriptome profiles of flounder spleen at three time points after V. anguillarum challenge. More than one billion high quality reads were obtained, approximately 80.70% of which were successfully mapped to the reference genome of flounder. A total of 6060, 4688 and 4235 differentially expressed genes (DEGs) were captured at 6, 12 and 24-h post-infection, respectively. The DEGs exhibited dynamic changes in expression and were assigned into four different profiles based on expression trend. GO and KEGG analysis showed that the DEGs were enriched in various immune-related terms, including response to stimulation, immune system and pathways of cytokine-cytokine receptor interaction, Jak-STAT signaling and Toll-like receptor signaling. Furthermore, a network of highly interactive DEGs involved in 11 immune-related pathways was detected by utilizing the weighted co-expressing network analysis (WGCNA). Accordingly, 26 hub genes were discovered that constituted an elaborate immune regulatory network and functioned mainly in pathogen recognition, antigen processing, and molecular signaling. The results of this study provided the first systematical transcriptome profile of flounder in association with V. anguillarum infection and can serve as a valuable resource of target genes for future studies on the molecular mechanisms underlying the immune defense of flounder against bacterial infection.

Comparative expression analysis in small cell lung carcinoma reveals neuroendocrine pattern change in primary tumor versus lymph node metastases

Background

Recent preclinical data suggest that neuroendocrine (NE) subtype of small cell lung cancer (SCLC) has strong therapeutic relevance. NE high tumors are associated with immune desert and NE low tumors are considered to have an immune oasis phenotype. Our aim was to investigate the NE phenotypes of surgically resected SCLC tumors according to inter-tumor heterogeneity.

Methods

Expression analysis for 2,560 genes was performed in 32 surgically resected SCLC patients' primary tumors and corresponding lymph node (LN) metastases. To analyze tumor heterogeneity, we examined the differences in the gene expression of primary tumors versus LN metastases. We performed cluster analysis and heat map to divide patients into NE high and low subtypes by using the top NE-associated genes described in preclinical studies.

Results

We found 6% (n=154) genes with significant differences and only 13.1% (n=336) of all genes in the panel had a strong correlation between the primary tumor and LN metastases. Cluster analysis clearly distinguished SCLC NE high versus low subtypes both in primary tumor (20 vs. 12, respectively) and LNs (23 vs. 9, respectively). As for inter-tumor heterogeneity, in case of five patients, a change in the NE pattern was observed. Specifically, we found significant downregulation of the NE-associated genes CAV1 (P=0.004), CAV2 (P=0.029) and ANXA3 (P=0.035) in their LN metastases compared to their primary tumor.

Conclusions

Our data confirm the results of preclinical studies and clearly distinguish NE low and high differentiation clusters in SCLC. Moreover, they highlight the gene expression discordance between primary tumors and corresponding LN metastases suggesting that the NE pattern of metastatic LNs might not reflect that of the primary tumor. Altogether, by shedding light on the diversity of SCLC, the current study might help to improve patient selection and treatment in this devastating disease.

Keywords

Small cell lung cancer (SCLC); neuroendocrine tumor; lymph node metastasis; tumor heterogeneity; RNA sequencing.

A Comparative Analysis of Edwardsiella tarda-Induced Transcriptome Profiles in RAW264.7 Cells Reveals New Insights into the Strategy of Bacterial Immune Evasion

Edwardsiella tarda is a Gram-negative bacterial pathogen with a broad host range, including fish, reptiles, and mammals. One prominent virulence feature of E. tarda is its ability to survive and replicate in host phagocytes, but the relevant molecular mechanism is largely unknown. In this study, we examined the transcriptome profiles of RAW264.7 cells, a murine macrophage cell line, infected with live E. tarda or stimulated with dead E. tarda for 4 h and 8 h. Eighteen libraries were constructed, and an average of 69 million clean reads per library were obtained, with ~81.63% of the reads being successfully mapped to the reference genome. In total, 208 and 232 differentially expressed genes (DEGs) were identified between live and dead E. tarda-treated cells at 4 h and 8 h post-infection, respectively. The DEGs were markedly enriched in the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with immunity. Live E. tarda differed strikingly from dead E. tarda in the regulation of immune related genes. Compared with dead E. tarda-treated cells, live E. tarda-treated cells exhibited marked and significant suppression in the induction of a large amount of immune genes, including RIG-I-like receptors, cytokines, and interferon-related genes. Furthermore, some of the immune genes highly regulated by live E. tarda formed complicated interaction networks with each other. Together, the results of this study revealed a transcriptome profile specifically induced by the active virulence elements of live E. tarda during the infection process, thus adding new insights into the intracellular infection mechanism of E. tarda. This study also provided a valuable set of target genes for further study of the immune evasion strategy of E. tarda.

Genome-wide identification and characterization of heat shock protein family 70 provides insight into its divergent functions on immune response and development of Paralichthys olivaceus

Flatfish undergo extreme morphological development and settle to a benthic in the adult stage, and are likely to be more susceptible to environmental stress. Heat shock proteins 70 (hsp70) are involved in embryonic development and stress response in metazoan animals. However, the evolutionary history and functions of hsp70 in flatfish are poorly understood. Here, we identified 15 hsp70 genes in the genome of Japanese flounder (Paralichthys olivaceus), a flatfish endemic to northwestern Pacific Ocean. Gene structure and motifs of the Japanese flounder hsp70 were conserved, and there were few structure variants compared to other fish species. We constructed a maximum likelihood tree to understand the evolutionary relationship of the hsp70 genes among surveyed fish. Selection pressure analysis suggested that four genes, hspa4l, hspa9, hspa13, and hyou1, showed signs of positive selection. We then extracted transcriptome data on the Japanese flounder with Edwardsiella tarda to induce stress, and found that hspa9, hspa12b, hspa4l, hspa13, and hyou1 were highly expressed, likely to protect cells from stress. Interestingly, expression patterns of hsp70 genes were divergent in different developmental stages of the Japanese flounder. We found that at least one hsp70 gene was always highly expressed at various stages of embryonic development of the Japanese flounder, thereby indicating that hsp70 genes were constitutively expressed in the Japanese flounder. Our findings provide basic and useful resources to better understand hsp70 genes in flatfish.

Functional differentiation of three phosphatidylinositol 3-kinase catalytic subunit alpha (PIK3CA) in response to Vibrio anguillarum infection in turbot (Scophthalmus maximus)

PIK3CA has been extensively investigated from its molecular mechanism perspective and association with its mutations in different types of cancers. However, little has been reported regarding the pathological significance of PIK3CA expression in teleost. Here, in our present study, three PIK3CA genes termed SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were firstly identified in the genome of turbot S. maximus. Although these three genes located in different chromosomes, all of them share the same five domains. Phylogenetic and synteny analysis indicated that SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were three paralogs that may originate from duplication of the same ancestral PIK3CA gene. Subcellular localization analysis confirmed the cytoplasm distribution of these three paralogs. All three SmPIK3CA were ubiquitously expressed in examined tissues in turbot, with the higher expression levels in immune-related tissues such as blood, spleen, kidney, gills and intestines. Upon Vibrio anguillarum challenge, SmPIK3CAa and SmPIK3CA-like transcripts were significantly induced in spleen, intestine and blood despite of differential expression levels and responsive time points. Additionally, individuals in resistant group showed significantly higher expression level of both two genes than in the susceptible group. Moreover, four SNPs (102, 2530, 3027 and 3060) and one haplotype (Hap2) located in exon region of SmPIK3CA-like were identified and confirmed to be associated with V. anguillarum resistance in turbot by association analysis in different populations. Taken together, these results suggested that functional differentiation occurred in three SmPIK3CA paralogs with Vibrio anguillarum resistance and SmPIK3CAa and SmPIK3CA-like probable play potential roles in innate immune response to pathogenic invasions in turbot.

Metabolites-Enabled Survival of Crucian Carps Infected by Edwardsiella tarda in High Water Temperature

Temperature is one of the major factors that affect the outbreak of infectious disease. Lines of evidences have shown that virulence factors can be controlled by thermo-sensors in bacterial pathogens. However, how temperature influences host's responses to the pathogen is still largely unexplored, and the study of this might pave the way to develop strategies to manage pathogenic bacterial infection. In the present study, we show that finfish Carassius carassius, the crucian carp that is tolerant to a wide range of temperatures, is less susceptible to bacterial infection when grown in 20°C than in 30°C. The different responses of C. carassius to bacterial infection could be partially explained by the distinct metabolisms under the specific temperatures: C. carassius shows elevated tricarboxylic acid cycle (TCA cycle) but decreased taurine and hypotaurine metabolism as well as lower biosynthesis of unsaturated fatty acids at 30°C. The decreased abundance of palmitate, threonine, and taurine represents the most characteristic metabolic feature. Consistently, exogenous palmitate, threonine, or taurine enhances the survival of C. carassius to bacterial infection at 30°C in a dose-dependent manner. This effect could be attributed to the inhibition on the TCA cycle by the three metabolites. This notion is further supported by the fact that low concentration of malonate, a succinate dehydrogenase inhibitor, increases the survival of C. carassius at 30°C as well. On the other hand, addition of the three metabolites rescued the decreased expression of pro-inflammatory cytokines including TNF-α1, TNF-α2, IL-1β1, IL-1β2, and lysozyme at 30°C. Taken together, our results revealed an unexpected relationship between temperature and metabolism that orchestrates the immune regulation against infection by bacterial pathogens. Thus, this study shed light on the modulation of finfish physiology to fight against bacterial infection through metabolism.

Edwardsiella tarda-induced miR-7a functions as a suppressor in PI3K/AKT/GSK3β signaling pathway by targeting insulin receptor substrate-2 (IRS2a and IRS2b) in Paralichthys olivaceus

To study the effect of Edwardsiella tarda infection on miRNAs expression profile in Japanese flounder, fish were injected intraperitoneally with E. tarda. The miRNAs involved in regulating immune responses were analyzed by high-throughput sequencing. A total of 164 mature miRNAs were identified, of which 17 miRNAs were differentially expressed (DE miRNAs) after E. tarda infection, indicating that they were immune-related miRNAs. To further examine the relationship between the miRNAs and their predicted target mRNAs, a total of 22 predicted target mRNAs, mainly related to endocytic signaling pathway, NF-κB signaling pathway, and p53 signaling pathway, were detected with miRNA mimics in HEK-293T cells by dual-luciferase reporter experiments. Finally, we confirmed that insulin receptor substrate-2 (IRS2a and IRS2b) were regulated by miR-7a. And the target sites of the 3' untranslated region (UTR) of IRS2a and IRS2b were verified by dual-luciferase reporter experiments. Furthermore, we found that the E. tarda and LPS significantly increased host miR-7a expression. In vivo and in vitro studies revealed that IRS2-mediated PI3K/AKT/GSK3β signaling pathway was suppressed. Taken together, these results implied that miR-7a might be a key regulator of PI3K/AKT/GSK3β signaling pathway via suppressing the IRS2a and IRS2b genes.

Comparative Study of Immune Reaction Against Bacterial Infection From Transcriptome Analysis

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

Transcriptome profiling reveals key roles of phagosome and NOD-like receptor pathway in spotting diseased Strongylocentrotus intermedius

Spotting disease is a common disease in the process of aquaculture and restocking of the sea urchin Strongylocentrotus intermedius and leads to mass mortality. To characterize the molecular processes and candidate genes related to spotting disease in S. intermedius, we conducted next-generation sequencing to assess the key genes/pathways in spotting diseased sea urchin (DUG) compared to healthy ones (HUG). A total of 321.1 million clean reads were obtained and assembled into 93,877 Unigenes with an N50 of 1185 bp, in which 86.48% of them matched to the genome sequence of the sea urchin S. purpuratus and 27,456 Unigenes mapped to Nr database. Salmon expression analysis revealed 1557 significantly differently expressed genes (DEGs) between DUG and HUG. These DEGs were enriched into 151 KEGG pathways including a core set of immune correlated pathways notably in phagosome and NOD-like receptor signaling. DUG displayed an obvious downregulation in these immune pathways. The expression patterns of six DEGs were confirmed by RT-qPCR, and the expressions were consistent with the results of RNA-seq. Furthermore, 15,990 SSRs were identified and a total of 235,249 and 295,567 candidate SNPs were identified from DUG and HUG, respectively. All these results provided basic information for our understanding of spotting disease outbreak in sea urchin.

Genome-wide association analysis of intra-specific QTL associated with the resistance for enteric septicemia of catfish

Disease resistance is one of the most important traits for aquaculture industry. For catfish industry, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the most severe disease, causing enormous economic losses every year. In this study, we used three channel catfish families with 900 individuals (300 fish per family) and the 690K catfish SNP array, and conducted a genome-wide association study to detect the quantitative trait loci (QTL) associated with ESC resistance. Three significant QTL, with two of located on LG1 and one on LG26, and three suggestive QTL located on LG1, LG3, and LG21, respectively, were identified to be associated with ESC resistance. With a well-assembled- and -annotated reference genome sequence, genes around the involved QTL regions were identified. Among these genes, 37 genes had known functions in immunity, which may be involved in ESC resistance. Notably, nlrc3 and nlrp12 identified here were also found in QTL regions of ESC resistance in the channel catfish × blue catfish interspecific hybrid system, suggesting this QTL was operating within both intra-specific channel catfish populations and interspecific hybrid backcross populations. Many of the genes of the Class I MHC pathway, for mediated antigen processing and presentation, were found in the QTL regions. The positional correlation found in this study and the expressional correlation found in previous studies indicated that Class I MHC pathway was significantly associated with ESC resistance. This study validated one QTL previously identified using the second and fourth generation of the interspecific hybrid backcross progenies, and identified five additional QTL among channel catfish families. Taken together, it appears that there are only a few major QTL for ESC disease resistance, making marker-assisted selection an effective approach for genetic improvements of ESC resistance.