Gene Set Enrichment Analysis Reveals That Fucoidan Induces Type I IFN Pathways in BMDC

Fucoidan, a sulfated polysaccharide extracted from brown seaweed, has been proposed to effectively treat and prevent various viral infections. However, the mechanisms behind its antiviral activity are not completely understood. We investigate here the global transcriptional changes in bone marrow-derived dendritic cells (BMDCs) using RNA-Seq technology. Through both analysis of differentially expressed genes (DEG) and gene set enrichment analysis (GSEA), we found that fucoidan-treated BMDCs were enriched in virus-specific response pathways, including that of SARS-CoV-2, as well as pathways associated with nucleic acid-sensing receptors (RLR, TLR, NLR, STING), and type I interferon (IFN) production. We show that these transcriptome changes are driven by well-known regulators of the inflammatory response against viruses, including IRF, NF-κB, and STAT family transcription factors. Furthermore, 435 of the 950 upregulated DEGs are classified as type I IFN-stimulated genes (ISGs). Flow cytometric analysis additionally showed that fucoidan increased MHCII, CD80, and CD40 surface markers in BMDCs, indicative of greater antigen presentation and co-stimulation functionality. Our current study suggests that fucoidan transcriptionally activates PRR signaling, type I IFN production and signaling, ISGs production, and DC maturation, highlighting a potential mechanism of fucoidan-induced antiviral activity.


Introduction
Fucoidan is a sulfated polysaccharide found in various species of brown seaweed. It is largely composed of L-fucose and sulfate ester groups, but can contain other monosaccharides (mannose, galactose, glucose, xylose, etc.), and uronic acids depending on the extraction process and seaweed species [1]. Fucoidan has been reported to possess a variety of biological activities, including antiviral, anti-inflammatory, anticoagulant, and antitumor effects, making it an attractive target for the development of nutraceuticals and functional foods [2][3][4][5][6]. While fucoidan has garnered particular interest for its antiviral properties, further investigation is required to understand its underlying molecular and cellular mechanisms.
Rapid and robust synthesis of type I IFNs constitute the primary protective response to most viral infections [29]. Specifically, IFNs activate the JAK-STAT pathway to drive the expression of hundreds of IFN-stimulated genes (ISGs), which in turn directly and indirectly exert antiviral effects, as well as promote antigen presentation and adaptive immunity activation [6,30]. Accumulating evidence indicates that SARS-CoV-2 targets the host's type I interferon (IFN) response at multiple steps, including viral sensing and IFN signaling [31][32][33]. Viruses are sensed by pattern recognition receptors (PRRs) such as toll-like receptors (TLRs), nucleotide-binding oligomerization domain-like receptors (NLRs), retinoic acid-inducible gene-I-like receptors (RLRs), and C-type lectin receptors (CLRs) [34,35]. By inhibiting PRR signaling, viruses such as SARS-CoV-2 may evade detection and attenuate IFN production [32,36,37]. They may also directly inhibit IFN production and signaling, as suggested by patients with mild coronavirus disease 2019 (COVID- 19) demonstrating early and strong IFN responses whereas those with severe disease showing a deficiency in IFN synthesis [38,39].
One of the primary cells indicated in early viral detection and IFN production are dendritic cells [40,41]. Plasmacytoid DCs (pDCs) in particular are able to produce large amounts of type I IFN in response to viral infection, far exceeding those produced by other cells [42]. Indeed, patients with severe COVID-19 displayed low levels of circulating pDCs with a pro-inflammatory phenotype, suggesting a diminished ability to produce type I IFN [42][43][44][45][46][47]. In addition, DCs bridge the innate and adaptive antiviral immune responses via their role as professional antigen presenting cells (APCs). Specifically, conventional DCs (cDCs) capture viral pathogens for presentation to naïve T cells, inducing cytotoxic T-cell activity and antibody production that are critical to virus eradication [48,49]. cDCs were similarly reported to express decreased levels of MHC-II in COVID-19 patients, which would lead to impaired adaptive immune system activation [50,51].
In the present study, we analyzed the effects of fucoidan treatment on the transcriptome of bone-marrow derived DCs (BMDCs) given the importance of DCs in the antiviral response. Our findings illustrate that fucoidan promotes viral sensing, IFN signaling, and adaptive immune activation in BMDCs, highlighting its potential for further investigation against SARS-CoV-2 in in vivo studies.

Bone Marrow-Derived Dendritic Cell (BMDC) Culture
Bone marrow cells were obtained from the femur and tibia of 6-to 10-week-old female C57BL/6 mice (Damulscience, Korea) and differentiated to dendritic cells as previously described [52]. Briefly, bone marrow cells were cultured in RPMI medium containing 10% FBS, 1% penicillin/streptomycin, 50 µM 2-Mercaptoethanol, and 20 ng/mL GM-CSF for 3 days. The media was replaced with fresh supplemented media and cultured for an additional 3 days. Nonadherent immature dendritic cells were harvested on day 6, treated with 400 µg /mL fucoidan    [53,54] or vehicle (0.04% DMSO) in complete media for 24 h, and subjected to RNA extraction for RNA-seq analysis. For flow cytometry analysis, immature dendritic cells were treated with 100 to 400 µg /mL fucoidan or vehicle in complete media for 24 h.

RNA Sequencing
Total RNA was isolated using the Trizol reagent (Thermo) and RNeasy kit (Qiagen) and the purified mRNAs were used to synthesize double-stranded cDNAs using a SuperScript Double-Stranded cDNA Synthesis Kit (Thermo Fisher Scientific). Sequencing libraries were constructed using a NEBNext Ultra II DNA Library Prep Kit for Illumina (NEB) and sequenced on the Illumina HiSeq 2500. Raw fastq files were quality-checked with FastQC (v0.11.9), and subjected to trimming with 'Trim galore!' (v0.6.4) of the Python tool (FastQC score Q30 over) to remove adapters and low-quality reads. For subsequent mapping, mouse reference genome (mm10) data were obtained from Illumina iGenomes (https://sapac. support.illumina.com/sequencing/sequencing_software/igenome.html). Sequences were then aligned to the mouse reference genome (mm10) with STAR (v2.7.3a). For the mapped data, low read counts and duplicate reads were removed, and the total read count, counts per million (cpm), and log FC were obtained.

Analysis of Differentially Expressed Genes in IDEP
Gene expression read counts were exported and analyzed through the integrated Differential Expression and Pathway (iDEP) tool (http://bioinformatics.sdstate.edu/idep9 3/). An MA plot was created to identify differentially expressed genes (DEGs) using the DESeq2 method with a threshold false discovery rate (FDR) < 0.1 and fold-change > 2. To characterize the biological function of the DEGs, KEGG pathway enrichment analysis was performed and visualized on a hierarchical tree and network. In addition, TF motifs enriched in gene promoters defined as 300 bp upstream of DEGs were identified and visualized as a hierarchical tree and network.

GSEA
The Broad Institute's Gene Set Enrichment Analysis (GSEA) software (version 4.1.0) was used with the C2 gene sets in the Molecular Signatures Database (MSigDB). A normalized enrichment score (NES) was also calculated by GSEA in which differences in pathway size (i.e., gene set size) were considered, allowing for comparisons between pathways within the analysis. This ES reflects the degree to which a gene set is overrepresented at the top or bottom of a ranked list of genes. Positive and negative NES values represent enrichment at the top and bottom of the list, respectively. The analysis involved 1000 gene set permutations with gene sets limited to 15-500. Initially, a FDR of 25% was used for all bioinformatics analyses.
Leading-edge analysis was performed after each GSEA to determine the core genes that have the highest impact on the biological process of a given gene-set. GSEA was run on a database of C2, and all canonical pathways and the enriched gene sets in KEGG, Reactome, and WikiPathways (WP) (p < 0.05; FDR < 0.25) were selected for leading-edge analysis.
To understand the role of TFs in pathway perturbations after fucoidan treatment, TF enrichment analysis was performed with GSEA using the TFT_Legacy subset of C3: regulatory target gene sets (p < 0.05; FDR < 0.25). In addition, GSEA leading-edge analysis of the enriched gene sets of potential TF targets were performed.

Functional Enrichment Visualization and Annotation of GSEA Results
GSEA was performed using the 'CP: Canonical Pathway' gene sets in MSigDB and the GSEA results were visualized using the EnrichmentMap (version 3.3.2) module for Cytoscape (version 3.8.2). The enrichment results were mapped as a network of gene-sets (nodes) where the nodes represent statistically significant terms (Node cutoff, Q-value 0.08) and the links (edges) represent the degree of gene-set similarity (Edge cutoff, 0.25). For gene-sets annotations, the enriched gene-sets were grouped by AutoAnnotate software (version 1.3.4, http://apps.cytoscape.org/apps/autoannotate) according to the Markov cluster (MCL) algorithm based on the edge weights of similarity coefficients and were automatically annotated using the WordCloud algorithm with a maximum of three words per label. Nodes are colored by enrichment score (ES), and edges are sized on the basis of the number of genes shared by the connected pathways.

Statistical Analysis
Unless otherwise specified, data are shown as mean ± SEM and each experiment was repeated three times. Data were analyzed by the two tailed unpaired t-test or one-way ANOVA with Dunnett's post hoc analysis using GraphPad Prism (v7.02, GraphPad).

Fucoidan Treatment Upregulates Gene Sets Associated with Immune Responses to Exogenous Stimuli
To investigate whether fucoidan can induce changes to the antigen-presenting cell transcriptome, mouse BMDCs were grown in GM-CSF-contained medium for 6 days and treated with fucoidan or vehicle (0.04% DMSO) for 24 h. RNA was isolated from the harvested cells (n = 3) and cDNA libraries were generated and sequenced for transcriptome analysis ( Figure 1). After raw RNA-seq data from BMDC samples were processed to define a ranked gene list, DNA sequence reads were quality filtered by trimming to remove lowquality bases and mapped to a genome-wide reference set of transcripts to enable counting of reads per transcript. Read counts were then aggregated at the gene level (counts per gene) and counts per million (cpm), and log FC values were obtained for bioinformatics analyses.
We began our analysis by generating a list of differentially expressed genes (DEGs) between fucoidan-and control-treated BMDCs using iDEP, a web-based tool designed for the analysis and interpretation of omics data [55]. With the DESeq2 package, we identified 950 upregulated and 537 downregulated genes using a false discovery rate (FDR) < 0.1 and fold-change > 2 ( Figure 2A). The MA plot demonstrates that fucoidan treatment dramatically upregulated several genes involved in immune function, including Irf7, Il12b, H2-K1, Hmox1, and Ifit3 ( Figure 2B). We began our analysis by generating a list of differentially expressed genes (DEGs) between fucoidan-and control-treated BMDCs using iDEP, a web-based tool designed for the analysis and interpretation of omics data [55]. With the DESeq2 package, we identified 950 upregulated and 537 downregulated genes using a false discovery rate (FDR) < 0.1 and fold-change > 2 ( Figure 2A). The MA plot demonstrates that fucoidan treatment dramatically upregulated several genes involved in immune function, including Irf7, Il12b, H2-K1, Hmox1, and Ifit3 ( Figure 2B). To investigate biological themes among our list of DEGs, we used the DEG2 function in iDEP to identify enriched KEGG pathways (Table 1). We found that downregulated genes were strongly involved in cell cycle and DNA replication, suggesting that fucoidan treatment inhibited cell proliferation. In contrast, upregulated genes were related to pathogen recognition and signaling pathways, including C-type lectin receptor (CLR), To investigate biological themes among our list of DEGs, we used the DEG2 function in iDEP to identify enriched KEGG pathways (Table 1). We found that downregulated genes were strongly involved in cell cycle and DNA replication, suggesting that fucoidan treatment inhibited cell proliferation. In contrast, upregulated genes were related to pathogen recognition and signaling pathways, including C-type lectin receptor (CLR), Toll-like receptor (TLR), and NOD-like receptor (NLR) pathways. As expected, downstream cytokine and IFN signaling pathways such cytokine-cytokine receptor interaction, TNF signaling, JAK-STAT signaling, and virus-specific responses were also overrepresented. Other enriched pathways included NF-κB signaling, cancer, and antigen processing and presentation. We then constructed a hierarchical clustering tree ( Figure S1) and network ( Figure 2C) in iDEP to summarize the overlap among enriched KEGG pathways [55,56]. Within upregulated genes, viral infection pathways were closely associated, including that of measles, influenza A, Epstein-Barr virus, and Kaposi sarcoma-associated herpes virus. These pathways were further associated with pathogen recognition and response pathways related to CLR, NLR, TLR, and JAK-STAT signaling. Interestingly, pathways related to adaptive immune responses such as antigen processing and presentation and type I diabetes mellitus pathways formed a separate cluster.
We then examined whether gene expression changes in fucoidan-treated BMDCs were driven by specific transcription factors (TFs). Target genes of specific TFs that were overrepresented in our list of DEGs (Table 2 and Table S1) were visualized as a hierarchical clustering tree and network ( Figure S2 (Table S1). On the other hand, upregulated DEGs were enriched with target genes of the IRF, STAT and NF-κB TF families (Table 2). These target genes tended to be known regulators of immune function, including cytokines, chemokines, PI3K/AKT/mTOR pathway, cytoplasmic RNA sensors, costimulatory molecules, MHC, and antigen presentation [57][58][59].

Gene Set Enrichment Analysis (GSEA) Showed That Fucoidan Strongly Induces Interferon Signaling Pathways in BMDCs
Given our above findings suggesting that fucoidan induces transcriptome changes to promote immune surveillance and signaling in BMDCs, we sought to expand our analysis and generate a ranked list of genes without prior filtering by p-value and fold change thresholds. Mapped RNA-seq data were normalized as counts per million (cpm) using the edgeR package in R. The resulting expression data were analyzed with the Broad Institute's gene set enrichment analysis (GSEA) software using the C2 CP (canonical pathways) sets from the MSigDB gene database [60]. GSEA of normalized counts of all expressed genes showed that 728 gene sets were significantly enriched (FDR < 25%) by fucoidan treatment, while 551 gene sets were significantly enriched in control BMDC samples. Even under fairly stringent criteria (p-value < 0.01), there were 402 gene sets deemed significantly enriched in fucoidan-treated BMDCs and 359 gene sets in control BMDCs.
We then applied EnrichmentMap clustering and annotation to collapse redundant gene sets and generate a network visualization (Figure 3). Like the network function in iDEP, the EnrichmentMap application groups similar gene sets (nodes) based on the number of overlapping genes [61]. We found that the largest cluster in fucoidan-treated BMDCs was 'induction of interferon signaling' (52 gene sets) followed by 'ECM collagen crosslinking' (20), 'activation of PI3K FGFR2' (20), 'NRF2 transcriptional activity' (4), 'tetrasaccharide linker sequence' (4), 'therapy induced HIF1 pathway' (3), 'chemokine and chemokine receptors' (2), and 'scavenger scavenging ligands' (2). Downregulated gene sets included that of DNA strand repair, citrate cycle, and oxidative phosphorylation.  Given our EnrichmentMap clusters, we looked for specific pathways associated with interferon signaling such as downstream antiviral effectors and upstream pathogen recognition. Similar to our findings with iDEP, fucoidan-treated BMDCs were enriched in virus-specific response pathways including that of SARS-CoV-2, respiratory syncytial virus, parainfluenza virus 3, influenza A, HCMV, KSHV, oncolytic virus, adenovirus, Given our EnrichmentMap clusters, we looked for specific pathways associated with interferon signaling such as downstream antiviral effectors and upstream pathogen recognition. Similar to our findings with iDEP, fucoidan-treated BMDCs were enriched in virus-specific response pathways including that of SARS-CoV-2, respiratory syncytial virus, parainfluenza virus 3, influenza A, HCMV, KSHV, oncolytic virus, adenovirus, reovirus, MERS-CoV, and Ebola virus (Table S2). Interestingly, RLR, TLR, NLR, and STING pathways, all of which involve nucleic acid-sensing PRRs, were also enriched (Table S3). Other upregulated PRR pathways include that of dectin-1, class A scavenger receptors (SRs), and complement receptor 3 (CR3), which are consistent with previous reports of fucoidandependent signaling [62]. Altogether, fucoidan treatment stimulates the expression of gene sets related to nucleic acid-sensing, IFN induction, and virus-specific immunity in BMDCs, suggestive of an antiviral phenotype.

GSEA Leading-Edge Analysis Identified a Cluster of Type I IFN Signaling-Centered Host Defense Mechanisms
Within an enriched gene set, there is often a core group of genes that drives the enrichment signal and thus is thought to have the greatest phenotypic impact. To determine these gene subsets in fucoidan-treated BMDCs, C2 canonical pathway gene sets derived from the KEGG, Reactome, and WikiPathways (WP) databases were used to run a GSEA (p < 0.05; FDR < 0.25) for a subsequent leading-edge analysis. The resulting set-to-set diagrams are shown in Figure 4, where the intensity of a cell color demonstrates the number of shared leading-edge genes between the two gene sets. Gene sets were manually grouped into B1-B9 based on the presence of overlapping leading genes (as indicated by cell color intensity) and signaling themes ( Figure 4A). The B1 group consists of gene sets related to TLR, type I IFN, IL1R, RLR, and NFκB signaling, while the B2 group includes type I IFN, IFNs, and SARS-CoV-2 infection signaling. Other notable leading genes include NRF2 and scavenger receptor signaling pathways.
In terms of color intensity and biologic function, groups B1 and B2 may also be grouped together as one large cluster A1. Within cluster A1, various PRR signaling pathways such as RIG-I-like receptors (RLRs), Toll-like receptors (TLRs), NOD-like receptors (NLRs), and STING pathways are overrepresented (Table 3). Type I IFN, IFNs, IL1, and NFκB signaling pathways are also included in this cluster, as well as several SARS coronavirus-relevant innate immune pathways. These findings further support that fucoidan treatment may activate host defense pathways against viral pathogens, in which type I IFN signaling plays a critical role.
Among the top ten enriched gene sets in cluster A1, representative enrichment plots are shown in Figure 4B with their leading-edge genes listed in Table S4. Notably, the enrichment score of the interferon a/b signaling Reactome gene set ( Figure 4B, left) appeared to be driven by genes related to RNA-sensing pathways, including endosomal TLRs, MDA5, RIG-I, and PKR. Similarly, gene sets specific to anti-SARS-CoV-2 activity ( Figure 4B, middle and right) appeared to be driven by type I IFNs, IFN receptors, JAK-STAT signaling, and several interferon-stimulated genes (ISGs). Additional leading-edge analysis of cluster A1 generated a list of shared leading-edge genes based on the number of pathways that share this gene (Table 4). Consistent with previous findings, type I IFN pathways (Ifnb1, Ifna4, Ifnar2, Ifnar1), JAK-STAT signaling (Stat1, Stat2, Jak1), and ISG genes were highly shared among enriched pathways, as well as those of NF-κB, MAPK, and PI3K-AKT pathways. Genes of DC costimulatory molecules (Cd86, Cd80, Cd40), inflammatory cytokines, and chemokines were also observed. Set-to-set and heatmap diagrams are shown in Figure  S3. These results further support that viral recognition and type I interferon induction pathways are strongly induced in fucoidan-treated BMDCs.   Type I  Type II   119  316  37 Type III ( 0 ) TF targets   GSEA was then performed using the regulatory target gene (TFT) set in C3 to identify TFs specifically linked to fucoidan-induced gene expression changes. Unsupervised ranking for confidence of enrichment identified IRF, NF-κB and STAT targets as among the most significant TF-related changes between samples ( Figure 4C). As expected, the E2F TF family was enriched in control BMDCs (Table S5).

A1
Based on our results suggesting the importance of IFN-stimulation, we analyzed DEGs previously generated in iDEP via Interferome v2.01 [61], a database of IFN-regulated genes. Interestingly, 435 of the 950 upregulated genes belonged to type I IFN-stimulated ISGs ( Figure 4D), further supporting that fucoidan treatment leads to upregulation of ISGs in BMDCs.

Fucoidan Stimulation Promotes BMDC Maturation and Activation In Vitro
Type I IFN signaling is known to promote DC maturation and thereby stimulate activation of the adaptive immune response [63]. Consistent with our findings that fucoidan induces strong type I IFN induction in BMDCs, GSEA of the same samples showed enrichment of Lindstedt dendritic cell maturation A (NES 2.16, FDR q-value < 0.0001), B (NES 2.20, FDR q-value < 0.0001), and C (NES 1.80, FDR q-value = 0.026) gene sets ( Figure 5A). Furthermore, the leading-edge genes were those critical for adaptive immune activation, such as co-stimulatory molecules (Table S6). Analysis of normalized RNA-seq data (cpm) had also shown up-regulation of co-stimulatory molecules (e.g., CD80, CD86, and CD40) and CD83 in fucoidan-treated BMDCs ( Figure 5B). To complement our transcriptome analyses, we assessed the capacity of fucoidan to induce maturation of BMDCs by measuring the surface expression of DC maturation and T-cell co-stimulation markers by flow cytometry. After 24 h of treatment, MHCII, CD40, and CD80 were significantly upregulated ( Figure 5C). Statistical difference by using the unpaired t-test relative to vehicle (B) or, one-way ANOVA with Dunnett's post-test analysis compared to stained control(0 μg/mL) (C). The bar above the data displays ±SEM. Asterisks depict levels of significance as follows: NS, not significant (p ≥ 0.05); * p < 0.05; ** p < 0.01; *** p < 0.001.

Discussion
In the midst of a viral pandemic, fucoidan has been rapidly garnering interest for its antiviral activity [24,26]. Here, we analyzed global transcriptomic changes in fucoidantreated BMDCs to better understand fucoidan's antiviral properties. We found that fucoidan treatment upregulated viral-sensing receptors, their downstream signaling molecules, and subsequent effectors, effectively establishing an antiviral state in BMDCs.
Dendritic cells play a central role in the innate immune system to clear viral infections [29,34,64]. DCs recognize viral PAMPs and activate multiple signaling cascades, ultimately leading to the transcriptional induction of type I IFNs [6,29]. Ligand binding of type I IFN receptors will in turn activate the JAK-STAT pathway, driving the expression of hundreds of IFN-stimulated genes (ISGs) that directly and indirectly exert antiviral activities [6,65]. Both the DEGs and enriched gene sets in fucoidan-treated BMDCs showed strong enrichment of nucleic acid-sensing PRR signaling pathways such as that of RLR, TLR, NLR, and STING. Other enriched gene sets included that of inflammatory cytokines and chemokines, type I IFN production, and their respective receptors. Statistical difference by using the unpaired t-test relative to vehicle (B) or, one-way ANOVA with Dunnett's post-test analysis compared to stained control (0 µg/mL) (C). The bar above the data displays ±SEM. Asterisks depict levels of significance as follows: NS, not significant (p ≥ 0.05); * p < 0.05; ** p < 0.01; *** p < 0.001.

Discussion
In the midst of a viral pandemic, fucoidan has been rapidly garnering interest for its antiviral activity [24,26]. Here, we analyzed global transcriptomic changes in fucoidantreated BMDCs to better understand fucoidan's antiviral properties. We found that fucoidan treatment upregulated viral-sensing receptors, their downstream signaling molecules, and subsequent effectors, effectively establishing an antiviral state in BMDCs.
Dendritic cells play a central role in the innate immune system to clear viral infections [29,34,64]. DCs recognize viral PAMPs and activate multiple signaling cascades, ultimately leading to the transcriptional induction of type I IFNs [6,29]. Ligand binding of type I IFN receptors will in turn activate the JAK-STAT pathway, driving the expression of hundreds of IFN-stimulated genes (ISGs) that directly and indirectly exert antiviral activities [6,65]. Both the DEGs and enriched gene sets in fucoidan-treated BMDCs showed strong enrichment of nucleic acid-sensing PRR signaling pathways such as that of RLR, TLR, NLR, and STING. Other enriched gene sets included that of inflammatory cytokines and chemokines, type I IFN production, and their respective receptors. Furthermore, fucoidan stimulation led to a strong induction of hundreds of ISGs as demonstrated by our Interferome analysis. These data are consistent with the general understanding that multiple PRR signaling pathways converge to activate inflammatory transcription factors such as NF-κB, IFN regulatory factor (IRF) 3, and IRF7 to subsequently induce cytokine, chemokine, and interferon (IFN) expression. As such, fucoidan may induce a self-enhancing, antiviral program in BMDCs, highlighting its clinical potential for suppressing virus activity.
In antiviral defense, DCs also act to activate and amplify the adaptive immune response [40][41][42]. DCs are professional antigen presenting cells (APCs) that capture viral pathogens for degradation and antigen presentation [41,49]. Naïve T cells recognize viral antigens only after processing and presentation by DCs in association with MHC molecules, leading to differentiation of antigen-specific Th1 CD4+ and CTL CD8+ T cells to eradicate the virus [48]. Our analysis of DEGs in fucoidan-treated BMDCs demonstrated strong induction of IL12, which is known to be a potent inducer of IFN-gamma production in NK cells to produce a Th1-skewed response [66,67]. Furthermore, GSEA and counts from our RNA-seq data showed that fucoidan treatment results in strong enrichment of dendritic cell maturation genes and co-stimulatory molecules, such CD80, CD86, CD40, and CD83, ( Figure 5B). DEGs in fucoidan-treated BMDCs also included those related to DC activation, such as MHC I gene H2-K1, MHCI antigen presentation-related genes Tap1 and B2m, and MHC II genes H2-M3 and H2-Q4. These findings were further confirmed by our flow cytometry analysis showing that fucoidan upregulates surface expression of MHCII, CD80, and CD40 in BMDC cultures. Altogether, fucoidan may also enhance antigen presentation and co-stimulation functionality in BMDCs.
Given its ability to enhance both the innate and adaptive roles of DCs, fucoidan may effectively prevent and control a wide range of viral infectious diseases. While clinical observations during the recent outbreaks of coronaviruses have highlighted the importance of type I IFNs as a potent antiviral molecule, a state of hyper-inflammation has been linked to disease severity and poor outcomes following infection [68]. To appropriately manage the symptoms of COVID-19, it is crucial to balance adaptive immune activation for viral clearance, without leading to uncontrolled inflammation [68]. Interestingly, our results also suggest that fucoidan induces NRF2-based anti-inflammatory signaling pathways and heme oxygenase 1 (HO-1). It was recently demonstrated that several NRF2 activators suppress the inflammatory response to SARS-CoV-2 in human cells [69]. Upregulation of the NRF2-transcriptional target HO-1 (human gene name HMOX1) has also been linked to antiviral activity by catalyzing the degradation of heme into biliverdin, Fe 2+ , and CO, all of which have putative anti-SARS-CoV-2 activity [70]. As a potent inducer of both type I IFN signaling and NRF2 anti-inflammatory signaling, fucoidan-based nutraceuticals or functional foods may act as an effective, balanced strategy against COVID-19 [71].
The specific receptors responsible for fucoidan's antiviral activity are not well known. As hydrophilic macromolecules, polysaccharides generally bind scavenger receptors (SRs), TLRs, complement receptor 3 (CR3), and C-type lectin receptors (CLRs) [62]. Sulfated polysaccharides in particular may bind multiple receptors due to their polyvalency. In our study, both CLR and TLR signaling pathways were among the KEGG pathways overrepresented in DEGs in fucoidan-treated BMDCs (Table 1). These were confirmed in our the subsequent GSEA, which additionally showed the enrichment of RLRs, CR3, class A SRs, and cytosolic DNA sensing molecules (Supp. Table 3). These findings are notable given that CLRs are known to route viral antigens such as Dectin1 to MHC class II compartments for presentation to CD4 + T cells [71]. In addition, SR-A has been shown to be a primary endocytic receptor that will internalize fucoidan and activate TLR3 and TLR9, thereby promoting type I IFN signaling [62,72]. It has been further reported that fucoidan may directly activate TLR-2 and TLR-4 [73]. Our findings are consistent with previous literature suggesting that TLRs, SRs, and CLRs may function as fucoidan receptors and subsequently lead to type I IFN signaling and ISG production.
In summary, we have demonstrated that the sulfated polysaccharide fucoidan induces global transcriptomic changes in BMDCs that may promote antiviral activity. Fucoidantreated cells showed overrepresentation of virus-specific response pathways, including that of SARS-CoV-2, as well as gene sets related to nucleic acid-sensing PRR signaling, type I interferon (IFN) signaling, and ISG production. Our findings may serve as a key foundation for further research on the development of fucoidan-based strategies for the prevention and treatment of COVID-19 and other human viral diseases. In the future, critical observations obtained in the bioinformatic analyses need to be verified using other methods and also fucoidan will need to be investigated in animal models as an antiviral nutraceutical for boosting the type 1 interferon response to RNA viruses, such as influenza and coronaviruses.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/nu14112242/s1, Figure S1: Visualization of the relationship among enriched KEGG pathways in DEGs using tree; Table S1: Enriched TF target RegNetwork; Figure S2: Visualization of the relationship among enriched TF target gene sets using (A) tree and (B) network; Table S2: Virusspecific immune pathways enriched in fucoidan-treated BMDCs; Table S3: PRR signaling pathways enriched in fucoidan-treated BMDCs; Table S4: Leading-edge genes within representative enrichment plots shown in Figure 4B; Figure S3: Leading-edge analysis of cluster A1 enriched genes; Table S5: GSEA of TF targets in fucoidan vs control BMDCs; Table S6: Leading-edge genes within LINDSTEDT DC maturation gene set enrichment plots in Figure 5A.

Informed Consent Statement: Not applicable.
Data Availability Statement: All raw and processed sequencing data generated in this study have been submitted to the NCBI Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) under accession number GSE200944 and the Korean Nucleotide Archive (KoNA; https://kobic.re.kr/ kona) under KoNA access ID PRJKA220176.