- 1Department of Ocean Sciences, Memorial University of Newfoundland, St. John’s, NL, Canada
- 2Centre for Marine Applied Research, Dartmouth, NS, Canada
- 3Aquatic Research Cluster, Core Research Equipment and Instrument Training (CREAIT) Network, Ocean Sciences Centre, Memorial University of Newfoundland, St. John’s, NL, Canada
- 4Faculty of Medicine, Memorial University of Newfoundland, St. John’s, NL, Canada
- 5Marine Microbial Pathogenesis and Vaccinology Laboratory, Department of Ocean Sciences, Memorial University of Newfoundland, St. John’s, NL, Canada
The economic importance of lumpfish (Cyclopterus lumpus) is increasing, but several aspects of its immune responses are not well understood. To discover genes and mechanisms involved in the lumpfish antiviral response, fish were intraperitoneally injected with either the viral mimic polyinosinic:polycytidylic acid [poly(I:C)] or phosphate-buffered saline (PBS; vehicle control), and head kidneys were sampled 24 hours post-injection (hpi) for transcriptomic analyses. RNA sequencing (RNA-Seq) (adjusted p-value <0.05) identified 4,499 upregulated and 3,952 downregulated transcripts in the poly(I:C)-injected fish compared to the PBS-injected fish. Eighteen genes identified as differentially expressed by RNA-Seq were included in a qPCR study that confirmed the upregulation of genes encoding proteins with antiviral immune response functions (e.g., rsad2) and the downregulation of genes (e.g., jarid2b) with potential cellular process functions. In addition, transcript expression levels of 12 members of the interferon regulatory factor (IRF) family [seven of which were identified as poly(I:C)-responsive in this RNA-Seq study] were analyzed using qPCR. Levels of irf1a, irf1b, irf2, irf3, irf4b, irf7, irf8, irf9, and irf10 were significantly higher and levels of irf4a and irf5 were significantly lower in the poly(I:C)-injected fish compared to the PBS-injected fish. This research and associated new genomic resources enhance our understanding of the genes and molecular mechanisms underlying the lumpfish response to viral mimic stimulation and help identify possible therapeutic targets and biomarkers for viral infections in this species.
Introduction
Lumpfish (Cyclopterus lumpus) are commonly used as an environmentally friendly solution for sea lice control (e.g., Lepeophtheirus salmonis) in Atlantic salmon (Salmo salar) farms in the North Atlantic region (1–3). Sea lice infestations lead to decreased fish health, growth, and, consequently, market value (4). Additionally, lumpfish as a biological method to control sea lice help reduce the reliance on chemical treatments, which contribute to environmental pollution. However, lumpfish farming faces several challenges. For example, lumpfish are susceptible to several infectious diseases, which can be transferred to other aquatic hosts such as Atlantic salmon (2, 5). Atlantic salmon and lumpfish can both be infected with pathogens such as the bacterium Renibacterium salmoninarum and the viral hemorrhagic septicemia (VHS) virus (2, 6). While lumpfish may be infected by viral pathogens (e.g., Supplementary Table S1), the development of vaccines for farmed lumpfish has thus far focused on bacterial pathogens (7–12) rather than viruses. The development of vaccines for the protection of lumpfish against viral infection is a priority (7).
Polyinosinic:polycytidylic acid [poly(I:C)] is a synthetic analog of double-stranded RNA (dsRNA) that can mimic viral infections (i.e., elicit a potent antiviral-like response) in several species including teleosts and is commonly used as an immune stimulant in aquaculture research (13–15). Poly(I:C) was previously used in several studies to evaluate the antiviral response of zebrafish (Danio rerio), Chinook salmon (Oncorhynchus tshawytscha), seven-band grouper (Epinephelus septemfasciatus), and Atlantic salmon, as it mimics RNA viral pathogens of fish (15–19). For example, positive-sense single-stranded RNA (ssRNA) viruses such as nodaviruses and flaviviruses form dsRNA intermediates during their replication cycle (20–22). Also, it has been reported that lumpfish can be infected with dsRNA viruses (Supplementary Table S1), such as C. lumpus toti-like virus (CLuTLV) (23), and ssRNA viruses that likely produce dsRNA during their replication cycle including C. lumpus virus (CLuV) (24) and nervous necrosis virus (NNV) (25).
RNA sequencing (RNA-Seq) is a highly robust method for assessing global gene expression responses, as it generates accurate and reproducible results (26). Transcriptomic studies may be used to identify genes and pathways that respond to immune challenges. As examples, the transcriptomic responses to viral infection [e.g., infectious salmon anemia virus (ISAV) and infectious pancreatic necrosis virus (IPNV)] and/or poly(I:C) were previously explored in Atlantic salmon, rainbow trout (Oncorhynchus mykiss) (27), red-spotted grouper (Epinephelus akaara) (28), yellowhead catfish (Tachysurus fulvidraco) (29), ya-fish (Schizothorax prenanti) (30), and yellow catfish (Pelteobagrus fulvidraco) (31). To our knowledge, the lumpfish transcriptomic response to a viral pathogen has not been studied. While lumpfish primary leukocytes’ transcript expression responses to poly(I:C) have been studied recently (32, 33), the lumpfish systemic immune antiviral transcriptomic response [e.g., response of immune tissue/organ such as head kidney to in vivo stimulation with poly(I:C)] had not been characterized prior to the current study. Characterization of the lumpfish head kidney transcriptomic response to poly(I:C) will provide a foundation for understanding the mechanisms underlying the lumpfish immune response to viral infections, thereby aiding in the development of vaccines and the improvement of aquaculture practices.
In addition to the transcriptomic response to poly(I:C), several aspects of lumpfish antiviral mechanisms remain uncharacterized. Interferon regulatory factor (IRF) family members are key elements of fish immune responses (34). Transcript expression levels of several members of the IRF family (e.g., IRF3, 5, and 7) are upregulated following viral infections (34). IRFs are transcription factors, and their activation leads to the induced expression of interferons (IFNs) and IFN-stimulated genes (ISGs), which play crucial roles in antiviral responses. Additionally, irf expression is dysregulated in response to several stressors such as viral and bacterial infections, heat shock, and toxins (35–37). Moreover, while members of the IRF family are suggested to be highly conserved in their structure across vertebrate species, they also play some species-dependent roles (38–40). Several members of the IRF family were characterized in teleost species such as common carp (Cyprinus carpio L.) (41), Japanese flounder (Paralichthys olivaceus) (42), turbot (Scophthalmus maximus) (43, 44), Japanese seabass (Lateolabrax japonicus) and Atlantic cod (Gadus morhua) (45–47). However, the antiviral response of the lumpfish irf family members remained unknown prior to this study. In addition to transcriptomic profiling, in this study, we focused on evolutionary aspects and poly(I:C) responses of irf family members to investigate if these genes play conserved roles in lumpfish.
In the current study, we analyzed the transcriptomic response of lumpfish head kidney to intraperitoneal (IP) injection with poly(I:C) using RNA-Seq. Real-time quantitative polymerase chain reaction (qPCR) analyses were then utilized to assess expression levels of selected transcripts (e.g., representing hub genes in pathway analyses and well-known antiviral biomarkers) to confirm the results of the RNA-Seq analyses and, specifically, to elucidate the expression profiles of the 12 members of the lumpfish irf family in response to poly(I:C). Also, we investigated the molecular phylogeny of IRF members from lumpfish and three other teleost species representing different superorders to improve our understanding of the evolutionary history of the IRF family across Teleostei. The results of the current study enhance our knowledge of the genes and molecular pathways involved in the antiviral immune responses of fishes.
Materials and methods
Animals, experimental design, and sample collection
Juvenile lumpfish were raised at the Dr. Joe Brown Aquatic Research Building (JBARB), Ocean Sciences Centre, Memorial University of Newfoundland, Canada (48). At 300 days post-hatch, 18 fish (average weight ± standard deviation = 85.4 ± 14.6 g) were randomly selected and distributed into three 500-L tanks. The fish were held in 8–10°C filtered and UV-treated seawater, with a flow rate of 7.5 L/min, and the oxygen level was maintained at a saturation range of 95%–110%. The photoperiod was 12-h light and 12-h dark. The fish were fed a commercial diet (Marine grower diet, Zeigler Bros., Inc., Gardners, PA, USA) at 0.5% of the average body weight per day.
After 1 month, the poly(I:C) challenge study was conducted. Briefly, the fish were fasted for 24 h and then lightly anesthetized with MS-222 (50 mg/L, Syndel Laboratories, Vancouver, BC, Canada). Three individuals per tank were intraperitoneally injected with either poly(I:C) [2 µg/g fish, dissolved in phosphate-buffered saline (PBS)] or PBS (sham-injection/vehicle control). Immediately after injection, each fish was stitched with nylon surgical thread on one fin [the left pectoral fin for poly(I:C) or the right pectoral fin for PBS] and then returned to the same tank. At 24 hours post-injection (hpi), fish were euthanized using MS-222 (400 mg/L), and the head kidneys were collected. The samples were flash-frozen in liquid nitrogen and stored at −80°C until RNA extraction. All procedures in this experiment were performed following the Canadian Council of Animal Care guidelines (Memorial University of Newfoundland Animal Care Protocol, 17–03-RG and 18–01-MR) and in accordance with ARRIVE guidelines (https://arriveguidelines.org). Figure 1A depicts an overview of the workflow in this study.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g001.jpg)
Figure 1. (A) Overview of the methods used in the current study. This figure was generated using Biorender.com. (B) Heatmap and hierarchical clustering of polyinosinic:polycytidylic acid [poly(I:C)]-responsive transcripts [i.e., differentially expressed transcripts (DETs)]. The transcripts per million (TPM; log2 transformed) values for each DET are clustered along the vertical axis, the ID of the sample subjected to RNA-Seq is shown on the base of the figure, and the sample cluster is plotted on the top of the map. (C) Principal component analysis of DET of all the samples used in sequencing using TPMs.
Total RNA sample preparation
The RNA extractions were performed as described in Emam et al. (49). Briefly, the head kidney samples were homogenized in TRIzol Reagent (Invitrogen/Thermo Fisher Scientific, Burlington, ON, Canada) using stainless steel beads (5 mm; QIAGEN, Mississauga, ON, Canada) and a TissueLyser (QIAGEN), and the total RNA extractions were then completed following the manufacturer’s instructions. Total RNA samples (~25 µg) were treated with 6.8 Kunitz units of DNaseI (RNase-Free DNase Set, QIAGEN) for 10 min at room temperature and then purified using the RNeasy MinElute Cleanup kit (QIAGEN) following the manufacturer’s instructions. The integrity and the purity of the purified RNA were evaluated using 1% agarose gel electrophoresis and NanoDrop spectrophotometry (NSW-1000), respectively. The RNA samples used in this study showed high integrity (i.e., tight 28S and 18S ribosomal RNA bands at a ratio of ~2:1) and purity (i.e., A260/280 and A260/230 ratios above 1.9 and 2.0, respectively).
RNA−seq analysis
Library construction and RNA-Seq services were performed at the Centre d’expertise et de services (CES), Génome Québec, Montréal, QC, Canada. Prior to library construction, RNA quality was further evaluated using the Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA), and all samples were of high integrity (i.e., with RNA Integrity Numbers of 9.8–10). Ten libraries (i.e., for head kidney RNA samples from five of the nine poly(I:C)-injected fish and five of the nine PBS-injected fish) were constructed using the NEBNext mRNA Library Prep Reagent Set for Illumina (New England Biolabs, Whitby, ON, Canada) following the manufacturer’s instructions. Sequencing was performed on the Illumina NovaSeq 6000 [S2] with paired-end 100 bp (PE100) and at least 50 million reads per library. A summary of the RNA-Seq and library construction can be found in Supplementary Table S2.
RNA-seq data processing
The quality control, trimming, and filtering of low-quality reads were performed using FastQC version 0.11.8 and Trimmomatic version 0.3971. The sequencing reads were aligned to the lumpfish genome [C. lumpus (assembly fCycLum1.pri), ID: 86363] using HISAT2 version 2.1.0. StringTie version 2.0 was used to assemble and calculate the expression levels of all transcripts using reference gene models provided in the form of Gene Transfer Format (GTF) annotation files that are distributed with the lumpfish genome. StringTie was also used to assemble and quantify novel genes and transcripts. The accuracy of transcript assembly was evaluated using gffcompare version 0.11.2. The read count data used for differential expression analysis were obtained from the python script “prepDE.py” provided by the StringTie authors. The differential expression analysis was performed using DESeq2 version 1.28.1 in the Bioconductor package at an adjusted p-value (padj) threshold of <0.05, |log2 fold change| > 1, and filtering for transcripts with expression transcripts per million (TPMs) greater than 1 in at least two replicates in each group.
A heatmap with the TPM values (log2 transformed) of all differentially expressed transcripts (DETs) was generated using the “heatmap3” function of the “gplots” package in R (2023.03.1; Figure 1B). A principal component analysis (PCA; Figure 1C) was performed using the standardized TPMs of all of the identified DETs. A volcano plot was used to show significance (−log10 padj) against the expression log2 fold change (LFC). Gene Ontology (GO) term enrichment analysis was conducted using the ClueGO (50) plugin of Cytoscape (51) version 3.5 and using the C. lumpus ontology database. The enrichment analyses considered biological process, cellular component, and molecular function. The distribution of up- and downregulated transcripts within each leading GO term were plotted by LFC using density plots (Figure 2B). Leading GO terms were then manually classified based on related function into 1) intracellular processes and regulation of gene expression; 2) immune system, movement, cell structure, and apoptosis; and 3) cell signaling and response to stimuli (Figure 2). A circular bar plot was used to show GO terms with the highest percentage of DETs, after removing redundant GO terms, using library(tidyverse). The standard lumpfish gene symbols for all transcripts in our RNA-Seq dataset were used as references for the enrichment analyses, which were performed at the transcript level since ClueGO conducts the analyses using standard symbols. The heatmap and hierarchical clustering of DETs (using TPMs; log2 transformed) that participated in enriching “response to virus”, “regulation of retinoic acid receptor signaling pathway”, and “response to lipid” were generated using TBtools software (52).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g002.jpg)
Figure 2. (A) Gene Ontology (GO) term enrichment and pathway term network analysis of differentially expressed transcripts (DETs). The GO term enrichment analysis was performed using ClueGO plugin in Cytoscape. The p-value was adjusted at 0.05, kappa score level was ≥0.4 on ClueGO, and Benjamini–Hochberg correction was used. Biological process, cellular component, and molecular function were the selected ontologies on ClueGO. Nodes represent enriched GO terms. A complete list of the enriched GO terms is found in Supplementary Table S4, while the leading GO terms are also labeled in the figure. (B) Density plots of the fold change for the leading GO terms, showing the upregulated and downregulated genes in each GO term.
qPCR overview
To confirm the results of the RNA-Seq analyses, transcript expression levels of 18 genes (either hub genes or known antiviral biomarkers) that were identified as differentially expressed were also assessed using qPCR. In addition, expression levels of the 12 members of the lumpfish irf family [seven of which were also identified as poly(I:C)-responsive by RNA-Seq] were assessed to elucidate the response of this gene family to poly(I:C). Levels were assessed in all of the 18 study samples (i.e., head kidney samples from the poly(I:C)- and PBS-injected fish at 24 hpi; n = 9 per group).
cDNA synthesis and qPCR parameters
First-strand cDNA templates for qPCR were synthesized in 20-µL reactions from 1 µg of DNaseI-treated, column-purified total RNA using random primers (250 ng; Invitrogen), dNTPs (0.5 mM final concentration; Invitrogen), and M-MLV reverse transcriptase (200 U; Invitrogen) with the manufacturer’s first-strand buffer (1× final concentration) and DTT (10 mM final concentration) at 37°C for 50 min.
The qPCR amplifications were performed in 13 µL reactions containing 1× Power SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA, USA), 50 nM of both the forward and reverse primers, and the indicated cDNA quantity. All reactions were performed in triplicate and included no-template controls (NTCs). Amplifications were performed using either the ViiA 7 Real-Time PCR system or the QuantStudio 6 Flex Real-Time PCR system (384-well format) (Applied Biosystems). The qPCR analysis program included 1 cycle of 50°C for 2 min, 1 cycle of 95°C for 10 min, and 40 cycles of 95°C for 15 s and 60°C for 1 min, with fluorescence detection at the end of each 60°C step, and was followed by dissociation curve analysis.
Primer design and quality assurance testing
Previously published (2) and newly designed qPCR primers were used in this study. For the 25 qPCR-analyzed genes that were identified as differentially expressed in the RNA-Seq analyses, a BLASTn search of the non-redundant nucleotide (nr/nt) sequence database of the National Center for Biotechnology Information (NCBI) [C. lumpus (taxid: 8103) sequences only] was performed to identify annotated sequences corresponding to the transcript sequence generated herein [i.e., transcript of interest (TOI)]. This search also determined if gene paralogues/isoforms were present; the additional five irf family members were identified using BLASTn. A database of the sequences obtained from GenBank was compiled using Vector NTI (Vector NTI Advance 11.5.4, Life Technologies, Carlsbad, CA, USA). For a given gene, if paralogues/isoforms were present, multiple sequence alignments were performed using AlignX (Vector NTI Advance 11.5.4) to identify regions where paralogue/isoform-specific qPCR primers for the TOI could be designed (i.e., in an area with ≥3-bp difference between them). However, if sequences for transcript variants were present, primers were designed in a region that was conserved among the variants and generated identical amplicons. Most primers were designed using either PrimerQuest (www.idtdna.com/Primerquest/Home/Index) or Primer3 (53, 54) however, some were manually designed in paralogue/isoform-specific areas to ensure specificity. All primers had a melting temperature (Tm) of 60°C and were located in the CDS.
The qPCR primers utilized herein were subjected to quality control testing as described previously (55). All showed single-product amplification and the absence of primer dimer in the NTC using dissociation curve analysis. Amplification efficiencies were calculated for two cDNA pools generated from all individuals in the PBS group and from all individuals in the poly(I:C) group. Standard curves were generated using a 5-point 1:3-fold dilution series starting with cDNA representing 10 ng of input total RNA. The reported efficiencies are an average of the two values. The sequences, amplicon sizes, and efficiencies for all primer pairs used in the qPCR analyses are presented in Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g003.jpg)
Figure 3. Phylogenetic tree analysis of putative interferon regulatory factor (IRF) orthologues across the four teleost superorders. Putative IRF amino acid sequences from a fish representing each of the four teleost superorders, namely, Protacanthopterygii [Salmo salar (Atlantic salmon)], Acanthopterygii [Cyclopterus lumpus (lumpfish)], Paracanthopterygii [Gadus morhua (Atlantic cod)], and Ostariophysi [Danio rerio (zebrafish)] were collected from the National Center for Biotechnology Information (NCBI) non-redundant protein database (see Supplementary Table S5). The 59 amino acid sequences were aligned using ClustalW, and the tree was constructed using the Neighbor-Joining method with Poisson correction; the bootstrap test of phylogeny was performed with 10,000 replicates in the MEGA 11 (v.11.0.13) (60) software. The numbers at the branch points represent the bootstrap values, and branch lengths are proportional to calculated evolutionary distances. The scale represents the number of amino acid substitutions per site. The Atlantic salmon IRF paralogues were named as suggested in Clark et al. (2021), with the exceptions being IRF2–1 and IRF2–2, which were named as in Crossman et al. (2023); IRF1a (alias IRF11). The four subgroups—IRF1-G (IRF1 and IRF2), IRF3-G (IRF3 and IRF7), IRF4-G (IRF4, IRF8, IRF9, and IRF10), and IRF5-G (IRF5 and IRF6)—are shown in different colors. Corresponding proteins to transcripts explored using qPCR in the current study are marked by "*".
Endogenous control (normalizer) selection
Transcript expression levels of each gene of interest (GOI) were normalized to expression levels of two endogenous control genes. These endogenous controls were selected from five candidate normalizers [ribosomal protein lateral stalk subunit p1 (rplp1), ribosomal protein l32 (rpl32), poly(A) binding protein cytoplasmic 1b (papbc1b), eukaryotic translation initiation factor 3 subunit d (eif3d), eukaryotic translation elongation factor 1 alpha 2a (ef1a2a)]. Briefly, the fluorescence threshold cycle (CT) values for all 18 samples were measured for each of these transcripts using cDNA representing 5 ng of input total RNA and then analyzed using geNorm (qBase plus, Biogazelle NV, Zwijnaarde, Belgium) (56). This analysis identified rpl32 and ef1a2a as the most stably expressed normalizers, with geNorm M-values of 0.20 and 0.21, respectively.
Experimental qPCR and data analyses
The qPCR analyses were conducted according to MIQE guidelines (57). cDNA representing 5 ng of input total RNA was used as a template in the PCRs. The relative quantity (RQ) of each GOI in each of the 18 samples was then determined using the qBase relative quantification framework (58, 59). This was performed using the CT values measured for each GOI, with normalization to both rpl32 and ef1a2a and with the amplification efficiencies incorporated. For each GOI, the sample with the lowest normalized expression was used as the internal calibrator (i.e., assigned an RQ value = 1.0). The RQ values are presented as mean ± SE (Figures 4, 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g004.jpg)
Figure 4. qPCR results of selected transcripts that were identified as differentially expressed in response to polyinosinic:polycytidylic acid [poly(I:C)] in RNA-Seq analyses. Transcript levels are presented as mean ± SE relative quantity (RQ) values (i.e., values for the transcript of interest were normalized to both rpl32 and ef1a2a transcript levels and were calibrated to the individual with the lowest normalized expression level of that given transcript). For transcripts exhibiting homogeneity of variance across samples, significance was assessed using t-tests and is denoted with asterisks. For transcripts with unequal variance across samples, significance was assessed using the Mann–Whitney U test and is denoted with underlined asterisks. For both methods, significance levels are “**” for p ≤ 0.01, and “***” for p ≤ 0.001. FC, fold-change [mean RQ values for poly(I:C)/mean RQ values for phosphate-buffered saline (PBS)]. The plotted transcripts represent (A) Toll-like receptor 7 (tlr7), (B) ATP-dependent RNA helicase lgp2 (dhx58), (C) activating transcription factor 3 (atf3), (D) Interleukin-1 beta-like (il1b), (E) Interferon alpha/beta receptor 2-like (ifnar2), (F) Interferon-induced GTPbinding protein Mx-like (mx1), (G) Radical S-adenosyl methionine domain containing 2 (rsad2), (H) Cholesterol 25-hydroxylase like 3 (ch25hl3), E3 ubiquitin/ISG15 ligase, (I) TRIM25-like (trim25), (J) Tripartite motif-containing protein 16-like (trim16), (K) ADAM metallopeptidase domain 22 (adam22), (L) ADAM metallopeptidase with thrombospondin type 1 motif, 15a (adamts15a), (M) Sacsin (sacs), (N) Sacsin-like (sacs-like), (O) Adenosine monophosphate deaminase 2b (ampd2b), (P) Jumonji, AT rich interactive domain 2b (jarid2b), (Q) PHD finger protein 8 (phf8), (R) Vesicle-associated membrane protein 8 (endobrevin), vamp8.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g005.jpg)
Figure 5. qPCR analysis of the response of the 12 irf family members in lumpfish to intraperitoneal (IP) challenge with polyinosinic:polycytidylic acid [poly(I:C)]. Transcript levels are presented as mean ± SE relative quantity (RQ) values (i.e., values for the transcript of interest were normalized to both rpl32 and ef1a2a transcript levels and were calibrated to the individual with the lowest normalized expression level of that given transcript). Significance was assessed using t-tests and is denoted with asterisks (“**” for p ≤ 0.01, and “***” for p ≤ 0.001). FC, fold-change [mean RQ values for poly(I:C)/mean RQ values for phosphate-buffered saline (PBS)]. The irfs with underlined gene symbols were identified as differentially expressed in the RNA-Seq analysis and were validated using qPCR. The irfs with non-underlined symbols were not differentially expressed in the RNA-Seq analysis using the preidentified cutoff criteria. Both irf1a [named in National Center for Biotechnology Information (NCBI) as irf1-like] and irf4b (named in NCBI as irf4-like) were named based on the closest orthologues in the phylogenetic tree (Figure 3). (A) interferon regulatory factor 1a (irf1a), (B) Interferon regulatory factor 1b (irf1b), (C) Interferon regulatory factor 2 (irf2), (D) Interferon regulatory factor 3 (irf3), (E) Interferon regulatory factor 4a (irf4a), (F) Interferon regulatory factor 4b (irf4a), (G) Interferon regulatory factor 5 (irf5), (H) Interferon regulatory factor 6 (irf6), (I) Interferon regulatory factor 7 (irf7), (J) Interferon regulatory factor 8 (irf8), (K) Interferon regulatory factor 9 (irf9), (L) interferon regulatory factor 10. (irf10).
Phylogenetic tree analysis of putative IRF orthologues
Putative orthologous amino acid (AA) sequences for IRF family members from a fish representing each of the four teleost superorders, namely, Protacanthopterygii (Atlantic salmon), Acanthopterygii (lumpfish), Paracanthopterygii (Atlantic cod), and Ostariophysi (zebrafish; D. rerio), were collected from the NCBI GenBank non-redundant (nr) protein database (59 in total). The GenBank accession numbers and AA sequences are provided in Supplementary Table S5. The sequences were subjected to BLASTP analyses to help identify all of the IRF isoforms/paralogues for each species and to ensure that all sequences used in the tree were unique (i.e., the tree did not include transcript variants; if present, the best representative sequence was selected). Phylogenetic and molecular evolutionary analyses were conducted using MEGA 11 (v.11.0.13) (60). Briefly, a multiple sequence alignment was performed using the ClustalW algorithm. The phylogenetic tree was then constructed using the Neighbor-Joining method with the Poisson correction; the bootstrap test of phylogeny was performed with 10,000 replicates.
Statistical analysis
All of the residuals were tested for homoscedasticity and normality (i.e., Levene’s and Shapiro–Wilk tests). Significant (p < 0.05) differences in transcript expression levels between the PBS- and poly(I:C)-injected groups were assessed using either Student’s t-test or the Mann–Whitney U test (for genes that failed the normality test). These analyses were performed using SPSS (IBM SPSS Statistics, Version 25, Armonk, NY, USA). PCA (Figure 1C) was performed using PRIMER 7 (PRIMER-E Ltd., Auckland, New Zealand). The scatter plot for the LFC (Supplementary Figure S2) was generated using the function “ggscatter” for the library “ggpubr”.
Results
Lumpfish head kidney transcriptome assemblies and RNA-Seq analyses
Transcriptome sequencing and assemblies
In this study, RNA-Seq was used to profile the responses of lumpfish to viral mimic, poly(I:C). Supplementary Table S2 summarizes the RNA-Seq read quality control for the 10 samples. The average number of raw reads across all samples was 79.7 M (range, ~62M to 92M). On average, 98% of the read pairs (across all samples) survived the trimming process (range, 97.8% to 98.2%). The average percentage of reads that dropped during trimming was 0.4% (range, 0.3% to 0.5%). Overall, the results show that most reads were successfully trimmed and kept. Both uniquely and multi-mapped reads were used for transcript assembly. Overall, ~94%–96% of the processed reads in each sample were mapped to the lumpfish genome (Supplementary Table S2).
Poly(I:C)-responsive transcripts in lumpfish head kidney
We identified 4,499 upregulated and 3,952 downregulated transcripts by poly(I:C) in the head kidney of lumpfish (Supplementary Table S3, Figures 1B, 6). All of the samples belonging to a given group (i.e., PBS- or poly(I:C)-injected) clustered together based on the expression of all of the identified DETs. In addition, the PBS- and poly(I:C)-injected fish were clearly segregated in the PCA space. PC1 explained 38.7% of the variability, and PC2 explained 10.4% of the variability (Figure 1C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g006.jpg)
Figure 6. Volcano plot reporting −log10 adjusted p-values (padj) against log2 fold-changes. Dots represent differentially expressed transcripts. The colored dots represent the most significant transcripts from the upregulated and downregulated probes (red and green, respectively) (Supplementary Table S3). Dots are labeled with the gene symbols of the putative human orthologues, and one uncharacterized is labeled with the locus annotation.
There were 311 GO terms enriched in the poly(I:C)-induced transcript list including 248 biological processes, five cellular components, and seven molecular functions. The most enriched biological process GO terms were related to “defense response”, “cytokine-mediated signaling pathway”, and “response to other organisms” (Figure 2). Enriched molecular function GO terms included “immune receptor activity”, “cytokine receptor binding”, and “interleukin-1 binding”. Enriched cellular component GO terms included “chromosome”, “nuclear protein-containing complex”, and “lipid droplet” (Supplementary Table S4). GO terms with the highest percentage of DETs included “regulation of retinoic acid receptor signaling pathway” (Supplementary Figure S1).
The TPMs of the genes represented by the enriched GO term “response to virus” were used to generate the heatmap plotted in Figure 7. It showed upregulation of transcripts including lysosomal trafficking regulator (lyst), complement component 1, q subcomponent binding protein (c1qbp), DEAD-box helicase 3 X-linked a (ddx3xa), interleukin-12 subunit beta (il12b), stimulator of interferon response cGAMP Interactor 1 (sting1), toll-like receptor 7 (tlr7), cholesterol 25-hydroxylase A (ch25ha), mitochondrial antiviral signaling protein (mavs), BCL2 apoptosis regulator B (bcl2b), and irf2 with poly(I:C) challenge. Also, it showed downregulation of several genes, for example, vesicle-associated membrane protein 8 (vamp8), spondin2a, and scavenger receptor cysteine-rich type 1 protein M130-like.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g007.jpg)
Figure 7. Heatmap and hierarchical clustering of differentially expressed transcripts (DETs) associated with “response to virus” using transcripts per million (TPMs; log2 transformed).
To investigate dysregulation in metabolism during a viral mimic challenge, we plotted a heatmap for DETs representing the enriched GO terms “regulation of retinoic acid receptor signaling pathway” and “response to lipid” GO terms (Figures 8A, B). The GO term regulation of the retinoic acid receptor signaling pathway was one of the GO terms with the highest percentage of DETs. Figure 8A shows the upregulation of different transcripts including tripartite motif containing 16-like (trim16-like) and tripartite motif containing 25-like (trim25-like). Also, it shows downregulation of genes including dehydrogenase/reductase 3B (dhrs3b), kruppel-like factor 17 (klf17), c-terminal binding protein 2a (ctbp2a), and others representing mitogen-activated protein kinase 11-like, vesicle-associated membrane protein 8 and fibroleukin with poly(I:C) injection.
The GO term “response to lipid” was enriched with a mixture of upregulated and downregulated genes (Figure 8B). For example, mucosa-associated lymphoid tissue lymphoma translocation protein 1 (malt1), NLR family CARD domain containing 3-like (nlrc3l), and B-cell lymphoma 2 B (bcl2b) were upregulated with poly(I:C). While frizzled class receptor 4 (fzd4), annexin a2 receptor (anxa2r), prostaglandin E receptor 2 (ptger2a), peroxisome proliferator-activated receptor alpha B (pparab), and peroxisome proliferator-activated receptor gamma coactivator 1 alpha (ppargc1a) were downregulated with poly(I:C) stimulation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1439465/fimmu-15-1439465-HTML/image_m/fimmu-15-1439465-g008.jpg)
Figure 8. Heatmap and hierarchical clustering of differentially expressed transcripts (DETs) that participated in enriching “regulation of retinoic acid receptor signaling pathway” in panel (A) and “response to lipid” in panel (B) using transcripts per million (TPM; log2 transformed).
qPCR validation
The qPCR results confirmed the RNA-Seq results for all of the selected transcripts (Supplementary Figure S2). There was a significant correlation (p = 0.007; Pearson’s coefficient = 0.61) between the qPCR and the RNA-Seq results (Supplementary Figure S2).
The qPCR analysis confirmed the poly(I:C) induction of tlr7, ATP-dependent RNA helicase lgp2 (dhx58), activating transcription factor 3 (atf3), interleukin 1 beta (il1b), interferon alpha/beta receptor 2 (ifnar2), interferon-induced GTP-binding protein Mx-like (mx1), radical S-adenosyl methionine domain containing 2 (rsad2), ch25ha, trim25, trim16, ADAM metallopeptidase domain 22 (adam22), and spastic ataxia of Charlevoix-Saguenay (sacsin)-like (sacs-like) (Figure 4).
The qPCR analysis confirmed the poly(I:C) repression of ADAM metallopeptidase with thrombospondin type 1 motif 15 (adamts15a), sacs, adenosine monophosphate deaminase 2b (ampd2b), jumonji and at-rich interaction domain containing 2 (jarid2b), histone lysine demethylase phf8 (phf8), and vamp8 (Figure 4).
Lumpfish irf family qPCR and phylogenetic analysis
To characterize the response of the 12 irf family members in lumpfish stimulated with poly(I:C), their transcript expression levels were assessed using qPCR. Seven members, i.e., irf1b, irf2, irf3, irf7, irf8, irf9, and irf10, were identified as dysregulated with poly(I:C) challenge in the RNA-Seq analyses, whereas irf1a, irf4a, irf4b, irf5, and irf6 were not detected as differentially expressed using a cutoff level of LFC = 1.0. With the exception of irf6, all irf family members subjected to qPCR analysis were found to be poly(I:C)-responsive. The transcript levels of irf1a, irf1b, irf2, irf3, irf7, irf8, irf9, and irf10 were upregulated with poly(I:C) injection when compared with the PBS-injected group. The levels of irf4a and irf5 were suppressed in response to poly(I:C) injection. The levels of irf6 were not significantly different between the poly(I:C)- and PBS-injected groups (Figure 5).
A total of 59 IRF family members in species representing the four teleost superorders (lumpfish, Acanthopterygii; Atlantic salmon, Protacanthopterygii; zebrafish, Ostariophysi; and Atlantic cod Paracanthopterygii), identified in the GenBank nr protein database, were used to build a phylogenetic tree (Figure 3). As anticipated, the teleost IRF sequences cluster into the four previously defined subgroups: IRF1-G (IRF1 and IRF2), IRF3-G (IRF3 and IRF7), IRF4-G (IRF4, IRF8, IRF9, and IRF10), and IRF5-G (IRF5 and IRF6) (61). The phylogenetic tree shows that lumpfish IRF2, IRF4a, IRF4b, IRF5, IRF6, IRF7, IRF8, and IRF10 are evolutionarily more closely related to the Atlantic cod compared with the Atlantic salmon or zebrafish putative orthologues (Figure 3). However, the tree also reveals that lumpfish IRF3 is most closely related to Atlantic salmon IRF3 (i.e., sharing a branch point); likewise, lumpfish IRF9 is more closely related to Atlantic salmon IRF9 paralogues (IRF9–1 and IRF9–2, arising from duplication in the salmonid lineage) than to Atlantic cod or zebrafish IRF9 sequences (Figure 3). Finally, for IRF1a and IRF1b, zebrafish and Atlantic salmon sequences are more closely related to each other than either are to the lumpfish putative orthologous sequences.
Discussion
The current RNA-Seq results detected an extensive global gene expression response (i.e., 4,499 upregulated and 3,952 downregulated DETs) in the head kidney of poly(I:C)-challenged lumpfish, shedding light on the molecular mechanisms and immune pathways involved in response to this viral mimic. GO term analyses identified enriched immune-related GO terms consistent with an immune response to a viral challenge based on the current knowledge about molecular antiviral responses in teleost fish (13, 14). Several leading immune-related GO terms (e.g., “antigen processing and presentation of peptide antigen via MHC class I”, “leukocyte apoptotic process”, “pattern recognition receptor signaling pathway”, “response to virus”, “type I interferon signaling pathway”, and “NIK/NF-kappaB signaling”) were primarily represented by poly(I:C)-induced DETs compared with downregulated ones (Figure 2B). Many of the leading GO terms within cellular component and biological process categories (e.g., “RNA splicing, via transesterification reactions with bulged adenosine as nucleophile”, “chromosome”, “nuclear protein-containing complex”, and “cell migration”) were represented by a higher number of downregulated DETs compared with upregulated DETs. These results suggest a possible shift in cellular activity toward fighting infection, highlighting gene expression regulation patterns toward the induction of immune-relevant genes.
Our GO term enrichment results generally agreed with those of a recent in vitro study (32) investigating the effects of poly(I:C) on the transcriptome of lumpfish primary leukocytes after 24 h of exposure, especially in terms of general immune-relevant GO terms such as “cytokine receptor binding” and “response to virus”. However, the present study identified a more varied suite of enriched biological processes, including GO terms related to cell migration (e.g., “leukocyte migration”), apoptosis (e.g., “leukocyte apoptotic process”), and adaptive immunity (e.g., “antigen processing and presentation of peptide antigen via MHC class I”), which may be a consequence of analyzing transcriptome changes in head kidney samples as opposed to primary leukocytes. Notably, in Rao et al. (32), biological processes associated with the metabolism of nitrogen-containing compounds (e.g., “peptide metabolic process”) were dominant. Our analyses also identified enriched GO terms related to retinoic acid (e.g., “regulation of retinoic acid receptor signaling pathway”) like in Rao et al. (32) but also other terms associated with lipids (e.g., “lipid droplet” and “response to lipid”). We anticipate that single-cell RNA-Seq in lumpfish will likely allow the assignment of transcript expression changes to specific cells within the head kidney, elucidating GO terms enriched in each type of cell.
During viral infection, the host’s metabolism undergoes various changes to combat the virus and support the immune response. This may be evidenced in the current study by enrichment in GO terms relevant to metabolism (e.g., GO terms related to lipid and retinoic acid metabolism; Supplementary Table S3, Figures 2A, B). In mammals, lipid droplets (62, 63) and retinoic acid (64) are key players in the viral infection mechanisms and the inflammatory processes they trigger. The GO term “regulation of retinoic acid receptor signaling pathway” was one of the GO terms with the highest percentage of DETs (Supplementary Figure S1). The heatmap for “regulation of retinoic acid receptor signaling pathway” (Figure 8A) shows that several trim16-like transcripts were upregulated by poly(I:C), except for one downregulated transcript (accession number: XM_034555257.1). The heatmap also showed decreased transcript levels for ftr83 (tripartite motif-containing protein 11), dhrs3b, foxc1a (forkhead box C1a), klf17 (krüppel-like factor 17), ctbp2a (C-terminal binding protein 2a), and znf536 (zinc finger protein 536) in the poly(I:C)-injected fish (Figure 8A). These poly(I:C)-responsive genes encode proteins involved in regulating the retinoid metabolic process (dhrs3b), regulation of cell proliferation (foxc1a), differentiation (klf17), and response to external stress (e.g., ctbp2a and znf536) (65–68). Their regulation may be instrumental to mounting the immune response to the viral mimic challenge in lumpfish and highlights the role of vitamin A (retinoic acid) during viral infection in lumpfish. However, vitamin A modulation of the immune response during viral infection requires further research.
Several upregulated and downregulated genes by poly(I:C) identified herein contributed to enriching the GO term “response to lipid” (Figure 8B). For example, transcripts encoding receptors (e.g., fzd4, anxa2r, and ptger2a) and transcription factors involved in metabolism, energy homeostasis, and immunomodulation (69) (e.g., pparab and ppargc1a) were found downregulated by the viral mimic challenge (Figure 8B). Others were found to be upregulated, for example, transcripts encoding proteins involved in the activation of the transcription factor NF-κB (i.e., malt1), intracellular pattern recognition receptors previously suggested to regulated innate immune response (70) (e.g., nlrc3l), and apoptosis (e.g., bcl2b) (71). Overall, the enrichment of “response to lipid” and other lipid-related GO terms may highlight the role of lipids during viral infection (62, 72–74), emphasizing the potential of dietary lipids to modulate the antiviral fish immune response.
The enriched GO term “response to virus” was predominantly represented by upregulated transcripts involved in antiviral immune defense [e.g., mavs, and signal transducer and activator of transcription 1a (stat1a)] (75), inflammation (e.g., ch25h, tumor necrosis factor-alpha, and il6) (76), oxidative stress (e.g., c1qbp) (77), and apoptosis regulation (e.g., bcl2b) (78). Although to a lesser extent, this GO term was also represented by poly(I:C)-repressed genes related to various cellular process, like vesicle trafficking (e.g., vamp8) (79) and autophagy (e.g., autophagy related 7 homolog and atg7) (80) as well as deubiquitination and protein metabolism (e.g., mitochondrial E3 ubiquitin protein ligase 1b and mul1b) (81). The dysregulated transcripts representing this GO term indicate the involvement of genes with diverse functions in the antiviral immune response.
qPCR results of lumpfish response to poly(I:C)
In the current study, all of the 18 genes used in the qPCR validated the RNA-Seq, indicating the reliability of the RNA-Seq results (Supplementary Figure S1). In the qPCR study, two pattern recognition receptor (PRR)-encoding genes, tlr7 and dhx58, were found to be upregulated by poly(I:C). Tlr7 contributed to 59 enriched GO terms, including “NIK/NF-kappaB signaling”, “pattern recognition receptor signaling pathway”, “response to virus”, and “regulation of cytokine production”. Similarly, dhx58 contributed to enriching 83 GO terms, including “response to virus”, “positive regulation of intracellular signal transduction”, and “regulation of response to stimulus” (Supplementary Table S3). PRRs detect conserved pathogen-associated molecular patterns (PAMPs) and consequently activate the innate immune response (82). Upon recognition of viral single-stranded RNA molecules, TLR7 triggers the production of type I IFNs and pro-inflammatory cytokines, which are critical for the host’s antiviral response (83–85). DHX58 is a member of the RLR family and an ATP-dependent RNA helicase, also known as LGP2 or RIG-I-like receptor 1 (RLR1), and has several antiviral roles such as regulation of TLRs and RLRs (86). Dhx58 was also found to be upregulated in Atlantic cod spleen and brain after IP injection with poly(I:C) in the brains of nodavirus-positive Atlantic cod (87, 88) and lumpfish infected with R. salmoninarum (2). Dhx58 was upregulated in lumpfish larvae after oral immunization against Vibrio anguillarum (8); also, it was highly upregulated in lumpfish leukocytes after stimulation with poly(I:C) (32). The results of previous and current studies suggest that DHX58 plays a role in both the antibacterial and antiviral immune responses of lumpfish, whereas TLR7 dysregulation occurs following the activation of its antiviral responses.
In the current study, the transcript levels of atf3 were strongly upregulated (over 69-fold) in the head kidneys of lumpfish stimulated with poly(I:C), as it was previously reported for Atlantic salmon macrophages (13) and Atlantic cod spleen (89). Mammalian ATF3 is activated by cellular stress response pathways and plays a role in the host’s immune response to viral infections (89, 90). ATF3 upregulation during viral infections promotes a stronger immune response and increased resistance in various mammalian species including mice and humans (91–93). Additionally, human ATF3 can directly inhibit the replication of some viruses by suppressing their transcription (93) and regulating the expression of the host’s immune-related genes such as IFN-induced and pro-inflammatory cytokines (94). While our results indicate the involvement of atf3 in lumpfish’s antiviral immune response, its viral inhibitory and regulatory functions are yet to be investigated in this species.
Two lumpfish transcripts encoding proteins classified as cytokines and cytokine receptors (i.e., il1b and ifnar2) were found to be over fourfold upregulated by poly(I:C). Several ILs contributed to enriching key GO terms, including “response to interleukin-1”, “signaling receptor binding”, “leukocyte migration”, and “defense response to other organism” (Supplementary Table S3). IL1B is a pro-inflammatory cytokine mediating the immune response of fish to viral and bacterial infection (95, 96). The production of IL1B by immune cells such as macrophages is triggered by the detection of viral nucleic acid by PRRs, such as TLRs (97). Previous studies showed il1b induction in the kidney of Sockeye salmon (Oncorhynchus nerka) infected with infectious hematopoietic necrosis virus (IHNV) (95, 96). Type I interferons, including IFNA, are produced following viral detection (98) and can activate cellular antiviral immune mechanisms (e.g., the expression of IFN-stimulated genes) and the recruitment of immune cells such as natural killer cells and T cells (98). It was previously reported that IFNA inhibited Salmonid Alphavirus Subtype 3 replication in a salmon cell line (i.e., TO cells originated from head kidney leukocytes) (99). Additionally, ifna was found upregulated in the head kidney of Atlantic salmon with New Piscine Orthomyxovirus (POMV) infection (100). These results collectively emphasize the conserved roles of IL1B and IFNA in the antiviral responses of lumpfish as in other teleost fishes.
Antiviral markers mx1, rsad2 (alias viperin), and ch25hl3 were found upregulated by poly(I:C) stimulation. Rsad2 and ch25hl3 contributed to enriching several GO terms, e.g., “response to virus” (Figure 7) and “defense response to other organism” (Supplementary Table S3). Atlantic salmon mx and rsad2 showed strong upregulation in ISAV-infected TO cells (101), poly(I:C)-stimulated macrophages (13), and the head kidney of poly(I:C)-injected fish (55). MX1 plays a role in the salmon immune response to viral infections, notably myxoviruses such as ISAV (102), and spring viraemia of carp virus (SVCV) (103). The poly(I:C)-induced mx1 can enhance resistance to viral infections in salmonids (104). RLR-activated RSAD2 (105) regulates the RLR signaling pathway through phosphorylation of downstream targets, e.g., MAVS (mitochondrial antiviral-signaling protein) and IRF3 [interferon regulatory factor 3; one of the top upregulated transcripts in the current study], leading to amplified antiviral response (106). CH25H plays a role in the immune response of salmon to bacterial infections (e.g., Renibacterium salmoninarum and Piscirickettsia salmonis) (49, 107, 108). Viral infections can lead to increased expression of ch25h in salmon (109). Additionally, it has been found that fish CH25H can directly inhibit the replication of some viruses (110). Altogether, mx1, rsad2, and ch25hl3 responses seen herein reflect the activation of antiviral agents by poly(I:C) in lumpfish and suggest these transcripts as potential antiviral biomarkers for this species.
The mRNA levels of trim16 and trim25, which play roles as immune regulators, were upregulated by poly(I:C) in lumpfish. TRIM25 shared in enriching several GO terms, such as “signaling receptor binding”, “response to virus”, “regulation of cytokine production”, and “response to lipid” (Supplementary Table S3). TRIM16 and TRIM25 are E3 ubiquitin ligases that play key roles in the host’s immune response to viral infections (111). Lumpfish leukocytes stimulated with poly(I:C) showed higher levels of trim25 from 6 to 24 h post-challenge (32). In Atlantic salmon TO cells, trim16 and trim25 were strongly upregulated by ISAV infection (101). The literature regards TRIM16L as a negative regulator of IFN-mediated antiviral responses in fish. However, the role of this protein in antiviral immune responses may be cell type-dependent, based on the available gene expression regulation data from fish and human cells exposed to viral infection (65, 112, 113). Human TRIM25 is reported to be able to target and degrade viral proteins (e.g., influenza A virus), thereby inhibiting viral replication (114). The conserved induction of trim16 and trim25 found in the current study highlights the importance of these factors in the antiviral response of lumpfish and suggests their role in the regulatory mechanism by which lumpfish respond to viral pathogens.
Adam22 was found upregulated with poly(I:C) stimulation and contributed to enriching “organonitrogen compound metabolic process” GO term (Supplementary Table S3). In a similar direction to our results, transcript levels of adam22 were previously found upregulated in the brain of Atlantic cod injected with poly(I:C) (88). Human ADAM22 has been shown to mediate the entry of the human rhinovirus (HRV) (115). However, the role of ADAM22 during antiviral immune response remains to be elucidated in lumpfish.
Lumpfish sacs-like (accession number: XM_034562115) was upregulated by poly(I:C) stimulation, whereas its paralogue (XM_034549198) was poly(I:C)-suppressed (Figure 4, Supplementary Table S3). The levels of Atlantic cod sacs have been reported to increase in the brain of nodavirus carrier fish (88) and Atlantic cod macrophages stimulated with poly(I:C) (14). It was also found upregulated in the brain of sockeye salmon infected with IHNV (116). The current study results may indicate paralogue-specific functions for lumpfish sacs genes. Opposite transcriptional regulation was previously reported for some paralogues in zebrafish and salmon (117, 118). However, further research is needed to understand the implications of the different sacs regulation patterns in lumpfish antiviral immune responses.
We identified several immune-related genes downregulated in poly(I:C)-stimulated lumpfish. The transcript levels of ampd2b were found to be downregulated by poly(I:C) and contributed to enriching the “GTP metabolic process” GO term. AMPD2 is an enzyme involved in the regulation of cellular energy levels, i.e., purine metabolism by converting adenosine monophosphate (AMP) to inosine monophosphate (IMP) (119). Therefore, the downregulation of ampd2b may indicate a decrease in energy production or utilization in the head kidney cells of the lumpfish, which could be part of the response to the stress caused by the viral mimic. Transcripts encoding proteins that are involved in epigenetic regulation and histone modification, such as jarid2b (120) and phf8 (121), were suppressed (qPCR and RNA-Seq results) in the head kidney of the poly(I:C)-injected fish (Figure 4, Supplementary Table S3). Jarid2b contributed to enriching several GO terms, including “regulation of cellular protein metabolic process”, and phf8 contributed to enriching different GO terms including “immune system development” (Supplementary Table S4). Although histone modifications were found to play a protective role in the body’s defense against viral infections (122), the functions of jarid2b and phf8 in epigenetic interaction and immune response of teleost fish require further study. The transcript levels of vamp8 were downregulated in both the RNA-Seq and qPCR results (Figure 4, Supplementary Table S3). Also, they shared in enriching several GO terms, e.g., “regulation of cell activation” and “positive regulation of multicellular organismal process” (Supplementary Table S4). VAMP8 plays several roles in intracellular membrane trafficking and fusion (123). VAMP8 may also be implicated in the release of cytokine and the inhibition of phagocytosis (124). The downregulation of vamp8 may be part of the above-hypothesized inhibited cellular function or the immune response regulation.
Different members of irf family response to poly(I:C)
IRFs are a family of transcription factors that play a key role in the host’s immune response to viral infections, especially in regulating IFN and interferon-stimulated genes (125). In the current study, the qPCR results showed that irf1a, irf1b, irf2, irf3, irf4b, irf7, irf8, irf9, and irf10 were significantly and strongly upregulated (FC range, 3.2 for irf4b and irf7 to 23.7 for irf1b) by poly(I:C) injection (Figure 5). Similar to qPCR results, the RNA-Seq also identified significant induction of irf1b, irf2, irf3, irf7, irf8, irf9, and irf10 by poly(I:C) (Supplementary Table S3, Figure 5). Various irfs (e.g., irf1, irf2, irf3, irf4b, irf7, irf9, and irf10) (46) were previously found upregulated with poly(I:C) or viral infection in teleost species (13, 126–129). The irf family members (i.e., irf1, irf2, irf3, irf7, irf8, and irf10) detected in our transcriptomic analysis contributed to enriching various GO terms, for example, “pattern recognition receptor signaling pathway”, “response to type I interferon”, “regulation of signaling”, and “cytokine production”. In contrast, irf9 was only involved in enriching “cytokine production”. These IRFs have been shown to play several roles in the activation of innate immunity and the production of type I interferons, which are important in the defense against viral infections (46, 47). IRF1 is a negative regulator of cytokine-induced cell proliferation in mammals (130). It has been reported that irf2 positively regulates the antiviral responses of large yellow croaker (Larimichthys crocea) (129). IRF7 and IRF3, key family members involved in antiviral responses, are activated downstream of the RLR and RLR pathways and enhance the expression of several immune genes such as IFNs (131). Human IRF8 supports the rapid expansion of virus-specific natural killer cells by enhancing the expression of genes involved in the cell cycle (132). IRF9 mediates the type I interferon responses, resulting in the production of IFN-induced genes (133). Zebrafish IRF 2, 4b, and 10 were suggested to be negative regulators of IFN (134), indicating that their role in the host’s antiviral response may be different among species. The induction of irf genes by poly(I:C), alongside dysregulation of several genes involved in TLR, RLR, and IFN pathways seen herein, highlights the importance of these transcription factors in antiviral responses of lumpfish. However, despite conserved structure, IRF family members may have species-specific regulatory functions (135, 136); further studies are needed to functionally characterize the lumpfish IRFs.
Unlike other lumpfish irfs studied here, irf4a and irf5 were significantly downregulated (less than twofold) in response to poly(I:C). IRF4 was previously reported to play a role in the differentiation of immune cells and the regulation of the immune response (137–139). In agreement with these findings, seabream irf5 was found to be downregulated with NNV at 12 h post-infection (140). Poly(I:C)-dependent irf4a and irf5 downregulation seen herein suggest their potential role in the regulation of antiviral responses in lumpfish.
The phylogenetic analysis of IRF sequences from lumpfish, Atlantic salmon, Atlantic cod, and zebrafish was used to examine the evolutionary history of IRF family members in lumpfish in the current study. The majority of lumpfish IRF family members (e.g., IRF2, IRF4a, IRF4b, IRF5, IRF6, IRF7, and IRF10) were grouped with the corresponding Atlantic cod orthologues. Overall, the phylogenetic tree suggests a high degree of similarity and evolutionary conservation between the IRF family members in lumpfish and their orthologues in other species. The observed grouping supports the notion that these specific IRF genes have been conserved over evolutionary time in teleost fishes, highlighting their functional importance across species, e.g., lumpfish and Atlantic cod. This finding contributes to our understanding of the evolutionary relationships and conservation of IRF genes in lumpfish.
Conclusion
Our findings suggest that poly(I:C) injection dysregulated diverse pathways associated with the antiviral immune system, cellular differentiation, cytokine production and response, NF-κB signaling, response to retinoic acid and lipids, and cell migration in the lumpfish head kidney. The leading GO terms related to cellular processes were enriched with more downregulated transcripts than the upregulated ones (e.g., “chromosome”). In contrast, GO terms with immune-relevant enriched pathways were dominated by upregulated genes. Our qPCR results validated the upregulation of genes involved in innate immunity and antiviral defense mechanisms and the downregulation of those with putative roles in cellular processes (e.g., histone modification: jarid2b and phf8).
The regulation of several lumpfish irf family members with poly(I:C) injections suggests their involvement in the host’s antiviral response. However, the functional characterization of IRF family members in lumpfish requires additional investigation. The results of the current study provide valuable insight into the underlying mechanisms of the induction of the innate immune system using poly(I:C) and suggest potential targets for developing therapeutic strategies and evaluating vaccine efficacy in lumpfish.
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 below: PRJNA1082277 (SRA).
Ethics statement
All procedures in the present study were approved (Protocol numbers: 17-03-RG and 18-01-MR) by the Animal Care Committee of Memorial University, following the guidelines of the Canadian Council on Animal Care and in accordance with ARRIVE guidelines (https://arriveguidelines.org). The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
ME: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. SK: Formal analysis, Writing – review & editing. KE: Writing – original draft, Investigation, Writing – review & editing. AC-S: Investigation, Writing – review & editing. JH: Writing – original draft, Writing – review & editing. XX: Writing – original draft, Writing – review & editing. HP: Writing – original draft, Writing – review & editing. RG: Writing – original draft, Writing – review & editing. JS: Writing – original draft, Writing – review & editing. MR: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant to MLR (2020–04519) and an Ocean Frontier Institute Vitamin Research Fund award to RLG and HP.
Acknowledgments
We would like to thank the JBARB team for helping with rearing the fish. We want to thank the Atlantic Computational Excellence Network (ACENET) and Compute Canada for facilitating the computational resources for the RNA-Seq data bioinformatics analyses.
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.2024.1439465/full#supplementary-material
References
1. Gendron RL, Hyde T, Paradis H, Cao T, Machimbirike VI, Segovia C, et al. CD45 in ocular tissues during larval and juvenile stages and early stages of V. Anguillarum infection in young lumpfish (Cyclopterus lumpus). Fish Shellfish Immunol. (2022) 128:523–35. doi: 10.1016/j.fsi.2022.08.023
2. Gnanagobal H, Cao T, Hossain A, Dang M, Hall JR, Kumar S, et al. Lumpfish (Cyclopterus lumpus) is susceptible to Renibacterium Salmoninarum infection and induces cell-mediated immunity in the chronic stage. Front Immunol. (2021) 12:733266. doi: 10.3389/fimmu.2021.733266
3. Powell A, Treasurer JW, Pooley CL, Keay AJ, Lloyd R, Imsland AK, et al. Use of lumpfish for sea-lice control in salmon farming: challenges and opportunities. Rev Aquac. (2018) 10:683–702. doi: 10.1111/raq.12194
4. Torrissen O, Jones S, Asche F, Guttormsen A, Skilbrei OT, Nilsen F, et al. Salmon lice – impact on wild salmonids and salmon aquaculture. J Fish Dis. (2013) 36:171–94. doi: 10.1111/jfd.12061
5. Snieszko SF. “9 - NUTRITIONAL FISH DISEASES.” In: Halver JE, editor. Fish Nutrition. Academic Press (1972). p. 403–37. doi: 10.1016/B978-0-12-319650-7.50014-6
6. Guðmundsdóttir S, Vendramin N, Cuenca A, Sigurðardóttir H, Kristmundsson A, Iburg TM, et al. Outbreak of viral haemorrhagic septicaemia (VHS) in lumpfish (Cyclopterus lumpus) in Iceland caused by VHS virus genotype IV. J Fish Dis. (2019) 42:47–62. doi: 10.1111/jfd.12910
7. Tollaksvik T. (2023). The University of Bergen. Available online at: https://bora.uib.no/bora-xmlui/handle/11250/3085928 (Accessed December 1, 2023).
8. Dang M, Cao T, Vasquez I, Hossain A, Gnanagobal H, Kumar S, et al. Oral Immunization of Larvae and Juvenile of Lumpfish (Cyclopterus lumpus) against Vibrio Anguillarum Does Not Influence Systemic Immunity. Vaccines (Basel). (2021) 9:819. doi: 10.3390/vaccines9080819
9. Erkinharju T, Dalmo RA, Hansen M, Seternes T. Cleaner fish in aquaculture: review on diseases and vaccination. Rev Aquac. (2021) 13:189–237. doi: 10.1111/raq.12470
10. Erkinharju T, Strandskog G, Vågnes Ø, Hordvik I, Dalmo RA, Seternes T. Intramuscular vaccination of Atlantic lumpfish (Cyclopterus lumpus L.) induces inflammatory reactions and local immunoglobulin M production at the vaccine administration site. J Fish Dis. (2019) 42:1731–43. doi: 10.1111/jfd.13101
11. Erkinharju T, Lundberg MR, Isdal E, Hordvik I, Dalmo RA, Seternes T. Studies on the antibody response and side effects after intramuscular and intraperitoneal injection of Atlantic lumpfish (Cyclopterus lumpus L.) with different oil-based vaccines. J Fish Dis. (2017) 40:1805–13. doi: 10.1111/jfd.12649
12. Erkinharju T, Dalmo RA, Vågsnes Ø, Hordvik I, Seternes T. Vaccination of Atlantic lumpfish (Cyclopterus lumpus L.) at a low temperature leads to a low antibody response against Aeromonas salmonicida. J Fish Dis. (2018) 41:613–23. doi: 10.1111/jfd.12760
13. Eslamloo K, Xue X, Hall JR, Smith NC, Caballero-Solares A, Parrish CC, et al. Transcriptome profiling of antiviral immune and dietary fatty acid dependent responses of Atlantic salmon macrophage-like cells. BMC Genomics. (2017) 18:706. doi: 10.1186/s12864-017-4099-2
14. Eslamloo K, Xue X, Booman M, Smith NC, Rise ML. Transcriptome profiling of the antiviral immune response in Atlantic cod macrophages. Dev Comp Immunol. (2016) 63:187–205. doi: 10.1016/j.dci.2016.05.021
15. Ruyra A, Torrealba D, Morera D, Tort L, MacKenzie S, Roher N. Zebrafish liver (ZFL) cells are able to mount an anti-viral response after stimulation with Poly (I:C). Comp Biochem Physiol B Biochem Mol Biol. (2015) 182:55–63. doi: 10.1016/j.cbpb.2014.12.002
16. Lulijwa R, Alfaro AC, Merien F, Burdass M, Meyer J, Venter L, et al. Metabolic and immune responses of Chinook salmon (Oncorhynchus tshawytscha) smolts to a short-term poly (I:C) challenge. J Fish Biol. (2020) 96:731–46. doi: 10.1111/jfb.14266
17. Jensen I, Albuquerque A, Sommer A-I, Robertsen B. Effect of poly I:C on the expression of Mx proteins and resistance against infection by infectious salmon anaemia virus in Atlantic salmon. Fish Shellfish Immunol. (2002) 13:311–26. doi: 10.1006/fsim.2001.0406
18. Andresen AMS, Gjøen T. Chitosan nanoparticle formulation attenuates poly (I:C) induced innate immune responses against inactivated virus vaccine in Atlantic salmon (Salmo salar). Comp Biochem Physiol D Genomics Proteomics. (2021) 40:100915. doi: 10.1016/j.cbd.2021.100915
19. Nishizawa T, Takami I, Yoshimizu M, Oh M-J. Required dose of fish nervous necrosis virus (NNV) for Poly(I:C) immunization of sevenband grouper Epinephelus septemfasciatus. Aquaculture. (2011) 311:100–4. doi: 10.1016/j.aquaculture.2010.12.009
20. Nishikiori M, den Boon JA, Unchwaniwala N, Ahlquist P. Crowning touches in positive-strand RNA virus genome replication complex structure and function. Annu Rev Virol. (2022) 9:193–212. doi: 10.1146/annurev-virology-092920-021307
21. Wu SX, Ahlquist P, Kaesberg P. Active complete in vitro replication of nodavirus RNA requires glycerophospholipid. Proc Natl Acad Sci U.S.A. (1992) 89:11136–40. doi: 10.1073/pnas.89.23.11136
22. Bandín I, Souto S. Betanodavirus and VER disease: A 30-year research review. Pathogens. (2020) 9:106. doi: 10.3390/pathogens9020106
23. Sandlund L, Mor SK, Singh VK, Padhi SK, Phelps NBD, Nylund S, et al. Comparative molecular characterization of novel and known piscine toti-like viruses. Viruses. (2021) 13:1063. doi: 10.3390/v13061063
24. Edwards M, Bignell JP, Papadopoulou A, Trani E, Savage J, Joseph AW, et al. First detection of Cyclopterus lumpus virus in England, following a mortality event in farmed cleaner fish. Bull EAFP. (2022) 43:28–37. doi: 10.48045/001c.56559
25. Toffan A, De Salvador M, Scholz F, Pretto T, Buratin A, Rodger HD, et al. Lumpfish (Cyclopterus lumpus, Linnaeus) is susceptible to viral nervous necrosis: Result of an experimental infection with different genotypes of Betanodavirus. J Fish Dis. (2019) 42:1667–76. doi: 10.1111/jfd.13088
26. Kukurba KR, Montgomery SB. RNA sequencing and analysis. Cold Spring Harb Protoc. (2015) 2015:951–69. doi: 10.1101/pdb.top084970
27. Aedo JE, Aravena-Canales D, Dettleff P, Fuentes-Valenzuela M, Zuloaga R, Rivas-Aravena A, et al. RNA-seq analysis reveals the dynamic regulation of proteasomal and autophagic degradation systems of rainbow trout (Oncorhynchus mykiss) skeletal muscle challenged with infectious pancreatic necrosis virus (IPNV). Aquaculture. (2022) 552:738000. doi: 10.1016/j.aquaculture.2022.738000
28. Sandamalika WMG, Liyanage DS, Lim C, Yang H, Lee S, Jeong T, et al. Differential gene expression of red-spotted grouper (Epinephelus akaara) in response to lipopolysaccharide, poly I:C, and nervous necrosis virus revealed by RNA-seq data. Fish Shellfish Immunol. (2022) 131:939–44. doi: 10.1016/j.fsi.2022.11.006
29. Liu Q-N, Tang Y-Y, Zhou M-J, Luo S, Li Y-T, Wang G, et al. Differentially expressed genes involved in immune pathways from yellowhead catfish (Tachysurus fulvidraco) after poly (I:C) challenge. Int J Biol Macromolecules. (2021) 183:340–5. doi: 10.1016/j.ijbiomac.2021.04.167
30. Du X, Li Y, Li D, Lian F, Yang S, Wu J, et al. Transcriptome profiling of spleen provides insights into the antiviral mechanism in Schizothorax prenanti after poly (I: C) challenge. Fish Shellfish Immunol. (2017) 62:13–23. doi: 10.1016/j.fsi.2017.01.004
31. Liu Y, Xin Z-Z, Zhang D-Z, Wang Z-F, Zhu X-Y, Tang B-P, et al. Transcriptome analysis of yellow catfish (Pelteobagrus fulvidraco) liver challenged with polyriboinosinic polyribocytidylic acid (poly I:C). Fish Shellfish Immunol. (2017) 68:395–403. doi: 10.1016/j.fsi.2017.07.030
32. Rao SS, Lunde HS, Dolan DWP, Fond AK, Petersen K, Haugland GT. Transcriptome-wide analyses of early immune responses in lumpfish leukocytes upon stimulation with poly(I:C). Front Immunol. (2023) 14:1198211. doi: 10.3389/fimmu.2023.1198211
33. Rao SS, Nelson PA, Lunde HS, Haugland GT. Evolutionary, comparative, and functional analyses of STATs and regulation of the JAK-STAT pathway in lumpfish upon bacterial and poly(I:C) exposure. Front Cell Infect Microbiol. (2023) 13:1252744. doi: 10.3389/fcimb.2023.1252744
34. Jefferies CA. Regulating IRFs in IFN driven disease. Front Immunol. (2019) 10:325. doi: 10.3389/fimmu.2019.00325
35. Yan R, van Meurs M, Popa ER, Jongman RM, Zwiers PJ, Niemarkt AE, et al. Endothelial interferon regulatory factor 1 regulates lipopolysaccharide-induced VCAM-1 expression independent of NFκB. J Innate Immun. (2017) 9:546–60. doi: 10.1159/000477211
36. Lai CF, Wang T-Y, Yeh M-I, Chen T-Y. Characterization of orange-spotted grouper (Epinephelus coioides) interferon regulatory factor 4 regulated by heat shock factor 1 during heat stress in response to antiviral immunity. Fish Shellfish Immunol. (2020) 106:755–67. doi: 10.1016/j.fsi.2020.08.033
37. Han C, Huang W, Peng S, Zhou J, Zhan H, Li W, et al. Characterization and expression analysis of the interferon regulatory factor (IRF) gene family in zig-zag eel (Mastacembelus armatus) against Aeromonas veronii infection. Dev Comp Immunol. (2023) 140:104622. doi: 10.1016/j.dci.2022.104622
38. Yanai H, Negishi H, Taniguchi T. The IRF family of transcription factors. Oncoimmunology. (2012) 1:1376–86. doi: 10.4161/onci.22475
39. Li W, Zhao G, Jiao Z, Xiang C, Liang Y, Huang W, et al. Nuclear import of IRF11 via the importin α/β pathway is essential for its antiviral activity. Dev Comp Immunol. (2023) 141:104649. doi: 10.1016/j.dci.2023.104649
40. Langevin C, Aleksejeva E, Passoni G, Palha N, Levraud J-P, Boudinot P. The antiviral innate immune response in fish: evolution and conservation of the IFN system. Journal of Molecular Biology (2013) 425:4904–20. doi: 10.1016/j.jmb.2013.09.033
41. Zhu Y, Qi C, Shan S, Zhang F, Li H, An L, et al. Characterization of common carp (Cyprinus carpio L.) interferon regulatory factor 5 (IRF5) and its expression in response to viral and bacterial challenges. BMC Vet Res. (2016) 12:127. doi: 10.1186/s12917-016-0750-4
42. Hu G-B, Lou H-M, Dong X-Z, Liu Q-M, Zhang S-C. Characteristics of the interferon regulatory factor 5 (IRF5) and its expression in response to LCDV and poly I:C challenges in Japanese flounder, Paralichthys olivaceus. Dev Comp Immunol. (2012) 38:377–82. doi: 10.1016/j.dci.2012.06.001
43. 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
44. Hu G, Xia J, Lou H, Liu Q, Lin J, Yin X, et al. Cloning and expression analysis of interferon regulatory factor 7 (IRF-7) in turbot, Scophthalmus maximus. Dev Comp Immunol. (2011) 35:416–20. doi: 10.1016/j.dci.2010.12.004
45. Dong X, Xu H, Mai K, Xu W, Zhang Y, Ai Q. Cloning and characterization of SREBP-1 and PPAR-α in Japanese seabass Lateolabrax japonicus, and their gene expressions in response to different dietary fatty acid profiles. Comp Biochem Physiol B Biochem Mol Biol. (2015) 180:48–56. doi: 10.1016/j.cbpb.2014.10.001
46. Inkpen SM, Solbakken MH, Jentoft S, Eslamloo K, Rise ML. Full characterization and transcript expression profiling of the interferon regulatory factor (IRF) gene family in Atlantic cod (Gadus morhua). Dev Comp Immunol. (2019) 98:166–80. doi: 10.1016/j.dci.2019.03.015
47. Inkpen SM, Hori TS, Gamperl AK, Nash GW, Rise ML. Characterization and expression analyses of five interferon regulatory factor transcripts (Irf4a, Irf4b, Irf7, Irf8, Irf10) in Atlantic cod (Gadus morhua). Fish Shellfish Immunol. (2015) 44:365–81. doi: 10.1016/j.fsi.2015.02.032
48. Deering MJ, Paradis H, Ahmad R, Al-Mehiawi AS, Gendron RL. The role of dietary vitamin A in mechanisms of cataract development in the teleost lumpfish (Cyclopterus lumpus L). J Fish Dis. doi: 10.1111/jfd.13899
49. Emam M, Eslamloo K, Caballero-Solares A, Lorenz EK, Xue X, Umasuthan N, et al. Nutritional immunomodulation of Atlantic salmon response to Renibacterium salmoninarum bacterin. Front Mol Biosci. (2022) 9:931548. doi: 10.3389/fmolb.2022.931548
50. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics. (2009) 25:1091–3. doi: 10.1093/bioinformatics/btp101
51. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. (2003) 13:2498–504. doi: 10.1101/gr.1239303
52. Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. (2020) 13:1194–202. doi: 10.1016/j.molp.2020.06.009
53. Koressaar T, Remm M. Enhancements and modifications of primer design program Primer3. Bioinformatics. (2007) 23:1289–91. doi: 10.1093/bioinformatics/btm091
54. Kõressaar T, Lepamets M, Kaplinski L, Raime K, Andreson R, Remm M. Primer3_masker: integrating masking of template sequence with primer design software. Bioinformatics. (2018) 34:1937–8. doi: 10.1093/bioinformatics/bty036
55. Caballero-Solares A, Hall JR, Xue X, Eslamloo K, Taylor RG, Parrish CC, et al. The dietary replacement of marine ingredients by terrestrial animal and plant alternatives modulates the antiviral immune response of Atlantic salmon (Salmo salar). Fish Shellfish Immunol. (2017) 64:24–38. doi: 10.1016/j.fsi.2017.02.040
56. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. (2002) 3:research0034.1. doi: 10.1186/gb-2002-3-7-research0034
57. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. (2009) 55:611–22. doi: 10.1373/clinchem.2008.112797
58. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCT method. Methods. (2001) 25:402–8. doi: 10.1006/meth.2001.1262
59. Pfaffl MW. A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. (2001) 29:e45–5. doi: 10.1093/nar/29.9.e45
60. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. (2021) 38:3022–7. doi: 10.1093/molbev/msab120
61. Nehyba J, Hrdlicková R, Bose HR. Dynamic evolution of immune system regulators: the history of the interferon regulatory factor family. Mol Biol Evol. (2009) 26:2539–50. doi: 10.1093/molbev/msp167
62. Zhou Y, Pu J, Wu Y. The role of lipid metabolism in influenza A virus infection. Pathogens. (2021) 10:303. doi: 10.3390/pathogens10030303
63. Farías MA, Diethelm-Varela B, Navarro AJ, Kalergis AM, González PA. Interplay between lipid metabolism, lipid droplets, and DNA virus infections. Cells. (2022) 11:2224. doi: 10.3390/cells11142224
64. Sarohan AR. COVID-19: endogenous retinoic acid theory and retinoic acid depletion syndrome. Med Hypotheses. (2020) 144:110250. doi: 10.1016/j.mehy.2020.110250
65. Langevin C, Aleksejeva E, Houel A, Briolat V, Torhy C, Lunazzi A, et al. FTR83, a member of the large fish-specific finTRIM family, triggers IFN pathway and counters viral infection. Front Immunol. (2017) 8:617. doi: 10.3389/fimmu.2017.00617
66. Qin Z, Ren F, Xu X, Ren Y, Li H, Wang Y, et al. ZNF536, a novel zinc finger protein specifically expressed in the brain, negatively regulates neuron differentiation by repressing retinoic acid-induced gene transcription. Mol Cell Biol. (2009) 29:3633–43. doi: 10.1128/MCB.00362-09
67. Bin L, Deng L, Yang H, Zhu L, Wang X, Edwards MG, et al. Forkhead box C1 regulates human primary keratinocyte terminal differentiation. PloS One. (2016) 11:e0167392. doi: 10.1371/journal.pone.0167392
68. McConnell BB, Yang VW. Mammalian krüppel-like factors in health and diseases. Physiol Rev. (2010) 90:1337–81. doi: 10.1152/physrev.00058.2009
69. Yang Q, Xie Y, Alexson SEH, Dean Nelson B, DePierre JW. Involvement of the peroxisome proliferator-activated receptor alpha in the immunomodulation caused by peroxisome proliferators in mice. Biochem Pharmacol. (2002) 63:1893–900. doi: 10.1016/S0006-2952(02)00923-1
70. Li S, Chen X, Hao G, Geng X, Zhan W, Sun J. Identification and characterization of a novel NOD-like receptor family CARD domain containing 3 gene in response to extracellular ATP stimulation and its role in regulating LPS-induced innate immune response in Japanese flounder (Paralichthys olivaceus) head kidney macrophages. Fish Shellfish Immunol. (2016) 50:79–90. doi: 10.1016/j.fsi.2016.01.029
71. Deng H, Yue JK, Zusman BE, Nwachuku EL, Abou-Al-Shaar H, Upadhyayula PS, et al. B-cell lymphoma 2 (Bcl-2) and regulation of apoptosis after traumatic brain injury: A clinical perspective. Medicina (Kaunas). (2020) 56:300. doi: 10.3390/medicina56060300
72. Jármay K, Karácsony G, Nagy A, Schaff Z. Changes in lipid metabolism in chronic hepatitis C. World J Gastroenterol. (2005) 11:6422–8. doi: 10.3748/wjg.v11.i41.6422
73. Grunfeld C, Feingold KR. Regulation of lipid metabolism bycytokines during host defense. Nutrition. (1996) 12:S24–6. doi: 10.1016/0899-9007(95)00073-9
74. Abu-Farha M, Thanaraj TA, Qaddoumi MG, Hashem A, Abubaker J, Al-Mulla F. The role of lipid metabolism in COVID-19 virus infection and as a drug target. Int J Mol Sci. (2020) 21:3544. doi: 10.3390/ijms21103544
75. Ren Z, Ding T, Zuo Z, Xu Z, Deng J, Wei Z. Regulation of MAVS expression and signaling function in the antiviral innate immune response. Front Immunol. (2020) 11:1030. doi: 10.3389/fimmu.2020.01030
76. Magoro T, Dandekar A, Jennelle LT, Bajaj R, Lipkowitz G, Angelucci AR, et al. IL-1β/TNF-α/IL-6 inflammatory cytokines promote STAT1-dependent induction of CH25H in Zika virus–infected human macrophages. J Biol Chem. (2019) 294:14591–602. doi: 10.1074/jbc.RA119.007555
77. Wang J, Huang CL-H, Zhang Y. Complement C1q binding protein (C1QBP): physiological functions, mutation-associated mitochondrial cardiomyopathy and current disease models. Front Cardiovasc Med. (2022) 9:843853. doi: 10.3389/fcvm.2022.843853
78. Singh R, Letai A, Sarosiek K. Regulation of apoptosis in health and disease: the balancing act of BCL-2 family proteins. Nat Rev Mol Cell Biol. (2019) 20:175–93. doi: 10.1038/s41580-018-0089-8
79. Behrendorff N, Dolai S, Hong W, Gaisano HY, Thorn P. Vesicle-associated membrane protein 8 (VAMP8) is a SNARE (Soluble N-ethylmaleimide-sensitive factor attachment protein receptor) selectively required for sequential granule-to-granule fusion. J Biol Chem. (2011) 286:29627–34. doi: 10.1074/jbc.M111.265199
80. Sebti S, Prébois C, Pérez-Gracia E, Bauvy C, Desmots F, Pirot N, et al. BAT3 modulates p300-dependent acetylation of p53 and autophagy-related protein 7 (ATG7) during autophagy. Proc Natl Acad Sci USA. (2014) 111:4115–20. doi: 10.1073/pnas.1313618111
81. Calle X, Garrido-Moreno V, Lopez-Gallardo E, Norambuena-Soto I, Martínez D, Peñaloza-Otárola A, et al. Mitochondrial E3 ubiquitin ligase 1 (MUL1) as a novel therapeutic target for diseases associated with mitochondrial dysfunction. IUBMB Life. (2022) 74:850–65. doi: 10.1002/iub.2657
82. Carty M, Bowie AG. Recent insights into the role of Toll-like receptors in viral infection. Clin Exp Immunol. (2010) 161:397–406. doi: 10.1111/j.1365-2249.2010.04196.x
83. Ali S, Mann-Nüttel R, Schulze A, Richter L, Alferink J, Scheu S. Sources of type I interferons in infectious immunity: plasmacytoid dendritic cells not always in the driver’s seat. Front Immunol. (2019) 10:778. doi: 10.3389/fimmu.2019.00778
84. Kawai T, Sato S, Ishii KJ, Coban C, Hemmi H, Yamamoto M, et al. Interferon-α induction through Toll-like receptors involves a direct interaction of IRF7 with MyD88 and TRAF6. Nat Immunol. (2004) 5:1061–8. doi: 10.1038/ni1118
85. Blasius AL, Beutler B. Intracellular toll-like receptors. Immunity. (2010) 32:305–15. doi: 10.1016/j.immuni.2010.03.012
86. Su C, Tang Y, Zheng C. DExD/H-box helicases: multifunctional regulators in antiviral innate immunity. Cell Mol Life Sci. (2022) 79:2. doi: 10.1007/s00018-021-04072-6
87. Rise ML, Hall J, Rise M, Hori T, Gamperl AK, Kimball J, et al. Functional genomic analysis of the response of Atlantic cod (Gadus morhua) spleen to the viral mimic polyriboinosinic polyribocytidylic acid (pIC). Dev Comp Immunol. (2008) 32:916–31. doi: 10.1016/j.dci.2008.01.002
88. Rise ML, Hall JR, Rise M, Hori TS, Browne MJ, Gamperl AK, et al. Impact of asymptomatic nodavirus carrier state and intraperitoneal viral mimic injection on brain transcript expression in Atlantic cod (Gadus morhua). Physiol Genomics. (2010) 42:266–80. doi: 10.1152/physiolgenomics.00168.2009
89. Feng CY, Rise ML. Identification and molecular cloning of Atlantic cod (Gadus morhua) activating transcription factor 3 (ATF3) transcript and its induction in spleen following intraperitoneal polyriboinosinic polyribocytidylic acid injection. Fish Shellfish Immunol. (2011) 31:475–81. doi: 10.1016/j.fsi.2011.06.002
90. Shu M, Du T, Zhou G, Roizman B. Role of activating transcription factor 3 in the synthesis of latency-associated transcript and maintenance of herpes simplex virus 1 in latent state in ganglia. Proc Natl Acad Sci U.S.A. (2015) 112:E5420–5426. doi: 10.1073/pnas.1515369112
91. Sood V, Sharma KB, Gupta V, Saha D, Dhapola P, Sharma M, et al. ATF3 negatively regulates cellular antiviral signaling and autophagy in the absence of type I interferons. Sci Rep. (2017) 7:8789. doi: 10.1038/s41598-017-08584-9
92. Rosenberger CM, Clark AE, Treuting PM, Johnson CD, Aderem A. ATF3 regulates MCMV infection in mice by modulating IFN-γ expression in natural killer cells. Proc Natl Acad Sci. (2008) 105:2544–9. doi: 10.1073/pnas.0712182105
93. Hai T, Wolford CC, Chang Y-S. ATF3, a hub of the cellular adaptive-response network, in the pathogenesis of diseases: is modulation of inflammation a unifying component? Gene Expr. (2010) 15:1–11. doi: 10.3727/105221610X12819686555015
94. Labzin LI, Schmidt SV, Masters SL, Beyer M, Krebs W, Klee K, et al. ATF3 is a key regulator of macrophage IFN responses. J Immunol. (2015) 195:4446–55. doi: 10.4049/jimmunol.1500204
95. Polinski MP, Bradshaw JC, Rise ML, Johnson SC, Garver KA. Sockeye salmon demonstrate robust yet distinct transcriptomic kidney responses to rhabdovirus (IHNV) exposure and infection. Fish Shellfish Immunol. (2019) 94:525–38. doi: 10.1016/j.fsi.2019.09.042
96. Polinski MP, Bradshaw JC, Inkpen SM, Richard J, Fritsvold C, Poppe TT, et al. De novo assembly of Sockeye salmon kidney transcriptomes reveal a limited early response to piscine reovirus with or without infectious hematopoietic necrosis virus superinfection. BMC Genomics. (2016) 17:848. doi: 10.1186/s12864-016-3196-y
97. Takeuchi O, Akira S. Pattern recognition receptors and inflammation. Cell. (2010) 140:805–20. doi: 10.1016/j.cell.2010.01.022
98. McNab F, Mayer-Barber K, Sher A, Wack A, O’Garra A. Type I interferons in infectious disease. Nat Rev Immunol. (2015) 15:87–103. doi: 10.1038/nri3787
99. Xu C, Guo T-C, Mutoloki S, Haugland Ø, Marjara IS, Evensen Ø. Alpha interferon and not gamma interferon inhibits salmonid alphavirus subtype 3 replication in vitro. J Virol. (2010) 84:8903. doi: 10.1128/JVI.00851-10
100. Samsing F, Alexandre P, Rigby M, Taylor RS, Chong R, Wynne JW. Transcriptome response of atlantic salmon (Salmo salar) to a new piscine orthomyxovirus. Pathogens. (2020) 9:807. doi: 10.3390/pathogens9100807
101. Workenhe ST, Hori TS, Rise ML, Kibenge MJT, Kibenge FSB. Infectious salmon anaemia virus (ISAV) isolates induce distinct gene expression responses in the Atlantic salmon (Salmo salar) macrophage/dendritic-like cell line TO, assessed using genomic techniques. Mol Immunol. (2009) 46:2955–74. doi: 10.1016/j.molimm.2009.06.015
102. Valenzuela-Miranda D, Boltaña S, Cabrejos ME, Yáñez JM, Gallardo-Escárate C. High-throughput transcriptome analysis of ISAV-infected Atlantic salmon Salmo salar unravels divergent immune responses associated to head-kidney, liver and gills tissues. Fish Shellfish Immunol. (2015) 45:367–77. doi: 10.1016/j.fsi.2015.04.003
103. Wei X, Li XZ, Zheng X, Jia P, Wang J, Yang X, et al. Toll-like receptors and interferon associated immune factors responses to spring viraemia of carp virus infection in common carp (Cyprinus carpio). Fish Shellfish Immunol. (2016) 55:568–76. doi: 10.1016/j.fsi.2016.05.043
104. Saint-Jean SR, Pérez-Prieto SI. Effects of salmonid fish viruses on Mx gene expression and resistance to single or dual viral infections. Fish Shellfish Immunol. (2007) 23:390–400. doi: 10.1016/j.fsi.2006.11.012
105. Zhang X, Yang F, Li K, Cao W, Ru Y, Chen S, et al. The insufficient activation of RIG-I–like signaling pathway contributes to highly efficient replication of porcine picornaviruses in IBRS-2 cells. Mol Cell Proteomics : MCP. (2021) 20:100147. doi: 10.1016/j.mcpro.2021.100147
106. Rustagi A, Gale M. Innate antiviral immune signaling, viral evasion and modulation by HIV-1. J Mol Biol. (2014) 426:1161–77. doi: 10.1016/j.jmb.2013.12.003
107. Eslamloo K, Caballero-Solares A, Inkpen SM, Emam M, Kumar S, Bouniot C, et al. Transcriptomic profiling of the adaptive and innate immune responses of Atlantic salmon to Renibacterium salmoninarum infection. Front Immunol. (2020) 11:567838. doi: 10.3389/fimmu.2020.567838
108. Xue X, Caballero-Solares A, Hall JR, Umasuthan N, Kumar S, Jakob E, et al. Transcriptome profiling of Atlantic salmon (Salmo salar) parr with higher and lower pathogen loads following Piscirickettsia salmonis infection. Front Immunol. (2021) 12:789465. doi: 10.3389/fimmu.2021.789465
109. Adamek M, Davies J, Beck A, Jordan L, Becker AM, Mojzesz M, et al. Antiviral actions of 25-hydroxycholesterol in fish vary with the virus-host combination. Front Immunol. (2021) 12:581786. doi: 10.3389/fimmu.2021.581786
110. Zhang Y, Wang L, Huang X, Wang S, Huang Y, Qin Q. Fish cholesterol 25-hydroxylase inhibits virus replication via regulating interferon immune response or affecting virus entry. Front Immunol. (2019) 10:322. doi: 10.3389/fimmu.2019.00322
111. van Gent M, Sparrer KMJ, Gack MU. TRIM proteins and their roles in antiviral host defenses. Annu Rev Virol. (2018) 5:385–405. doi: 10.1146/annurev-virology-092917-043323
112. Yu Y, Huang X, Zhang J, Liu J, Hu Y, Yang Y, et al. Fish TRIM16L exerts negative regulation on antiviral immune response against grouper iridoviruses. Fish Shellfish Immunol. (2016) 59:256–67. doi: 10.1016/j.fsi.2016.10.044
113. Cho JY, Kim J, Kim J-W, Lee D, Kim D-G, Kim Y-S, et al. Characterization of TRIM16, a member of the fish-specific finTRIM family, in olive flounder Paralichthys olivaceus. Fish Shellfish Immunol. (2022) 127:666–71. doi: 10.1016/j.fsi.2022.07.003
114. Choudhury NR, Trus I, Heikel G, Wolczyk M, Szymanski J, Bolembach A, et al. TRIM25 inhibits influenza A virus infection, destabilizes viral mRNA, but is redundant for activating the RIG-I pathway. Nucleic Acids Res. (2022) 50:7097–114. doi: 10.1093/nar/gkac512
115. Marschall M, Strojan H, Kiener R, Wangen C, Sonntag E, Müller R, et al. Differential upregulation of host cell protein kinases by the replication of α-, β- and γ-herpesviruses provides a signature of virus-specific signalling. J Gen Virol. (2020) 101:284–9. doi: 10.1099/jgv.0.001370
116. Müller A, Sutherland BJG, Koop BF, Johnson SC, Garver KA. Infectious hematopoietic necrosis virus (IHNV) persistence in Sockeye Salmon: influence on brain transcriptome and subsequent response to the viral mimic poly(I:C). BMC Genomics. (2015) 16:634. doi: 10.1186/s12864-015-1759-y
117. Eslamloo K, Kumar S, Xue X, Parrish KS, Purcell SL, Fast MD, et al. Global gene expression responses of Atlantic salmon skin to Moritella viscosa. Sci Rep. (2022) 12:4622. doi: 10.1038/s41598-022-08341-7
118. Sommer F, Torraca V, Kamel SM, Lombardi A, Meijer AH. Frontline Science: Antagonism between regular and atypical Cxcr3 receptors regulates macrophage migration during infection and injury in zebrafish. J Leukocyte Biol. (2020) 107:185–203. doi: 10.1002/JLB.2HI0119-006R
119. Tomczyk M, Glaser T, Slominska EM, Ulrich H, Smolenski RT. Purine nucleotides metabolism and signaling in huntington’s disease: search for a target for novel therapies. Int J Mol Sci. (2021) 22:6545. doi: 10.3390/ijms22126545
120. Zhou M, Yao Z, Zhao M, Fang Q, Ji X, Chen H, et al. Molecular cloning and expression responses of jarid2b to high-temperature treatment in nile tilapia (Oreochromis niloticus). Genes (Basel). (2022) 13:1719. doi: 10.3390/genes13101719
121. Loenarz C, Ge W, Coleman ML, Rose NR, Cooper CDO, Klose RJ, et al. PHF8, a gene associated with cleft lip/palate and mental retardation, encodes for an Nϵ-dimethyl lysine demethylase. Hum Mol Genet. (2010) 19:217–22. doi: 10.1093/hmg/ddp480
122. Tsai K, Cullen BR. Epigenetic and epitranscriptomic regulation of viral replication. Nat Rev Microbiol. (2020) 18:559–70. doi: 10.1038/s41579-020-0382-3
123. Jahn R, Scheller RH. SNAREs — engines for membrane fusion. Nat Rev Mol Cell Biol. (2006) 7:631–43. doi: 10.1038/nrm2002
124. Verboogen DRJ, González Mancha N, ter Beest M, van den Bogaart G. Fluorescence Lifetime Imaging Microscopy reveals rerouting of SNARE trafficking driving dendritic cell activation. eLife. (2017) 6:e23525. doi: 10.7554/eLife.23525
125. Clark TC, Boudinot P, Collet B. Evolution of the IRF family in salmonids. Genes (Basel). (2021) 12:238. doi: 10.3390/genes12020238
126. Bergan V, Kileng Ø, Sun B, Robertsen B. Regulation and function of interferon regulatory factors of Atlantic salmon. Mol Immunol. (2010) 47:2005–14. doi: 10.1016/j.molimm.2010.04.015
127. Sun F, Zhang Y-B, Liu T-K, Gan L, Yu F-F, Liu Y, et al. Characterization of fish IRF3 as an IFN-inducible protein reveals evolving regulation of IFN response in vertebrates. J Immunol. (2010) 185:7573–82. doi: 10.4049/jimmunol.1002401
128. Ai K, Luo K, Xia L, Gao W, Hu W, Qi Z, et al. Functional characterization of interferon regulatory factor 5 and its role in the innate antiviral immune response. Fish Shellfish Immunol. (2018) 72:31–6. doi: 10.1016/j.fsi.2017.10.042
129. Chen X, Guan Y, Li K, Luo T, Mu Y, Chen X. IRF1 and IRF2 act as positive regulators in antiviral response of large yellow croaker (Larimichthys crocea) by induction of distinct subgroups of type I IFNs. Dev Comp Immunol. (2021) 118:103996. doi: 10.1016/j.dci.2021.103996
130. Romeo G, Fiorucci G, Chiantore MV, Percario ZA, Vannucchi S, Affabris E. Review: IRF-1 as a negative regulator of cell proliferation. J Interferon Cytokine Res. (2002) 22:39–47. doi: 10.1089/107999002753452647
131. Bonjardim CA, Ferreira PCP, Kroon EG. Interferons: signaling, antiviral and viral evasion. Immunol Lett. (2009) 122:1–11. doi: 10.1016/j.imlet.2008.11.002
132. Adams NM, Lau CM, Fan X, Rapp M, Geary CD, Weizman O-E, et al. Transcription factor IRF8 orchestrates the adaptive natural killer cell response. Immunity. (2018) 48:1172–1182.e6. doi: 10.1016/j.immuni.2018.04.018
133. Paul A, Tang TH, Ng SK. Interferon regulatory factor 9 structure and regulation. Front Immunol. (2018) 9:1831. doi: 10.3389/fimmu.2018.01831
134. An L-L, Zhao X, Gong X-Y, Li Y-L, Qu Z-L, Sun H-Y, et al. Promoter binding and nuclear retention features of zebrafish IRF family members in IFN response. Front Immunol. (2022) 13:861262. doi: 10.3389/fimmu.2022.861262
135. Huang B, Qi ZT, Xu Z, Nie P. Global characterization of interferon regulatory factor (IRF) genes in vertebrates: Glimpse of the diversification in evolution. BMC Immunol. (2010) 11:22. doi: 10.1186/1471-2172-11-22
136. Stein C, Caccamo M, Laird G, Leptin M. Conservation and divergence of gene families encoding components of innate immune response systems in zebrafish. Genome Biol. (2007) 8:R251. doi: 10.1186/gb-2007-8-11-r251
137. Tussiwand R, Everts B, Grajales-Reyes GE, Kretzer NM, Iwata A, Bagaitkar J, et al. Klf4 expression in conventional dendritic cells is required for T helper 2 cell responses. Immunity. (2015) 42:916–28. doi: 10.1016/j.immuni.2015.04.017
138. Rengarajan J, Mowen KA, McBride KD, Smith ED, Singh H, Glimcher LH. Interferon regulatory factor 4 (IRF4) interacts with NFATc2 to modulate interleukin 4 gene expression. J Exp Med. (2002) 195:1003–12. doi: 10.1084/jem.20011128
139. Klein U, Casola S, Cattoretti G, Shen Q, Lia M, Mo T, et al. Transcription factor IRF4 controls plasma cell differentiation and class-switch recombination. Nat Immunol. (2006) 7:773–82. doi: 10.1038/ni1357
Keywords: antiviral, lumpfish (Cyclopterus lumpus), IRF, RNA seq, qPCR
Citation: Emam M, Kumar S, Eslamloo K, Caballero-Solares A, Hall JR, Xue X, Paradis H, Gendron RL, Santander J and Rise ML (2024) Transcriptomic response of lumpfish (Cyclopterus lumpus) head kidney to viral mimic, with a focus on the interferon regulatory factor family. Front. Immunol. 15:1439465. doi: 10.3389/fimmu.2024.1439465
Received: 28 May 2024; Accepted: 08 July 2024;
Published: 15 August 2024.
Edited by:
Iddya Karunasagar, Nitte University, IndiaReviewed by:
Rahul Krishnan, Kerala University of Fisheries and Ocean Studies, IndiaXiuzhen Sheng, Ocean University of China, China
Copyright © 2024 Emam, Kumar, Eslamloo, Caballero-Solares, Hall, Xue, Paradis, Gendron, Santander and Rise. 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: Mohamed Emam, melsayedemam@mun.ca; Matthew L. Rise, mrise@mun.ca
†These authors have contributed equally to this work