Xiaoyan Wu
Jing Xing
Xiaoqian Tang
Xiuzhen Sheng
Heng Chi
Wenbin Zhan- Laboratory of Pathology and Immunology of Aquatic Animals, Key Laboratory of Mariculture, MOE, Ocean University of China, Qingdao, China
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.
Introduction
Vaccines induce protective immune responses in fish (1, 2). The immune protection involves the relative percentage survival (RPS), antibody production and T/B lymphocyte response against challenge etc (3–5). RPS has been widely used to assess vaccine protection (6, 7). Besides of the RPS, the production of antibodies is the strategy of vaccines. Antibodies specifically recognize and bind antigens, which promote phagocytosis and achieve clearance of pathogens (8, 9). T/B lymphocytes are essential component of the adaptive immune response for assessing vaccine protection (10, 11). Study in turbot (Scophthalmus maximus) immunized with inactivated bivalent vaccine (IVVah1) showed high RPS to Vibrio anguillarum and Vibrio harveyi infection were maintained from week 4 to week 8 post immunization. In addition to RPS, antibody levels showed a trend of increasing and then decreasing, peaking at weeks 2 and 3. After 8 weeks of immunization against the pathogen, antibodies increased within a week. This indicates that at the protein and individual level, turbot immunized with IVVah1 are protected against the pathogen (12). In the study of vaccine-immunized flounder (Paralichthys olivaceus) challenged with Edwardsiella tarda, antibodies increased significantly at week 3. The percentage of T lymphocytes peaked at day 7. The percentage of IgM+ B lymphocytes showed a trend of increasing and then decreasing, peaking at week 2. This suggested that inactivated vaccine enhanced protective humoral and cellular immune responses in fish after challenge (13). In addition to this, the immune protection of fish vaccines needs to be supported by fine grained networks and comprehensive data. This is urgently needed for the development and application of effective vaccines.
Flounder is valuable marine fish. In recent years, with the rapid development of intensive aquaculture, the incidence of disease outbreaks has been increasing, which seriously affects the culture of flounder (14). E. tarda is a gram-negative bacterium with intracellular parasitism (15). It causes high mortality and significant economic losses in flounder through invasion of epithelial cells, production of toxins, and evasion of phagocyte-mediated killing (16, 17). Safe and effective vaccines are essential for the control of bacterial diseases (18, 19), among which inactivated vaccines have become commercialized vaccine species (20). For elucidating the protective immunity mechanism of inactivated vaccines, it is essential to study the response of immune organs in pre-vaccinated flounder after infection (21). Similar to mammals, the spleen of fish is the main peripheral lymphoid organ (22). Here most of the antigens in the blood are captured and engulfed by macrophages. At the same time, it is also the site of aggregation of T and B lymphocytes for antigen presentation and initiation of adaptive immune responses. The spleen is an essential organ for resolving the immune response in fish (23–25).
Transcriptome sequencing analyzes gene expression dynamics at the level of individual transcripts, which can contribute to the resolution of immune-related genetic information and functional molecules (26, 27). RNA-Seq has been applied to analyze immune-related genes and signaling pathways in fish such as flounder, grass carp, large yellow croaker and turbot (28–31). The spleen red blood cells of flounder challenged with E. tarda were analyzed by transcriptome. 21 key genes were identified, mainly involved in antigen processing and presentation, pathogen recognition and inflammation (32). The phagosome pathway was activated in the head kidney of turbot inoculated with bivalent inactivated bacteria vaccine Aeromonas salmonicida and Vibrio scophthalmi. Antigenic peptide transport protein 1 (TAP1), complement fraction 3 (C3) and mannose receptor (MR) were significantly upregulated (33). However, it is often limited to the analysis of tissues such as head kidney, gill and blood at different time points after immunization or infection. This is still inadequate for supporting the molecular mechanisms of immune protection in fish. In addition, genes in organisms are functionally interconnected and together control the activities of the organism (34). For example, kidney from flounder infected with E. tarda was sequenced. Immune-related genes were used to construct protein-protein interaction networks (35). The identification of hub genes helps to understand the protective mechanism of the vaccine.
In this study, flounder were immunized with inactivated E. tarda vaccine for five weeks and then challenged with E. tarda. The infection status of flounder was evaluated. Spleens of flounder were sampled for transcriptome analysis at 7th post immunization day, 1st and 7th post challenge day. RNA sequencing (RNA-Seq) was performed using the Illumina Novaseq6000 platform. The results showed that an average of 46 million clean reads were obtained per library, with clean reads accounting for more than 99% of raw reads and Q30 greater than 92% in all cases; while the percentage of clean reads being mapped to the reference genome was more than 89% in all cases. 1422, 1210 and 1929 differentially expressed genes (DEGs) were identified at 7th post immunization, 1st and 7th post challenge day, respectively. Differential genes were annotated into Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional databases. These genes are involved in many immune processes, mainly Toll-like receptor signaling pathways, antigen processing and presentation, and Th1 and Th2 cell differentiation. At 1st and 7th post challenge day, respectively, 73 and 141 immune-related DEGs were used to construct protein-protein interaction network, predicting 11 and 30 hub genes involved in the immune response. This study provides the basis for further elucidation of the immune protection against bacterial infection in flounder immunized with inactivated vaccine.
Materials and methods
Experiment fish
Healthy flounder (P. olivaceus) (length range: 15-17 cm) were purchased from a farm in Rizhao, Shandong Province, China. The experimental fish were kept in the laboratory basement for two weeks, during which the water temperature was maintained at 21 ± 1°C in aerated seawater. The fish were fed with commercial pellets at 3% of body weight per day, and 1/3 of the seawater in the tanks was replaced. Before the experiment, flounders were randomly tested as free of pathogenic bacteria (36). Fish were anesthetized with tricaine methylate (MS-222, Sigma, USA) before tissue sampling. The treatment of fish in this study was approved by the Institutional Animal Care and Use Committee of the Ocean University of China (permit number: 20150101).
Inactivated Edwardsiella tarda vaccine
E. tarda HC01090721 strain was isolated from the ascites of diseased flounder by researchers in our laboratory and stored in brain heart infusion (BHI) medium containing 15% glycerol at -80°C (37). The inactivated E. tarda vaccine was prepared by referring to the previous process (38). First, the conserved strains were inoculated on BHI solid medium by continuous scribing and cultured at 28°C for 24 h. Single colonies were inoculated on BHI liquid medium in oscillating incubator at 28°C and 180 rpm for expansion. After the bacteria had grown to the logarithmic growth phase, bacterial precipitates were obtained by centrifugation at 8000 g for 5 min and resuspended in sterilized 0.01 M phosphate-buffered saline (PBS; pH=7.4). The bacterial concentration was adjusted to 10× 109 CFU/mL and inactivated by the adding 0.5% formalin (V/V) and shaking at 4°C for 72 h. The precipitate was collected by centrifugation at 8000 g for 5 min and resuspended with sterile PBS. 200 µL of inactivated bacteria were coated in BHI solid medium and incubated for 24 h at 28°C. The inactivation was proved to be successful if no colonies grew. The bacterial concentration was adjusted to 2.0×109 CFU/mL and stored in a refrigerator at 4°C until use.
Vaccination, challenge and sampling
Healthy flounder were randomly divided into two groups, the PBS group was injected intraperitoneally with 100 μL of PBS with complete Freund’s adjuvant (1:1) and the inactivated vaccine group was injected intraperitoneally with 100 μL of inactivated E. tarda with complete Freund’s adjuvant (1:1). After 5 weeks of immunization, 1.0 × 106 CFU of E. tarda was inoculated for the challenge experiment. At 7th day post immunization and 1st and 7th day post challenge, respectively, spleens from three fish in each group were randomly sampled and snap-frozen in liquid nitrogen for RNA extraction. Samples were marked: CS-7 (CS-7-1/CS-7-2/CS-7-3) for fish on the seventh day after immunization with PBS, IPS-7 (IPS -7-1/IPS -7-2/IPS-7-3) for fish on the seventh day after immunization with inactivated vaccine, ACS-1 (ACS-1-1/ACS-1-2/ACS-1-3) for fish on the first day after five weeks of PBS immunization with E. tarda infection, AIPS-1 (AIPS-1-1/AIPS-1-2/AIPS-1-3) for fish on the first day after five weeks of vaccine immunization with E. tarda infection, ACS-7 (ACS-7-1/ACS-7-2/ACS-7-3) for fish on the seventh day after five weeks of PBS immunization with E. tarda infection, AIPS-7 (AIPS-7-1/AIPS-7-2/AIPS-7-3) for fish on the seventh day after five weeks of vaccine immunization with E. tarda infection. To detect the infection status of sampled fish, remaining spleen tissue after challenge was immersed in RNALater (TaKaRa) for detection of bacterial load and also fully embedded in Tissue-Tek O.C.T. Compound (Sakura Finetek USA) and stored at -80°C for 4 h for immunofluorescence analysis. The specific experimental procedure is shown in Figure 1.
Figure 1 Experimental procedure of flounder after immunization with inactivated vaccine and subsequent E. tarda challenge.
RNA sample preparation, library construction and sequencing
For transcriptome analysis, 18 RNA-seq libraries were constructed using spleens. RNA was extracted from spleen using the Trizol according to the manufacturer’s instructions for RNA-seq analysis and subsequent RT-qPCR validation of transcriptomic data. RNA samples from spleens were subjected to RNase free agarose gel electrophoresis and Agilent 2100 Bioanalyzer for quality and integrity. The mRNA was enriched by Oligo(dT) beads and fragmented with fragmentation buffer. The mRNA was reverse transcribed into cDNA using the NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). cDNA was purified using 1.8X Agencourt AMPure XP Beads, followed by end repair, the addition of base A, and ligation of sequencing adapters. Size selection of ligated fragments was performed by agarose gel electrophoresis and polymerase chain reaction (PCR). Sequencing of 18 cDNA libraries from spleen samples was performed on the Illumina Novaseq6000 platform at Gene Denovo Biotechnology Co. (Guangzhou, China).
Data quality control, sequence alignment analysis and expression statistics
Raw reads of the spleen cDNA libraries contained low quality bases. The reads containing more than 10% unknown nucleotides and containing more than 50% low quality (Q-value ≤ 20) bases and containing adapter are removed according to FASTP (version 0.18.0) (39) to obtain high quality clean reads from the spleen cDNA libraries. The clean reads were aligned to the ribosomal RNA (rRNA) database using the short reads alignment tool bowtie2 (40), and the aligned ribosomal reads in the spleen cDNA libraries were removed. The retained clean reads were used for subsequent transcriptome analysis. The clean reads were mapped to the flounder reference genome (GenBank project accession: PRJNA344006) using HISAT2.2.4 (41). The mapped reads from 18 spleen cDNA libraries were assembled respectively by using StringTie v1.3.1 (42, 43). For each transcribed region, FPKM (fragment per kilobase of transcript per million mapped reads) values were calculated using RSEM (44) to quantify its expression abundance and variation.
Analysis of differentially expressed genes, GO and KEGG functional enrichment
Splenic differentially expressed genes (DEGs) between the PBS and inactivated vaccine groups were analyzed using DESeq2 (45), and the screening criteria for differential genes were p ≤ 0.05 and expression fold change > 1.5. To further understand the biological functions of the genes, DEGs were mapped to terms in the GO database (http://www.geneontology.org/), and the number of genes per term was calculated. GO terms that were significantly enriched to differential genes compared to background genes were identified. KEGG (http://www.genome.jp/kegg/) is the main database on Pathway. Pathways that were significantly enriched by differential genes compared to background genes were obtained. A hypergeometric test was performed to identify significantly enriched GO terms and KEGG pathways using p-value ≤ 0.05 as the threshold.
Protein-protein interaction network analysis
Venn diagram was used to demonstrate shared and unique DEGs for the three comparison groups after immunization and challenge. String (http://string-db.org) was applied to analyze protein-protein interactions that upregulate immune-related DEGs at 1st and 7th after challenge day, respectively. Cystoscape (V3.7.1) (https://cytoscape.org/) was used to visualize network file in which genes as nodes and interaction relationships as lines of the network. Studying the protein-protein interaction network helps to identify hub genes.
Quantitative real-time polymerase chain reaction
To analyze the bacterial load in the spleen of flounder post E. tarda infection, spleen DNA was extracted using the TIANamp Marine Animals DNA Kit (Tiangen, Beijing, China) according to the manufacturer’s instructions. Specific primers (F: TAGGGAGGAAGGTGTGAA; R: CTCTAGCTTGCCAGTCTT) were used for the amplification of E. tarda gene fragments. Each sample was taken in triplicate. The bacterial load in the spleen (log10 E. tarda cells/0.1 g) was quantified according to a previously established standard curve (46).
To verify the reliability of the transcriptome data, immune-related genes were selected for qRT-PCR. RNA was adjusted to the concentration of 1 μg/μL using Nanodrop 8000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). cDNA was obtained using HiScript III RT SuperMix for qPCR (Vazyme, Nanjing, China) for RNA reverse transcription. Gene specific primers were designed using Primer Premier 5.0 and the primer sequences are listed in Table 1. β-actin was used as internal control. The mRNA levels of spleens in the PBS group at each time were set to 1. qRT-PCR was performed using the LightCycler® 480 II Real Time System (Roche, Basel, Switzerland). Each reaction system contained 10 μL of 2× Universal SYBR Green Fast qPCR Mix, 2 μL of cDNA template, 0.4 μL each of forward and reverse primers, and 7.2 μL of DEPC water. The reaction procedure was: pre-denaturation at 95°C for 3 min, 40 cycles including annealing at 95°C for 5 s and extension at 60°C for 30 s. Reactions were performed in triplicate. Gene expression levels were analyzed by the 2-ΔΔCt method.
Table 1 List of primers for qRT-PCR validation of differentially expressed genes.
Indirect immunofluorescence assay
To analyze the infection of the spleen by E. tarda as described previously (47), in brief, sections obtained using Cryostats (Leica CM1900) were fixed in cold acetone for 15 min. Sections were washed three times with phosphate-buffered saline containing 0.05% Tween-20 (PBST) and then sealed with 5% BSA at 37°C for 45 min. Sections were then washed three times with PBST for 5 min each and incubated with rabbit anti-E. tarda polyclonal antibody (1:2000) as primary antibody for 1 h. After washing away the unconjugated primary antibody, sections were incubated with Alexa Flour 488-conjugated goat anti-rabbit IgG (1:1000, Thermo Fisher Scientific, Waltham, MA, USA) for 45 min in the dark at 37 °C. After three washes with PBST, the sections were counterstained with DAPI on the nuclei and incubated for 15 min at room temperature in the dark. After blocking with glycerol, the green positive signal was observed under fluorescent microscope (Olympus DP70, Tokyo, Japan). Rabbit negative serum was used as a negative control (1:1000 diluted in PBS).
Statistical analysis
All experiments were performed three times. Data were presented as mean ± standard deviation (SD). Data analysis was performed by SPSS 20.0 (IBM, Armonk, NY, USA). T-test was used to examine the differences between the PBS and inactivated vaccine groups. Differences were considered statistically significant when *p < 0.05. Graphs were plotted using GraphPad Prism 9 (Inc. San Diego, CA, USA).
Results
The load of E. tarda in the spleen of flounder after challenge
The presence of E. tarda in spleen samples was characterized quantitatively by qRT-PCR and qualitatively by indirect immunofluorescence (Figure 2). At 1st day post challenge, the bacterial loads in the spleen of PBS and inactivated vaccine groups showed small differences (p > 0.05). At 7th day post challenge, the spleen bacterial loads were significantly higher in the PBS group than in the inactivated vaccine group (p < 0.05). A large amount of specific green fluorescence was present in PBS group, whereas a small amount of green fluorescence was present in the vaccine group.
Figure 2 Bacterial load in the spleen of flounder sampled from the PBS and inactivated vaccine groups was determined by qRT-PCR at 1st and 7th post challenge day (A). Values are shown as mean ± SD (N = 3). Asterisks on the bars represent statistical significance (*p < 0.05). Distribution of E tarda in the spleen (B) observed by indirect immunofluorescence (Bar = 20 um). The corresponding enlargement of the diagram in the dotted box is shown in the upper right of the image (Bar = 5 um). The green fluorescent represents signal of E tarda. Cell nuclei were stained with DAPI in blue.
Transcriptome sequencing and read mapping of flounder
At 7th day post immunization, 1st and 7th day post challenge, respectively, 18 libraries of spleen RNA from PBS and inactivated vaccine group flounders were sequenced. To ensure data quality, raw reads were quality controlled. The six groups produced an average of 43875712, 44863902, 43172840, 55465904, 40599382, and 52051965 clean reads, and the ratio of clean reads to raw reads was above 99%. Among these clean reads, the rates of Q20 and Q30 were above 97% and 92%, respectively, which indicated that the sequencing results were of good quality. Clean reads were mapped to the flounder reference genome at a rate of over 89% in all cases, with an average of 20,542 genes annotated to each library (Table 2). The raw transcriptome sequencing data were submitted to the Sequence Read Archive (SRA) in NCBI. The accession number is PRJNA870695.
Table 2 Summary of transcriptome data from splenic samples.
Analysis of differentially expressed genes after immunization
The distribution of DEGs between the PBS and inactivated vaccine groups is shown using the volcano plot (Figure 3A). At 7th day post immunization, there were 1020 genes upregulated and 402 genes downregulated (Figure 3D). Compared to the PBS group, interleukin-21 receptor-like (IL21R), interleukin-1 beta-like (IL1β), C-C motif chemokine 25-like (CCL25), C-C motif chemokine 20-like (CCL20) and other interleukin and chemokine-related genes were significantly upregulated. Cytoplasmic dynein 1 intermediate chain 1 (DYNC1I1), thrombospondin-4-B-like (THBS4B), thrombospondin-2 isoform X1 (THBS2), tubulin beta chain isoform X2 (TUBB), tubulin beta-4B chain-like (TUBB4B), cathepsin L1-like (CTSS), CD209 antigen-like protein E (CD209E) and other phagosome-related genes were significantly upregulated. Contactin-1 (CNTN1A), claudin-23-like (CLDN23), claudin-3-like (CLDN3), neuroligin-3 (NLGN3), neuroligin-4, X-linked (NLGN4X), contactin-associated protein 1 (CNTNAP1) and other cell adhesion molecules-related genes were significantly upregulated. Phosphatidylinositol 4,5-bisphosphate 3-kinase catalytic subunit alpha isoform-like (PIK3CA), phosphatidylinositol 3-kinase regulatory subunit beta-like (PIK3R2) and other phosphatidylinositol 3’-kinase-related genes were significantly upregulated. Double-stranded RNA-specific adenosine deaminase (ADAR), cyclic GMP-AMP synthase (MB21D1) and other cytosolic DNA-sensing-related genes were Significantly downregulated (Table 3).
Figure 3 Volcano plots analysis of significantly different genes (DEGs) between samples from the PBS and inactivated vaccine groups at 7th post immunization day (A) and 1st (B) and 7th (C) post challenge day. The red (expression up-regulated) and blue (expression down-regulated) points indicate differential expression of genes, and black points are no differences. Statistical plot of DEGs for the PBS and inactivated vaccine groups at three time points (D).
Table 3 Summary of immune-related genes in the inactivated vaccine group compared to the PBS group at 7th day post immunization.
GO functional analysis showed that DEGs were involved in biological process (BP), molecular function (MF) and cellular component (CC) after immunization (Figure 4). Top 20 GO terms of DEGs enrichment were shown according to p < 0.05, with some GO terms associated with immunization. DEGs were mainly enriched in receptor binding and actin binding of MF; cytoskeleton, extrinsic component of membrane and intermediate filament cytoskeleton of CC; biological adhesion and multicellular organismal process of BP (Figure 5A, Table 4).
Figure 4 Gene ontology (GO) enrichment analysis of DEGs in the inactivated vaccine group compared to the PBS group at 7th post immunization day and 1st and 7th post challenge day. 3-level GO annotations are distributed in three categories (biological process, molecular function and cellular component). Yellow represents up-regulated expression and blue represents down-regulated expression.
Figure 5 GO enrichment circle plot of Top 20 GO terms at 7th post immunization day (A) and 1st (B) and 7th (C) post challenge day (p < 0.05). The first circle is the ID of the GO term, yellow represents the biological process, pink represents the molecular function, and green represents the cellular component. The second circle is the number of the GO term in the background genes and the p value. The third circle shows the proportion of up-regulated (deep purple) and down-regulated (blue) genes. The fourth circle shows the Rich Factor value of each GO term (Rich Factor value is the number of differential genes divided by the number of background genes in the GO term).
Table 4 ID, description, class and specific P-value of Top 20 GO terms at 7th post immunization day.
DEGs were enriched to six branches in KEGG (Figure 6), and bubble plots were used to demonstrate the top 20 signaling pathways (p < 0.05). At 7th post immunization day, DEGs were significantly enriched to Phagosome, Cell adhesion molecules (CAMs), PI3K-Akt signaling pathway and NF-kappa B signaling pathway (Figure 7A).
Figure 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs at 7th post immunization day and 1st and 7th post challenge day. The enriched pathways were counted in terms of level 1 and 2 classifications.
Figure 7 KEGG enrichment bubble plot of Top 20 pathways between PBS and inactivated vaccine groups at 7th post immunization day (A) and 1st (B) and 7th (C) post challenge day (p < 0.05).
Early transcriptomic responses in the spleen of flounder challenged with E. tarda
On the first day after five weeks of immunization with E. tarda infection, 1210 DEGs (620 up-regulated genes and 590 down-regulated genes) were identified in the inactivated vaccine group compared with PBS group (Figures 3B, D). Toll-like receptor 3 (TLR3), toll-like receptor 7 (TLR7), toll-like receptor 8 (TLR8), toll/interleukin-1 receptor domain-containing adapter protein isoform X1 (TIRAP) and other toll-like receptor-related genes were significantly upregulated. Junctional adhesion molecule B-like (JAM2), cadherin-4-like isoform X1 (CDH4), neural cell adhesion molecule 1-like isoform X1 (NCAM1), contactin-1a-like (CNTN1A) and other cell adhesion molecules-related genes were significantly upregulated. Inhibin beta B chain-like (INHBB), C-X-C chemokine receptor type 4-like (CXCR4), interleukin-21 receptor-like (IL21R), bone morphogenetic protein receptor type-1B isoform X1 (BMPR1B) and other Cytokine-cytokine receptor-related genes were significantly upregulated. Rano class II histocompatibility antigen, A beta chain-like (RT1-B), H-2 class II histocompatibility antigen, A-U alpha chain-like (RT1-Ba) and other Antigen processing and presentation-related genes were significantly upregulated. DNA replication ATP-dependent helicase/nuclease DNA2 isoform X1 (DNA2), replication factor C subunit 5 (RFC5), replication factor C subunit 2 (RFC2) and other cell cycle-related genes were significantly downregulated. G2 and S phase-expressed protein 1 isoform X1 (GTSE1), serine/threonine-protein kinase Chk1 (CHEK1), cyclin-G2-like (CCNG2) and other p53-related genes were significantly downregulated (Table 5).
Table 5 Summary of immune-related genes in the inactivated vaccine group compared to the PBS group at 1st post challenge day.
GO functional analysis showed that at 1st day post challenge, DEGs were mainly enriched to motor activity and tubulin binding of MF; cytoskeleton and microtubule cytoskeleton of CC; DNA metabolic process, cell cycle, and cell cycle process of BP (p < 0.05) (Figure 5B, Table 6).
Table 6 ID, description, class and specific P-value of Top 20 GO terms at 1st post challenge day.
KEGG functional analysis showed that DEGs were significantly enriched to Antigen processing and presentation, Cell adhesion molecules (CAMs), p53 signaling pathway (p < 0.05) (Figure 7B).
Late transcriptomic responses in the spleen of flounder challenged with E. tarda
On the seventh day after five weeks of immunization with E. tarda infection, compared with the PBS group, 1106 genes were up-regulated and 823 genes were down-regulated (Figures 3C, D).
Toll-like receptor 7 (TLR7), signal transducer and activator of transcription 1-alpha/beta-like isoform X1 (STAT1), toll-like receptor 5 (TLR5) and other Toll-like receptor-related genes were significantly upregulated. Interleukin-12 receptor subunit beta-2-like (ILL2RB2), interleukin-12 subunit beta-like (IL12B), protein jagged-2-like (JAG2), interleukin-6 (IL6), interleukin-6 receptor subunit alpha-like (IL6R), heat shock protein HSP 90-alpha (HSP90A.1) and other T cell differentiation-related genes were significantly upregulated. Fibroblast growth factor 1 (FGF1), lysophosphatidic acid receptor 3 (LPAR3), lysophosphatidic acid receptor 4 (LPAR4), laminin subunit beta-1(LAMB1), integrin beta-4 isoform X1 (ITGB4) and other PI3K-Akt signaling pathway-related genes were significantly upregulated. Claudin-3-like (CLDN3), E-selectin (SELE), L-selectin-like (SELL) and other Cell adhesion molecules-related genes were significantly upregulated. Mitogen-activated protein kinase 14A-like isoform X2 (MAPK14A), mitogen-activated protein kinase 8-like isoform X1 (MAPK8) and other MAPK-related genes were significantly downregulated (Table 7).
Table 7 Summary of immune-related genes in the inactivated vaccine group compared to the PBS group at 7th post challenge day.
GO functional analysis showed that DEGs were mainly enriched to NADH dehydrogenase activity, G-protein coupled nucleotide receptor activity and iron ion binding of MF, glycerol-3-phosphate metabolic process and endothelial cell differentiation of BP (Figure 5C, Table 8).
Table 8 ID, description, class and specific P-value of Top 20 GO terms at 7th post challenge day.
KEGG functional analysis showed that DEGs were significantly enriched to Th1 and Th2 cell differentiation, Th17 cell differentiation, IL-17 signaling pathway, Antigen processing and presentation, Toll-like receptor signaling pathway and RIG-I-like receptor signaling pathway, which are closely related to immunity (Figure 7C).
Compared with immunization, the protection of vaccine on flounder after challenge
To explore the protection of vaccine after challenge compared to immunization. Venn diagram was used to show the common and specific profiles of upregulated DEGs at three time points. A specific 512 DEGs were upregulated at 1st day post challenge compared to 7th day post immunization (Figure 8B). Immune-related 73 DEGs were used to construct protein-protein interaction networks (Table S1). Their expression levels at three time points are represented by heat map (Figure 8D). 57 DEGs showed interactive network relationships. 11 hub genes (TLR7, TLR3, CXCR4, TLR8 etc) were identified in the network according to the number of node connections, most of which belong to the Toll-like receptor signaling pathway (Figure 9, Table 9). A specific 1042 DEGs were upregulated at 7th day post challenge compared to 7th day post immunization (Figure 8C). Immune-related 141 DEGs were used to construct protein-protein interaction networks (Table S2). Their expression levels at three time points are also represented by heat map (Figure 8E). 127 DEGs showed interactive network relationships. 30 hub genes (IL6, STAT1, HSP90A.1, TLR7, IL12B etc) were identified in the network according to the number of node connections (Figure 10, Table 10). In addition, a total of 9 DEGs were upregulated at the three time points. The immune-related DEGs were IL21R, Lpar4 (Figure 8A).
Figure 8 Venn diagram of upregulated genes between the PBS and inactivated vaccine groups at the three time points (A). Venn diagram of genes upregulated between the seventh day after immunization and the first day after challenge (B) and heat map of 73 immune-related genes specifically upregulated on the first day after challenge between the PBS and inactivated vaccine groups at the three time points (D). Venn diagram of genes upregulated between day 7 after immunization and day 7 after challenge (C) and heat map of 141 immune-associated genes specifically upregulated on day 7 after challenge in the PBS and inactivated vaccine groups at three time points (E). The color shades in the heat map represent the gene expression levels. Closer to red indicates greater expression and closer to blue less expression.
Figure 9 Protein interaction network analysis of immune-related upregulated genes on the first day after challenge.
Table 9 Description and degree of connectivity of the 11 immune-associated upregulated hub genes on the first day after challenge.
Figure 10 Protein interaction network analysis of immune-related upregulated genes on day 7 after challenge.
Table 10 Description and degree of connectivity of the 30 immune-associated upregulated hub genes on day 7 after challenge.
qRT-PCR validation of transcriptomic immune-related genes
The expression of ten immune-related genes at each time point was randomly examined using qRT-PCR. At 7th day post immunization, the expression of IL1β, CXCL12, IL21R, RT1-B and LPAR4 were upregulated, while the expression of TLR7, CCL19, IRF5 and DHX58 were down-regulated (Figure 11A). At 1st day post challenge, the expressions of TLR7, IL21R, MAPK8, CXCL14, LPAR4 and TLR8 were up-regulated, while the expressions of GADD45β, PPAP2B, CCL25 and CCL20 expressions were down-regulated (Figure 11B). At 7th day post challenge, the expressions of IRF5, HSP90α.1, IL21R, TLR7, STAT1 and LPAR4 were up-regulated, while the expressions of IL5Rα, CXCL8, MAPK14A and JUN were down-regulated (Figure 11C). The expression levels of immune-related genes were basically consistent with the transcriptome results at the three time points, indicating that the sequencing results were accurate and reliable.
Figure 11 Comparison of gene expression levels between quantitative real-time polymerase chain reaction (qRT-PCR) and transcriptome sequencing at 7th post immunization day (A) and 1st (B) and 7th (C) post challenge day.
Discussion
Immune response induced by vaccines is essential factor in protecting fish from pathogens (12, 48). In previous studies, flounder immunization with inactivated vaccines was effective against the challenge of E. tarda in terms of RPS, innate and specific immunity (13). The transcriptome at the genetic level can provide rich information on the mechanisms of the immune response in fish after vaccination and the status of the pathogen after infection. The spleen is the secondary immune organ with an important role in hematopoiesis and immunity. Here, flounder spleens were collected for RNA-Seq at 7th day post immunization, 1st and 7th day post E. tarda challenge, respectively. DEGs were functionally annotated to explore immune response.
Vaccines trigger the body to produce the immune response in order to provide protection in case of pathogenic challenge. When vaccination is administered, the immune system makes a complex series of responses (3, 5, 49). First, many cytokines and chemokines are induced, causing an inflammatory response (50, 51). On the seventh day of marbled sleepy goby vaccination with inactivated iridovirus and rhabdovirus bivalent vaccine, pro-inflammatory cytokines such as IL-1β, IFN-γ, and IL-2R are upregulated to activate the innate immune response and, more importantly, to trigger a specific immune response (52). In study of IL-1β, IL-8, TNF-α and G-CSF as adjuvants for OmpV of E. tarda subunit vaccine, IL-1β and IL-8 were reported to significantly enhance serum antibodies and sIg+B lymphocytes, and the expression of genes (CD4-1, CD4-2, MHCIIα and IgM) related to cellular and humoral immunity (53, 54). Similarly, in this study, inflammatory and chemotaxis-related genes such as IL21R, IL1b, CCL25, and CCL20 were significantly upregulated at day 7 after immunization, activating cell recruitment and laying the foundation for triggering adaptive immunity. Phagocytosis is the process of host defense against pathogens. In innate immunity, macrophages, dendritic cells and neutrophils take up antigens through endocytosis. The antigen is then digested by lysosomes and the antigen signal is presented to the specific immune system (55). In adaptive immunity, B lymphocytes have been demonstrated to have phagocytic effect in fish. A previous study showed that B lymphocytes of dental flounder could phagocytose inactivated Lactococcus lactis (L. lactis). Transcriptome sequencing analysis of B lymphocytes after L. lactis stimulation showed that many DEGs were enriched to the phagocytic pathway. Further studies revealed the key role of Fc receptor (FcR) in regulating phagocytosis and bactericidal activity of B lymphocytes (56, 57). After 2 and 4 weeks of turbot (Scophthalmus maximus) inoculation with the bivalent inactivated bacteria vaccine Aeromonas salmonicida and Vibrio scophthalmi, Phage-associated genes such as Calreticulin (CALR), Antigen peptide transporter 1 (TAP1), and Integrin beta-3 (αVβ3) C-type mannose receptor 2 (MRC2) were upregulated in the kidney (33). Integrin beta-5-like (ITGB5), cyptoplasmic dynein 1 intermediate chain 1 (DYNCLI1), and thrombospondin-4-B-like (THBS4B) of the Phagosome pathway were also significantly enriched in this study, indicating that the early immune response underwent antigen processing. Cell adhesion molecules are involved in recognition between cells by means of ligand and receptor binding. In addition, in the immune response, it transmits signals for the antigen delivery process (58). In muscles around the injection site of flounder (Paralichthys olivaceus) vaccinated with VAA DNA vaccine, cell adhesion molecules enhance the local immune response by mediating the recruitment of immune cells to the site of inflammation (59). In addition, the PI3K-Akt signaling pathway is important node for signaling. It promotes cell proliferation, differentiation and anti-apoptosis after receiving extracellular signals (60). PIK3CD and PIK3R2 were identified as hub genes involved in the immune response in flounder gill infected with E. tarda. It was found that the expression of PIK3CD decreased continuously, while the expression of PIK3R2 increased and then decreased during infection. Combined with its transduction of antigenic signals, it affects the specific immune response, which is inhibited by E. tarda (61). In the present study, genes related to CAM (CNTN1A, CLDN23, CLDN3, NLGN3, NLGN4X) and PI3K family (PIK3R2, PIK3CA) were significantly upregulated after immunization.
At 1st day post challenge, pattern recognition receptors (PRRs), such as Toll-like receptors, elicit the body’s immune response by binding to pathogen-associated molecular patterns (PAMPs) (62). Eleven TLR family members were identified in flounder (63). TLR3, TLR7, TLR8 and TIRAP were significantly upregulated in this study. In mammals, TLR3, TLR7 and TLR8, located in intracellular vesicles, are critical receptors for the recognition of viral nucleic acids in the antiviral response (64). TLR7 was upregulated in both head kidney and spleen within 48 h after tongue sole infection with Pseudomonas fluorescence, and knockdown of TLR7 resulted in significantly higher bacterial load in tissues than in control (65). In the intestine of black rockfish Sebastes schlegelii infected with E. tarda, TLR3 showed upregulation at 2h, 6h, 12h and 24h (66). The TLR family of bony fish has a more complex immune response to pathogens such as viruses and bacteria than that of mammals. In the protective response of vaccines, antigen delivery is the initiation of adaptive immunity. Antigen signals are presented to T cells, which can specifically bind and kill target cells or release cytokines to stimulate B cell proliferation and differentiation (67). Among them, MHC molecules expressed on antigen-presenting cells such as dendritic cells and B cells are the markers of presentation (68). MHC class I presents antigenic fragments to CD8+ T cells for killing of target cells by releasing cytotoxic particles (perforin and granzyme) (69, 70). MHC class II presents antigenic fragments to CD4+ T cells to achieve expanded and increased immune effector functions through synthesis and release of cytokines (71). The expression of MHC Iα, MHC IIα, CD4-1 and CD8α was significantly upregulated in immunized tissues when flounder was immunized with inactivated E. tarda vaccine. In addition, these genes were also significantly elevated in the spleen and head kidney after five weeks of immunization with E. tarda infection. This indicates that immunization activates both cellular and humoral immune responses (13, 47). Vaccinated Arctic Charr showed significant expression of TLR7 after infection with Aeromonas salmonicida, activating B cells and DC to produce IFN-α and triggering Th1 and CD8+ T cell responses to demonstrate the effect of vaccination (72). Antigen processing and presentation-related genes (TAP1, TAP2, ABCB9 and PSME2) were identified in flounder spleen erythrocytes infected with E. tarda for 24h (32). Full-length transcriptome sequencing was performed on several tissues (liver, kidney, intestine, skin, gill) involved in immune and metabolic processes in black rockfish (Sebastes schlegelii). Four immune-related genes annotated as H-2 class II histocompatibility antigens were mined (73). In the present study, RT1-B and RT1-Ba were identified. In addition, Cytokine-cytokine receptor interaction as well as adhesion molecules (CAMs) remained functional after E. tarda infection. Toll-like receptors form the initial barrier against pathogens by specifically recognizing pathogens. In addition, TLR signaling induces DCs to produce IL-1β, IL-6, IL-12, and chemotactic cytokines that regulate antigen-specific Th1 and Th2 cell differentiation, linking innate and adaptive immunity (74–77). Marbled sleepy goby vaccinated with inactivated iridovirus and rhabdovirus bivalent vaccine showed consistent upregulation of MHC I, CD8, TCR, MHC II, CD4 and IgM expression after 2 days of Oxyeleotris marmoratus rhabdovirus challenge, indicating rapid induction of cellular and humoral immunity (52). At 7th day post challenge, TLR5, TLR7, il6, il1b, and IL12B were still highly expressed, suggesting that TLR may play a role in both early and late stages of infection. In mammals, stimulated by different antigens, CD4+ T lymphocytes differentiate into different cell subtypes (Th1, Th2, Th17 and Treg cells). This mechanism has also been demonstrated in teleost fish (78). Flounder were immunized with the NADP-dependent isocitrate dehydrogenase (IDH) subunit vaccine of E. tarda. The expression of Th1 and Th2 immune-related genes (IL-1β, TNF-α, IL-8, IL-6, NKEF, IFN-γ) was significantly increased. After infection with E. tarda, the bacterial load in the tissues was significantly reduced and the RPS reached 73.3%, which provided good protection against edwardsiellosis (46). In the present study, Th1 and Th2 cell differentiation, Th17 cell differentiation were significantly enriched, indicating that activation of cellular immunity plays an important role in vaccine protection against bacterial infection.
According to previous experiments, the immune response was initially activated at 7th day post immunization (36, 79). The response of the organism to the pathogen is strong and rapid after infection, with 24 h being important time point. 73 immune-related genes were activated on day 1 after E. tarda challenge compared to day 7 after immunization. 11 hub gene were identified, in which TLR family members were more enriched. TLR7 acts in organelles such as endoplasmic reticulum and lysosomes, which recognize viral single-stranded RNA and induce IFN-α, cytokine and chemokine production (80). In addition, the Toll-IL-1 receptor domain interacts with the junctional protein MyD88 in the antimicrobial immune response, activating the downstream NF-κB signaling pathway and producing pro-inflammatory cytokines. The bacterial load in the tissues of tongue sole (Cynoglossus semilaevis) infected with Pseudomonas fluorescence TSS was significantly enhanced after knockdown of CsTLR7. The results suggest that CsTLR7 has a positive role in the clearance of bacterial pathogens (65). In mammals, TLR3 specifically recognizes viral double-stranded RNA (dsRNA) and triggers toll interleukin-I receptor domain (TIR) through the myeloid differentiation factor 88 (MyD88) non-dependent pathway, activating downstream type I interferon gene expression and the NF-κB signaling pathway to induce an antiviral response in the organism (81). High expression of TLR3 was detected in the spleen and head kidney of channel catfish infected with virulent Edwardsiella ictalurid (82). TLR3 showed upregulation in the intestine of black rockfish Sebastes schlegelii for 24 h after Edwardsiella tarda infection (66). In the present study, TLR3 expression was upregulated on the first day after infection, suggesting that TLR3 in fish also plays a role in resistance to bacterial infection. TLR8, which is highly homologous to TLR7, functions in the lysosome to recognize bacterial or viral single-stranded RNA (83). TLR8 expression was detected in mucosal tissues (skin, gill and intestine) of turbot after infection with Vibrio anguillarum and Streptococcus iniae (84). Based on their central position in the interaction network, TLRs may play the important function in the anti-infection response of the fish spleen.
141 immune-related genes were activated on day 7 after E. tarda challenge compared to day 7 after immunization. Inflammatory cytokines (IL6, IL12B, IL6R, IL12RB2, IL17D), transcription factors (STAT1, IRF3) showed the strong interaction according to the functional classification. IL-12, which is involved in Th1 differentiation, is produced by DCs. IL12B (p40) and IL-12A (p35) together encode IL12, which acts by inducing the production of IFN-γ. IL-12Rβ1 and IL-12Rβ2 form the IL-12R complex, which is mainly expressed by activated T cells and natural killer cell. IL-12 and IL-12R binding activates the JAK2/STAT4 pathway to increase IFN-γ production as well as induce shift of T cells to Th1 phenotype (85). Th1 cells secrete IFN-γ, TNF-β, IL-2, etc. to mediate cellular immunity, which effectively defends against infection by intracellular pathogens (86, 87). A previous study showed that T-bet is a transcription factor involved in the immune response of Th1 cells. The expression was significantly increased after E. tarda infection, which laterally corroborates the important role of Th1 cells in responding to pathogenic infections. Moreover, IFN-γ and IL-2 were able to upregulate T-bet expression and contribute to the differentiation of Th0 cells to Th1 cell type (88). On day 14 of immunization with inactivated iridovirus and rhabdovirus bivalent vaccine and on day 7 of Oxyeleotris marmoratus iridovirus infection, the expression of inflammatory cytokines such as IL-12 and IL-1β was upregulated in the spleen of marbled sleepy goby, realizing the inflammatory of the innate immune response and leading to the development of adaptive immunity (52). In the present study, IL-6 was the hub gene of the protective immune network at 7th post challenge day. Previous studies have shown that significant upregulation of IL-6 was also detected after immunization of flounder with rIDH vaccine (46). IL-6 is pro-inflammatory cytokine produced mainly by macrophages and Th2 cells. IL-6 forms a complex with IL-6R, which binds to the membrane protein gp130, which activates intracellular signal transduction to function. It can promote the proliferation activation of T cells and the expression of IL-2 receptor on the surface of T cells, which further assists the proliferation and differentiation of B cells and the production of antibodies to participate in the humoral immune response (89, 90). A previous study showed that cyclosporine A (CsA) inhibited T lymphocyte expression by blocking activation of the transcription factor NFAT, and then inhibited B lymphocyte expression and antibody production. This suggests that T lymphocytes have an important regulatory role on B lymphocytes in the immune response (91, 92). In response to cytokines (IL-6, IFNγ) and growth factors (epidermal growth factor), STAT1 forms dimers that are transported to the nucleus to regulate apoptosis and the cell cycle (93). After P. olivaceus was infected with E. tarda for 8h and 48h, STAT1 was identified as hub gene in blood, gill, and kidney expression profiling (35, 61, 94). In the present study, IL21R was significantly upregulated at all three time points, corroborating its protective role in immunity and infection. IL21R is cytokine receptor for IL21. It is expressed on activated NK cells and belongs to the type I cytokine receptor. The receptor binds to IL21, leading to the activation of several downstream molecules (JAK1, JAK3, STAT1 and STAT3), while inducing the proliferation and differentiation of T cells, B cells and natural killer (NK) cells (95, 96).
The spleen, as an important lymphoid organ, is the main site of the immune response. Previous studies have found strong immune response in the spleen after flounder immunization or infection (13, 47, 97, 98). In the study of flounder immunized with inactivated E. tarda vaccine, the uptake of antigen and antigen presentation-related immune genes (MHC Iα, MHC IIα, CD4-1 and CD8α etc) were significantly elevated in the spleen. The immunohistochemical results of the spleen showed that CD4+ and CD8+ T lymphocyte were distributed around the melanocyte macrophage center (MMC). This fully demonstrates the important role of the spleen in capturing antigens, aggregating macrophages and lymphocyte populations, and presenting antigens to lymphocytes to activate the adaptive immunity (47, 97). Also in the study of immunized flounder infected with E. tarda, the immune response in the spleen varied significantly at the tissue, T/B lymphocyte and genetic levels and was the main tissue for monitoring the protective effect (13, 98). Keeping in line with previous studies, in this work, the spleen was sampled for transcriptome sequencing analysis after immunization with inactivated vaccine and infection with E. tarda. In terms of immune and protective responses, the strong immune response (cytokines, T- and B-cell related factors, etc.) in the spleen was also confirmed. In addition, in the immune system of fish, there are temporal differences in the production of cytokines, T/B lymphocytes, and antibodies, which are closely related to immune protection (12, 13). The overall immune response requires a multitemporal analysis, which is important for a comprehensive and detailed elucidation of vaccine immunization mechanism (99, 100). For example, the humoral immune response and the production of antibodies require a longer response time to be effective. IgM was significantly upregulated in marbled sleepy goby immunized with inactivated iridovirus and rhabdovirus bivalent vaccine for two weeks, indicating that the humoral immune response was activated (52). In the vaccine evaluation, antibodies were produced in fish during the immunization phase of two weeks to four weeks (12, 101, 102). The study of multiple time points of immune response is also worth doing to fully reveal the mechanism of fish vaccines.
Conclusions
This study investigated the mechanism of immune protection in flounder inoculated with inactivated E. tarda vaccine. Transcriptome sequencing analysis of the flounder spleen was performed on the seventh day after immunization and on the first and seventh day after five weeks of immunization with E. tarda infection. 1422, 1210 and 1929 DEGs were identified, respectively, which were significantly enriched in immune-related pathways such as Toll-like receptor signaling pathway and Th1 and Th2 cell differentiation. In addition, immune-related hub genes were identified after E. tarda infection compared to after immunization, in which TLR family members (TLR3, TLR5, TLR7, TLR8), pro-inflammatory cytokines and their receptors (IL6, IL6R, IL12B, IL12RB2, IL17D) were the main regulators that exerted immune protection. These genes are closely associated with the recognition and presentation of pathogens and the activation of cellular immunity. This study analyzed the transcriptional profiles of flounder spleen after immunization and infection, providing basis for further elucidation of the immune protection mechanisms in flounder immunized with inactivated E. tarda vaccine.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The animal study was reviewed and approved by the Institutional Animal Care and Use Committee of the Ocean University of China (permit number: 20150101).
Author contributions
XW and JX contributed to the conception and design of this experiment, performed most of the experiments and statistical analysis, drafted and revised the manuscript. XT, HC, and XS participated in the design of the study, helped analyzed experiments and data. JX and WZ design the study, provided reagents, instruments and experiment space. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the National Natural Science Foundation of China (32173005; 31730101; 31672684; 31672685), Shandong Provincial Natural Science Foundation (ZR2020KC025), the National Key Research and Development Program of China (2018YFD0900503).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.1058599/full#supplementary-material
References
1. Weisberg SP, Ural BB, Farber DL. Tissue-specific immunity for a changing world. Cell (2021) 184:1517–29. doi: 10.1016/j.cell.2021.01.042
2. Castellino F, Galli G, Del Giudice G, Rappuoli R. Generating memory with vaccination. Eur J Immunol (2009) 39:2100–5. doi: 10.1002/eji.200939550
3. Munang’andu HM, Evensen Ø. Correlates of protective immunity for fish vaccines. Fish Shellfish Immunol (2019) 85:132–40. doi: 10.1016/j.fsi.2018.03.060
4. Dadar M, Dhama K, Vakharia VN, Hoseinifar SH, Karthik K, Tiwari R, et al. Advances in aquaculture vaccines against fish pathogens: Global status and current trends. Rev Fish Sci Aquac (2017) 25:184–217. doi: 10.1080/23308249.2016.1261277
5. Sommerset I, Krossøy B, Biering E, Frost P. Vaccines for fish in aquaculture. Expert Rev Vaccines (2005) 4:89–101. doi: 10.1586/14760584.4.1.89
6. Hwang JY, Kwon MG, Seo JS, Hwang SD, Jeong JM, Lee JH, et al. Current use and management of commercial fish vaccines in Korea. Fish Shellfish Immunol (2020) 102:20–7. doi: 10.1016/j.fsi.2020.04.004
7. Flores-Kossack C, Montero R, Köllner B, Maisey K. Chilean Aquaculture and the new challenges: Pathogens, immune response, vaccination and fish diversification. Fish Shellfish Immunol (2020) 98:52–67. doi: 10.1016/j.fsi.2019.12.093
8. Romstad AB, Reitan LJ, Midtlyng P, Gravningen K, Evensen Ø. Antibody responses correlate with antigen dose and in vivo protection for oil-adjuvanted, experimental furunculosis (Aeromonas salmonicida subsp. salmonicida) vaccines in Atlantic salmon (Salmo salar l.) and can be used for batch potency testing of vaccin. Vaccine (2013) 31:791–6. doi: 10.1016/j.vaccine.2012.11.069
9. Romstad AB, Reitan LJ, Midtlyng P, Gravningen K, Evensen Ø. Development of an antibody ELISA for potency testing of furunculosis (Aeromonas salmonicida subsp salmonicida) vaccines in Atlantic salmon (Salmo salar l). Biologicals (2012) 40:67–71. doi: 10.1016/j.biologicals.2011.09.011
10. Yamasaki M, Araki K, Nakanishi T, Nakayasu C, Yamamoto A. Role of CD4+ and CD8α+ T cells in protective immunity against edwardsiella tarda infection of ginbuna crucian carp, carassius auratus langsdorfii. Fish Shellfish Immunol (2014) 36:299–304. doi: 10.1016/j.fsi.2013.11.016
11. Soveral L de F, de Almeida PA, Kreutz Y, Ribeiro VA, Frandoloso R, Kreutz LC. Modulation of expression of proinflammatory genes and humoral immune response following immunization or infection with aeromonas hydrophila in silver catfish (Rhamdia quelen). Fish Shellfish Immunol Rep (2022) 3:100053. doi: 10.1016/j.fsirep.2022.100053
12. Zhang J, Hu Y, Sun Q, Li X, Sun L. An inactivated bivalent vaccine effectively protects turbot (Scophthalmus maximus) against vibrio anguillarum and vibrio harveyi infection. Aquaculture (2021) 544:737158. doi: 10.1016/j.aquaculture.2021.737158
13. Wu X, Xing J, Tang X, Sheng X, Chi H, Zhan W. Protective cellular and humoral immune responses to edwardsiella tarda in flounder (paralichthys olivaceus) immunized by an inactivated vaccine. Mol Immunol (2022) 149:77–86. doi: 10.1016/j.molimm.2022.06.008
14. Seikai T. Flounder culture and its challenges in Asia. Rev Fish Sci (2002) 10:421–32. doi: 10.1080/20026491051721
15. Mohanty BR, Sahoo PK. Edwardsiellosis in fish: A brief review. J Biosci (2007) 32:1331–44. doi: 10.1007/s12038-007-0143-8
16. Hu Yh, Zhou Hz, Jin Qw, Zhang J. The serine protease autotransporter tsh contributes to the virulence of edwardsiella tarda. Vet Microbiol (2016) 189:68–74. doi: 10.1016/j.vetmic.2016.04.021
17. Yu JE, Yoo AY, Choi KH, Cha J, Kwak I, Kang HY. Identification of antigenic edwardsiella tarda surface proteins and their role in pathogenesis. Fish Shellfish Immunol (2013) 34:673–82. doi: 10.1016/j.fsi.2012.11.019
18. Kramarz P, Lopalco PL, Huitric E, Pastore Celentano L. Vaccine-preventable diseases: The role of the European centre for disease prevention and control. Clin Microbiol Infect (2014) 20:2–6. doi: 10.1111/1469-0691.12430
19. Jansen KU, Anderson AS. The role of vaccines in fighting antimicrobial resistance (AMR). Hum Vaccines Immunother (2018) 14:2142–9. doi: 10.1080/21645515.2018.1476814
20. Dong Y, Dai T, Wei Y, Zhang L, Zheng M, Zhou F. A systematic review of SARS-CoV-2 vaccine candidates. Signal Transduct Target Ther (2020) 5:1–14. doi: 10.1038/s41392-020-00352-y
21. Wall DJ, Patel MM, Chung JR, Lee B, Dawood FS. Antibody response and protection after receipt of inactivated influenza vaccine: A systematic review. Pediatrics (2021) 147(6):e2020019901. doi: 10.1542/peds.2020-019901
22. Soulliere C, Dixon B. Immune system organs of bony fishes. Ref Modul Life Sci (2017) Chapter 12179. doi: 10.1016/B978-0-12-809633-8.12179-X
23. Liu Y, Zhang S, Jiang G, Yang D, Lian J, Yang Y. The development of the lymphoid organs of flounder, paralichthys olivaceus, from hatching to 13 months. Fish Shellfish Immunol (2004) 16:621–32. doi: 10.1016/j.fsi.2003.10.001
24. Aslam R, Kapur R, Segel GB, Guo L, Zufferey A, Ni H, et al. The spleen dictates platelet destruction, anti-platelet antibody production, and lymphocyte distribution patterns in a murine model of immune thrombocytopenia. Exp Hematol (2016) 44:924–930.e1. doi: 10.1016/j.exphem.2016.07.004
25. Raida MK, Buchmann K. Temperature-dependent expression of immune-relevant genes in rainbow trout following yersinia ruckeri vaccination. Dis Aquat Organ (2007) 77:41–52. doi: 10.3354/dao01808
26. Chu Y, Corey DR. RNA Sequencing: Platform selection, experimental design, and data interpretation. Nucleic Acid Ther (2012) 22:271–4. doi: 10.1089/nat.2012.0367
27. Oshlack A, Robinson MD, Young MD. From RNA-seq reads to differential expression results. Genome Biol (2010) 11:1–10. doi: 10.1186/gb-2010-11-12-220
28. Zhang Y, Lu L, Li C, Shao G, Chen X. Transcriptome analysis revealed multiple immune processes and energy metabolism pathways involved in the defense response of the large yellow croaker larimichthys crocea against pseudomonas plecoglossicida. Comp Biochem Physiol - Part D Genomics Proteomics (2021) 40:100886. doi: 10.1016/j.cbd.2021.100886
29. Wang H, Tang X, Sheng X, Xing J, Chi H, Zhan W. Transcriptome analysis reveals temperature-dependent early immune response in flounder (Paralichthys olivaceus) after hirame novirhabdovirus (HIRRV) infection. Fish Shellfish Immunol (2020) 107:367–78. doi: 10.1016/j.fsi.2020.10.020
30. Zhang J, Sun L. Transcriptome analysis reveals temperature-regulated antiviral response in turbot scophthalmus maximus. Fish Shellfish Immunol (2017) 68:359–67. doi: 10.1016/j.fsi.2017.07.038
31. Lu X, Chen HM, Qian XQ, Gui JF. Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish. Comp Biochem Physiol - Part D Genomics Proteomics (2020) 35:100688. doi: 10.1016/j.cbd.2020.100688
32. Sun B, Li X, Ning X, Sun L. Transcriptome analysis of paralichthys olivaceus erythrocytes reveals profound immune responses induced by edwardsiella tarda infection. Int J Mol Sci (2020) 21(9):3094. doi: 10.3390/ijms21093094
33. Xiu Y, Su L, Guo B, Yi J, Yang Z, Zhou S. Transcriptomic profiling reveals the protection mechanism of bivalent inactivated bacteria vaccine aeromonas salmonicida and vibrio scophthalmi in turbot (Scophthalmus maximus). Aquaculture (2022) 560:738479. doi: 10.1016/j.aquaculture.2022.738479
34. Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: Protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res (2015) 43:D447–52. doi: 10.1093/nar/gku1003
35. Liu X, Li Z, Wu W, Liu Y, Liu J, He Y, et al. Sequencing-based network analysis provides a core set of gene resource for understanding kidney immune response against edwardsiella tarda infection in Japanese flounder. Fish Shellfish Immunol (2017) 67:643–54. doi: 10.1016/j.fsi.2017.06.051
36. Gao YL, Tang XQ, Sheng XZ, Xing J, Zhan WB. Immune responses of flounder paralichthys olivaceus vaccinated by immersion of formalin-inactivated edwardsiella tarda following hyperosmotic treatment. Dis Aquat Organ (2015) 116:111–20. doi: 10.3354/dao02909
37. Tang X, Zhan W, Sheng X, Chi H. Immune response of Japanese flounder paralichthys olivaceus to outer membrane protein of edwardsiella tarda. Fish Shellfish Immunol (2010) 28:333–43. doi: 10.1016/J.FSI.2009.11.015
38. Gao Y, Tang X, Sheng X, Xing J, Zhan W. Antigen uptake and expression of antigen presentation-related immune genes in flounder (Paralichthys olivaceus) after vaccination with an inactivated edwardsiella tarda immersion vaccine, following hyperosmotic treatment. Fish Shellfish Immunol (2016) 55:274–80. doi: 10.1016/j.fsi.2016.05.042
39. Chen S, Zhou Y, Chen Y, Gu J. Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics (2018) 34:i884–90. doi: 10.1093/bioinformatics/bty560
40. Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods (2012) 9:357–9. doi: 10.1038/nmeth.1923
41. Kim D, Langmead B, Salzberg1 SL. HISAT: a fast spliced aligner with low memory requirements daehwan HHS public access. Nat Methods (2015) 12:357–60. doi: 10.1038/nmeth.3317.HISAT
42. Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and ballgown. Nat Protoc (2016) 11:1650–67. doi: 10.1038/nprot.2016.095
43. Pertea M, Pertea GM, Antonescu CM, Chang T-C, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol (2015) 33:290–5. doi: 10.1038/nbt.3122
44. Li B, Dewey CN. RSEM: Accurate transcript quantification from RNA-seq data with or without a reference genome. BMC Bioinf (2011) 12:323. doi: 10.1186/1471-2105-12-323
45. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15:1–21. doi: 10.1186/s13059-014-0550-8
46. Tang X, Liu F, Sheng X, Xing J, Zhan W. Recombinant NADP-dependent isocitrate dehydrogenase of edwardsiella tarda induces both Th1 and Th2 type immune responses and evokes protective efficacy against edwardsiellosis. Vaccine (2018) 36:2337–45. doi: 10.1016/j.vaccine.2018.03.022
47. Tang X, Gong J, Zeng C, Sheng X, Xing J, Zhan W. Dynamic distribution of formalin-inactivated edwardsiella tarda in flounder (Paralichthys olivaceus) post intraperitoneal vaccination. Fish Shellfish Immunol (2019) 89:393–402. doi: 10.1016/j.fsi.2019.04.022
48. Lauscher A, Krossøy B, Frost P, Grove S, König M, Bohlin J, et al. Immune responses in Atlantic salmon (Salmo salar) following protective vaccination against infectious salmon anemia (ISA) and subsequent ISA virus infection. Vaccine (2011) 29:6392–401. doi: 10.1016/j.vaccine.2011.04.074
49. Adams A. Progress, challenges and opportunities in fish vaccine development. Fish Shellfish Immunol (2019) 90:210–4. doi: 10.1016/J.FSI.2019.04.066
50. Ge H, Farris CM, Tong M, Maina A, Richards AL. Transcriptional profiles of cytokines and chemokines reveal important pro-inflammatory response from endothelial cells during orientia tsutsugamushi infection. Microbes Infect (2019) 21:313–20. doi: 10.1016/j.micinf.2019.01.002
51. Fang S, Xu T, Xiong M, Zhou X, Wang Y, Haydu LE, et al. Role of immune response, inflammation, and tumor immune response–related Cytokines/Chemokines in melanoma progression. J Invest Dermatol (2019) 139:2352–2358.e3. doi: 10.1016/j.jid.2019.03.1158
52. Guo X, Zhou Y, Fu X, Lin Q, Liu L, Liang H, et al. Transcriptomic profiles reveal that inactivated iridovirus and rhabdovirus bivalent vaccine elicits robust adaptive immune responses against lethal challenge in marbled sleepy goby. Fish Shellfish Immunol (2020) 98:429–37. doi: 10.1016/j.fsi.2020.01.047
53. Guo M, Tang X, Sheng X, Xing J, Zhan W. Comparative study of the adjuvant potential of four Th0 cytokines of flounder (Paralichthys olivaceus) on an e. tarda subunit vaccine. Dev Comp Immunol (2018) 86:147–55. doi: 10.1016/j.dci.2018.05.001
54. Guo M, Tang X, Sheng X, Xing J, Zhan W. The effects of IL-1β, IL-8, G-CSF and TNF-α as molecular adjuvant on the immune response to an e. tarda subunit vaccine in flounder (Paralichthys olivaceus). Fish Shellfish Immunol (2018) 77:374–84. doi: 10.1016/j.fsi.2018.04.009
55. Wu L, Li L, Gao A, Ye J, Li J. Antimicrobial roles of phagocytosis in teleost fish: Phagocytic b cells vs professional phagocytes. Aquac Fish (2022). doi: 10.1016/j.aaf.2021.12.008
56. Hao Y, Tang X, Xing J, Sheng X, Chi H, Zhan W. Regulatory role of fc receptor in mIgM+ b lymphocyte phagocytosis in flounder (Paralichthys olivaceus). Front Immunol (2021) 12:804244. doi: 10.3389/fimmu.2021.804244
57. Tang X, Yang S, Sheng X, Xing J, Zhan W. Transcriptome analysis of immune response of mIgM+ b lymphocytes in Japanese flounder (Paralichthys olivaceus) to lactococcus lactis in vitro revealed that IFN I-3 could enhance their phagocytosis. Front Immunol (2019) 10:1622. doi: 10.3389/fimmu.2019.01622
58. Liu Y, Zhang Y, Ding Y, Zhuang R. Platelet-mediated tumor metastasis mechanism and the role of cell adhesion molecules. Crit Rev Oncol Hematol (2021) 167:103502. doi: 10.1016/j.critrevonc.2021.103502
59. Xing J, Jiang X, Xu H, Sheng X, Tang X, Chi H, et al. Local immune responses to VAA DNA vaccine against listonella anguillarum in flounder (Paralichthys olivaceus). Mol Immunol (2021) 134:141–9. doi: 10.1016/j.molimm.2021.03.012
60. Juntilla MM, Koretzky GA. Critical roles of the PI3K/Akt signaling pathway in T cell development. Immunol Lett (2008) 116:104–10. doi: 10.1016/j.imlet.2007.12.008
61. Li Z, Liu X, Cheng J, He Y, Wang X, Wang Z, et al. Transcriptome profiling provides gene resources for understanding gill immune responses in Japanese flounder (Paralichthys olivaceus) challenged with edwardsiella tarda. Fish Shellfish Immunol (2018) 72:593–603. doi: 10.1016/j.fsi.2017.11.041
62. Magnadottir B. Immunological control of fish diseases. Mar Biotechnol (2010) 12:361–79. doi: 10.1007/s10126-010-9279-x
63. Zhu Y, Li S, Su B, Xue T, Cao M, Li C. Genome-wide identification, characterization, and expression of the toll-like receptors in Japanese flounder (Paralichthys olivaceus). Aquaculture (2021) 545:737127. doi: 10.1016/j.aquaculture.2021.737127
64. Kaur A, Baldwin J, Brar D, Salunke DB, Petrovsky N. ScienceDirect toll-like receptor ( TLR ) agonists as a driving force behind next-generation vaccine adjuvants and cancer therapeutics. Curr Opin Chem Biol (2022) 70:102172. doi: 10.1016/j.cbpa.2022.102172
65. Li XP, Sun L. TLR7 is required for optimal immune defense against bacterial infection in tongue sole (Cynoglossus semilaevis). Fish Shellfish Immunol (2015) 47:93–9. doi: 10.1016/j.fsi.2015.08.025
66. Cao M, Yan X, Yang N, Fu Q, Xue T, Zhao S, et al. Genome-wide characterization of toll-like receptors in black rockfish sebastes schlegelii: Evolution and response mechanisms following edwardsiella tarda infection. Int J Biol Macromol (2020) 164:949–62. doi: 10.1016/j.ijbiomac.2020.07.111
67. Horton R, Wilming L, Rand V, Lovering RC, Bruford EA, Khodiyar VK, et al. Gene map of the extended human MHC. Nat Rev Genet (2004) 5:889–99. doi: 10.1038/nrg1489
68. Abualrous ET, Sticht J, Freund C. Major histocompatibility complex (MHC) class I and class II proteins: impact of polymorphism on antigen presentation. Curr Opin Immunol (2021) 70:95–104. doi: 10.1016/j.coi.2021.04.009
69. Gao GF, Jakobsen BK. Molecular interactions of coreceptor CD8 and MHC class I: The molecular basis for functional coordination with the T-cell receptor. Immunol Today (2000) 21:630–6. doi: 10.1016/S0167-5699(00)01750-3
70. Corr M, Slanetz AE, Boyd LF, Jelonek MT, Khilko S, Al-ramadi BK, et al. T Cell Receptor-MHC Class I Peptide Interactions. Affinity, Kinetics, and Specificity Science (1994) 265:946–9. doi: 10.1126/science.8052850
71. Farber DL, Lugman M, Acuto O, Bottomly K. Control of memory CD4 T cell activation: MHC class II molecules on APCs and CD4 ligation inhibit memory but not naive CD4 T cells. Immunity (1995) 2:249–59. doi: 10.1016/1074-7613(95)90049-7
72. Braden LM, Whyte SK, Brown ABJ, Van Iderstine C, Letendre C, Groman D, et al. Vaccine-induced protection against furunculosis involves pre-emptive priming of humoral immunity in arctic charr. Front Immunol (2019) 10:120. doi: 10.3389/fimmu.2019.00120
73. Cao M, Zhang M, Yang N, Fu Q, Su B, Zhang X, et al. Full length transcriptome profiling reveals novel immune-related genes in black rockfish (Sebastes schlegelii). Fish Shellfish Immunol (2020) 106:1078–86. doi: 10.1016/j.fsi.2020.09.015
74. Iwasaki A, Medzhitov R. Toll-like receptor control of the adaptive immune responses. Nat Immunol (2004) 5:987–95. doi: 10.1038/ni1112
75. Schnare M, Barton GM, Holt AC, Takeda K, Akira S, Medzhitov R. Toll-like receptors control activation of adaptive immune responses. Nat Immunol (2001) 2:947–50. doi: 10.1038/ni712
76. Duan T, Du Y, Xing C, Wang HY, Wang RF. Toll-like receptor signaling and its role in cell-mediated immunity. Front Immunol (2022) 13:812774. doi: 10.3389/fimmu.2022.812774
77. Palliser D, Ploegh H, Boes M. Myeloid differentiation factor 88 is required for cross-priming In vivo. J Immunol (2004) 172:3415–21. doi: 10.4049/jimmunol.172.6.3415
78. Tian HF, Xing J, Tang XQ, Chi H, Sheng XZ, Zhan WB. Cluster of differentiation antigens: essential roles in the identification of teleost fish T lymphocytes. Mar Life Sci Technol (2022) 4:303–16. doi: 10.1007/s42995-022-00136-z
79. Du Y, Tang X, Sheng X, Xing J, Zhan W. Immune response of flounder (Paralichthys olivaceus) was associated with the concentration of inactivated edwardsiella tarda and immersion time. Vet Immunol Immunopathol (2015) 167:44–50. doi: 10.1016/j.vetimm.2015.07.002
80. Gorden KB, Gorski KS, Gibson SJ, Kedl RM, Kieper WC, Qiu X, et al. Synthetic TLR agonists reveal functional differences between human TLR7 and TLR8. J Immunol (2005) 174:1259–68. doi: 10.4049/jimmunol.174.3.1259
81. Chen Y, Lin J, Zhao Y, Ma X, Yi H. Toll-like receptor 3 (TLR3) regulation mechanisms and roles in antiviral innate immune responses. J Zhejiang Univ Sci B (2021) 22:609–32. doi: 10.1631/jzus.B2000808
82. Bilodeau AL, Waldbieser GC. Activation of TLR3 and TLR5 in channel catfish exposed to virulent edwardsiella ictaluri. Dev Comp Immunol (2005) 29:713–21. doi: 10.1016/j.dci.2004.12.002
83. Zhang S, Hu Z, Tanji H, Jiang S, Das N, Li J, et al. Small-molecule inhibition of TLR8 through stabilization of its resting state. Nat Chem Biol (2018) 14:58–64. doi: 10.1038/nchembio.2518
84. Dong X, Su B, Zhou S, Shang M, Yan H, Liu F, et al. Identification and expression analysis of toll-like receptor genes (TLR8 and TLR9) in mucosal tissues of turbot (Scophthalmus maximus l.) following bacterial challenge. Fish Shellfish Immunol (2016) 58:309–17. doi: 10.1016/j.fsi.2016.09.021
85. Schwinge D, von Haxthausen F, Quaas A, Carambia A, Otto B, Glaser F, et al. Dysfunction of hepatic regulatory T cells in experimental sclerosing cholangitis is related to IL-12 signaling. J Hepatol (2017) 66:798–805. doi: 10.1016/j.jhep.2016.12.001
86. Liu Q, Hua M, Zhang C, Wang R, Liu J, Yang X, et al. NLRP3-activated bone marrow dendritic cells play antileukemic roles via IL-1β/Th1/IFN-γ in acute myeloid leukemia. Cancer Lett (2021) 520:109–20. doi: 10.1016/j.canlet.2021.06.014
87. Liu W, Li M, Zhang X, Zhou Z, Shen Z, Shen X. Association of polymorphisms in Th1/Th2-related cytokines (IFN-γ, TGFβ1, IL-1β, IL-2, IL-4, IL-18) with oral lichen planus: A pooled analysis of case–control studies. J Dent Sci (2022) 18. doi: 10.1016/j.jds.2022.08.032
88. Tian H, Xing J, Tang X, Chi H, Sheng X, Zhan W. Identification and characterization of a master transcription factor of Th1 cells, T-bet, within flounder (Paralichthys olivaceus). Front Immunol (2021) 12:704324. doi: 10.3389/fimmu.2021.704324
89. Jin RM, Huang HZ, Zhou Y, Wang YY, Fu HC, Li Z, et al. Characterization of mandarin fish (Siniperca chuatsi) IL-6 and IL-6 signal transducer and the association between their SNPs and resistance to ISKNV disease. Fish Shellfish Immunol (2021) 113:139–47. doi: 10.1016/j.fsi.2021.04.003
90. Luo P, Wang Y, Zhao C, Guo J, Shi W, Ma H, et al. Bazedoxifene exhibits anti-inflammation and anti-atherosclerotic effects via inhibition of IL-6/IL-6R/STAT3 signaling. Eur J Pharmacol (2021) 893:173822. doi: 10.1016/j.ejphar.2020.173822
91. Xing J, Luo K, Xiao Y, Tang X, Zhan W. Influence of CD4-1+, CD4-2+ and CD8+ T lymphocytes subpopulations on the immune response of b lymphocytes in flounder (Paralichthys olivaceus) immunized with thymus-dependent or thymus-independent antigen. Fish Shellfish Immunol (2019) 84:979–86. doi: 10.1016/j.fsi.2018.11.004
92. Xing J, Xiao Y, Tang X, Sheng X, Zhan W. Inhibition of cyclosporine a or rapamycin on T lymphocyte counts and the influence on the immune responses of b lymphocytes in flounder (Paralichthys olivaceus). Fish Shellfish Immunol (2017) 66:78–85. doi: 10.1016/j.fsi.2017.05.017
93. O’Brown ZK, Van Nostrand EL, Higgins JP, Kim SK. The inflammatory transcription factors NFκB, STAT1 and STAT3 drive age-associated transcriptional changes in the human kidney. PloS Genet (2015) 11:1–28. doi: 10.1371/journal.pgen.1005734
94. Li Z, Liu X, Liu J, Zhang K, Yu H, He Y, et al. 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). Dev Comp Immunol (2018) 78:100–13. doi: 10.1016/j.dci.2017.09.013
95. Wang LN, Cui YX, Ruge F, Jiang WG. Interleukin 21 and its receptor play a role in proliferation, migration and invasion of breast cancer cells. Cancer Genomics Proteomics (2015) 12:211–22.
96. Australia R, Klein C, Leonard WJ, Wan C-K. IL-21 signaling in immunity [version 1; referees: 3 approved]. F1000research (2016) 5:1–10. doi: 10.12688/f1000research.7634.1
97. Xing J, Ma J, Tang X, Sheng X, Zhan W. Characterizations of CD4-1, CD4-2 and CD8β T cell subpopulations in peripheral blood leucocytes, spleen and head kidney of Japanese flounder (Paralichthys olivaceus). Mol Immunol (2017) 85:155–65. doi: 10.1016/J.MOLIMM.2017.02.015
98. Xing J, Zhang Z, Luo K, Tang X, Sheng X, Zhan W. T And b lymphocytes immune responses in flounder (Paralichthys olivaceus) induced by two forms of outer membrane protein K from vibrio anguillarum: Subunit vaccine and DNA vaccine. Mol Immunol (2020) 118:40–51. doi: 10.1016/j.molimm.2019.12.002
99. Rauta PR, Nayak B, Das S. Immune system and immune responses in fish and their role in comparative immunity study: A model for higher organisms. Immunol Lett (2012) 148:23–33. doi: 10.1016/j.imlet.2012.08.003
100. Secombes CJ, Wang T. The innate and adaptive immune system of fish. Infectious Disease in Aquaculture (2012) p. 3–68. doi: 10.1533/9780857095732.1.3
101. Feng L, Wang Q, Shan C, Yang C, Feng Y, Wu J, et al. An adenovirus-vectored COVID-19 vaccine confers protection from SARS-COV-2 challenge in rhesus macaques. Nat Commun (2020) 11:1–11. doi: 10.1038/s41467-020-18077-5
Keywords: immune protection, RNA-seq, spleen, Paralichthys olivaceus, Edwardsiella tarda
Citation: Wu X, Xing J, Tang X, Sheng X, Chi H and Zhan W (2022) Splenic protection network revealed by transcriptome analysis in inactivated vaccine-immunized flounder (Paralichthys olivaceus) against Edwardsiella tarda infection. Front. Immunol. 13:1058599. doi: 10.3389/fimmu.2022.1058599
Received: 30 September 2022; Accepted: 24 October 2022;
Published: 09 November 2022.
Edited by:
Hai-peng Liu, Xiamen University, ChinaCopyright © 2022 Wu, Xing, Tang, Sheng, Chi and Zhan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing Xing, xingjing@ouc.edu.cn; Wenbin Zhan, wbzhan@ouc.edu.cn