Skip to main content

Profiling the microRNA signature of the peripheral sensory ganglia in experimental autoimmune encephalomyelitis (EAE)



Multiple sclerosis is an autoimmune disease with a distinct female bias, as well as a high prevalence of neuropathic pain in both sexes. The dorsal root ganglia (DRG) contain the primary sensory neurons that give rise to pain, and damage to these neurons may lead to neuropathic pain. Here, we investigate the sex differences of the DRG transcriptome in a mouse model of MS.


Next-generation sequencing was used to establish RNA and microRNA profiles from the DRG of mice with MOG35–55-induced EAE, a model of CNS inflammation that mimics aspects of MS. Differential expression and multiple meta-analytic approaches were used to compare expression profiles in immunized female and male mice. Differential expression of relevant genes and microRNAs were confirmed by qPCR.


Three thousand five hundred twenty genes and 29 microRNAs were differentially expressed in the DRG of female mice with MOG35–55-EAE, while only 189 genes and 3 microRNAs were differentially expressed in males with MOG35–55-EAE. Genes related to the immune system were uniquely regulated in immunized female mice. Direct comparison of sex within disease indicates significant differences in interferon and phagosomal pathways between the sexes. miR-21a-5p is the primary dysregulated microRNA in both sexes, with females having additional dysregulated microRNAs, including miR-122-5p.


This study provides evidence that females are uniquely affected by MOG35–55-EAE and that this difference may result from additional signaling not present in the male. The altered transcriptome of females correlates with other studies finding hyperactivity of pain-sensing neurons and suggests underlying sex-specific pathways for neuropathic pain.


Multiple sclerosis (MS) is an autoimmune disease of the central nervous system, characterized by demyelination, disturbances in neuronal function, and progressive neurodegeneration [1]. Experimental autoimmune encephalomyelitis (EAE) is a mouse model used to study the pathophysiology of MS, where an autoimmune response is artificially induced against the myelin of the central nervous system (CNS). The model produces immune-related demyelination of the CNS [2,3,4]. People with MS experience altered motor function but recent estimates report that 63% of people living with MS are additionally negatively impacted by disturbances to sensation such as chronic, neuropathic pain [5]. Neuropathic pain arises from damage to the neurons involved in transmitting the “pain” signal to the brain. This damage can be localized to any part comprising the pain system, from the pain-sensing periphery to the pain-perceiving central nervous system [6].

The dorsal root ganglia (DRG) residing in the peripheral nervous system (PNS) contain primary sensory neurons including pain-sensing neurons. The precise contribution of DRG neurons and the mechanisms mediating neuropathic pain in MS and EAE remain elusive. While transcriptome dysregulation has been studied in the CNS and immune cell populations in the context of EAE [7,8,9], less research has been done on the PNS, a potential source of neuropathic pain in the disease. Recently, our laboratory has described differences in the development of pain states between female and male mice in the EAE model. These differences include differences in immune and neurodegenerative histochemical and oxidative stress profiles [2, 4]. These sex differences extend beyond the CNS, affecting the DRG [3]. Knowing that sex differences exist in EAE and that there are functional alterations correlated to pain, we sought to describe the transcriptional profile of the DRG in EAE from female and male mice with the disease.

MicroRNAs (miRNAs) are ~ 22 nucleotide long, single-stranded RNA molecules that repress translation of RNA species with complementary sequences [10]. miRNAs have risen to the forefront of research into the regulation of the transcriptome in pathophysiological states. It is estimated that the majority of all transcriptionally related pathways are under miRNA regulation [11]. Increasing amounts of evidence suggest that “miRNA-signatures” in various tissues can be correlated with specific disease states. As miRNAs can function as one-to-many vectors by regulating the expression of multiple genes, we were also interested to see if a “miRNA-signature” existed at the level of the DRG in addition to a functional analysis of the “transcriptome signature.” This information could yield insights into the mechanisms of pain in EAE and MS more generally.

Here, we identify EAE-related miRNAs that are differentially expressed in both male and female DRGs. We compare the two sex-specific signatures at the transcriptional level for mice with EAE. We describe the analysis of RNASeq and miRSeq datasets, and functional clustering of the differentially expressed genes and miRNAs into disease-relevant signatures. Our findings highlight extensive sex differences that must be considered with animal models of disease and identify potential targets for pain-modifying therapies in the disease.


MOG35–55-EAE generation

Female and male 6–8-week-old C57BL/6 mice were acquired from Charles River Laboratories, Canada, and habituated for 2 weeks in the housing facility with baseline handling and behavioral testing. After this period, MOG35–55-EAE, (hereon referred to as “EAE”) was induced by 50 μg subcutaneous injections of myelin oligodendrocyte glycoprotein (sp) (MOG35–55; Peptide Synthesis Facility, University of Calgary) emulsified in complete Freund’s adjuvant (CFA) at a concentration of 1.5 mg/ml. An additional intraperitoneal injection of 300 ng of pertussis toxin in 0.2 ml saline was given on the day of induction as well as 48 h later. Control (termed CFA hereon after) mice received identical injection protocols except for the lack of emulsified MOG. CFA is known to produce transient alterations in immune activity [12] but is still an effective tool for EAE induction and an appropriate control. Mice were assessed daily for clinical progression of EAE following a four-point scale: grade 0, normal mouse; grade 1, flaccid tail (disease onset); grade 2, mild hindlimb weakness with quick righting reflex; grade 3, severe hindlimb weakness with slow righting reflex; grade 4, hindlimb paralysis in one hindlimb or both. Tissue was extracted on the first day a mouse presented with clinical signs (see below) that is denoted as “onset.”

Tissue collection and total RNA extraction

Animals were euthanized by Euthanyl® (sodium pentobarbital) overdose injected intrapertioneally. After injection, animals were monitored for level of consciousness and dissections did not proceed until no response to toe pinch or corneal contact was observed. Cardiac punctures were performed to confirm euthanization and animals were perfused with 10 mL of saline. Whole DRGs were immediately extracted and suspended in Qiazol (Qiagen Biosystems) and stored at − 80 °C for later processing. Upon collection of the full tissue set, total RNA was extracted using miRNeasy® Mini kits (Qiagen Biosystems).

Next-generation sequencing

Total RNA was supplied to ArrayStar Inc. for sequencing. Total RNA was enriched by oligoDT magnetic beads, and library preparation was completed using KAPA Stranded RNA-Seq Library Prep Kit. Transcriptome data was acquired by sequencing on an Illumina HiSeq 4000 machine according to manufacturer’s instructions. MicroRNA data was acquired by 3′-adapter ligation with T4 RNA ligase 2 (truncated), 5′-adapter ligation with T4 RNA ligase, cDNA synthesis with RT primer, PCR amplification, and extraction and purification of ~ 130–150 bp PCR amplified fragments and sequencing on an Illumina HiSeq 2000 machine according to manufacturer’s instructions.

Analysis of differential expression

Both RNASeq and miRSeq differential expression was investigated through the R [13] package DESeq2. Briefly, raw transcript sequences from NGS were aligned to the UCSC “mm10” database of known genes using the Genomic Features package function “exonsBy.” miR read counts were directly uploaded into R. The resulting data structures were processed through the DESeq2 package following the author’s instructions found in their published vignette [14,15,16]. The results of each contrasting analysis (e.g., M EAE vs. M CFA) were published in a tabular format for usage in subsequent analyses. Genes were considered significantly differentially expressed using the threshold of q <  0.1 after Benjamini-Hochberg correction. Clustergram and PCA analyses were performed using MATLAB (MATLAB and Statistics Toolbox Release 2017b, The MathWorks, Inc., Natick, MA, USA).

Ontological analysis with GOrilla/REViGo

Single gene lists ranked by p-adjusted values were analyzed with the GOrilla [17] web software, using a p value threshold of 10E−9. Significant GO terms from each category (process, function, component) were then processed through REViGO [18] to screen for redundant terminology.

Functional analysis with ingenuity pathway analysis

Activation patterns of canonical pathways were analyzed through the use of IPA (QIAGEN Inc., [19]. Briefly, gene lists ranked by p-adjusted values were analyzed using the “Core Analysis” feature using default parameters and a p-adjusted cutoff of 0.1.

qPCR validation

For validation of DEmiRs, 50 ng of total RNA was reverse transcribed using the miScript II RT Kit (Qiagen Biosystems). miRs were detected with the Mm_miR-21_2 and Mm_miR-122a_1 miScript Primer Assays (Qiagen Biosystems) using the miScript SYBR Green PCR Kit (Qiagen Biosystems). Values detected were standardized to sno72. For validation of DEGs, 100 ng of total RNA was processed with DNase I (Invitrogen) followed by reverse transcription using Superscript III (Invitrogen). Transcripts were detected using the following primers: Ppia: RT2 qPCR Primer Assay for Mouse Ppia (Qiagen Biosystems), C3: Forward, GAGGCACATTGTCGGTGGTG, Reverse, CCAGGATGGACATAGTGGCG, C5ar1: Forward, ATGGACCCCATAGATAACAGCA, Reverse, GAGTAGATGATAAGGGCTGCAAC, Stat1: RT2 qPCR Primer Assay for Mouse Stat1 (Qiagen Biosystems) RT2 qPCR Primer Assay for Mouse Tlr7 (Qiagen Biosystems), RT2 qPCR Primer Assay for Mouse Tlr8 (Qiagen Biosystems), and RT2 qPCR Primer Assay for Mouse Trem2 (Qiagen Biosystems). Transcripts detected were standardized to Ppia values. All qPCR reactions were 20 uL volumes utilizing a StepOnePlus machine (Applied Biosystems).

qPCR statistical analysis

Statistical analysis was performed using Graphpad Prism 6.01. qPCR analysis utilized two-way ANOVA testing with Sidak’s correction for multiple testing. qPCR statistics were performed on untransformed values but plotted as linearized values for ease of visualization. Alpha was set to 0.05 for all testing.



To detect a transcriptional signature within the DRG of EAE mice, we performed next-generation sequencing (NGS) of DRGs taken from mice at the onset of the disease when pain hypersensitivity is evident [2, 3, 20]. Roughly 20,000 genes were detected by sequencing. We observed gross differential expression of genes (DEGs) in mice subjected to EAE both in female (3520 DEGs) and male (189 DEGs) mice relative to their sex-specific CFA controls (Fig. 1a, Additional file 1: Table S1 and Additional file 2: Table S2) (Wald test; q <  0.1). To visualize variation in animal profiles and expression patterns, we performed a principal component analysis of all animals by all genes. This unbiased approach projects numerous complex variables into simpler dimensions, allowing for appreciation of relatedness between samples based on their proximity. The first two components captured almost half of the total variation within the entire sequencing set, with 29 and 19%, respectively (Fig. 1b). Cluster analysis (MATLAB clustergram) of all animals by the sum of female and male DEGs showed high specificity for EAE over CFA control, but minimal specificity for sex (Fig. 1c).

Fig. 1
figure 1

EAE primarily affects female DRGs over males at the transcriptional level. a Joined volcano plots of DEGs; female DEGs = 3520, male DEGs = 189. Vertical dotted lines represent log10(q = 0.1) as a cutoff threshold for significance. b Principal component analysis suggests separation of the EAE and CFA DRGs when analyzed by all mRNA sequencing data. Top two components are plotted. c Heatmaps of all DEGs in EAE further visualize the clustering of EAE. Values are Reads/Totals Reads normalized to the average of the CFA condition. DEG = differentially expressed gene

Detection of EAE-related sex-specific signatures

To contextualize the significant genes into functional groups from known ontologies, we enlisted the DAVID [21], REACTOME [22], GOrilla [17], REViGO [18], and PANTHER [23,24,25] algorithms. Analysis of DEGs generated by DESeq analysis for significantly enriched GOslim (Gene Ontology) terms for Cell Compartment (Fig. 2a) showed a similar pattern of gene regulation for females and males. The Cell Compartment analysis references genes known to localize in specific regions of the cell. While many terms were identical in enrichment (i.e., “membrane”), further analysis of GO terms revealed some unique variations in Cell Compartment, as males included enrichment in mitochondrial terms whereas females were enriched in vesicular related terms (Tables 1 and 2). Remarkably, analysis of DEGs by REACTOME pathways to cluster molecular events revealed a striking sex-specific signature (Fig. 2b). The majority of pathways enriched in females such as immune system-related pathways were completely absent in males. Male EAE mice contained significant enrichment of a single unique term: “SLC-mediated transport,” with all others sharing overlap with females.

Fig. 2
figure 2

Female and male DRGs exhibit unique functional RNA signatures. a Both sexes show similar PANTHER Pathways GO cell compartment terms. b However, REACTOME pathways of female and male DEGs show a distinct profile, suggesting that the function of each set of DEGs may be different. Color legends represent grouped terms under the REACTOME hierarchy. Wedges are scaled to # DEGs in pathway/# genes defined in the respective pathway

Table 1 GOrilla/REViGO analysis of differentially expressed genes in female EAE. Top 10 non-dispensable terms are reported
Table 2 GOrilla/REViGO analysis of differentially expressed genes in male EAE. Top 10 non-dispensable terms are reported

Validation of RNASeq

To confirm the differential expression of transcripts indicated by RNASeq analysis and of special interest through subsequent meta-analyses, we performed qPCR on relevant target genes of interest (Fig. 3). Validation targets were chosen based on criteria of differential expression between EAE and CFA condition, differential expression between female and male EAE conditions, and inclusion in significant canonical pathways of Ingenuity Pathway Analysis. (Additional file 1: Table S1, Additional file 2: Table S2, Additional file 5: Table S5 and Additional file 6: Table S6). In females, we confirmed the upregulation of C3, C5ar1, Stat1, Trem2, Tlr7, and Tlr8, in concordance with RNASeq predications. Interestingly, we also detected significant downregulation of C3, Stat1, Trem2, and Tlr8 in male EAE compared to their control (Table 3).

Fig. 3
figure 3

Validation of DEGs by qPCR. Bar graphs of a C3, b C5ar1, c Stat1, d Trem2, e Tlr7, and f Tlr8. *p <  0.05, two-way ANOVAs with Sidak’s multiple comparison test. Bars indicate geometric mean with 95% confidence interval

Table 3 qPCR statistics for DEGs and DEmiRs


In addition to sequencing for gene coding fragments, we acquired NGS data for almost 2000 microRNAs (miRNAs). Differential expression was determined using an identical methodology to the RNASeq data above. We detected 29 differentially expressed miRNAs (DEmiRs) in the DRG of EAE females relative to CFA female controls, but only 3 DEmiRs in EAE males (Fig. 4a, Additional file 3: Table S3 and Additional file 4: Table S4) (Wald test; q <  0.1). All DEmiRs detected in males (miR-21a-5p, miR-21c, and miR-142a-5p) were similarly detected in females and displayed similar patterns of upregulation. Multidimensional variation was visualized by PCA analysis of all expressed miRNAs (Fig. 4b), and cluster analysis of DEmiR expression values was carried out in a similar manner to the RNAseq data set (Fig. 4c). We did not, however, identify a separation by disease by miRSeq variation unlike the RNASeq suggesting that any disease-related alterations may be more subtle.

Fig. 4
figure 4

miR-21a-5p is the candidate miR for EAE in females and males. a Joined volcano plots of DEmiRs; female DEmiRs = 29, male DEmiRs = 3. Vertical dotted lines represent log10(q = 0.1) as a cutoff threshold for significance. b Principal component analysis shows minimal separation of EAE and CFA animals by all miR sequencing data. Top two components are plotted. c Heatmaps of all DEmiRs in EAE. Values are Reads / Totals Reads normalized to the average of the CFA condition. d qPCR validation of miR-21a-5p in female and male DRGs. e qPCR validation of miR-122-5p in female and male DRGs. *p <  0.05, two-way ANOVAs with Sidak’s multiple comparison test. Bars indicate geometric mean with 95% confidence interval. DEmiR = differentially expressed miRs

Differences in transcriptional signatures by sex within EAE

We next wanted to directly compare females and males with EAE to determine the effect of sex on the transcriptional profiles in the DRG. To accomplish this, we first compared the groups of “Female EAE” and “Male EAE” to find ~ 6000 DEGs. To account for genes differentially expressed as a result of sex (such as Y-linked genes), we removed DEGs from an “all Female” vs. “all Male” comparison, leaving 869 DEGs (Additional file 5: Table S5). We were interested in the functional significance of the DEGs and relied upon Ingenuity Pathway Analysis to investigate hypothetical activation of known pathways. Direct comparison yielded terminology related to metabolism (i.e., LXR/RXR activation), phagocytosis (i.e., acute phase response), and interferon signaling (Table 4 and Additional file 6: Table S6). Negative log(p values) correspond to increasing levels of significance (−log(≤ 0.05) ≥ 1.3), ratios represent the proportion of DEGs/# of genes in canonical pathway, and z-scores represent the predicted activation state of the pathway, with positive and negative being activated and inhibited, respectively. The heightened activation state of these pathways suggests an increased amount of immune-related inflammation in the DRGs of females compared to males with EAE.

Table 4 Canonical pathways differentially expressed between female vs. male EAE

Validation of miRSeq

We next wanted to validate the differential expression of miRNAs in both females and males with EAE. We identified miR-21a-5p as the most significantly dysregulated miRNA in both female and male EAE relative to respective CFA control, as well as miR-122-5p as being significantly dysregulated in female EAE alone. We validated that miR-21a-5p (hereon referred to as miR-21) was increased in both female and male EAE mice relative to CFA controls by 3.78 and 4.76-fold, respectively (Fig. 5a, b). There was no significant difference in the level of miR-21 upregulation between the sexes (2-way ANOVA, disease main effect F1,34 = 79.21, p <  0.0001, followed by Sidak’s multiple comparison test for disease within females, t = 6.794, q <  0.0001, and males, t = 6.021, q <  0.0001). miR-122-5p (hereon referred to as miR-122) was confirmed to be significantly upregulated in females only relative to CFA female controls by 4.36-fold, with a significant difference between the sexes (2-way ANOVA, sex main effect F1,32 = 9.803, p = 0.0037, followed by Sidak’s multiple comparison test for disease within females, t = 2.675, q = 0.0232, and males, t = 1.940, q = 0.1188) (Fig. 5c, d).

Fig. 5
figure 5

Schematic of a hypothetical DRG in EAE. Although we observe increased miR-21 in both sexes, only females display the increased inflammatory signaling (see reference)


EAE is a complex disease that primarily targets the CNS. Less is known about how the PNS is affected in EAE. We have shown that a disease-related transcriptional signature is detectable in the DRG of EAE mice and that sex is an important factor in the characterization of this signature. To accomplish this, we utilized next-generation sequencing for mRNA and miRNAs from the DRG at the clinical onset of EAE and profiled the differentially expressed genes and miRNAs to acquire a functional signature for each sex. The identification of functional signatures for complex disease states or comorbidities such as pain is relevant for prompting strategic experimental design and stratifying future therapies.

Detecting the signature

We first wanted to describe the transcriptome of our animals, considering sex as a factor. Sex differences have been previously described in immunity [26,27,28] and inflammatory responses [29, 30], and both of these processes play a key role in our disease model. Indeed, we found extensive differential expression of the transcriptome in female EAE mice but only small amounts of differentially expressed genes (DEGs) in male EAE mice. Large changes in gene expression have been reported in the DRG after peripheral nerve injury, and these changes have been associated with neuropathic pain behaviors [31, 32]. This supports previous evidence from our laboratory, indicating that female but not male DRG neurons are hyper-responsive at the onset of EAE [3]. In addition to these differences in the quantity of DEGs, we also found major differences in the functional signature of this altered transcriptional profile. Our data indicates that there are extensive immune and cell-signaling-related processes disturbed in the female DRGs, a potential by-product of immune cells transiently localized there [33].

Sex-specific disease effects

Comparing the sexes directly, we noted the appearance of phagosomal activation pathways, complement, and interferon signaling as being differentially affected in female EAE versus their CFA controls (Table 3). The phagosome-related terms contained references to DEGs such as Trem2, Csf1r, Tlr7, and most of the complement system. Taken together, this data indicates sex-specific dysfunction of inflammatory activity in phagocytic cells. Laffont et al. have reported increased interferon responses of estrogen activated human plasma dendritic cells to TLR7 ligands, suggesting that this mechanism may be sex-hormone linked [34]. Tlr7 is an X-linked gene and this may account for the sex-bias we observed, as X chromosome genes may escape silencing.

Interferon signaling is an approved target for MS therapy, as beta-interferons have been utilized in relapsing-remitting MS since the late 1990s [35]. Although the effects of interferon inhibitors are modest, all clinical trials conducted to date have been with mixed female and male populations. The possibility of female specific interferon dysregulation raises the possibility of higher efficacy in a subset of the clinical population. We validated the female-EAE-specific increased expression of Stat1—a downstream effector of types I and II interferons—suggesting that interferon signaling is activated in these animals. In male EAE mice however, Stat1 levels were significantly reduced suggestive of a silencing of interferon signaling. As beta-interferon therapy would activate these pathways, it follows that their effect may have differential effects based on sex.

Activation of interferon signaling through Stat1 typically results in increased inflammation and enhanced immune activity [36]. While these pathways do seem to be activated in the diseased group of females, it is apparent that it is not the case for the males. The high activation of interferon-signaling and Stat1 in females is complemented by a similar increase in transcripts associated with inflammatory macrophages, including most members of the complement system. Complement proteins are involved with the recruitment and phagocytic capacity of innate immune cells [37, 38]. Based on our transcriptomic data, females may be poised to specifically respond to interferon stimulation, resulting in heightened phagocytic activity in the DRG in response to EAE.

Accounting for the signature

To account for the transcriptomic changes observed in EAE, we also sequenced the miR-transcriptome. We observed a small amount of overlap between the detected DEmiRs in each sex, corresponding to miRNAs previously linked to EAE and the MS literature: miRs-21 [39, 40] and − 142 [41]. The most significant miRNA for both sexes, miR-21, was both highly expressed and confirmed to be significantly upregulated in both female and male EAE DRGs. MS brain lesions have been shown to be enriched with miR-21 in patients with MS but not neuromyelitis optica (NMO) [42]. NMO is an autoimmune disease with similar symptomatology, potentially indicating that this signature is specific to MS/EAE pathology. Therefore, miR-21 may be suitable as a marker of EAE. miR-21 has been well studied in the context of cancer, giving rise to its label of an “oncomiR” [43,44,45]. Interest in miR-21 as a regenerative marker is also emerging with evidence that it is pro-regenerative [46] and immunomodulatory [47]. Macrophages uptake miR-21 laden exosomes after peripheral nerve injuries [48] where it can act as a direct ligand for Toll-like receptors 7 and 8 [49, 50]. Interestingly, miR-21 is affected by sex hormones, containing a direct response element for androgens [51] and directly affected by estradiol signaling [52, 53]. As miR-21 was similarly upregulated in both sexes, we additionally screened for miRNAs unique to sex. miR-122 was only upregulated in the diseased female DRG. miR-122 is primarily described for its expression in the liver; however, it is known to increase in the circulation after injury [54, 55] and may interact with type I interferon signaling [56]. Together with the transcriptomic signatures, our data indicates that the female DRGs may be uniquely affected by the epigenetic regulation of miRs-21 and -122. Although EAE is canonically not thought of as a PNS disease, the presence of the typical “MS-miR-signature” indicates that female DRGs may be affected concordantly with CNS lesions.


To our knowledge, this is the first presentation of transcriptome and miR-transcriptome dysregulation in DRGs using the EAE mouse model. We have shown that the miR-signature is primarily defined by miR-21, bearing similarity to previously described findings in the spinal cord and serum of mice with EAE, and MS lesions. We present evidence of a female-biased dysregulation of both “omic” analyses, the presence of distinct functional transcriptomic signatures for each sex, and that interferon signaling and phagosome function may be indicated for future analysis. Acknowledging fundamental differences between the sexes is not only important for experimental design, but for the development of potential sex-specific therapies. Our study provides a platform to study sex-specific changes in the DRG of EAE and the role that miRNAs may have on their transcriptome.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.



Complete Freund’s adjuvant


Central nervous system


Differential expression of genes


Dorsal root ganglia


Experimental autoimmune encephalomyelitis


Gene Ontology




MicroRNA sequencing


Myelin oligodendrocyte glyocoproteinn


Multiple sclerosis


Next-generation sequencing


Peripheral nervous system


RNA sequencing


  1. Lassmann H, Brück W, Lucchinetti CF. The immunopathology of multiple sclerosis: an overview. Brain Pathol. 2007;17(2):210–8.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Catuneanu A, Paylor JW, Winship I, Colbourne F, Kerr BJ. Sex differences in central nervous system plasticity and pain in experimental autoimmune encephalomyelitis. Pain. 2019;160(5):1037–49.

    Article  PubMed  Google Scholar 

  3. Mifflin KA, Yousuf MS, Thorburn KC, Huang J, Pérez-Muñoz ME, Tenorio G, et al. Voluntary wheel running reveals sex-specific nociceptive factors in murine experimental autoimmune encephalomyelitis. Pain. 2019;160(4):870–81.

    Article  CAS  PubMed  Google Scholar 

  4. Mifflin KA, Frieser E, Benson C, Baker G, Kerr BJ. Voluntary wheel running differentially affects disease outcomes in male and female mice with experimental autoimmune encephalomyelitis. J Neuroimmunol. 2017;305:135–44.

    Article  CAS  PubMed  Google Scholar 

  5. Foley PL, Vesterinen HM, Laird BJ, Sena ES, Colvin LA, Chandran S, et al. Prevalence and natural history of pain in adults with multiple sclerosis: systematic review and meta-analysis. Pain. 2013;154(5):632–42.

    Article  PubMed  Google Scholar 

  6. Woolf CJ, Mannion RJ. Neuropathic pain: aetiology, symptoms, mechanisms, and management. Lancet. 1999;353(9168):1959–64.

    Article  CAS  PubMed  Google Scholar 

  7. Itoh N, Itoh Y, Tassoni A, Ren E, Kaito M, Ohno A, et al. Cell-specific and region-specific transcriptomics in the multiple sclerosis model: focus on astrocytes. Proc Natl Acad Sci U S A. 2018;115(2):E302–9.

    Article  CAS  PubMed  Google Scholar 

  8. Gillett A, Maratou K, Fewings C, Harris RA, Jagodic M, Aitman T, et al. Alternative splicing and transcriptome profiling of experimental autoimmune encephalomyelitis using genome-wide exon arrays. PLoS One. 2009;4(11):e7773.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  9. Juźwik CA, Drake S, Lécuyer M-A, Johnson RM, Morquette B, Zhang Y, et al. Neuronal microRNA regulation in experimental autoimmune encephalomyelitis. Sci Rep. 2018;8(1):13437.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Hammond SM. An overview of microRNAs. Adv Drug Deliv Rev. 2015;87:3–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Vidigal JA, Ventura A. The biological functions of miRNAs: lessons from in vivo studies. Trends Cell Biol. 2015;25(3):137–47.

    Article  CAS  PubMed  Google Scholar 

  12. Billiau A, Matthys P. Modes of action of Freund’s adjuvants in experimental models of autoimmune diseases. J Leukoc Biol. 2001;70(6):849–60.

    CAS  PubMed  Google Scholar 

  13. Team R. R: a language and environment for statistical computing; 2013.

    Google Scholar 

  14. Lawrence M, Huber W, Pagès H, Aboyoun P, Carlson M, Gentleman R, et al. Software for computing and annotating genomic ranges. PLoS Comput Biol. 2013;9(8):e1003118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Morgan M, Pagès H, Obenchain VHN. Rsamtools: binary alignment (BAM), FASTA, variant call (BCF), and tabix file import; 2018.

    Google Scholar 

  16. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. Eden E, Lipson D, Yogev S, Yakhini Z. Discovering motifs in ranked lists of DNA sequences. PLoS Comput Biol. 2007;3(3):e39.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6(7):e21800.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Krämer A, Green J, Pollard J Jr, Tugendreich S. Causal analysis approaches in ingenuity pathway analysis. Bioinformatics. 2014;30(4):523.

    Article  PubMed  CAS  Google Scholar 

  20. Thorburn KC, Paylor JW, Webber CA, Winship IR, Kerr BJ. Facial hypersensitivity and trigeminal pathology in mice with experimental autoimmune encephalomyelitis. Pain. 2016;157(3):627–42.

    Article  PubMed  Google Scholar 

  21. Huang DW, Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, et al. The DAVID gene functional classification tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007;8(9):R183.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Fabregat A, Jupe S, Matthews L, Sidiropoulos K, Gillespie M, Garapati P, et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018;46(D1):D649–55.

    Article  CAS  PubMed  Google Scholar 

  23. Mi H, Muruganujan A, Huang X, Ebert D, Mills C, Guo X, et al. Protocol update for large-scale genome and gene function analysis with the PANTHER classification system (v.14.0). Nat Protoc. 2019;14(3):703–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Mi H, Thomas P. PANTHER pathway: an ontology-based pathway database coupled with data analysis tools; 2009. p. 123–40.

    Google Scholar 

  25. Mi H, Muruganujan A, Ebert D, Huang X, Thomas PD. PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res. 2019;47(D1):D419–26.

    Article  CAS  PubMed  Google Scholar 

  26. Klein SL, Flanagan KL. Sex differences in immune responses. Nat Rev Immunol. 2016;16(10):626–38.

    Article  CAS  PubMed  Google Scholar 

  27. Whitacre CC. Sex differences in autoimmune disease. Nat Immunol. 2001;2(9):777–80.

    Article  CAS  PubMed  Google Scholar 

  28. Rosen S, Ham B, Mogil JS. Sex differences in neuroimmunity and pain. J Neurosci Res. 2017;95(1–2):500–8.

    Article  CAS  PubMed  Google Scholar 

  29. Fairweather D. Sex differences in inflammation during atherosclerosis. Clin Med Insights Cardiol. 2014;8(Suppl 3):49–59.

    PubMed  Google Scholar 

  30. Berkley KJ, Zalcman SS, Simon VR. Sex and gender differences in pain and inflammation: a rapidly maturing field. Am J Physiol Integr Comp Physiol. 2006;291(2):R241–4.

    Article  CAS  Google Scholar 

  31. Wu S, Marie Lutz B, Miao X, Liang L, Mo K, Chang Y-J, et al. Dorsal root ganglion transcriptome analysis following peripheral nerve injury in mice. Mol Pain. 2016;12:174480691662904.

    Article  CAS  Google Scholar 

  32. Gong L, Wu J, Zhou S, Wang Y, Qin J, Yu B, et al. Global analysis of transcriptome in dorsal root ganglia following peripheral nerve injury in rats. Biochem Biophys Res Commun. 2016;478(1):206–12.

    Article  CAS  PubMed  Google Scholar 

  33. Yousuf MS, Noh M-C, Friedman TN, Zubkow K, Johnson JC, Tenorio G, et al. Sensory neurons of the dorsal root ganglia become hyperexcitable in a T-cell-mediated MOG-EAE model of multiple sclerosis. eNeuro. 2019;6(2) ENEURO.0024-19.2019.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Laffont S, Rouquié N, Azar P, Seillet C, Plumas J, Aspord C, et al. X-chromosome complement and estrogen receptor signaling independently contribute to the enhanced TLR7-mediated IFN-α production of plasmacytoid dendritic cells from women. J Immunol. 2014 Dec 1;193(11):5444–52.

    Article  CAS  PubMed  Google Scholar 

  35. Filippi M, Bar-Or A, Piehl F, Preziosa P, Solari A, Vukusic S, et al. Multiple sclerosis. Nat Rev Dis Prim. 2018 Dec 8;4(1):43.

    Article  PubMed  Google Scholar 

  36. Ivashkiv LB. IFNγ: signalling, epigenetics and roles in immunity, metabolism, disease and cancer immunotherapy. Nat Rev Immunol. 2018 Sep 19;18(9):545–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Noris M, Remuzzi G. Overview of complement activation and regulation. Semin Nephrol. 2013 Nov;33(6):479–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Sarma JV, Ward PA. The complement system. Cell Tissue Res. 2011 Jan;343(1):227–35.

    Article  CAS  PubMed  Google Scholar 

  39. Junker A, Krumbholz M, Eisele S, Mohan H, Augstein F, Bittner R, et al. MicroRNA profiling of multiple sclerosis lesions identifies modulators of the regulatory protein CD47. Brain. 2009;132(12):3342–52.

    Article  PubMed  Google Scholar 

  40. Lescher J, Paap F, Schultz V, Redenbach L, Scheidt U, Rosewich H, et al. MicroRNA regulation in experimental autoimmune encephalomyelitis in mice and marmosets resembles regulation in human multiple sclerosis lesions. J Neuroimmunol. 2012;246(1–2):27–33.

    Article  CAS  PubMed  Google Scholar 

  41. Talebi F, Ghorbani S, Chan WF, Boghozian R, Masoumi F, Ghasemi S, et al. MicroRNA-142 regulates inflammation and T cell differentiation in an animal model of multiple sclerosis. J Neuroinflammation. 2017;14(1):55.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Vaknin-Dembinsky A, Charbit H, Brill L, Abramsky O, Gur-Wahnon D, Ben-Dov IZ, et al. Circulating microRNAs as biomarkers for rituximab therapy, in neuromyelitis optica (NMO). J Neuroinflammation. 2016 Dec 8;13(1):179.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Folini M, Gandellini P, Longoni N, Profumo V, Callari M, Pennati M, et al. miR-21: an oncomir on strike in prostate cancer. Mol Cancer. 2010;9:12.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Krichevsky AM, Gabriely G. miR-21: a small multi-faceted RNA. J Cell Mol Med. 2009;13(1):39–53.

    Article  CAS  PubMed  Google Scholar 

  45. Medina PP, Nolde M, Slack FJ. OncomiR addiction in an in vivo model of microRNA-21-induced pre-B-cell lymphoma. Nature. 2010;467(7311):86–90.

    Article  CAS  PubMed  Google Scholar 

  46. Strickland IT, Richards L, Holmes FE, Wynick D, Uney JB, Wong L-F. Axotomy-induced miR-21 promotes axon growth in adult dorsal root ganglion neurons. PLoS One. 2011;6(8):e23423.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Murugaiyan G, da Cunha AP, Ajay AK, Joller N, Garo LP, Kumaradevan S, et al. MicroRNA-21 promotes Th17 differentiation and mediates experimental autoimmune encephalomyelitis. J Clin Invest. 2015;125(3):1069–80.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Simeoli R, Montague K, Jones HR, Castaldi L, Chambers D, Kelleher JH, et al. Exosomal cargo including microRNA regulates sensory neuron to macrophage communication after nerve trauma. Nat Commun. 2017;8(1):1778.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. Zhang Z-J, Guo J-S, Li S-S, Wu X-B, Cao D-L, Jiang B-C, et al. TLR8 and its endogenous ligand miR-21 contribute to neuropathic pain in murine DRG. J Exp Med. 2018;215(12):3019–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Fabbri M, Paone A, Calore F, Galli R, Gaudio E, Santhanam R, et al. MicroRNAs bind to Toll-like receptors to induce prometastatic inflammatory response. Proc Natl Acad Sci U S A. 2012;109(31):E2110–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Ribas J, Ni X, Haffner M, Wentzel EA, Salmasi AH, Chowdhury WH, et al. miR-21: an androgen receptor-regulated microRNA that promotes hormone-dependent and hormone-independent prostate cancer growth. Cancer Res. 2009;69(18):7165–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Queirós AM, Eschen C, Fliegner D, Kararigas G, Dworatzek E, Westphal C, et al. Sex- and estrogen-dependent regulation of a miRNA network in the healthy and hypertrophied heart. Int J Cardiol. 2013;169(5):331–8.

    Article  PubMed  Google Scholar 

  53. Wickramasinghe NS, Manavalan TT, Dougherty SM, Riggs KA, Li Y, Klinge CM. Estradiol downregulates miR-21 expression and increases miR-21 target gene expression in MCF-7 breast cancer cells. Nucleic Acids Res. 2009;37(8):2584–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Wu Y, Gao C, Cai S, Xia M, Liao G, Zhang X, et al. Circulating miR-122 is a predictor for virological response in CHB patients with high viral load treated with nucleos(t)ide analogs. Front Genet. 2019;10:243.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Park H-K, Jo W, Choi H-J, Jang S, Ryu J-E, Lee H-J, et al. Time-course changes in the expression levels of miR-122, −155, and −21 as markers of liver cell damage, inflammation, and regeneration in acetaminophen-induced liver injury in rats. J Vet Sci. 2016;17(1):45–51.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Yoshikawa T, Takata A, Otsuka M, Kishikawa T, Kojima K, Yoshida H, et al. Silencing of microRNA-122 enhances interferon-α signaling in the liver through regulating SOCS3 promoter methylation. Sci Rep. 2012;2(1):637.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references


The authors wish to thank G. Tenorio for assistance in making the figures.


Funding for this project was provided by operating grants from the MS Society of Canada (MSSC), the Alberta MS Collaboration, and the University Hospital Fund (University of Alberta). TNF and AC were supported by studentships from the MS Society of Canada. MSY was supported by the Alexander Graham Bell Canada Graduate Scholarship from NSERC.

Author information

Authors and Affiliations



AC conducted the EAE experiment. TNF, MSY, and MD performed the tissue extraction and qPCR validation. TNF analyzed and interpreted the sequencing data. TNF, CAJ, AEF, and BJK were involved in experimental design and manuscript preparation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bradley J. Kerr.

Ethics declarations

Ethics approval and consent to participate

All animal studies complied with the Canadian Council on Animal Care Guidelines and Policies and were approved by the Animal Care and Use Committee: Health Sciences for the University of Alberta.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Spreadsheet of RNASeq differential expression analysis between EAE and CFA condition in females.

Additional file 2: Table S2.

Spreadsheet of RNASeq differential expression analysis between EAE and CFA condition in males.

Additional file 3: Table S3.

Spreadsheet of miRSeq differential expression analysis between EAE and CFA condition in females.

Additional file 4: Table S4.

Spreadsheet of miRSeq differential expression analysis between EAE and CFA condition in males.

Additional file 5: Table S5.

Spreadsheet of RNASeq differential expression analysis between EAE conditions of females and males.

Additional file 6: Table S6.

Spreadsheet of Ingenuity Pathway Analysis of differentially expressed genes between female and male EAE conditions.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Friedman, T.N., Yousuf, M.S., Catuneanu, A. et al. Profiling the microRNA signature of the peripheral sensory ganglia in experimental autoimmune encephalomyelitis (EAE). J Neuroinflammation 16, 223 (2019).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: