Next Article in Journal
Peptide LCGA-17 Attenuates Behavioral and Neurochemical Deficits in Rodent Models of PTSD and Depression
Next Article in Special Issue
miR-182-5p Regulates Nogo-A Expression and Promotes Neurite Outgrowth of Hippocampal Neurons In Vitro
Previous Article in Journal
Serotonin 1A Receptor Pharmacotherapy and Neuroplasticity in Spinal Cord Injury
Previous Article in Special Issue
Association between Proton Pump Inhibitor Use and Parkinson’s Disease in a Korean Population
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis of Fumarate Compounds Identifies Unique Effects of Isosorbide Di-(Methyl Fumarate) on NRF2, NF-kappaB and IRF1 Pathway Genes

by
William R. Swindell
1,*,
Krzysztof Bojanowski
2,3 and
Ratan K. Chaudhuri
3,4
1
Department of Internal Medicine, The Jewish Hospital, Cincinnati, OH 45236, USA
2
Sunny BioDiscovery Inc., Santa Paula, CA 93060, USA
3
Symbionyx Pharmaceuticals Inc., Boonton, NJ 07005, USA
4
Sytheon Ltd., Boonton, NJ 07005, USA
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2022, 15(4), 461; https://doi.org/10.3390/ph15040461
Submission received: 27 February 2022 / Revised: 9 April 2022 / Accepted: 10 April 2022 / Published: 11 April 2022
(This article belongs to the Special Issue Therapeutic Agents for Neurological Disorders 2022)

Abstract

:
Dimethyl fumarate (DMF) has emerged as a first-line therapy for relapsing-remitting multiple sclerosis (RRMS). This treatment, however, has been limited by adverse effects, which has prompted development of novel derivatives with improved tolerability. We compared the effects of fumarates on gene expression in astrocytes. Our analysis included diroximel fumarate (DRF) and its metabolite monomethyl fumarate (MMF), along with a novel compound isosorbide di-(methyl fumarate) (IDMF). Treatment with IDMF resulted in the largest number of differentially expressed genes. The effects of DRF and MMF were consistent with NRF2 activation and NF-κB inhibition, respectively. IDMF responses, however, were concordant with both NRF2 activation and NF-κB inhibition, and we confirmed IDMF-mediated NF-κB inhibition using a reporter assay. IDMF also down-regulated IRF1 expression and IDMF-decreased gene promoters were enriched with IRF1 recognition sequences. Genes altered by each fumarate overlapped significantly with those near loci from MS genetic association studies, but IDMF had the strongest overall effect on MS-associated genes. These results show that next-generation fumarates, such as DRF and IDMF, have effects differing from those of the MMF metabolite. Our findings support a model in which IDMF attenuates oxidative stress via NRF2 activation, with suppression of NF-κB and IRF1 contributing to mitigation of inflammation and pyroptosis.

1. Introduction

Multiple sclerosis (MS) is a chronic immune-mediated demyelinating disease affecting the central nervous system [1]. The estimated worldwide number of MS patients is 2.8 million and prevalence has increased in every region for which estimates are available [2]. MS is a significant source of mortality accounting for approximately 19,000 deaths worldwide in 2016 [3] with disease burden greatest at higher latitudes [4]. The pathophysiology is initiated by breakdown of basement membrane, endothelial cells, and tight junctions that comprise the blood–brain barrier (BBB) [5]. This leads to entry of T-cells into the CNS [6], which then undergo activation to release cytokines, such as TNF, IL-1, IL-6, IL-17, and IL-22 [7]. These cytokines have direct toxic effects on myelin surrounding axons [8], but additionally can enhance BBB degradation to promote an influx of B cells and macrophages [9]. Oligodendrocytes counter this neuroinflammatory cascade, by actively promoting remyelination, but over time disease progression occurs due to loss of oligodendrocytes and worsening inflammation [10]. This progression is seen on MRI by development of white matter plaques [11], along with symptoms, such as dysarthria, nystagmus, tremor, and autonomic dysfunction [12]. In the most common disease subtype, relapse-remitting multiple sclerosis (RRMS), alternating periods of demyelination, and remyelination occur with exacerbation and remission of symptoms [13]. There is no definitive MS cure, but steroids can treat acute exacerbations, and disease-modifying therapies, such as beta-interferon, natalizumab, glatiramer acetate, and ocrelizumab, may prolong remission [14].
Fumaric acid esters have an immunomodulatory mechanism of action with demonstrated efficacy for treatment of autoimmune diseases, such as psoriasis and multiple sclerosis [15]. The effectiveness of dimethyl fumarate (DMF) in the setting of multiple sclerosis (RRMS) was initially demonstrated by two phase III trials (DEFINE and CONFIRM) [16,17]. These trials demonstrated reduced rates of disease relapse among adults (ages 18–55) with recent disease activity. In both trials, flushing was the most reported adverse effect (24–38%), although gastrointestinal (GI) adverse effects occurred in 7 to 19% of patients. It has subsequently been reported that 36% of patients in these trials experienced adverse GI effects after 1 year [18], and in real-world settings, GI side effects led to DMF discontinuation in 5–19% of patients [19,20]. To address this issue, a novel agent diroximel fumarate (DRF) has been developed, which undergoes esterase cleavage in the GI tract to generate monomethyl fumarate (MMF), 2-hydroxyethyl succinimide (HES), RDC-8439, and methanol [21]. MMF is believed to be the active metabolite of DMF and DRF 462 mg and DMF 240 mg given twice-daily (BID) appear bioequivalent in terms of metabolic MMF generation. A 96-week study (EVOLVE-MS-1) showed that 30.9% of patients given DRF 462 mg BID still developed adverse GI effects, but such effects led to discontinuation in only 6.3% of patients [18]. Consistent with this, a 5-week study showed that the number of days with at least moderate GI symptoms is reduced by 46% in patients given DRF 462 mg BID as compared to DMF 240 mg BID [22].
The improved tolerability of DRF appears to be related to its unique chemical structure, which is mostly (>90%) metabolized to HES and MMF, with RDC-8439 and methanol being generated as minor products (<10%) [21]. The generation of high methanol concentrations in the small intestine is likely responsible for adverse GI effects seen with DMF and DRF treatment [23]. This has raised the possibility that other fumarate esters with alternative chemical structures may also, like DRF, demonstrate improved tolerability based upon differences in their metabolic byproducts. We recently developed a novel fumarate compound, isosorbide di-(methyl fumarate) (IDMF), consisting of two methyl fumarate groups attached to a central isosorbide moiety [24,25]. IDMF is predicted to be less reactive than DMF due to its chemical structure and higher activation free energy barrier [24]. When topically applied to skin, IDMF is less harsh compared to DMF, without evidence of irritation or sensitization in multiple animal models (i.e., rat, rabbit, guinea pig) [24]. A prior transcriptome study showed that treatment of astrocytes with IDMF decreased expression of genes associated with a proliferative reactive astrocyte phenotype as well as NF-κB target genes [25]. No prior study, however, has compared transcriptional effects of IDMF to those of MMF or DRF and, thus, it remains unclear whether these fumarate compounds differentially impact key disease pathways associated with MS progression.
This study used microarray profiling to evaluate the effects of three fumarate compounds on gene expression in cultured astrocytes (MMF, DRF, and IDMF). Astrocytes are critical for maintaining BBB integrity and interact with oligodendrocytes to influence remyelination and neuroinflammation [26]. We identify individual genes and co-expression modules regulated by each fumarate and evaluate their effects on genes linked to the nuclear factor-erythroid factor 2-related factor 2 (NRF2), nuclear factor kappa B (NF-κB), and hydroxycarboxylic acid receptor 2 (HCAR2/GPR109A) pathways [27]. Expression responses relevant to MS are established by comparisons to literature-based gene–disease association databases [28,29,30,31,32], loci previously associated with MS in genome-wide association studies (GWAS) [33], and expression profiles of astrocytes isolated from MS patient white matter using laser capture microdissection [34]. Our results uncover unique effects of IDMF not replicated by other fumarates and provide insights into cellular mechanisms by which such compounds may influence MS disease activity.

2. Results

2.1. The MMF Expression Profile Differs from That of IDMF and DRF

We treated human fetal astrocytes with sterile distilled water (CTL), MMF, DRF, or IDMF (2.5 µM each) for 24 h and evaluated gene expression using in situ oligonucleotide microarrays (see Methods). The three fumarate compounds (Figure 1A–C) had disparate effects on gene expression. This could be discerned from evaluation of differential expression vectors in a two-dimensional principal component space, which showed that DRF and IDMF shifted expression in a direction opposite from that of MMF (Figure 1D). Consistent with this, we identified groups of genes for which the estimated FC (MMF/CTL) differed from those calculated for DRF and IDMF (Figure 1E). Among all astrocyte-expressed genes, DRF/CTL and IDMF/CTL FC estimates were highly correlated with each other (rs = 0.65), but weakly correlated with MMF/CTL FC estimates (rs ≤ 0.18; Figure 1F). These trends were further supported by inspection of self-organizing maps color-coded based upon FC estimates, which revealed similar patterns for DRF and IDMF differing from that of MMF.

2.2. IDMF Leads to Extensive Down-Regulation of Genes Associated with Metabolism and Post-Transcriptional Regulation

Differentially expressed genes (DEGs) were identified using generalized least square linear models with empirical Bayes moderated t-statistics (see Methods) [35]. We identified 450, 379, and 1576 DEGs with respect to the MMF vs. CTL, DRF vs. CTL and IDMF vs. CTL comparisons, respectively (p < 0.05 with FC > 1.25 or FC < 0.80; Figure S1K–P; Tables S1–S6, Supplementary Materials). Most of the 30 genes most strongly increased or decreased by MMF were not correspondingly regulated by DRF and IDMF (Figure 2A). The strongest overall effect of MMF (i.e., lowest p-value) was up-regulation of nuclear factor of activated T cells 2 (NFATC2) (Figure 2D), although DEGs uniquely up-regulated by MMF included C1orf54, BSCL2, and LSAMP, and DEGs uniquely down-regulated by MMF included ZC3H4, GPR157 and SLC20A2 (Figure 2A). The complete set of MMF-increased DEGs was most strongly associated with oxidation–reduction, FC-epsilon receptor and amino acid metabolism (Figure 2G), whereas MMF-decreased DEGs were associated with mast cell degranulation/activation, positive regulation of TNF production, and phospholipid metabolism (Figure 2H).
Most DRF-increased DEGs were not correspondingly altered by MMF or IDMF, although in contrast DRF-decreased genes were often correspondingly decreased by IDMF (Figure 2B). Examples of genes uniquely up-regulated by DRF include ABI3BP, GSAP, and ZNF180, while genes uniquely down-regulated by DRF include METTL8, COA7, and ZNF567 (Figure 2B). Overall, down-regulation of SAC1 like phosphatidylinositide phosphatase (SACM1L) was the strongest transcriptional effect of DRF (Figure 2E). The set of DRF-increased DEGs was associated with circadian behavior, purine nucleobase metabolism, and cardiocyte differentiation (Figure 2I), while DRF-decreased DEGs were associated with macromolecule modification, p53 regulation, and regulation of glycosylation (Figure 2J).
The majority of genes strongly regulated by IDMF were not significantly altered by MMF or DRF (Figure 2C). Unique effects of IDMF include up-regulation of OIT3, C9orf24, and ZDHHC11B, and down-regulation of TNFAIP3, KPNB1, and COPA (Figure 2C). The strongest overall effect of IDMF was down-regulation of actin related protein 3 (ACTR3) (Figure 2F). The complete set of IDMF-increased DEGs was most strongly associated with olfactory detection, phospholipid remodeling, and antifungal innate immunity (Figure 2K). The set of IDMF-decreased DEGs included approximately 300 genes associated with cellular metabolism and protein modification, along with 50 genes associated with post-transcriptional gene regulation (Figure 2L).
Enrichment of motifs in 5 kb sequences upstream of DEGs was assessed using semiparametric generalized additive logistic models (Figure S2, Supplementary Materials) [36]. IDMF-increased DEG promoters were enriched with 5-CATG/CATG-3 and 5-ATGGG/CCCAT-3 elements (Figure S2D). MMF-increased DEGs and DRF/IDMF-decreased DEGs were enriched for AT-rich motifs (Figure S2B,C,G). Both DRF- and IDMF-decreased DEG promoters were significantly enriched for a 5-TAATT/AATTA-3 element recognized by paired related homeobox 1 and 2 (PRRX1 and PRRX2) (Figure S2G).

2.3. IDMF Up-Regulates Co-Expressed Tight Junction Genes and Represses Groups of Co-Expressed Genes Localized to Mitochondria

We next identified 234 gene expression modules (25–196 genes/module) by clustering genes based upon their expression across an independent set of astrocyte microarray samples (see Methods). Genes within each module had a stronger tendency to be altered in the same direction by IDMF as compared to DRF or MMF (Figure 3A,B). Consequently, we identified 141 modules differentially expressed by IDMF, whereas only 82 and 25 were identified by DRF and MMF, respectively (Figure 3C). Modules most strongly altered by IDMF were similarly altered but to a lesser degree by DRF, although usually MMF had opposite effects (Figure 3B,D,E). The module DCTN1-74 was most strongly up-regulated by IDMF (p = 9.77 × 10−19; FDR = 2.29 × 10−16; Figure 3D,E). Genes belonging to this module included KLHL33, CARD14, and SCAF1 (Figure 3F), and such genes were often localized to tight junctions (Figure 3H) and associated with inhibition of translation and positive regulation of proliferation (Figure 3G). Likewise, module SLC8A1-109 was most strongly down-regulated by IDMF (p = 1.33 × 10−14; FDR = 1.55 × 10−12; Figure 3D,E). Genes belonging to this module included MRPS28, TIMM9, and GRSF1 (Figure 3I), and such genes were frequently localized to the mitochondrial membrane (Figure 3K) and associated with mitochondrial translation and gene expression (Figure 3J).

2.4. The Effects of DRF and IDMF on Gene Expression Are Consistent with NRF2 Pathway Activation

MMF significantly increased expression of the NRF2 target gene NAD(P)H quinone dehydrogenase 1 (NQO1) (Figure 4A) but did not significantly increase expression of other NRF2 targets, such as HMOX1, G6PD, and TXNRD1 (Figure 4B–D). Each of these NRF2 targets, however, was significantly increased by DRF and/or IDMF (Figure 4A–D). To evaluate NRF2 regulation more systematically, we evaluated 143 genes recently identified as belonging to a NRF2 expression signature, including 68 genes up-regulated with NRF2 activation (Figure 4E) and 75 genes down-regulated with NRF2 activation (Figure 4F) [37]. There was no significant trend towards increased or decreased expression of signature genes with MMF treatment (Figure 4G–J). However, among 56 astrocyte-expressed genes up-regulated by NRF2 activation, 73% were correspondingly increased by DRF (p < 0.001; Figure 4K) and such genes were enriched among DRF-increased genes (p < 0.001; Figure 4M). On the other hand, there was no significant trend to indicate DRF-decreased expression of genes down-regulated with NRF2 activation (Figure 4L,N). Expression patterns associated with IDMF treatment were slightly more consistent with NRF2 activation. Most genes up-regulated by NRF2 activation were IDMF-increased (61%; p = 0.18; Figure 4O) and such genes were enriched among IDMF-increased genes (p = 0.014; Figure 4Q). Likewise, most genes down-regulated by NRF2 activation were IDMF-decreased (56%, p = 0.018; Figure 4P) and such genes were enriched among IDMF-decreased genes (p = 0.0186; Figure 4R).

2.5. The Effects of MMF and IDMF on Gene Expression Are Consistent with NF-κB Suppression

An initial analysis of four well-validated NF-κB targets (NFKBIA, TNFAIP3, BCL2L1, ICAM1) showed that each was more strongly down-regulated by IDMF compared to MMF and DRF (Figure 5A–D). We next focused on the expression of 389 NF-κB (RELA) targets recently identified from the genomic analysis of a B-cell lymphoma cell line (BJAB) [38]. This gene set included 304 RELA-activated direct targets induced by NF-κB-stimulating treatments, associated with RELA binding (ChIP-Seq), and down-regulated by the dominant negative NF-κB inhibitor alpha (dnIκBα) (Figure 5E) [38]. The gene set included additional 85 RELA-suppressed targets up-regulated by dnIκBα expression (Figure 5F) [38].
There was no significant evidence to support inhibition of NF-κB signaling by DRF (p ≥ 0.184; Figure 5K–N). However, most RELA-activated genes were down-regulated by MMF and IDMF (p ≤ 0.031; Figure 5G–O), and consistent with this such genes were enriched among genes decreased by MMF and IDMF (p ≤ 0.024; Figure 5I–Q). There was no significant evidence to demonstrate up-regulation of RELA-suppressed genes by MMF or IDMF (Figure 5H,P,J,R), although the majority of such genes were up-regulated by IDMF (p = 0.17; Figure 5P).
We did not observe significant changes in the expression of genes encoding NF-κB subunit proteins (i.e., NFKB1, NFKB2, REL, RELA, RELB; Figure S3D–I, Supplementary Materials). Of six different NF-κB binding site motifs examined (Figure S3A), there was one motif significantly enriched in sequences upstream of IDMF-increased genes (p < 0.05; Figure S3B), but such motifs were not enriched in sequences upstream of genes decreased by any fumarate compound (Figure S3C).

2.6. MMF Significantly Up-Regulates the Expression of Genes Induced via the Niacin-HCAR2 Pathway

The response to DMF is also mediated by the niacin-activated HCAR2 (GPR109A) receptor [27]. We did not observe significant changes in the expression of HCAR2 (Figure S4A, Supplementary Materials) or suspected niacin-HCAR2 pathway targets (Figure S4B–D) [27,39,40]. We thus focused on the expression of genes with altered expression in monocytes treated with niacin (100 μM) for 6 h (GSE103381) [41]. MMF was the only fumarate for which niacin-increased genes were biased towards up-regulated expression, consistent with activation of the niacin-HCAR2 pathway (Figure S4G–R). Nonetheless, we identified some individual niacin-responsive genes significantly altered by each fumarate compound (Figure S4E–F).

2.7. IDMF Down-Regulates IRF1 Expression and IRF1 Binding Sites Are Enriched in Sequences Upstream of IDMF-Decreased Genes

Interferon regulatory factor 1 (IRF1) expression was significantly down-regulated by IDMF (Figure S5D, Supplementary Materials) although genes encoding other IRF factors were not significantly altered by any fumarate (Figure S5E,F,H–J). We evaluated enrichment of 19 IRF family motifs (Figure S5A) within 5 kb upstream regions of genes significantly altered by each fumarate (Figure S5B,C). An IRF1 motif was significantly enriched (FDR < 0.05) in regions upstream of MMF-increased DEGs, and there was marginal enrichment (p < 0.05) of IRF4 motifs in regions upstream of genes increased by MMF or IDMF (Figure S5B). Promoters of IDMF-decreased genes were strongly enriched (FDR < 0.05) with motifs linked to IRF1, IRF4, IRF6, and IRF7, and likewise promoters of DRF-increased genes were enriched for motifs associated with IRF1 and IRF7 (Figure S5C). Predicted IRF1 targets among IDMF-decreased genes included MID1, PSMB5, and GTF2H3 (Figure S5G). Putative IRF1 binding sites identified near such genes were concentrated approximately 500 kb upstream from each gene’s transcription start site (Figure S5K).

2.8. IDMF Decreases the Expression of A1 Reactive Astrocyte Marker Genes Associated with Cell Death, Catabolism, and Lysosomes

Astrocytes develop a reactive phenotype in response to CNS injury or disease [42]. The A1 and A2 sub-phenotypes are two reactive astrocyte subtypes described previously, with the A1 phenotype viewed as neurotoxic and the A2 phenotype viewed as neuroprotective [43]. The expression of two A1 marker genes, complement C3 (C3) and guanylate binding protein 2 (GBP2), were not significantly altered by MMF, DRF, or IDMF (Figure S6A,B, Supplementary Materials). Likewise, expression of A2 markers S100 calcium binding protein A10 (S100A10) and transforming growth factor beta 1 (TGFB1) were not significantly altered by any fumarate (Figure S6C,D).
We next evaluated the expression of gene sets associated with the A1 phenotype (n = 161 genes), A2 phenotype (n = 75), and pan-reactive astrogliosis (A1 + A2) [25,44]. MMF led to a significant (p = 0.03) but slight (3%) increase in the average expression of A2 genes (Figure S6E) but only altered expression of some A1 and no A2 marker genes (Figure S6F,G). DRF altered the expression of some A1 and A2 marker genes (Figure S6I,J) but there was no significant systemic pattern among all marker genes (Figure S6H). IDMF decreased expression of A2 marker genes by only 1% on average although the median FC differed significantly from other astrocyte-expressed genes (p = 0.04; Figure S6K). Notably, we identified 19 A1 and 6 A2 marker genes significantly increased or decreased by IDMF (p < 0.05; Figure S6L,M). Of the A1 marker genes, 11 were IDMF-decreased and these were associated with catalytic activity, positive regulation of programmed cell death, and lysosomes (Figure S6N–P).

2.9. IDMF Alters the Expression of More MS-Associated Genes Compared to MMF and DRF

MS-associated genes were identified from seven database sources (Figure 6A). Only two genes were linked to MS by all seven sources (HLA-DRB1 and HLA-DQB1) while six genes were linked to MS by six of the seven sources (Figure 6B). Only 6% of the MS-associated genes were linked to MS by three or more database sources (i.e., 331 of 5160 genes; Figure 6C). Several genes linked to MS by four or more sources were significantly altered by IDMF but not by MMF or DRF (e.g., TNFAIP3, CXCL8, VCAM1; Figure 6D). On average, MMF and DRF altered expression of MS-associated genes (4+ sources) by 9% and 10%, respectively, which was not significantly different compared to other genes with detectable expression (p ≥ 0.759; Figure 6E,G). However, IDMF altered the expression of MS-associated genes by 20% on average, which was significantly greater than that seen for other expressed genes (p = 0.002; Figure 6I). A larger number of MS-associated genes were also significantly altered by IDMF compared to MMF and DRF (p < 0.05; Figure 6F,H,I). For example, IDMF significantly increased expression of HMOX1, HLA-DRB5, and BTNL2, and decreased expression of TNFAIP3, IFNAR1, and ICAM1 (p < 0.05; Figure 6J). As a group, IDMF-increased genes (p < 0.05, FC > 1.25) associated with MS (3+ sources) were associated with response to LPS, myeloid leukocyte activation, and negative regulation of inflammatory response (Figure 6K). Likewise, IDMF-decreased genes (p < 0.05, FC < 0.80) associated with MS (3+ sources) were associated with bacterial response, positive regulation of immune process, and IL-1 production (Figure 6L).

2.10. Genes Most Strongly Altered by MMF, DRF, and IDMF Overlap Significantly with Genes near MS GWAS Loci

We next evaluated overlap between DEGs and those genes near MS-associated loci from the NHGRI-EBI GWAS catalog. Genes significantly altered by MMF, DRF, and IDMF (increased and decreased) did not overlap significantly with MS-associated genes from GWA studies, although such genes were approximately four times more frequent among genes altered by IDMF compared to MMF or DRF (Figure 7A,E,I). However, among the 30 genes most strongly altered by MMF, DRF, or IDMF (i.e., lowest p-value), there was significant overlap between such DEGs and genes located varying distances from MS GWAS loci (Figure 7B,F,J). For IDMF, the top 30 increased and decreased genes overlapped significantly with genes near GWAS loci (Figure 7J), whereas for other compounds overlap was only significant for decreased (MMF) or increased (DRF) genes (Figure 7B,F). As an alternative approach, we examined the average distance between the top 30 DEGs and MS GWAS loci (Figure 7C,D,G,H,K,L). This average distance was marginally less than expected (p < 0.10) in the case of MMF-decreased, DRF-increased, and IDMF-increased genes (Figure 7D,G,K). We further identified several IDMF-regulated genes that were overlapping with or near MS-associated GWAS loci (Figure 7M,N), including ALDH1L1, MLANA, SLC2A4RG, ZBTB38, and ZNF433 (Figure 7O–R).

2.11. IDMF Stimulates the Expression of Transcriptional Activation Pathways Correspondingly Up-Regulated in MS Patient Astrocytes

We next evaluated overlap between DEGs and genes with altered expression in astrocytes obtained from normal appearing white matter within brain tissue sections of MS patients (GSE83670) [34]. The effects of MMF on gene expression were negatively correlated with those seen in MS astrocytes (r = −0.03), with a significant overabundance (33.5%) of genes increased by MMF and decreased in MS astrocytes (Figure S7A, Supplementary Materials). Consistent with this, there was no significant overlap between genes altered by MMF and genes correspondingly altered in MS astrocytes (Figure S7B–D). Gene expression shifts in DRF- and IDMF-treated astrocytes were positively correlated with those seen in MS astrocytes (Figure S7E,I). There was significant overlap between DRF-decreased DEGs and MS-decreased genes (Figure S7G), but not with respect to DRF-increased genes or genes altered by IDMF (Figure S7F,J,K). GSEA analysis, however, showed that the top 100 genes increased by IDMF overlapped significantly with those genes increased in MS astrocytes, while the top 100 genes decreased by IDMF overlapped significantly with genes decreased in MS astrocytes (Figure S7L). This significant pattern was otherwise only seen with respect to DRF-decreased genes (not DRF-increased genes; Figure S7H). Genes increased by both IDMF and in MS astrocytes included MYO15B, TP53INP2, and RASSF2 (Figure S7M), and such genes were associated with positive regulation of transcription, developmental process, and cell differentiation (Figure S7O). Genes decreased by both IDMF and in MS astrocytes included SPPL2A, ABCE1, and PRKCI, and such genes were associated with macromolecule methylation, ephrin receptor signaling, and response to wounding (Figure S7P).

2.12. IDMF Represses NF-κB Induction in TNF-Stimulated Astrocytes

RT-PCR analyses were done to confirm microarray findings for selected IDMF-regulated genes (Figure 8B–H). IDMF increased expression of OSGIN1 by 60% based upon PCR, consistent with array findings (p < 0.05, Fisher’s LSD; Figure 8B). IDMF repressed expression of ICAM1 and MALT1, although this effect did not reach statistical significance (p > 0.05, Figure 8D,E). However, IDMF significantly decreased expression of IL6, TNFAIP3, IRF1 and CXCL8 (p < 0.05 each; Figure 8C,F,G,H). In the case of TNFAIP3 and CXCL8, expression in IDMF-treated cells was less than 10% of that in the CTL treatment.
Repression of TNFAIP3 (A20) and CXCL8 (IL8) expression by IDMF may be due to inhibition of the NF-κB pathway [45,46], which was suggested by our microarray analysis of NF-κB-responsive genes (Figure 5). To address this possibility, we transfected human fetal astrocytes with NF-κB reporter assay firefly/Renilla luciferase constructs and treated these cells with TNF-α alone (n = 6) or TNF-α + IDMF (n = 3) for 6 h. NF-κB induction was monitored using a dual-luciferase reporter assay system (i.e., luciferase/Renilla ratios) (Figure 8A). The NF-κB activation signal was reduced by approximately 50% in astrocytes treated with TNF + IDMF as compared to cells treated with TNF alone (p = 0.048, Wilcoxon rank sum test) (Figure 8A). IDMF thus repressed NF-κB induction in TNF-stimulated astrocytes.

2.13. IDMF Undergoes In Vitro Hydrolysis by Carboxylesterase to Generate MMF and IMMF

Carboxylesterase 2 (CES-2) is an endoplasmic reticulum serine esterase mediating first-pass hydrolysis of drugs in the small intestine [47]. We performed an in vitro study to assess fatty acid ester hydrolysis of IDMF by carboxylesterase-2. HPLC analyses showed that IDMF yielded MMF and isosorbide mono-(methyfumarate) (IMMF) as hydrolyzed products (Figure S8A, Supplementary Materials). IDMF was completely converted to MMF and IMMF after 15 min under assay conditions (Figure S8B–E). IDMF can therefore act as an MMF prodrug, similar to DRF, although it also generates a structurally similar compound (IMMF) as a secondary hydrolysate.

3. Discussion

This study used a transcriptomics approach to evaluate direct effects of MMF, DRF, and a novel fumarate (IDMF) on gene expression in cultured astrocytes. The effects of DRF and IDMF differed from those of MMF, although we could identify unique responses to each fumarate, with IDMF having the largest overall effect in terms of the number of differentially expressed genes and modules identified (Table 1). Gene expression responses to IDMF were consistent with NRF2 activation as well as inhibition of the NF-κB and IRF1 pathways (Figure 9 and Table 1). Genes most strongly altered by each fumarate overlapped significantly with genes near loci implicated in MS by genetic association studies, although genes altered by IDMF overlapped more strongly with such genes (Figure 7), and our analysis identified several MS-associated genes uniquely altered by IDMF alone (e.g., down-regulation of TNFAIP3, CXCL8, IL6, and IRF1) (Figure 6). These results demonstrate ways in which fumarate compounds lack bioequivalence in terms of their direct effect on CNS glial cells. This may have implications regarding therapeutic efficacy and potential for adverse effects when different fumarate compounds are used in vivo for treatment of MS or other autoimmune diseases.
The DMF mechanism of action has frequently been understood in terms of NRF2-dependent and NRF2-independent pathways, with the former postulated to mediate a neuroprotective effect on CNS cell types. NRF2 activation by DMF leads to transduction via the ERK1/2 MAPK pathway [48] with binding of NRF2 to antioxidant response elements (ARE) and induction of cytoprotective antioxidant genes having direct neuroprotective effects [49]. The importance of this mechanism was emphasized by early studies revealing that beneficial effects of DMF in a mouse model of chronic experimental autoimmune encephalitis (EAE) are abolished on the Nrf2(−/−) background [50]. DMF and MMF also lead to NRF2-dependent improvement in astrocyte viability following oxidative stress [49]. These effects of the NRF2 pathway are mediated by NRF2 target genes [51], but the current study showed that such targets respond differently to each fumarate compound. This was demonstrated using an externally validated NRF2 gene signature [37]. IDMF was the only fumarate to significantly increase expression of NRF2-increased genes and decrease expression of NRF2-decreased genes (Figure 4). IDMF also led to a potentially synergistic increase in G6PD and GSR expression (Figure 4). Glucose-6-phosphate dehydrogenase (G6PD) encodes a pentose phosphate pathway enzyme responsible for generation of NADPH [52], which is then required by glutathione disulfide reductase (GSR) to generate reduced glutathione [53]. Joint upregulation of both genes by IDMF may thus enhance GSR-dependent glutathione recycling, which is proposed to be a mechanism underlying the antioxidant activity of DMF [49,54,55]. IDMF also uniquely upregulated expression of the NRF2 target heme oxygenase-1 (HMOX1) (also HO-1) (Figure 4B). This gene encodes another NADPH-dependent enzyme, which catalyzes heme degradation to generate the antioxidants biliverdin and bilirubin [56] and further inhibits nitric oxide synthase to limit reactive nitrogen species production [51] (Figure 9). Fibroblasts lacking HMOX1 thus demonstrate hypersensitivity to oxidative stress agents [57] and induction of HMOX1 in astrocytes appears to mediate antioxidant effects of compounds such as resveratrol [58]. Fumarates we tested may thus have differing effects on redox homeostasis via such mechanisms.
The NRF2-independent effects of fumarates comprise a second set of mechanisms associated with inhibition of neuroinflammation with immunomodulatory effects on peripheral immune cells. These mechanisms are highlighted by the ability of DMF to ameliorate symptoms in an acute EAE mouse model on either a wild-type or Nrf2(−/−) genetic background [59]. Inhibition of NF-κB is an effect of DMF that is at least partly NRF2-independent, and we showed that an independently defined and validated group of NF-κB-activated genes [38] is down-regulated by MMF and IDMF but not by DRF (Figure 5). Several of such NF-κB target genes were down-regulated only by IDMF, such as TNFAIP3, NFKBIA, and ICAM1 (Figure 5A,B,D). Down-regulation of NF-κB targets by IDMF is consistent with findings from our earlier study also performed in astrocytes [25]. Both TNFAIP3 and NFKBIA encode inhibitory proteins that suppress the NF-κB pathway, and their down-regulation by IDMF is likely a compensatory response to loss of NF-κB signaling. TNFAIP3 (A20) encodes a deubiquitinating enzyme previously linked to MS by genetic association studies [60]. TNFAIP3 protein abundance is increased in white matter and cortical lesion astrocytes from MS patient brains [61], although mRNA expression is decreased in peripheral monocytes from MS patients [62]. Similarly, NFKBIA (IκBα) encodes an NF-κB inhibitor for which promoter polymorphisms have been associated with MS [63,64], and abundance of this protein is prominent in macrophage nuclei of MS plaques undergoing active demyelination [65]. The unique down-regulation of ICAM1 (CD54) by IDMF in this study may be a direct anti-inflammatory effect not seen with other fumarates [66] (Figure 9 and Table 1). ICAM1 encodes an endothelial adhesion molecule that promotes leucocyte extravasation, and its expression in astrocytes is regulated by NF-κB downstream of stimulation by cytokines such as TNF-α and IL-1β [67,68]. These findings support NF-κB inhibition as a NRF2-independent mechanism by which IDMF may have an anti-inflammatory effect (Figure 9).
Pro-inflammatory cytokines are downstream mediators within the NF-κB signaling pathway inhibited by DMF in a NRF2-independent fashion [69]. Consistent with this, a unique effect of IDMF in this study was down-regulation of the gene encoding IL-6 (Figure 9 and Table 1), which has now emerged as a therapeutic target in multiple autoimmune diseases [70]. In MS patients, IL-6 within the central nervous system is localized to astrocytes and concentrated in demyelinated regions [71,72]. This localized expression of IL-6 in the CNS may be synergistically enhanced by elevated IL-6 production in circulating blood cells [73,74], which can be inhibited by DMF in some cell types (e.g., B cells) [75]. Consistent with these results, a neutralizing IL-6 antibody delayed disease progression in EAE mice by inhibiting microglia activation, CNS inflammation and myelin loss [76], and astrocyte-specific knockout of IL-6 in such mice reduced inflammatory cell infiltration and demyelination [77]. Toxic effects of astrocyte-generated IL-6 have also been supported by mouse models with transgenic over-expression of IL-6 in astrocytes, which exhibit a neurologically diseased phenotype characterized by tremors, ataxia, and seizures, with loss of neuronal subpopulations and increased reactive astrocytosis and neovascularization [78]. These pathological effects of local IL-6 persist in the absence of IL-6 production by circulating blood cells [79]. In this study, neither MMF nor DRF inhibited IL6 mRNA levels in cultured astrocytes (Figure 6). In prior work, MMF was shown to increase IL6 expression in spinal cord from mice following ventral root crush lesion [80], although MMF was also reported to decrease IL6 expression in hippocampus [81], B cells [75], cardiomyocytes [82], liver [83], and intervertebral disk nucleus pulposus cells [84]. Effects of MMF and other fumarates on IL-6 production may therefore vary across contexts and cell types, although our work supports IL6 mRNA down-regulation as a unique effect of IDMF in astrocytes.
Genetic IRF1 variants have been associated with MS [85] and some treatments such as glatiramer acetate block cytokine-stimulated IRF1 up-regulation [86]. Multiple studies have reported resistance of IRF1(−/−) mice to EAE, although the precise cell type mediating this effect in unclear. Prior work has suggested that EAE progression in mice is mediated by the activity of IRF1 in spinal cord [87], myelin-specific T cells [88], and glial cells within the CNS, such as oligodendrocytes [89,90,91]. In this study, IDMF down-regulated IRF1 expression and multiple motifs recognized by IRF1 were enriched in upstream regions of genes down-regulated by DRF and IDMF (Figure S5). Down-regulation of IRF1 may attenuate cytotoxicity in oligodendroglial progenitor cells resulting from IFN-gamma pathway signaling [92,93]. Repression of the IRF1 pathway may also limit pyroptosis in oligodendrocytes by preventing pro-inflammatory caspase protein activation [94]. Along these lines, CASP4 expression was significantly down-regulated by IDMF but not by MMF or DRF. Finally, IRF1 repression is expected to down-regulate BAFF expression, which may attenuate non-canonical NF-κB activation [95] and induce apoptosis of autoreactive T cells [96]. Down-regulation of IRF1 signaling by IDMF may therefore have anti-inflammatory effects at multiple CNS and peripheral immune system sites (Figure 9).
Reactive astrocytosis is a stress-triggered process leading to a spectrum of astrocyte phenotypes with different functional capacities [26]. This spectrum has been conceptually characterized as polarization between inflammation-induced neurotoxic A1 and ischemia-induced neuroprotective A2 phenotypes [97]. Treatments that improve the clinical phenotype of EAE mouse models appear to promote A1 to A2 polarization [98,99,100], which in certain contexts may be mediated by NF-κB inhibition [98,101,102,103,104]. In this study, the best-supported A1 and A2 marker genes were not significantly altered by fumarate compounds (Figure S6A–D) and we did not see consistent directional changes among other genes linked to the A1 and A2 astrocyte phenotypes (Figure S6E,H,K). Indeed, IDMF both increased and decreased the expression of genes linked to the A1 phenotype (Figure S6L) and several such IDMF-decreased genes were localized to the lysosome and associated with catalytic activity (Figure S6N–P). Fumarate compounds thus appear to alter the expression of genes linked to reactive astrogliosis, although the overall pattern of change may not fit neatly within the A1 vs. A2 paradigm. This paradigm may be revised in future work if novel astrocyte sub-phenotypes can be established from studies of the neurodegenerative process [105].
The clinical use of DMF in MS patients has been limited by gastrointestinal (GI) adverse effects, which has prompted development of delayed-release formulations and novel agents such as DRF [22]. These adverse effects are predominantly due to the generation of methanol following first-pass hydrolysis in the small intestine [21]. Whereas orally administered DMF is hydrolyzed to MMF and methanol, DRF metabolism yields additional metabolites including metabolite 2-hydroxyethyl succinimide and RDC-8439 [23]. The generation of methanol following DRF hydrolysis appears to be a minor product, which may explain the improved GI tolerability of DRF as compared to MMF [23]. Our in vitro analysis shows that IDMF is hydrolyzed by carboxylesterase-2 to generate MMF and IMMF (Figure S8). This latter compound is structurally similar to IDMF. IDMF can therefore act as an MMF prodrug, similar to DRF, but may generate additional metabolites that may have IDMF-like bioactivity. Future in vivo studies are needed to evaluate systemic absorption of orally administered IDMF and the pharmacokinetics of its metabolites. This work can guide development of optimal delivery vehicles and further elucidate effects of IDMF and its metabolites on the CNS and peripheral immune cells.

4. Materials and Methods

4.1. Astrocyte Cell Cultures

Experiments were performed using human fetal astrocytes derived from cerebral cortex (cat. no. 882A-05f; Cell Applications, San Diego, CA, USA). Cells were grown to subconfluence in human astrocyte growth medium (cat. no. 820-500, Cell Applications). Stock solutions of MMF, DRF, and IDMF were prepared in DMSO (50 mM) and subsequent dilutions were made using distilled water. The experiment was performed using 12 replicate cultures, each of which was assigned to one of four experimental treatments (CTL, n = 3; MMF, n = 3, DRF, n = 3; IDMF, n = 3). CTL cells were treated with sterile distilled water. Test materials were added to subconfluent cultures at a concentration of 2.5 µM and cells were incubated for 24 h.

4.2. RNA Processing

RNA was extracted following this incubation period using the RNAeasy Mini Plus kit (Qiagen, Hilden, Germany) and robotic Qiacube Connect station (Qiagen). Samples were shipped on dry ice to the Thermo Fisher Microarray Research Services Laboratory (Santa Clara, CA, USA) for transcriptome analysis. RNA quality analysis was performed upon sample receipt using the NanoDrop Lite (Thermo Fisher Scientific, Waltham, MA, USA). All samples had 260/280 nm absorbance ratios higher than 1.90, consistent with high-purity RNA (Figure S8M).

4.3. Microarray Normalization and Quality Control

Gene expression profiling was performed using Clariom S Assays (n = 12 samples). Microarray pseudo-images generated from raw CEL files did not demonstrate evidence for spatial artifact with respect to 11 microarray samples (Figure S9B–L), although there was slight spatial artifact seen with respect to sample CTL-1 (Figure S9A, Supplementary Materials). However, eukaryotic hybridization spike-in controls showed expected trends for each array (Figure S9N), and likewise, 3 out of 4 poly-A RNA labeling controls spiked into RNA samples at varying concentrations showed expected trends (Figure S9O). Area under the curve (AUC) statistics evaluating the degree to which signals differed between exons and introns were close to 1.00 for all samples (≥0.93), consistent with good separation between these two probe set groups (Figure S9P). The distribution of probe-level model residuals was consistent among samples (Figure S9Q). The normalized unscaled standard error (NUSE) [106] median and interquartile range was higher for CTL-1 compared to other samples (Figure S9R,S). The relative log expression (RLE) [106] median was lower for CTL-1 compared to other samples (Figure S9T), whereas the RLE interquartile range was higher (Figure S9U).
The robust multichip average (RMA) algorithm (R package: oligo; R function: rma) [107] was used to generate normalized expression intensities for 27,189 probes. Of these, 19,937 probes were annotated with a protein coding gene, although in some cases more than one probe was annotated with the same human gene symbol. In such cases, we selected the single probe for each gene having the highest average expression across the 12 samples. This yielded a final set of 18,088 probes included in further analyses, with each probe associated with a unique protein-coding human gene symbol. Hierarchical cluster analysis of the 12 samples was performed using expression estimates from these 18,088 genes. This showed that MMF samples grouped together, and likewise IDMF samples were similar to one another (Figure S9V). Cluster analysis did suggest that CTL-1 was a possible outlier (Figure S9V). However, CTL-1 did not emerge as a strong outlier when samples were plotted with respect to the first two principal component axes (Figure S9W). Consistent with this, sample CTL-1 was not a statistically significant outlier with respect to the first principal component axis (p = 0.18, maximum normalized residual test, i.e., Grubb’s test; Figure S9X). Differential expression testing was thus performed using all 12 samples, with exclusion of any genes meeting outlier criteria with respect to the CTL-1 sample (see below).

4.4. Differential Expression Analysis

Tests for differential expression were performed with respect to three two-group comparisons (MMF vs. CTL, DRF vs. CTL, and IDMF vs. CTL). Differential expression tests were performed only for genes with detectable expression in at least 2 of the 6 samples involved in each two-sample comparison, with genes having an expression intensity above the 20th percentile on a given array considered to have detectable expression. Genes were also excluded if their standard deviation across all 12 samples was below the 5th percentile. Finally, given that sample CTL-1 had emerged as a borderline outlier in preliminary analyses (Figure S9), we excluded 132 genes for which expression estimates met outlier criteria with respect to the CTL-1 sample (i.e., absolute z-score greater than 2.50 with a value lower or higher than all other samples). After applying these filters, there remained 14,109, 14,847, and 14,118 genes that were included in differential expression analyses for the MMF vs. CTL, DRF vs. CTL, and IDMF vs. CTL comparisons, respectively.
Differential expression analysis was performed using generalized least square linear models with empirical Bayes moderated t-statistics (R package: limma; R functions: lmFit and eBayes) [35]. To control the false discovery rate, raw p-values were adjusted for multiple hypothesis testing using the Benjamini–Hochberg method [108]. The distributions of raw p-values for the MMF vs. CTL and DRF vs. CTL comparisons were right-skewed, suggesting weak treatment effects (Figure S1A,B). In contrast, the raw p-value distribution for the IDMF vs. CTL comparison was left-skewed, consistent with a stronger treatment effect (Figure S1C). The stronger treatment effect of IDMF compared to MMF and DRF could also be discerned from QQ plots of moderated T statistics (Figure S1D–F) and from raw p-value empirical CDFs (Figure S1G).
Applying a stringent significant threshold (FDR < 0.10 with FC > 1.25 or FC < 0.80), we identified 23 DEGs with respect to the IDMF vs. CTL comparison but did not identify any DEGs with respect to the MMF vs. CTL or DRF vs. CTL comparisons. In further analyses, therefore, we identified DEGs based upon a less stringent threshold (p < 0.05 with FC > 1.25 or FC < 0.80) (Figure S1H–J). We identified approximately twice as many MMF-increased DEGs compared to MMF-decreased DEGs (Figure S1H,K,N), but the numbers of increased and decreased DEGs were similar for both other comparisons (Figure S1I,J,L,M,O,P). For all comparisons, fold-change estimates were not biased among genes having low or high expression (Figure S1N–P).

4.5. Gene Expression Module Analysis

Gene expression modules were evaluated as an alternative to focusing on individual genes. This allowed us to identify differentially expressed modules (DEMs), defined as groups of co-expressed genes consistently altered in the same direction by each fumarate compound [109]. Expression modules were identified based upon an independent set of 24 microarray samples generated from prior studies of iPSC-derived astrocytes treated with TNF-α or α-synuclein proteins (GSE166768 and GSE166769) [110]. This study had utilized the same microarray platform design as in our current study (i.e., Clariom S Assays). Module analyses included only matching probes associated with 12,986 genes having detectable expression such that they were included in each differential expression comparison (MMF vs. CTL, DRF vs. CTL, IDMF vs. CTL).
Raw CEL files were normalized as above using the RMA algorithm [107] and the normalized expression intensity of each gene was centered (by subtracting its average expression across the 24 samples). The resulting normalized/centered expression values were then scaled to have the same mean (0) and standard deviation (1) across all array samples. The 12,986 genes were then clustered using a hierarchical approach with Euclidean distance and average linkage (R function: hclust). The resulting dendrogram was analyzed using a variable height branch pruning technique to detect 233 gene modules (R package: dynamicTreeCut; function: cutreeDynamicTree) [111]. The size of gene modules ranged from 25 to 196 genes. Approximately 7% of genes (929) were not assigned to any module and were thus excluded from further analyses. The 233 modules were labeled using the module size and the single member gene associated with the largest number of Gene Ontology terms. For example, module DCTN1-74 contained 74 genes and of these DCTN1 was associated with the most Gene Ontology terms. To assess differential module expression, a two-sample Wilcoxon rank sum test was used to assess whether FC estimates for module genes differed significantly from all other detectable genes [109]. To correct for multiple hypothesis testing among the 233 modules, raw p-values were adjusted using the Benjamini–Hochberg approach [108].

4.6. DEG Analysis and Annotation

DEGs analysis was performed using enrichment analyses with annotations drawn from several existing databases, including Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) [32,112]. Annotation-based gene set enrichment analysis was performed using a conditional hypergeometric test (R package: GOstats; function: hyperGTest) [113]. We used pre-compiled sets of target genes to evaluate the effects of fumarate compounds on the NRF2 and NF-kB pathways [37,38]. Effects of fumarates on the niacin-HCAR2 pathway were evaluated using data from a prior microarray study in which monocytes were treated with niacin (100 μM) for 6 h (GSE103381) [41]. Enrichment of motifs in DEG upstream regions was assessed using semiparametric generalized additive logistic models [36] with a pre-compiled set of 2935 empirically-determined binding sites (position weight matrix models) associated with mouse or human transcription factors or unconventional DNA binding proteins (uDBPs) [114]. Genes linked to MS via GWA studies were identified using the NHGRI-EBI GWAS catalog [33]. A more comprehensive set of MS-associated genes was identified using an aggregative approach described previously [25], which incorporated MS-associated genes included in 7 database sources (i.e., NHGRI-EBI GWAS catalog, MeSH, disease ontology, DisGeNet, MalaCards, and KEGG) [28,29,30,31,32,33]. DEGs identified in this study were further compared with those with altered expression in astrocytes from normal appearing white matter within brains obtained from MS patients (GSE83670) [34]. Finally, DEGs were compared against those previously identified as markers of A1 and A2 astrogliosis [42], which were identified as the human orthologues of genes previously identified from microarray analyses of astrocytes isolated from brains of mice subjected to inflammatory or ischemic injury, respectively [44].

4.7. Real-Time Quantitative PCR (RT-PCR)

The expression of selected genes was further evaluated using RT-PCR. cDNA was prepared using the AzuraQuant cDNA kit (Azura Genomics, Raynham, MA, USA) and PCR reactions were performed using the BioRad iCycler iQ Detection System with Fast Green qPCR Master Mix—Fluor (Azura Genomics). PCR primers were purchased from Realtimeprimers (Elkins Park, PA, USA). The ΔΔCt method [115] was used to estimate fold changes with cycle threshold values normalized based upon the expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH).

4.8. NF-κB Reporter Assay

Experiments were performed using human fetal astrocytes (HA) derived from cerebral cortex (cat no. 882A-05f; Cell Applications, San Diego, CA, USA) grown to subconfluence for 48 h in HA growth medium (cat. no. 820-500, Cell Applications) in a 96-well black wall tissue culture plate. HA growth medium was replaced by Opti-MEM (Thermo Fisher Scientific, Waltham, MA, USA) on the day of transfection and cells were transfected using NF-κB reporter assay firefly/Renilla luciferase constructs (cat. no. 60616; BPS Bioscience, San Diego, CA, USA) combined with plasmid transfection reagent (Santa Cruz Biotechnology, Dallas, TX, USA) for 18 h. The medium was changed back to HA growth medium and TNF (20 ng/mL) (n = 6 replicates) or TNF (20 ng/mL) + IDMF (2.5 µM) (n = 3) was added. NF-κB induction was quantified 6 h later using a dual-luciferase reporter assay system (cat. no. E1960; Promega, Madison, WI, USA). Signal quantification was obtained using a Thermo Fisher Scientific Luminoskan Ascent Microplate Luminometer. The instrument had passed DLReady™ validation at Promega Corporation and was certified as dual-luciferase assay-ready. NF-κB inhibition was quantified based upon the ratio between luciferase and Renilla signals, with the latter signal providing an internal control for transfection efficiency and cell viability.

4.9. Fatty Acid Ester Hydrolysis of IDMF with Carboxylesterase-2

IDMF (1 mM) was dissolved in DMSO (2%) and ethylene glycol (2%) in 20 mM HEPES (pH 7.4). We added 20 μL of enzyme (1 mg/mL) and the solution was incubated at 37 °C for 15, 30, 60, and 120 min. The reaction was stopped using Methanol. IDMF and MMF were quantified using a UV detector at a wavelength of 210 nm. HPLC was performed with 0.1% formic acid in water as the weak solvent and methanol as the strong solvent. The analysis was performed using a Luna C18 4.6 mm × 100 mm column with a flow rate of 1.0 mL/min (Phenomenex, Torrance, CA, USA).
The molecular weight of IMMF was determined using the API3000 triple-quadrupole mass spectrometer (AB Sciex, Framingham, MA, USA) and electrospray ionization (ESI). The Q1 scan was performed from 100–1500 Da with negative mode. MS analysis conditions were as follows: curtain gas (10 arbitrary units), nebulizer gas (14 arbitrary units), ion spray voltage (−4500 V), channel electron multiplier (CEM) at 2800 V. Declustering and entrance potentials were −130 V and −10 V, respectively.

5. Conclusions

Fumarate esters have been approved as treatments for RRMS and psoriasis but work is now underway to explore their utility as treatments for diverse other conditions, such as chronic pain, Parkinson’s disease, Alzheimer’s disease, stroke, glioblastoma multiforme, chronic lymphocytic leukemia, cutaneous T cell lymphoma, obstructive sleep apnea, rheumatoid arthritis, and atherosclerosis [27,116,117]. Development of next-generation fumarate compounds with improved properties is also ongoing although little work has been done so far to provide a side-by-side comparison of their effects. This study used transcriptomics to show that fumarates have differing effects on the antioxidant and inflammatory genes acting downstream of the NRF2, NF-κB, IRF1, and HCAR2 (GPR109A) pathways (Figure 9 and Table 1). Likewise, fumarates had differing effects on the expression of key genes that have been linked to multiple sclerosis (e.g., TNFAIP3, CXCL8, IL6, and IRF1). These findings highlight distinguishing features of fumarate compounds that may be consequential for their use as treatments for multiple sclerosis or other autoimmune diseases. This side-by-side comparison helps to define functional properties of different fumarate compounds with the goal of developing new therapies to meet the needs of MS patients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ph15040461/s1, Figure S1: Differential expression analyses; Figure S2: Motifs enriched in 5 kb sequences upstream of DEGs; Figure S3: NF-κB motif enrichment; Figure S4: Niacin-HCAR2 pathway genes; Figure S5: Interferon regulatory factor (IRF) motif enrichment; Figure S6: Genes associated with A1/A2 reactive astrocyte polarization; Figure S7: Genes altered in MS patient astrocytes (GSE83670) and their overlap with DEGs; Figure S8: In vitro kinetics of IDMF hydrolysis by carboxylesterase-2; Figure S9: Microarray quality control; Table S1. MMF-increased genes (p < 0.05 with FC > 1.25); Table S2. MMF-decreased genes (p < 0.05 with FC < 0.80); Table S3. DRF-increased genes (p < 0.05 with FC > 1.25); Table S4. DRF-decreased genes (p < 0.05 with FC < 0.80); Table S5. IDMF-increased genes (p < 0.05 with FC > 1.25); Table S6. IDMF-decreased genes (p < 0.05 with FC < 0.80).

Author Contributions

Conceptualization, K.B. and R.K.C.; validation, K.B. and R.K.C.; formal Analysis, W.R.S. and K.B.; investigation, W.R.S. and K.B.; resources, K.B. and R.K.C.; data curation, W.R.S.; writing—original draft preparation, W.R.S.; writing—review and editing, W.R.S., K.B. and R.K.C.; visualization, W.R.S.; supervision, R.K.C.; project administration, R.K.C.; funding acquisition, R.K.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Sytheon, Ltd. (Boonton, NJ, USA) and Symbionyx Pharmaceuticals (Boonton, NJ, USA).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The whole genome microarray data have been submitted to Gene Expression Omnibus (GEO) and are available under the accession GSE200539.

Acknowledgments

We thank Geovani Quijas, B.S. for technical assistance with the RT-PCR experiments.

Conflicts of Interest

WRS has received consulting reimbursement from Sytheon, Ltd., which provided funding for this study. WRS is a Symbionyx Pharmaceuticals shareholder and member of the Symbionyx advisory board. KB is CEO of Sunny BioDiscovery, Inc. (Santa Paula, CA, USA), has received consulting reimbursement from Sytheon, Ltd., and is Chief Scientific Officer of Symbionyx Pharmaceuticals. RKC is CEO and founder of Sytheon, Ltd., and Symbionyx Pharmaceuticals, with ownership interest in both companies. KB and RKC are listed as inventors of IDMF on composition of matter and application patents (US 10,597,402 and EP 3503866).

References

  1. Reich, D.S.; Lucchinetti, C.F.; Calabresi, P.A. Multiple Sclerosis. N. Engl. J. Med. 2018, 378, 169–180. [Google Scholar] [CrossRef]
  2. Walton, C.; King, R.; Rechtman, L.; Kaye, W.; Leray, E.; Marrie, R.A.; Robertson, N.; La Rocca, N.; Uitdehaag, B.; van der Mei, I.; et al. Rising prevalence of multiple sclerosis worldwide: Insights from the Atlas of MS, third edition. Mult. Scler. J. 2020, 26, 1816–1821. [Google Scholar] [CrossRef] [PubMed]
  3. Wallin, M.T.; Culpepper, W.J.; Nichols, E.; Bhutta, Z.A.; Gebrehiwot, T.T.; Hay, S.I.; Khalil, I.A.; Krohn, K.J.; Liang, X.; Naghavi, M.; et al. Global, regional, and national burden of multiple sclerosis 1990–2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet. Neurol. 2019, 18, 269–285. [Google Scholar] [CrossRef] [Green Version]
  4. Simpson, S., Jr.; Wang, W.; Otahal, P.; Blizzard, L.; van der Mei, I.A.F.; Taylor, B.V. Latitude continues to be significantly associated with the prevalence of multiple sclerosis: An updated meta-analysis. J. Neurol. Neurosurg. Psychiatry 2019, 90, 1193–1200. [Google Scholar] [CrossRef] [PubMed]
  5. Ortiz, G.G.; Pacheco-Moisés, F.P.; Macías-Islas, M.; Flores-Alvarado, L.J.; Mireles-Ramírez, M.A.; González-Renovato, E.D.; Hernández-Navarro, V.E.; Sánchez-López, A.L.; Alatorre-Jiménez, M.A. Role of the blood-brain barrier in multiple sclerosis. Arch. Med. Res. 2014, 45, 687–697. [Google Scholar] [CrossRef] [PubMed]
  6. Severson, C.; Hafler, D.A. T-cells in multiple sclerosis. Results Probl. Cell Differ. 2010, 51, 75–98. [Google Scholar] [CrossRef]
  7. Wang, K.; Song, F.; Fernandez-Escobar, A.; Luo, G.; Wang, J.H.; Sun, Y. The Properties of Cytokines in Multiple Sclerosis: Pros and Cons. Am. J. Med. Sci. 2018, 356, 552–560. [Google Scholar] [CrossRef]
  8. Hartung, H.P.; Jung, S.; Stoll, G.; Zielasek, J.; Schmidt, B.; Archelos, J.J.; Toyka, K.V. Inflammatory mediators in demyelinating disorders of the CNS and PNS. J. Neuroimmunol. 1992, 40, 197–210. [Google Scholar] [CrossRef]
  9. Varatharaj, A.; Galea, I. The blood-brain barrier in systemic inflammation. Brain Behav. Immun. 2017, 60, 1–12. [Google Scholar] [CrossRef] [Green Version]
  10. Yeung, M.S.Y.; Djelloul, M.; Steiner, E.; Bernard, S.; Salehpour, M.; Possnert, G.; Brundin, L.; Frisén, J. Dynamics of oligodendrocyte generation in multiple sclerosis. Nature 2019, 566, 538–542. [Google Scholar] [CrossRef]
  11. Tommasin, S.; Giannì, C.; De Giglio, L.; Pantano, P. Neuroimaging Techniques to Assess Inflammation in Multiple Sclerosis. Neuroscience 2019, 403, 4–16. [Google Scholar] [CrossRef] [PubMed]
  12. Saguil, A.; Kane, S.; Farnell, E. Multiple sclerosis: A primary care perspective. Am. Fam. Physician 2014, 90, 644–652. [Google Scholar] [PubMed]
  13. Hurwitz, B.J. The diagnosis of multiple sclerosis and the clinical subtypes. Ann. Indian Acad. Neurol. 2009, 12, 226–230. [Google Scholar] [CrossRef] [PubMed]
  14. Wiendl, H.; Gold, R.; Berger, T.; Derfuss, T.; Linker, R.; Mäurer, M.; Aktas, O.; Baum, K.; Berghoff, M.; Bittner, S.; et al. Multiple Sclerosis Therapy Consensus Group (MSTCG): Position statement on disease-modifying therapies for multiple sclerosis (white paper). Ther. Adv. Neurol. Disord. 2021, 14, 17562864211039648. [Google Scholar] [CrossRef]
  15. Zecca, C.; Caporro, M.; Adami, M.; Mainetti, C.; Gobbi, C. Fumaric acid esters in psoriasis and multiple sclerosis. Clin. Exp. Dermatol. 2014, 39, 488–491. [Google Scholar] [CrossRef]
  16. Gold, R.; Kappos, L.; Arnold, D.L.; Bar-Or, A.; Giovannoni, G.; Selmaj, K.; Tornatore, C.; Sweetser, M.T.; Yang, M.; Sheikh, S.I.; et al. Placebo-controlled phase 3 study of oral BG-12 for relapsing multiple sclerosis. N. Engl. J. Med. 2012, 367, 1098–1107. [Google Scholar] [CrossRef] [Green Version]
  17. Fox, R.J.; Miller, D.H.; Phillips, J.T.; Hutchinson, M.; Havrdova, E.; Kita, M.; Yang, M.; Raghupathi, K.; Novas, M.; Sweetser, M.T.; et al. Placebo-controlled phase 3 study of oral BG-12 or glatiramer in multiple sclerosis. N. Engl. J. Med. 2012, 367, 1087–1097. [Google Scholar] [CrossRef] [Green Version]
  18. Naismith, R.T.; Wolinsky, J.S.; Wundes, A.; LaGanke, C.; Arnold, D.L.; Obradovic, D.; Freedman, M.S.; Gudesblatt, M.; Ziemssen, T.; Kandinov, B.; et al. Diroximel fumarate (DRF) in patients with relapsing-remitting multiple sclerosis: Interim safety and efficacy results from the phase 3 EVOLVE-MS-1 study. Mult. Scler. J. 2020, 26, 1729–1739. [Google Scholar] [CrossRef]
  19. Min, J.; Cohan, S.; Alvarez, E.; Sloane, J.; Phillips, J.T.; van der Walt, A.; Koulinska, I.; Fang, F.; Miller, C.; Chan, A. Real-World Characterization of Dimethyl Fumarate-Related Gastrointestinal Events in Multiple Sclerosis: Management Strategies to Improve Persistence on Treatment and Patient Outcomes. Neurol. Ther. 2019, 8, 109–119. [Google Scholar] [CrossRef] [Green Version]
  20. Vollmer, B.; Nair, K.V.; Sillau, S.H.; Corboy, J.; Vollmer, T.; Alvarez, E. Comparison of fingolimod and dimethyl fumarate in the treatment of multiple sclerosis: Two-year experience. Mult. Scler. J.-Exp. Transl. Clin. 2017, 3, 2055217317725102. [Google Scholar] [CrossRef] [Green Version]
  21. Jonasson, E.; Sejbaek, T. Diroximel fumarate in the treatment of multiple sclerosis. Neurodegener. Dis. Manag. 2020, 10, 267–276. [Google Scholar] [CrossRef] [PubMed]
  22. Naismith, R.T.; Wundes, A.; Ziemssen, T.; Jasinska, E.; Freedman, M.S.; Lembo, A.J.; Selmaj, K.; Bidollari, I.; Chen, H.; Hanna, J.; et al. Diroximel Fumarate Demonstrates an Improved Gastrointestinal Tolerability Profile Compared with Dimethyl Fumarate in Patients with Relapsing-Remitting Multiple Sclerosis: Results from the Randomized, Double-Blind, Phase III EVOLVE-MS-2 Study. CNS Drugs 2020, 34, 185–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Palte, M.J.; Wehr, A.; Tawa, M.; Perkin, K.; Leigh-Pemberton, R.; Hanna, J.; Miller, C.; Penner, N. Improving the Gastrointestinal Tolerability of Fumaric Acid Esters: Early Findings on Gastrointestinal Events with Diroximel Fumarate in Patients with Relapsing-Remitting Multiple Sclerosis from the Phase 3, Open-Label EVOLVE-MS-1 Study. Adv. Ther. 2019, 36, 3154–3165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Bojanowski, K.; Ibeji, C.U.; Singh, P.; Swindell, W.R.; Chaudhuri, R.K.J.J.I. A Sensitization-Free Dimethyl Fumarate Prodrug, Isosorbide Di-(Methyl Fumarate), Provides a Topical Treatment Candidate for Psoriasis. JID Innov. 2021, 1, 100040. [Google Scholar] [CrossRef]
  25. Swindell, W.R.; Bojanowski, K.; Chaudhuri, R.K. A novel fumarate, isosorbide di-(methyl fumarate) (IDMF), replicates astrocyte transcriptome responses to dimethyl fumarate (DMF) but specifically down-regulates genes linked to a reactive phenotype. Biochem. Biophys. Res. Commun. 2020, 532, 475–481. [Google Scholar] [CrossRef]
  26. Nutma, E.; van Gent, D.; Amor, S.; Peferoen, L.A.N. Astrocyte and Oligodendrocyte Cross-Talk in the Central Nervous System. Cells 2020, 9, 600. [Google Scholar] [CrossRef] [Green Version]
  27. Hoogendoorn, A.; Avery, T.D.; Li, J.; Bursill, C.; Abell, A.; Grace, P.M. Emerging Therapeutic Applications for Fumarates. Trends Pharmacol. Sci. 2021, 42, 239–254. [Google Scholar] [CrossRef]
  28. Tsuyuzaki, K.; Morota, G.; Ishii, M.; Nakazato, T.; Miyazaki, S.; Nikaido, I. MeSH ORA framework: R/Bioconductor packages to support MeSH over-representation analysis. BMC Bioinform. 2015, 16, 45. [Google Scholar] [CrossRef] [Green Version]
  29. Kibbe, W.A.; Arze, C.; Felix, V.; Mitraka, E.; Bolton, E.; Fu, G.; Mungall, C.J.; Binder, J.X.; Malone, J.; Vasant, D.; et al. Disease Ontology 2015 update: An expanded and updated database of human diseases for linking biomedical knowledge through disease data. Nucleic Acids Res. 2015, 43, D1071–D1078. [Google Scholar] [CrossRef] [Green Version]
  30. Pinero, J.; Queralt-Rosinach, N.; Bravo, A.; Deu-Pons, J.; Bauer-Mehren, A.; Baron, M.; Sanz, F.; Furlong, L.I. DisGeNET: A discovery platform for the dynamical exploration of human diseases and their genes. Database J. Biol. Databases Curation 2015, 2015, bav028. [Google Scholar] [CrossRef]
  31. Rappaport, N.; Twik, M.; Plaschkes, I.; Nudel, R.; Iny Stein, T.; Levitt, J.; Gershoni, M.; Morrey, C.P.; Safran, M.; Lancet, D. MalaCards: An amalgamated human disease compendium with diverse clinical and genetic annotation and structured search. Nucleic Acids Res. 2017, 45, D877–D887. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016, 44, D457–D462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Welter, D.; MacArthur, J.; Morales, J.; Burdett, T.; Hall, P.; Junkins, H.; Klemm, A.; Flicek, P.; Manolio, T.; Hindorff, L.; et al. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014, 42, D1001–D1006. [Google Scholar] [CrossRef] [PubMed]
  34. Waller, R.; Woodroofe, M.N.; Wharton, S.B.; Ince, P.G.; Francese, S.; Heath, P.R.; Cudzich-Madry, A.; Thomas, R.H.; Rounding, N.; Sharrack, B.; et al. Gene expression profiling of the astrocyte transcriptome in multiple sclerosis normal appearing white matter reveals a neuroprotective role. J. Neuroimmunol. 2016, 299, 139–146. [Google Scholar] [CrossRef] [PubMed]
  35. Smyth, G.K. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 2004, 3, 1–25. [Google Scholar] [CrossRef]
  36. Swindell, W.R.; Johnston, A.; Xing, X.; Little, A.; Robichaud, P.; Voorhees, J.J.; Fisher, G.; Gudjonsson, J.E. Robust shifts in S100a9 expression with aging: A novel mechanism for chronic inflammation. Sci. Rep. 2013, 3, 1215. [Google Scholar] [CrossRef] [Green Version]
  37. Rooney, J.P.; Chorley, B.; Hiemstra, S.; Wink, S.; Wang, X.; Bell, D.A.; van de Water, B.; Corton, J.C. Mining a human transcriptome database for chemical modulators of NRF2. PLoS ONE 2020, 15, e0239367. [Google Scholar] [CrossRef] [PubMed]
  38. Zhao, M.; Joy, J.; Zhou, W.; De, S.; Wood, W.H., 3rd; Becker, K.G.; Ji, H.; Sen, R. Transcriptional outcomes and kinetic patterning of gene expression in response to NF-κB activation. PLoS Biol. 2018, 16, e2006347. [Google Scholar] [CrossRef] [Green Version]
  39. Maciejewski-Lenoir, D.; Richman, J.G.; Hakak, Y.; Gaidarov, I.; Behan, D.P.; Connolly, D.T. Langerhans cells release prostaglandin D2 in response to nicotinic acid. J. Investig. Dermatol. 2006, 126, 2637–2646. [Google Scholar] [CrossRef] [Green Version]
  40. Geisler, C.E.; Miller, K.E.; Ghimire, S.; Renquist, B.J. The Role of GPR109a Signaling in Niacin Induced Effects on Fed and Fasted Hepatic Metabolism. Int. J. Mol. Sci. 2021, 22, 4001. [Google Scholar] [CrossRef]
  41. Sarkar, S.; Yang, R.; Mirzaei, R.; Rawji, K.; Poon, C.; Mishra, M.K.; Zemp, F.J.; Bose, P.; Kelly, J.; Dunn, J.F.; et al. Control of brain tumor growth by reactivating myeloid cells with niacin. Sci. Transl. Med. 2020, 12, eaay9924. [Google Scholar] [CrossRef] [PubMed]
  42. Ding, Z.B.; Song, L.J.; Wang, Q.; Kumar, G.; Yan, Y.Q.; Ma, C.G. Astrocytes: A double-edged sword in neurodegenerative diseases. Neural Regen. Res. 2021, 16, 1702–1710. [Google Scholar] [CrossRef] [PubMed]
  43. Liddelow, S.A.; Guttenplan, K.A.; Clarke, L.E.; Bennett, F.C.; Bohlen, C.J.; Schirmer, L.; Bennett, M.L.; Münch, A.E.; Chung, W.S.; Peterson, T.C.; et al. Neurotoxic reactive astrocytes are induced by activated microglia. Nature 2017, 541, 481–487. [Google Scholar] [CrossRef] [PubMed]
  44. Zamanian, J.L.; Xu, L.; Foo, L.C.; Nouri, N.; Zhou, L.; Giffard, R.G.; Barres, B.A. Genomic analysis of reactive astrogliosis. J. Neurosci. Off. J. Soc. Neurosci. 2012, 32, 6391–6410. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Lai, T.Y.; Wu, S.D.; Tsai, M.H.; Chuang, E.Y.; Chuang, L.L.; Hsu, L.C.; Lai, L.C. Transcription of Tnfaip3 is regulated by NF-κB and p38 via C/EBPβ in activated macrophages. PLoS ONE 2013, 8, e73153. [Google Scholar] [CrossRef] [PubMed]
  46. Mukaida, N.; Okamoto, S.; Ishikawa, Y.; Matsushima, K. Molecular mechanism of interleukin-8 gene expression. J. Leukoc. Biol. 1994, 56, 554–558. [Google Scholar] [CrossRef] [PubMed]
  47. Imai, T. Human carboxylesterase isozymes: Catalytic properties and rational drug design. Drug Metab. Pharmacokinet. 2006, 21, 173–185. [Google Scholar] [CrossRef]
  48. Wang, Q.; Chuikov, S.; Taitano, S.; Wu, Q.; Rastogi, A.; Tuck, S.J.; Corey, J.M.; Lundy, S.K.; Mao-Draayer, Y. Dimethyl Fumarate Protects Neural Stem/Progenitor Cells and Neurons from Oxidative Damage through Nrf2-ERK1/2 MAPK Pathway. Int. J. Mol. Sci. 2015, 16, 13885–13907. [Google Scholar] [CrossRef] [Green Version]
  49. Scannevin, R.H.; Chollate, S.; Jung, M.Y.; Shackett, M.; Patel, H.; Bista, P.; Zeng, W.; Ryan, S.; Yamamoto, M.; Lukashev, M.; et al. Fumarates promote cytoprotection of central nervous system cells against oxidative stress via the nuclear factor (erythroid-derived 2)-like 2 pathway. J. Pharmacol. Exp. Ther. 2012, 341, 274–284. [Google Scholar] [CrossRef] [Green Version]
  50. Linker, R.A.; Lee, D.H.; Ryan, S.; van Dam, A.M.; Conrad, R.; Bista, P.; Zeng, W.; Hronowsky, X.; Buko, A.; Chollate, S.; et al. Fumaric acid esters exert neuroprotective effects in neuroinflammation via activation of the Nrf2 antioxidant pathway. Brain A J. Neurol. 2011, 134, 678–692. [Google Scholar] [CrossRef] [Green Version]
  51. Wakabayashi, N.; Slocum, S.L.; Skoko, J.J.; Shin, S.; Kensler, T.W. When NRF2 talks, who’s listening? Antioxid. Redox Signal. 2010, 13, 1649–1663. [Google Scholar] [CrossRef] [Green Version]
  52. Chen, L.; Zhang, C.; Wang, Y.; Li, Y.; Han, Q.; Yang, H.; Zhu, Y. Data mining and pathway analysis of glucose-6-phosphate dehydrogenase with natural language processing. Mol. Med. Rep. 2017, 16, 1900–1910. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ferreira, B.; Mendes, F.; Osório, N.; Caseiro, A.; Gabriel, A.; Valado, A. Glutathione in multiple sclerosis. Br. J. Biomed. Sci. 2013, 70, 75–79. [Google Scholar] [CrossRef] [PubMed]
  54. Albrecht, P.; Bouchachia, I.; Goebels, N.; Henke, N.; Hofstetter, H.H.; Issberner, A.; Kovacs, Z.; Lewerenz, J.; Lisak, D.; Maher, P.; et al. Effects of dimethyl fumarate on neuroprotection and immunomodulation. J. Neuroinflammation 2012, 9, 163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Hoffmann, C.; Dietrich, M.; Herrmann, A.K.; Schacht, T.; Albrecht, P.; Methner, A. Dimethyl Fumarate Induces Glutathione Recycling by Upregulation of Glutathione Reductase. Oxidative Med. Cell. Longev. 2017, 2017, 6093903. [Google Scholar] [CrossRef]
  56. Loboda, A.; Damulewicz, M.; Pyza, E.; Jozkowicz, A.; Dulak, J. Role of Nrf2/HO-1 system in development, oxidative stress response and diseases: An evolutionarily conserved mechanism. Cell. Mol. Life Sci. CMLS 2016, 73, 3221–3247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Poss, K.D.; Tonegawa, S. Reduced stress defense in heme oxygenase 1-deficient cells. Proc. Natl. Acad. Sci. USA 1997, 94, 10925–10930. [Google Scholar] [CrossRef] [Green Version]
  58. Quincozes-Santos, A.; Bobermin, L.D.; Latini, A.; Wajner, M.; Souza, D.O.; Gonçalves, C.A.; Gottfried, C. Resveratrol protects C6 astrocyte cell line against hydrogen peroxide-induced oxidative stress through heme oxygenase 1. PLoS ONE 2013, 8, e64372. [Google Scholar] [CrossRef] [Green Version]
  59. Schulze-Topphoff, U.; Varrin-Doyer, M.; Pekarek, K.; Spencer, C.M.; Shetty, A.; Sagan, S.A.; Cree, B.A.; Sobel, R.A.; Wipke, B.T.; Steinman, L.; et al. Dimethyl fumarate treatment induces adaptive and innate immune modulation independent of Nrf2. Proc. Natl. Acad. Sci. USA 2016, 113, 4777–4782. [Google Scholar] [CrossRef] [Green Version]
  60. Hoffjan, S.; Okur, A.; Epplen, J.T.; Wieczorek, S.; Chan, A.; Akkad, D.A. Association of TNFAIP3 and TNFRSF1A variation with multiple sclerosis in a German case-control cohort. Int. J. Immunogenet. 2015, 42, 106–110. [Google Scholar] [CrossRef]
  61. Perga, S.; Montarolo, F.; Martire, S.; Bonaldo, B.; Bono, G.; Bertolo, J.; Magliozzi, R.; Bertolotto, A. Overexpression of the ubiquitin-editing enzyme A20 in the brain lesions of Multiple Sclerosis patients: Moving from systemic to central nervous system inflammation. Brain Pathol. 2021, 31, 283–296. [Google Scholar] [CrossRef] [PubMed]
  62. Navone, N.D.; Perga, S.; Martire, S.; Berchialla, P.; Malucchi, S.; Bertolotto, A. Monocytes and CD4+ T cells contribution to the under-expression of NR4A2 and TNFAIP3 genes in patients with multiple sclerosis. J. Neuroimmunol. 2014, 272, 99–102. [Google Scholar] [CrossRef] [PubMed]
  63. Balood, M.; Mesbah-Namin, S.A.; Sanati, M.H.; Zahednasab, H.; Sahraian, M.A.; Ataei, M. Inhibitor IκBα promoter functional polymorphisms in patients with multiple sclerosis. Mol. Biol. Rep. 2014, 41, 613–616. [Google Scholar] [CrossRef]
  64. Miterski, B.; Böhringer, S.; Klein, W.; Sindern, E.; Haupts, M.; Schimrigk, S.; Epplen, J.T. Inhibitors in the NFkappaB cascade comprise prime candidate genes predisposing to multiple sclerosis, especially in selected combinations. Genes Immun. 2002, 3, 211–219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Gveric, D.; Kaltschmidt, C.; Cuzner, M.L.; Newcombe, J. Transcription factor NF-kappaB and inhibitor I kappaBalpha are localized in macrophages in active multiple sclerosis lesions. J. Neuropathol. Exp. Neurol. 1998, 57, 168–178. [Google Scholar] [CrossRef] [Green Version]
  66. Haarmann, A.; Nehen, M.; Deiß, A.; Buttmann, M. Fumaric Acid Esters Do Not Reduce Inflammatory NF-κB/p65 Nuclear Translocation, ICAM-1 Expression and T-Cell Adhesiveness of Human Brain Microvascular Endothelial Cells. Int. J. Mol. Sci. 2015, 16, 19086–19095. [Google Scholar] [CrossRef] [Green Version]
  67. Lee, J.K.; Kim, J.K.; Park, S.H.; Sim, Y.B.; Jung, J.S.; Suh, H.W. Lactosylceramide Mediates the Expression of Adhesion Molecules in TNF-α and IFNγ-stimulated Primary Cultured Astrocytes. Korean J. Physiol. Pharmacol. 2011, 15, 251–258. [Google Scholar] [CrossRef] [Green Version]
  68. Lee, S.J.; Park, J.Y.; Hou, J.; Benveniste, E.N. Transcriptional regulation of the intercellular adhesion molecule-1 gene by proinflammatory cytokines in human astrocytes. Glia 1999, 25, 21–32. [Google Scholar] [CrossRef]
  69. McGuire, V.A.; Ruiz-Zorrilla Diez, T.; Emmerich, C.H.; Strickson, S.; Ritorto, M.S.; Sutavani, R.V.; Weiβ, A.; Houslay, K.F.; Knebel, A.; Meakin, P.J.; et al. Dimethyl fumarate blocks pro-inflammatory cytokine production via inhibition of TLR induced M1 and K63 ubiquitin chain formation. Sci. Rep. 2016, 6, 31159. [Google Scholar] [CrossRef] [Green Version]
  70. Smolen, J.S.; Beaulieu, A.; Rubbert-Roth, A.; Ramos-Remus, C.; Rovensky, J.; Alecock, E.; Woodworth, T.; Alten, R. Effect of interleukin-6 receptor inhibition with tocilizumab in patients with rheumatoid arthritis (OPTION study): A double-blind, placebo-controlled, randomised trial. Lancet 2008, 371, 987–997. [Google Scholar] [CrossRef]
  71. Maimone, D.; Guazzi, G.C.; Annunziata, P. IL-6 detection in multiple sclerosis brain. J. Neurol. Sci. 1997, 146, 59–65. [Google Scholar] [CrossRef]
  72. Schönrock, L.M.; Gawlowski, G.; Brück, W. Interleukin-6 expression in human multiple sclerosis lesions. Neurosci. Lett. 2000, 294, 45–48. [Google Scholar] [CrossRef]
  73. Graber, J.J.; Allie, S.R.; Mullen, K.M.; Jones, M.V.; Wang, T.; Krishnan, C.; Kaplin, A.I.; Nath, A.; Kerr, D.A.; Calabresi, P.A. Interleukin-17 in transverse myelitis and multiple sclerosis. J. Neuroimmunol. 2008, 196, 124–132. [Google Scholar] [CrossRef]
  74. Li, R.; Rezk, A.; Miyazaki, Y.; Hilgenberg, E.; Touil, H.; Shen, P.; Moore, C.S.; Michel, L.; Althekair, F.; Rajasekharan, S.; et al. Proinflammatory GM-CSF-producing B cells in multiple sclerosis and B cell depletion therapy. Sci. Transl. Med. 2015, 7, 310ra166. [Google Scholar] [CrossRef]
  75. Li, R.; Rezk, A.; Ghadiri, M.; Luessi, F.; Zipp, F.; Li, H.; Giacomini, P.S.; Antel, J.; Bar-Or, A. Dimethyl Fumarate Treatment Mediates an Anti-Inflammatory Shift in B Cell Subsets of Patients with Multiple Sclerosis. J. Immunol. 2017, 198, 691–698. [Google Scholar] [CrossRef] [PubMed]
  76. Savarin, C.; Hinton, D.R.; Valentin-Torres, A.; Chen, Z.; Trapp, B.D.; Bergmann, C.C.; Stohlman, S.A. Astrocyte response to IFN-γ limits IL-6-mediated microglia activation and progressive autoimmune encephalomyelitis. J. Neuroinflammation 2015, 12, 79. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Erta, M.; Giralt, M.; Jiménez, S.; Molinero, A.; Comes, G.; Hidalgo, J. Astrocytic IL-6 Influences the Clinical Symptoms of EAE in Mice. Brain Sci. 2016, 6, 15. [Google Scholar] [CrossRef] [Green Version]
  78. Campbell, I.L.; Abraham, C.R.; Masliah, E.; Kemper, P.; Inglis, J.D.; Oldstone, M.B.; Mucke, L. Neurologic disease induced in transgenic mice by cerebral overexpression of interleukin 6. Proc. Natl. Acad. Sci. USA 1993, 90, 10061–10065. [Google Scholar] [CrossRef] [Green Version]
  79. Giralt, M.; Ramos, R.; Quintana, A.; Ferrer, B.; Erta, M.; Castro-Freire, M.; Comes, G.; Sanz, E.; Unzeta, M.; Pifarré, P.; et al. Induction of atypical EAE mediated by transgenic production of IL-6 in astrocytes in the absence of systemic IL-6. Glia 2013, 61, 587–600. [Google Scholar] [CrossRef]
  80. Carvalho, N.Z.M.; Chiarotto, G.B.; Bernardes, D.; Kempe, P.R.G.; Oliveira, A.L.R. Neuroprotection by dimethyl fumarate following ventral root crush in C57BL/6J mice. Brain Res. Bull. 2020, 164, 184–197. [Google Scholar] [CrossRef]
  81. Singh, N.; Saha, L.; Kumari, P.; Singh, J.; Bhatia, A.; Banerjee, D.; Chakrabarti, A. Effect of dimethyl fumarate on neuroinflammation and apoptosis in pentylenetetrazol kindling model in rats. Brain Res. Bull. 2019, 144, 233–245. [Google Scholar] [CrossRef] [PubMed]
  82. Zhao, J.; Cheng, Z.; Quan, X.; Xie, Z.; Zhang, L.; Ding, Z. Dimethyl fumarate protects cardiomyocytes against oxygen-glucose deprivation/reperfusion (OGD/R)-induced inflammatory response and damages via inhibition of Egr-1. Int. Immunopharmacol. 2020, 86, 106733. [Google Scholar] [CrossRef] [PubMed]
  83. Dwivedi, D.K.; Jena, G.; Kumar, V. Dimethyl fumarate protects thioacetamide-induced liver damage in rats: Studies on Nrf2, NLRP3, and NF-κB. J. Biochem. Mol. Toxicol. 2020, 34, e22476. [Google Scholar] [CrossRef] [PubMed]
  84. Zhu, H.; Chen, G.; Wang, Y.; Lin, X.; Zhou, J.; Wang, Z.; Suo, N. Dimethyl fumarate protects nucleus pulposus cells from inflammation and oxidative stress and delays the intervertebral disc degeneration. Exp. Ther. Med. 2020, 20, 269. [Google Scholar] [CrossRef]
  85. Fortunato, G.; Calcagno, G.; Bresciamorra, V.; Salvatore, E.; Filla, A.; Capone, S.; Liguori, R.; Borelli, S.; Gentile, I.; Borrelli, F.; et al. Multiple sclerosis and hepatitis C virus infection are associated with single nucleotide polymorphisms in interferon pathway genes. J. Interferon Cytokine Res. 2008, 28, 141–152. [Google Scholar] [CrossRef]
  86. Lu, H.; Zeng, C.; Zhao, H.; Lian, L.; Dai, Y. Glatiramer acetate inhibits degradation of collagen II by suppressing the activity of interferon regulatory factor-1. Biochem. Biophys. Res. Commun. 2014, 448, 323–328. [Google Scholar] [CrossRef]
  87. Tada, Y.; Ho, A.; Matsuyama, T.; Mak, T.W. Reduced incidence and severity of antigen-induced autoimmune diseases in mice lacking interferon regulatory factor-1. J. Exp. Med. 1997, 185, 231–238. [Google Scholar] [CrossRef] [Green Version]
  88. Buch, T.; Uthoff-Hachenberg, C.; Waisman, A. Protection from autoimmune brain inflammation in mice lacking IFN-regulatory factor-1 is associated with Th2-type cytokines. Int. Immunol. 2003, 15, 855–859. [Google Scholar] [CrossRef]
  89. Ren, Z.; Wang, Y.; Tao, D.; Liebenson, D.; Liggett, T.; Goswami, R.; Stefoski, D.; Balabanov, R. Central nervous system expression of interferon regulatory factor 1 regulates experimental autoimmune encephalomyelitis. J. Neuroimmune Pharmacol. 2010, 5, 260–265. [Google Scholar] [CrossRef]
  90. Ren, Z.; Wang, Y.; Liebenson, D.; Liggett, T.; Goswami, R.; Stefoski, D.; Balabanov, R. IRF-1 signaling in central nervous system glial cells regulates inflammatory demyelination. J. Neuroimmunol. 2011, 233, 147–159. [Google Scholar] [CrossRef]
  91. Ren, Z.; Wang, Y.; Tao, D.; Liebenson, D.; Liggett, T.; Goswami, R.; Clarke, R.; Stefoski, D.; Balabanov, R. Overexpression of the dominant-negative form of interferon regulatory factor 1 in oligodendrocytes protects against experimental autoimmune encephalomyelitis. J. Neurosci. 2011, 31, 8329–8341. [Google Scholar] [CrossRef] [Green Version]
  92. Itoh, T.; Horiuchi, M.; Wakayama, K.; Itoh, A.; Pleasure, D. Demyelination and oligodendrocytes. Rinsho Shinkeigaku 2011, 51, 1036–1039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Wang, Y.; Ren, Z.; Tao, D.; Tilwalli, S.; Goswami, R.; Balabanov, R. STAT1/IRF-1 signaling pathway mediates the injurious effect of interferon-gamma on oligodendrocyte progenitor cells. Glia 2010, 58, 195–208. [Google Scholar] [CrossRef] [PubMed]
  94. Loda, E.; Balabanov, R. Interferon regulatory factor 1 regulation of oligodendrocyte injury and inflammatory demyelination. Rev. Neurosci. 2012, 23, 145–152. [Google Scholar] [CrossRef]
  95. Sjöstrand, M.; Johansson, A.; Aqrawi, L.; Olsson, T.; Wahren-Herlenius, M.; Espinosa, A. The Expression of BAFF Is Controlled by IRF Transcription Factors. J. Immunol. 2016, 196, 91–96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Ma, N.; He, Y.; Xiao, H.; Han, G.; Chen, G.; Wang, Y.; Wang, K.; Hou, C.; Yang, X.; Marrero, B.; et al. BAFF maintains T-cell survival by inducing OPN expression in B cells. Mol. Immunol. 2014, 57, 129–137. [Google Scholar] [CrossRef] [PubMed]
  97. Liddelow, S.A.; Barres, B.A. Reactive Astrocytes: Production, Function, and Therapeutic Potential. Immunity 2017, 46, 957–967. [Google Scholar] [CrossRef] [Green Version]
  98. Brambilla, R.; Persaud, T.; Hu, X.; Karmally, S.; Shestopalov, V.I.; Dvoriantchikova, G.; Ivanov, D.; Nathanson, L.; Barnum, S.R.; Bethea, J.R. Transgenic inhibition of astroglial NF-kappa B improves functional outcome in experimental autoimmune encephalomyelitis by suppressing chronic central nervous system inflammation. J. Immunol. 2009, 182, 2628–2640. [Google Scholar] [CrossRef]
  99. Jing, Y.; Ma, R.; Chu, Y.; Dou, M.; Wang, M.; Li, X.; Zhu, L. Matrine treatment induced an A2 astrocyte phenotype and protected the blood-brain barrier in CNS autoimmunity. J. Chem. Neuroanat. 2021, 117, 102004. [Google Scholar] [CrossRef]
  100. Une, H.; Yamasaki, R.; Nagata, S.; Yamaguchi, H.; Nakamuta, Y.; Indiasari, U.C.; Cui, Y.; Shinoda, K.; Masaki, K.; Götz, M.; et al. Brain gray matter astroglia-specific connexin 43 ablation attenuates spinal cord inflammatory demyelination. J. Neuroinflammation 2021, 18, 126. [Google Scholar] [CrossRef]
  101. Xu, X.; Zhang, A.; Zhu, Y.; He, W.; Di, W.; Fang, Y.; Shi, X. MFG-E8 reverses microglial-induced neurotoxic astrocyte (A1) via NF-κB and PI3K-Akt pathways. J. Cell. Physiol. 2018, 234, 904–914. [Google Scholar] [CrossRef] [Green Version]
  102. Guo, H.; Fan, Z.; Wang, S.; Ma, L.; Wang, J.; Yu, D.; Zhang, Z.; Wu, L.; Peng, Z.; Liu, W.; et al. Astrocytic A1/A2 paradigm participates in glycogen mobilization mediated neuroprotection on reperfusion injury after ischemic stroke. J. Neuroinflammation 2021, 18, 230. [Google Scholar] [CrossRef] [PubMed]
  103. Guo, M.F.; Zhang, H.Y.; Li, Y.H.; Gu, Q.F.; Wei, W.Y.; Wang, Y.Y.; Zhang, X.J.; Liu, X.Q.; Song, L.J.; Chai, Z.; et al. Fasudil inhibits the activation of microglia and astrocytes of transgenic Alzheimer’s disease mice via the downregulation of TLR4/Myd88/NF-κB pathway. J. Neuroimmunol. 2020, 346, 577284. [Google Scholar] [CrossRef] [PubMed]
  104. Wang, J.; Liu, B.; Xu, Y.; Luan, H.; Wang, C.; Yang, M.; Zhao, R.; Song, M.; Liu, J.; Sun, L.; et al. Thioperamide attenuates neuroinflammation and cognitive impairments in Alzheimer’s disease via inhibiting gliosis. Exp. Neurol. 2021, 347, 113870. [Google Scholar] [CrossRef]
  105. King, A.; Szekely, B.; Calapkulu, E.; Ali, H.; Rios, F.; Jones, S.; Troakes, C. The Increased Densities, But Different Distributions, of Both C3 and S100A10 Immunopositive Astrocyte-Like Cells in Alzheimer’s Disease Brains Suggest Possible Roles for Both A1 and A2 Astrocytes in the Disease Pathogenesis. Brain Sci. 2020, 10, 503. [Google Scholar] [CrossRef] [PubMed]
  106. Bolstad, B.M.; Collin, F.; Brettschneider, J.; Simpson, K.; Cope, L.; Irizarry, R.A.; Speed, T.P. Quality assessment of Affymetrix GeneChip Data. In Bioinformatics and Computational Biology Solutions Using R and Bioconductor; Gentleman, R., Carey, V., Huber, W., Irizarry, R.A., Dudoit, S., Eds.; Springer: New York, NY, USA, 2005. [Google Scholar]
  107. Irizarry, R.A.; Bolstad, B.M.; Collin, F.; Cope, L.M.; Hobbs, B.; Speed, T.P. Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31, e15. [Google Scholar] [CrossRef] [PubMed]
  108. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: A powerful and practical approach to multiple testing. J. R. Stat. Soc. B 1995, 57, 289–300. [Google Scholar] [CrossRef]
  109. Swindell, W.R.; Johnston, A.; Xing, X.; Voorhees, J.J.; Elder, J.T.; Gudjonsson, J.E. Modulation of epidermal transcription circuits in psoriasis: New links between inflammation and hyperproliferation. PLoS ONE 2013, 8, e79253. [Google Scholar] [CrossRef]
  110. Russ, K.; Teku, G.; Bousset, L.; Redeker, V.; Piel, S.; Savchenko, E.; Pomeshchik, Y.; Savistchenko, J.; Stummann, T.C.; Azevedo, C.; et al. TNF-α and α-synuclein fibrils differently regulate human astrocyte immune reactivity and impair mitochondrial respiration. Cell Rep. 2021, 34, 108895. [Google Scholar] [CrossRef]
  111. Langfelder, P.; Zhang, B.; Horvath, S. Defining clusters from a hierarchical cluster tree: The Dynamic Tree Cut package for R. Bioinformatics 2008, 24, 719–720. [Google Scholar] [CrossRef]
  112. Blake, J.A.; Dolan, M.; Drabkin, H.; Hill, D.P.; Li, N.; Sitnikov, D.; Bridges, S.; Burgess, S.; Buza, T.; McCarthy, F.; et al. Gene Ontology annotations and resources. Nucleic Acids Res. 2013, 41, D530–D535. [Google Scholar] [CrossRef] [Green Version]
  113. Falcon, S.; Gentleman, R. Using GOstats to test gene lists for GO term association. Bioinformatics 2007, 23, 257–258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  114. Swindell, W.R.; Sarkar, M.K.; Stuart, P.E.; Voorhees, J.J.; Elder, J.T.; Johnston, A.; Gudjonsson, J.E. Psoriasis drug development and GWAS interpretation through in silico analysis of transcription factor binding sites. Clin. Transl. Med. 2015, 4, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  115. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  116. Ahmadi-Beni, R.; Najafi, A.; Savar, S.M.; Mohebbi, N.; Khoshnevisan, A. Role of dimethyl fumarate in the treatment of glioblastoma multiforme: A review article. Iran. J. Neurol. 2019, 18, 127–133. [Google Scholar] [CrossRef]
  117. Saidu, N.E.B.; Kavian, N.; Leroy, K.; Jacob, C.; Nicco, C.; Batteux, F.; Alexandre, J. Dimethyl fumarate, a two-edged drug: Current status and future directions. Med. Res. Rev. 2019, 39, 1923–1952. [Google Scholar] [CrossRef]
Figure 1. Differential expression comparison (MMF vs. DRF vs. IDMF). (AC) Fumarate compound molecular structures. (D) PC differential expression vectors. Arrows start at the CTL treatment (bivariate mean) and terminate at the bivariate mean for each treatment (MMF, DRF, or IDMF). Longer arrows correspond to a stronger treatment effect. (E) Heatmap. FC change estimates are shown for the top-ranked 50% of genes (n = 6493) having the largest absolute FC estimate (|FC|) among the three differential expression comparison. (F) Scatterplots. FC estimates are plotted to compare effects of MMF, DRF, and IDMF. The density of genes in each region is indicated (see color scale) and the Spearman correlation coefficient and p-value are shown (top margin). (G) Self-organizing maps (SOMs). Genes were assigned to regions within an SOM and the average fold-change of genes within each region is shown (see color scale). (H) SOM surface plots. The three-dimensional surface indicates the average FC among genes assigned to each SOM region. Plots are shown with varying rotations for viewing at multiple angles (45, 90, 135, 180, and 225 degrees).
Figure 1. Differential expression comparison (MMF vs. DRF vs. IDMF). (AC) Fumarate compound molecular structures. (D) PC differential expression vectors. Arrows start at the CTL treatment (bivariate mean) and terminate at the bivariate mean for each treatment (MMF, DRF, or IDMF). Longer arrows correspond to a stronger treatment effect. (E) Heatmap. FC change estimates are shown for the top-ranked 50% of genes (n = 6493) having the largest absolute FC estimate (|FC|) among the three differential expression comparison. (F) Scatterplots. FC estimates are plotted to compare effects of MMF, DRF, and IDMF. The density of genes in each region is indicated (see color scale) and the Spearman correlation coefficient and p-value are shown (top margin). (G) Self-organizing maps (SOMs). Genes were assigned to regions within an SOM and the average fold-change of genes within each region is shown (see color scale). (H) SOM surface plots. The three-dimensional surface indicates the average FC among genes assigned to each SOM region. Plots are shown with varying rotations for viewing at multiple angles (45, 90, 135, 180, and 225 degrees).
Pharmaceuticals 15 00461 g001
Figure 2. Top-ranked DEGs. (AC) Top 30 DEGs most strongly altered by (A) MMF, (B) DRF, and (C) IDMF. The top-ranked 30 increased and decreased genes are shown for each compound (i.e., lowest p-values with FC > 1.25 or FC < 0.80). (D) NFATC2. (E) SACM1L. (F) ACTR3. In (DF), average expression (log2 scale) is shown for each treatment (±1 standard error; * p < 0.05, comparison to CTL treatment). Expression is normalized to the average value in the CTL treatment for each gene. (GL) GO BP terms enriched among (G,I,K) increased and (H,J,L) decreased DEGs. The number of DEGs associated with each GO term is given in parentheses (left margin) and example DEGs are listed within each figure.
Figure 2. Top-ranked DEGs. (AC) Top 30 DEGs most strongly altered by (A) MMF, (B) DRF, and (C) IDMF. The top-ranked 30 increased and decreased genes are shown for each compound (i.e., lowest p-values with FC > 1.25 or FC < 0.80). (D) NFATC2. (E) SACM1L. (F) ACTR3. In (DF), average expression (log2 scale) is shown for each treatment (±1 standard error; * p < 0.05, comparison to CTL treatment). Expression is normalized to the average value in the CTL treatment for each gene. (GL) GO BP terms enriched among (G,I,K) increased and (H,J,L) decreased DEGs. The number of DEGs associated with each GO term is given in parentheses (left margin) and example DEGs are listed within each figure.
Pharmaceuticals 15 00461 g002
Figure 3. Gene expression module analysis. (A) Module co-expression networks. Each vertex represents one of 234 modules. Connections are drawn between modules for which medoids are correlated (rs ≥ 0.80). Vertex color reflects the average FC of module genes. Edge color reflects the average FC of the two modules joined. (B) Module cluster analysis. The 234 module medoids were clustered based upon the Euclidean distance (normalized to the [0, 1] interval). Average FC for each module is shown and module IDs are listed (right margin). (C) Number of differentially expressed modules (DEMs) for each comparison (MMF vs. CTL, DRF cv. CTL, IDMF vs. CTL). (D) IDMF-responsive modules (cluster analysis). The 10 modules most strongly altered by IDMF are shown. Cluster analysis was performed using Euclidean distance (normalized to the [0, 1] interval). Boxplots outline the middle 50% of FC estimates for each module (whiskers: 10th to 90th percentiles). (E) IDMF-responsive modules (network). The 10 modules from (D) are shown (vertex labels are abbreviated). Network parameters and color-coding are as described for part (A) above. (F) Module DCTN1-74 genes. The 16 genes most strongly increased by IDMF are shown (i.e., lowest p-values). (G) GO BP terms enriched among module DCTN1-74 genes. (H) GO CC terms enriched among module DCTN1-74 genes. (I) Module SLC8A1-109 genes. The 16 genes most strongly decreased by IDMF are shown (i.e., lowest p-values). (J) GO BP terms enriched among module SLC8A1-109 genes. (K) GO CC terms enriched among module SLC8A1-109 genes. In (G,H,J,K), the number of genes associated with each term is listed in parentheses (left margin) and example genes are listed within the figure.
Figure 3. Gene expression module analysis. (A) Module co-expression networks. Each vertex represents one of 234 modules. Connections are drawn between modules for which medoids are correlated (rs ≥ 0.80). Vertex color reflects the average FC of module genes. Edge color reflects the average FC of the two modules joined. (B) Module cluster analysis. The 234 module medoids were clustered based upon the Euclidean distance (normalized to the [0, 1] interval). Average FC for each module is shown and module IDs are listed (right margin). (C) Number of differentially expressed modules (DEMs) for each comparison (MMF vs. CTL, DRF cv. CTL, IDMF vs. CTL). (D) IDMF-responsive modules (cluster analysis). The 10 modules most strongly altered by IDMF are shown. Cluster analysis was performed using Euclidean distance (normalized to the [0, 1] interval). Boxplots outline the middle 50% of FC estimates for each module (whiskers: 10th to 90th percentiles). (E) IDMF-responsive modules (network). The 10 modules from (D) are shown (vertex labels are abbreviated). Network parameters and color-coding are as described for part (A) above. (F) Module DCTN1-74 genes. The 16 genes most strongly increased by IDMF are shown (i.e., lowest p-values). (G) GO BP terms enriched among module DCTN1-74 genes. (H) GO CC terms enriched among module DCTN1-74 genes. (I) Module SLC8A1-109 genes. The 16 genes most strongly decreased by IDMF are shown (i.e., lowest p-values). (J) GO BP terms enriched among module SLC8A1-109 genes. (K) GO CC terms enriched among module SLC8A1-109 genes. In (G,H,J,K), the number of genes associated with each term is listed in parentheses (left margin) and example genes are listed within the figure.
Pharmaceuticals 15 00461 g003
Figure 4. NRF2 target genes. (A) NQO1. (B) HMOX1. (C) G6PD. (D) TXNRD1. In (AD), average expression (log2 scale) is shown for each treatment (± 1 standard error; * p < 0.05, comparison to CTL treatment). Expression is normalized to the average value in the CTL treatment for each gene. (E) NRF2-increased genes. (F) NRF2-decreased genes. In (E,F), the value of |FC|max was calculated for each gene, where |FC|max is defined as max [abs(FCMMF, FCDRF, FCIDMF)], and the 23 NRF2-increased or NRF2-decreased genes with highest value of |FC|max are shown. Non-black gene labels are used if there is significant differential expression (p < 0.05) for any of the three comparisons, in which case the label color matches the comparison associated with the lowest differential expression p-value. (G,H,K,L,O,P) FC estimates for NRF2-responsive genes. FC estimates are plotted and the proportion of fumarate-increased (red) fumarate-decreased (blue) genes is shown (p-value: Fisher’s exact test). The percentage of increased/decreased genes is also indicated (see legend). (I,J,M,N,Q,R) The cumulate overlap is shown between NRF2-regulated genes and genes ranked based upon their response to each fumarate compound. A positive area statistic denotes enrichment among fumarate-increased genes, and a negative area statistic indicates enrichment among fumarate-decreased genes (p-value: Wilcoxon rank sum test).
Figure 4. NRF2 target genes. (A) NQO1. (B) HMOX1. (C) G6PD. (D) TXNRD1. In (AD), average expression (log2 scale) is shown for each treatment (± 1 standard error; * p < 0.05, comparison to CTL treatment). Expression is normalized to the average value in the CTL treatment for each gene. (E) NRF2-increased genes. (F) NRF2-decreased genes. In (E,F), the value of |FC|max was calculated for each gene, where |FC|max is defined as max [abs(FCMMF, FCDRF, FCIDMF)], and the 23 NRF2-increased or NRF2-decreased genes with highest value of |FC|max are shown. Non-black gene labels are used if there is significant differential expression (p < 0.05) for any of the three comparisons, in which case the label color matches the comparison associated with the lowest differential expression p-value. (G,H,K,L,O,P) FC estimates for NRF2-responsive genes. FC estimates are plotted and the proportion of fumarate-increased (red) fumarate-decreased (blue) genes is shown (p-value: Fisher’s exact test). The percentage of increased/decreased genes is also indicated (see legend). (I,J,M,N,Q,R) The cumulate overlap is shown between NRF2-regulated genes and genes ranked based upon their response to each fumarate compound. A positive area statistic denotes enrichment among fumarate-increased genes, and a negative area statistic indicates enrichment among fumarate-decreased genes (p-value: Wilcoxon rank sum test).
Pharmaceuticals 15 00461 g004
Figure 5. NF-kB (RELA) target genes. (A) NFKBIA. (B) TNFAIP3. (C) BCL2L1. (D) ICAM1. In (AD), average expression (log2 scale) is shown for each treatment (±1 standard error; * p < 0.05, comparison to CTL treatment). Expression is normalized to the average value in the CTL treatment for each gene. (E) RELA-activated genes. (F) RELA-suppressed genes. In (E,F), the value of |FC|max was calculated for each gene, where |FC|max is defined as max [abs(FCMMF, FCDRF, FCIDMF)], and the 23 RELA-activated or RELA-suppressed genes with highest value of |FC|max are shown. Non-black gene labels are used if there is significant differential expression (p < 0.05) for any of the three comparisons, in which case the label color matches the comparison associated with the lowest differential expression p-value. (G,H,K,L,O,P) FC estimates for RELA-activated genes. FC estimates are plotted and the proportion of fumarate-increased (red) fumarate-decreased (blue) genes is shown (p-value: Fisher’s exact test). The percentage of increased/decreased genes is also indicated (see legend). (I,J,M,N,Q,R) The cumulate overlap is shown between RELA-activated genes and genes ranked based upon their response to each fumarate compound. A positive area statistic denotes enrichment among fumarate-increased genes, and a negative area statistic indicates enrichment among fumarate-decreased genes (p-value: Wilcoxon rank sum test).
Figure 5. NF-kB (RELA) target genes. (A) NFKBIA. (B) TNFAIP3. (C) BCL2L1. (D) ICAM1. In (AD), average expression (log2 scale) is shown for each treatment (±1 standard error; * p < 0.05, comparison to CTL treatment). Expression is normalized to the average value in the CTL treatment for each gene. (E) RELA-activated genes. (F) RELA-suppressed genes. In (E,F), the value of |FC|max was calculated for each gene, where |FC|max is defined as max [abs(FCMMF, FCDRF, FCIDMF)], and the 23 RELA-activated or RELA-suppressed genes with highest value of |FC|max are shown. Non-black gene labels are used if there is significant differential expression (p < 0.05) for any of the three comparisons, in which case the label color matches the comparison associated with the lowest differential expression p-value. (G,H,K,L,O,P) FC estimates for RELA-activated genes. FC estimates are plotted and the proportion of fumarate-increased (red) fumarate-decreased (blue) genes is shown (p-value: Fisher’s exact test). The percentage of increased/decreased genes is also indicated (see legend). (I,J,M,N,Q,R) The cumulate overlap is shown between RELA-activated genes and genes ranked based upon their response to each fumarate compound. A positive area statistic denotes enrichment among fumarate-increased genes, and a negative area statistic indicates enrichment among fumarate-decreased genes (p-value: Wilcoxon rank sum test).
Pharmaceuticals 15 00461 g005
Figure 6. MS-associated genes. (A) MS gene databases. (B) Top MS-associated genes (six or more database sources). (C) Number of sources for MS-associated genes. (D) MS-associated genes (four or more sources) and their response to MMF, DRF, and IDMF. The 40 genes most strongly altered by one of the three compounds are shown (i.e., lowest p-value). (E,G,I) Average absolute FC (|log2(FC)|) of MS-associated genes (four or more sources). The null distribution shown was obtained by calculating the absolute FC in randomly sampled gene sets of the same size (10,000 simulation trials). (F,H,J) MS-associated genes most strongly altered by each compound (four or more sources). The number of database sources linking each to MS is shown in parentheses. Significantly altered genes (p < 0.05) are shown in red (increased) or blue (decreased) font. (K,L) GO BP terms enriched among MS-associated genes (three or more sources) altered by IDMF (p < 0.05 with FC > 1.25 or FC < 0.80). The number of genes associated with each GO BP term is indicated in parentheses (left margin) and example genes are listed within each figure.
Figure 6. MS-associated genes. (A) MS gene databases. (B) Top MS-associated genes (six or more database sources). (C) Number of sources for MS-associated genes. (D) MS-associated genes (four or more sources) and their response to MMF, DRF, and IDMF. The 40 genes most strongly altered by one of the three compounds are shown (i.e., lowest p-value). (E,G,I) Average absolute FC (|log2(FC)|) of MS-associated genes (four or more sources). The null distribution shown was obtained by calculating the absolute FC in randomly sampled gene sets of the same size (10,000 simulation trials). (F,H,J) MS-associated genes most strongly altered by each compound (four or more sources). The number of database sources linking each to MS is shown in parentheses. Significantly altered genes (p < 0.05) are shown in red (increased) or blue (decreased) font. (K,L) GO BP terms enriched among MS-associated genes (three or more sources) altered by IDMF (p < 0.05 with FC > 1.25 or FC < 0.80). The number of genes associated with each GO BP term is indicated in parentheses (left margin) and example genes are listed within each figure.
Pharmaceuticals 15 00461 g006
Figure 7. MS-associated genes from GWAS studies and their overlap with DEGs. (A,E,I) Overlap between DEGs and MS-associated genes from the NHGRI-EBI GWAS Catalog (reported, mapped, upstream, or downstream). DEGs were identified based on the same threshold in each analysis (p < 0.05, FC > 1.25 or FC < 0.80; bottom p-value: Fisher’s exact test for overlap). (B,F,J) Genes at varying distances from MS GWAS loci (horizontal axis) and their overlap with the top 30 genes most strongly increased/decreased by MMF, DRF, or IDMF (black line: overlap with respect to all other genes). Cases of significant overlap are indicated in the top margin (p < 0.05, Fisher’s exact test). (C,D,G,H,K,L) Average distance between the top 30 genes most strongly altered by MMF, DRF, or IDMF and the nearest MS GWAS locus (arrow). The distribution is obtained from 1000 simulation trials in which 30 genes were selected at random, with the average distance to the nearest MS GWAS locus calculated in each trial. (M,N) IDMF-increased/decreased genes nearest to an MS GWAS locus. For each gene, the distance to the nearest MS GWAS locus is shown (left) along with estimated fold-change (IDMF/CTL) (right). (OR) Selected IDMF-responsive genes near MS GWAS loci. Average expression (±1 standard error) is shown for each treatment (* p < 0.05, comparison to CTL).
Figure 7. MS-associated genes from GWAS studies and their overlap with DEGs. (A,E,I) Overlap between DEGs and MS-associated genes from the NHGRI-EBI GWAS Catalog (reported, mapped, upstream, or downstream). DEGs were identified based on the same threshold in each analysis (p < 0.05, FC > 1.25 or FC < 0.80; bottom p-value: Fisher’s exact test for overlap). (B,F,J) Genes at varying distances from MS GWAS loci (horizontal axis) and their overlap with the top 30 genes most strongly increased/decreased by MMF, DRF, or IDMF (black line: overlap with respect to all other genes). Cases of significant overlap are indicated in the top margin (p < 0.05, Fisher’s exact test). (C,D,G,H,K,L) Average distance between the top 30 genes most strongly altered by MMF, DRF, or IDMF and the nearest MS GWAS locus (arrow). The distribution is obtained from 1000 simulation trials in which 30 genes were selected at random, with the average distance to the nearest MS GWAS locus calculated in each trial. (M,N) IDMF-increased/decreased genes nearest to an MS GWAS locus. For each gene, the distance to the nearest MS GWAS locus is shown (left) along with estimated fold-change (IDMF/CTL) (right). (OR) Selected IDMF-responsive genes near MS GWAS loci. Average expression (±1 standard error) is shown for each treatment (* p < 0.05, comparison to CTL).
Pharmaceuticals 15 00461 g007
Figure 8. RT-PCR analyses and NF-κB luciferase reporter assay. (AG) RT-PCR analyses (OSGIN1, IL6, ICAM1, MALT1, TNFAIP3, IRF1, and CXCL8). Average relative expression is shown for each gene (±1 standard error; n = 3 replicates per treatment). Relative expression was calculated using the 2−∆∆Ct method and ∆Ct values were further normalized to the CTL treatment. Average relative expression is shown for each treatment (bottom margin). Treatments with different letters have significantly different average expression (p < 0.05, Fisher’s least significant difference). (H) NF-κB reporter assay. Human fetal astrocytes were transfected with firefly/Renilla luciferase constructs and treated with TNF-α (20 ng/ml) alone (n = 6) or TNF-α (20 ng/mL) + IDMF (2.5 µM) (n = 3) for 6 h. The ratio of luciferase and Renilla (internal control) signals was used to monitor NF-κB induction. The average ratio value is shown (±1 standard error) for each treatment (p-value: two-sample Wilcoxon rank sum test).
Figure 8. RT-PCR analyses and NF-κB luciferase reporter assay. (AG) RT-PCR analyses (OSGIN1, IL6, ICAM1, MALT1, TNFAIP3, IRF1, and CXCL8). Average relative expression is shown for each gene (±1 standard error; n = 3 replicates per treatment). Relative expression was calculated using the 2−∆∆Ct method and ∆Ct values were further normalized to the CTL treatment. Average relative expression is shown for each treatment (bottom margin). Treatments with different letters have significantly different average expression (p < 0.05, Fisher’s least significant difference). (H) NF-κB reporter assay. Human fetal astrocytes were transfected with firefly/Renilla luciferase constructs and treated with TNF-α (20 ng/ml) alone (n = 6) or TNF-α (20 ng/mL) + IDMF (2.5 µM) (n = 3) for 6 h. The ratio of luciferase and Renilla (internal control) signals was used to monitor NF-κB induction. The average ratio value is shown (±1 standard error) for each treatment (p-value: two-sample Wilcoxon rank sum test).
Pharmaceuticals 15 00461 g008
Figure 9. IDMF hypothesized mechanisms of therapeutic effect. IDMF activates NRF2-dependent signaling (red arrows) while inhibiting NF-κB activity and transcription of downstream genes (blue). IDMF triggers dissociation of NRF2 from KEAP1 to allow nuclear translocation of NRF2. Nuclear NRF2 dimerizes with musculoaponeurotic fibrosarcoma (Maf) and binds antioxidant response elements (ARE) to induce the transcription of target genes, such as HMOX1. The HMOX1 protein then acts as an enzyme to catalyze the formation of antioxidants such as bilirubin and biliverdin, while also inhibiting nitric oxide synthase (NOS). Together, these effects of HMOX1 reduce accumulation of reactive oxygen and reactive nitrogen species (ROS and RNS) to attenuate oxidative stress. The inhibition of NF-κB by IDMF prevents nuclear translocation and up-regulation of target genes such as ICAM1, CXCL8, IL6, and IRF1. Reduced expression of these targets leads to reduced inflammation via attenuation of ICAM1-dependent leukocyte adhesion, cytokine production with cytokine storm, and IRF1/CASP4-dependent inflammasome activation. Overall dampening of these inflammatory processes inhibits pyroptosis and, thus, limits progression of demyelination. See text for details. (Figure created with BioRender.com).
Figure 9. IDMF hypothesized mechanisms of therapeutic effect. IDMF activates NRF2-dependent signaling (red arrows) while inhibiting NF-κB activity and transcription of downstream genes (blue). IDMF triggers dissociation of NRF2 from KEAP1 to allow nuclear translocation of NRF2. Nuclear NRF2 dimerizes with musculoaponeurotic fibrosarcoma (Maf) and binds antioxidant response elements (ARE) to induce the transcription of target genes, such as HMOX1. The HMOX1 protein then acts as an enzyme to catalyze the formation of antioxidants such as bilirubin and biliverdin, while also inhibiting nitric oxide synthase (NOS). Together, these effects of HMOX1 reduce accumulation of reactive oxygen and reactive nitrogen species (ROS and RNS) to attenuate oxidative stress. The inhibition of NF-κB by IDMF prevents nuclear translocation and up-regulation of target genes such as ICAM1, CXCL8, IL6, and IRF1. Reduced expression of these targets leads to reduced inflammation via attenuation of ICAM1-dependent leukocyte adhesion, cytokine production with cytokine storm, and IRF1/CASP4-dependent inflammasome activation. Overall dampening of these inflammatory processes inhibits pyroptosis and, thus, limits progression of demyelination. See text for details. (Figure created with BioRender.com).
Pharmaceuticals 15 00461 g009
Table 1. Summary table. The table lists the number of DEGs identified (p < 0.05 with FC > 1.25 or FC < 0.80), top biological process (BP) terms most strongly enriched among DEGs, top sequence motifs enriched in DEG 5 kb upstream regions, DEGs correspondingly regulated by NRF2 and niacin, DEGs inversely regulated by NF-κB, IRF factors associated with enriched motifs in DEG upstream regions, A1 astrocyte marker down-regulated DEGs, MS-associated DEGs, and DEGs near MS GWAS loci. Results are summarized for increased DEGs (▲), decreased DEGs (▼), and the combined set of increased and decreased DEGs (▲ + ▼). See footnotes for further details.
Table 1. Summary table. The table lists the number of DEGs identified (p < 0.05 with FC > 1.25 or FC < 0.80), top biological process (BP) terms most strongly enriched among DEGs, top sequence motifs enriched in DEG 5 kb upstream regions, DEGs correspondingly regulated by NRF2 and niacin, DEGs inversely regulated by NF-κB, IRF factors associated with enriched motifs in DEG upstream regions, A1 astrocyte marker down-regulated DEGs, MS-associated DEGs, and DEGs near MS GWAS loci. Results are summarized for increased DEGs (▲), decreased DEGs (▼), and the combined set of increased and decreased DEGs (▲ + ▼). See footnotes for further details.
DEG GroupMMFDRFIDMF
No. DEGs identified 1308181673
142198903
▲ + ▼4503791576
No. DEMs identified 2144071
114270
▲ + ▼2582141
Top GO BP term 3Oxidation–reductionCircadian behaviorSmell stimulus detection
Mast cell degranulationMacromolecule modificationCellular metabolism
Top motif 45-AAATT/AATTT-35-CTAGCA/TGCTAG-35-CATG/CATG-3
5-GGGCG/CGCCC-35-TAATT/AATTA-35-TAAT/ATTA-3
NRF2-increased DEGs 5NQO1, SLC3A2, POP4TXNRD1, GCLC, PTGR1HMOX1, GSR, SLC7A11, SLC6A6
NRF2-decreased DEGs 5NoneNoneLTBP1, HNRNPU, PCYOX1
NF-kB-increased DEGs 6KLHL18ZZZ3, RIOK1, RHOBTB3NFKBIA, TNFAIP3, IRF1, NEDD4L
NF-kB-decreased DEGs 6FOSB, SLC3A2NoneEGR1
Niacin-increased DEGs 7CCNL1NoneSLC22A14
Niacin-decreased DEGs 7NoneCNIH1, ALG8HIGD1A, DHX9, AKT2
IRF motif enrichment 8NoneIRF1, IRF6IRF1, IRF4, IRF6, IRF7
A1 astrocyte DEGs 9CD82, CYP27A1TXNIPTERF1, GOT1, BPGM, SORT1
MS-associated DEGs 10CD47, SLC1A3HLA-DRB5, CNR2, PPARAHMOX1, HLA-DRB5, BTNL2
MPZGLO1, IFNAR2TNFAIP3, IFNAR2, ZNF433
DEGs near MS GWAS loci 11POPDC3, COG6, VWA8ALDH1L1, ACTRT3, ELAC1ALDH1L1, MLANA, SLC2A4RG
PRR5L, PLEC, ZC3H4PHGDH, PFDN4, EPPK1ZBTB38, ZNF433, STAT3
1 DEG criteria: p < 0.05 with FC > 1.25 or FC < 0.80. 2 DEM criteria: FDR < 0.05 (two-sample Wilcoxon rank sum test) (Figure 3). 3 GO BP term most strongly enriched among DEGs (Figure 2). 4 Motif most strongly enriched in regions 5 KB upstream of DEGs (Figure S2). 5 NRF2-regulated genes correspondingly altered by each fumarate (Figure 4). 6 NF-κB-regulated genes inversely altered by each fumarate (Figure 5). 7 Niacin-regulated genes correspondingly altered by each fumarate (Figure S4). 8 IRF transcription factors associated with recognition sequences enriched in sequences upstream of DEGs (Figure S5). 9 A1 neurotoxic astrocyte marker genes down-regulated by each fumarate (Figure S6). 10 MS-associated DEGs associated with each fumarate (Figure 6). 11 DEGs overlapping with or nearest to MS-associated GWAS loci (Figure 7).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Swindell, W.R.; Bojanowski, K.; Chaudhuri, R.K. Transcriptomic Analysis of Fumarate Compounds Identifies Unique Effects of Isosorbide Di-(Methyl Fumarate) on NRF2, NF-kappaB and IRF1 Pathway Genes. Pharmaceuticals 2022, 15, 461. https://doi.org/10.3390/ph15040461

AMA Style

Swindell WR, Bojanowski K, Chaudhuri RK. Transcriptomic Analysis of Fumarate Compounds Identifies Unique Effects of Isosorbide Di-(Methyl Fumarate) on NRF2, NF-kappaB and IRF1 Pathway Genes. Pharmaceuticals. 2022; 15(4):461. https://doi.org/10.3390/ph15040461

Chicago/Turabian Style

Swindell, William R., Krzysztof Bojanowski, and Ratan K. Chaudhuri. 2022. "Transcriptomic Analysis of Fumarate Compounds Identifies Unique Effects of Isosorbide Di-(Methyl Fumarate) on NRF2, NF-kappaB and IRF1 Pathway Genes" Pharmaceuticals 15, no. 4: 461. https://doi.org/10.3390/ph15040461

APA Style

Swindell, W. R., Bojanowski, K., & Chaudhuri, R. K. (2022). Transcriptomic Analysis of Fumarate Compounds Identifies Unique Effects of Isosorbide Di-(Methyl Fumarate) on NRF2, NF-kappaB and IRF1 Pathway Genes. Pharmaceuticals, 15(4), 461. https://doi.org/10.3390/ph15040461

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop