Identification of potential regulatory long non-coding RNA-associated competing endogenous RNA axes in periplaque regions in multiple sclerosis

Slow-burning inflammation at the lesion rim is connected to the expansion of chronic multiple sclerosis (MS) lesions. However, the underlying processes causing expansion are not clearly realized. In this context, the current study used a bioinformatics approach to identify the expression profiles and related lncRNA-associated ceRNA regulatory axes in the periplaque region in MS patients. Expression data (GSE52139) from periplaque regions in the secondary progressive MS spinal cord and controls were downloaded from the Gene Expression Omnibus database (GEO), which has details on mRNAs and lncRNAs. Using the R software’s limma package, the differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) were found. The RNA interactions were also found using the DIANA-LncBase, miRTarBase, and HMDD databases. The Pearson correlation coefficient was used to determine whether there were any positive correlations between DEmRNAs and DElncRNAs in the ceRNA network. Finally, lncRNA-associated ceRNA axes were created based on co-expression and connections between DElncRNA, miRNA, and DEmRNA. We used the Enrichr tool to enrich the biological process, molecular function, and pathways for DEmRNAs and DElncRNAs. A network of DEmRNAs’ protein-protein interactions was developed, and the top five hub genes were found using Cytoscape and STRING. The current study indicates that 15 DEmRNAs, including FOS, GJA1, NTRK2, CTNND1, and SP3, are connected to the MS ceRNA network. Additionally, four DElncRNAs (such as TUG1, ASB16-AS1, and LINC01094) that regulated the aforementioned mRNAs by sponging 14 MS-related miRNAs (e.g., hsa-miR-145-5p, hsa-miR-200a-3p, hsa-miR-20a-5p, hsa-miR-22-3p, hsa-miR-23a-3p, hsa-miR-27a-3p, hsa-miR-29b-3p, hsa-miR-29c-3p, hsa-miR-34a-5p) were found. In addition, the analysis of pathway enrichment revealed that DEmRNAs were enriched in the pathways for the “MAPK signaling pathway”, “Kaposi sarcoma-associated herpesvirus infection”, “Human immunodeficiency virus one infection”, “Lipid and atherosclerosis”, and “Amphetamine addiction”. Even though the function of these ceRNA axes needs to be investigated further, this study provides research targets for studying ceRNA-mediated molecular mechanisms related to periplaque demyelination in MS.

Introduction

Multiple sclerosis (MS) is a neurodegenerative, inflammatory disease affecting the central nervous system (CNS). In MS, it is now understood that axonal loss is the primary cause of permanent neurological impairment. Axonal dissection followed by axonal degeneration are thought to be primarily caused by acute inflammatory demyelination, which generally appears as new white matter lesions. Recent research suggests that slow-burning inflammatory demyelination at the edge of some chronic lesions is a key contributor to disease progression, including progressive axonal loss as well as worsening impairment (Prineas et al., 2001; Frischer et al., 2015; Dal-Bianco et al., 2017; Klistorner et al., 2021). These lesions, commonly referred to as smoldering or chronically active lesions (van der Valk and De Groot, 2000), are distinguished by a fully demyelinated, inactive hypocellular center and low-grade active demyelination in periplaque white matter, which has activated microglia/macrophages with myelin breakdown products (Lieury et al., 2014; Frischer et al., 2015). Despite the fact that the pathogenesis of chronic smoldering lesions is widely understood (Prineas et al., 2001), the periplaque region has received little attention up to this point. Although it is generally acknowledged that a combination of genetic and environmental factors may influence a person’s predisposition for MS, the underlying pathophysiological mechanisms and the sequence of events that contribute to lesion expansion remain poorly understood. The major histocompatibility complex (MHC) predominately determines disease occurrence in the genetically susceptible population, whereas modifiable environmental factors like smoking, Epstein-Barr virus (EBV) infection, increased body mass index (BMI) during adolescence, and low vitamin D levels may influence whether a person develops MS (Amato et al., 2018). The human leukocyte antigen (HLA) locus, which is known to encode molecules involved in essential immune functions, is the first genetic risk factor that was found decades ago (Wang et al., 2022). New research suggests that non-coding RNAs (ncRNAs), particularly long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), have the potential to regulate gene expression and offer novel insights into the development of MS (Sheng et al., 2015; Yang et al., 2018; Ghafouri-Fard and Taheri, 2020).

The competing endogenous RNA (ceRNA) hypothesis states that ncRNAs, which serve as miRNA “sponges,” compete with target mRNAs for binding to miRNAs as a novel regulatory mechanism (Salmena et al., 2011). This theory proposes that cross-talk between RNAs, including coding RNAs and ncRNAs, through miRNA complementary sequences known as miRNA response elements (MREs) generates a large-scale regulatory network throughout the transcriptome. If two RNA transcripts control each other via a ceRNA-mediated mechanism, then, according to the ceRNA hypothesis, the expression levels of these two RNA transcripts would be negatively correlated with the levels of target miRNAs and positively correlated with each other (Salmena et al., 2011). A few studies have suggested that MS may be caused by altered ceRNA-mediated gene regulation (Bian et al., 2020; Ding et al., 2021; Hao et al., 2022; Karimi et al., 2022; Wang et al., 2022). For instance, it has been demonstrated that the lncRNA Gm15575 is aberrantly expressed in MS patients and affects the functionality of Th17 in MS via ceRNA patterns (Bian et al., 2020). In a separate research, it was found that HOTAIR functions as a sponge for miR-136-5p, increasing AKT2-mediated NF-kB activation and hence favoring the microglial shift towards a proinflammatory M1-like phenotype (Duan et al., 2018), which is harmful in MS. Furthermore, lncRNA GAS5 has been postulated as a ceRNA for miR-137 involved in demyelination (Senousy et al., 2020). It is important to do more research on the associated expression patterns and processes in MS due to the fact that the so-called study on the ceRNA axes implicated in MS, particularly in the periplaque region, has received little attention to date. RNAs have greater “druggability” than proteins because their targeting primarily depends on sequence complementarity. Due to their unique properties, these molecules can be easily and inexpensively manufactured into novel targets for RNA-related drugs (Yu et al., 2020). CeRNA interaction networks involve a variety of factors, making it possible to study complex conditions like MS focusing on only one of the potential treatment targets (i.e., an instant variation in the levels of multiple disease-related RNAs) (Moreno-García et al., 2020).

This study aimed to identify the expression profiles and associated lncRNA-associated ceRNA regulatory axes in the periplaque region of MS patients using a bioinformatics approach.

Methods

In the current study, we used a system biology method for mining expression data related to periplaque regions in the spinal cord of MS with the accession number GSE52139. Our purposes were recognition of the differentially expressed lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) and to develop lncRNA-associated ceRNA regulatory axes by using previous bioinformatics approaches (Sabaie et al., 2021a; Sabaie et al., 2021b).

Data preparation of the gene expression profile

The gene expression profiles were downloaded from the NCBI Gene Expression Omnibus dataset (GEO, https://www.ncbi.nlm.nih.gov/geo/). The platform of GPL570 (HG-U133_Plus_2) Affymetrix Human Genome U133 and 2.0 Array (chip-based) was used for the mentioned database, in which both information of lncRNA and mRNA were included in the GPL570. The GSE52139 includes sixteen spinal cord samples. Eight periplaque samples were compared to normal-appearing white matter from identical patients (Lieury et al., 2014).

Data preprocessing and identifying DEmRNAs and DElncRNAs

Background correction and quantile normalization were done through the Robust Multichip Average (RMA) for the total files of data (Irizarry et al., 2003). The AgiMicroRna Bioconductor package version 2.46.0 was used for quality control (Romero and AgiMicroRna, 2020). The differential expression gene analysis (DEGA) was done by the limma package of R (version 3.52.2) (Ritchie et al., 2015) in Bioconductor (https://www.bioconductor.org/) (Huber et al., 2015). The linear model was accomplished via limma’s lmFit () function doing on the main factor “group” (i.e., periplaque against normal-appearing white matter). Subsequently, limma’s eBayes () function was applied to calculate differentially expressed genes from the linear fit model. Next, the whole list of lncRNAs’ genes was obtained (https://www.genen ames. org/) through approved HUGO Gene Nomenclature Committee (HGNC) symbols (Braschi et al., 2019). The dataset’s gene symbols were compared with the list of lncRNA genes, and the overlapped genes were chosen. According to this approach, we captured 1,070 lncRNAs. The paired student’s t-test was used to identify statistically significant mRNAs and lncRNAs. The following cut-offs were established for aberrantly expressed genes: (Prineas et al., 2001): |log2 fold change (log2FC) ≥ 1, and q < 0.01 (q is the False Discovery Rate (FDR) corrected p-value) for DEmRNAs and (Dal-Bianco et al., 2017) (log2FC) ≥ 0.585, and q < 0.01 for DElncRNAs. The implementation of less stringent selection criteria was prompted by the lower expression level of lncRNAs relative to mRNAs. Finally, using the Pheatmap (version 1.0.12) and the Enhanced Volcano (version1.14.0) R packages, the heat map and the volcano plot of DEGs were depicted.

RNA interaction pairs prediction

In this study, the DIANA-LncBase v3 (Karagkouni et al., 2020) was utilized to predict the interaction of the lncRNAs and miRNAs, which were experimentally confirmed. The criteria utilized in the DIANA-LncBase query were High “miRNA Confidence Levels” and Homo Sapiens “Species”. Additionally, the MS-related miRNAs were collected from the Human microRNA Disease Database (HMDD) version 3.2 (Huang et al., 2019). Also, interactions between the target mRNAs and collected miRNAs (used the MS-related miRNAs) were attained from miRTarBase (including persistent experimental shreds of evidence) (Huang et al., 2020). Then, these newly attained mRNAs were compared to previous mRNAs attained, and consequently, the duplicate mRNAs were used to construct the lncRNA-miRNA-mRNA ceRNA axes.

Analysis of correlation between DEmRNAs and DElncRNAs, and lncRNA-associated ceRNA axes construction

The Pearson correlation analysis was done to study the positive correlations in the regulatory axes of ceRNA among DEmRNAs and DElncRNAs. In the reverse expression pattern among DElncRNAs and the targeted DEmRNAs, interacted miRNAs targeted DEmRNAs, and DElncRNAs were removed from the network of ceRNA. For correlations’ visualization and calculation, both packages of Hmisc (version 4.7-1) and corrplot (version 0.92) were used. The Pearson correlation coefficients above 0.5 and p < 0.05 were considered inclusion criteria. Finally, for the construction of ceRNA regulatory axes, we applied the Cytoscape software (version 3.8.0) (Shannon et al., 2003).

DEmRNAs and DElncRNAs enrichment analysis

Biological process (BP), molecular function (MF) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of the DEmRNAs in the ceRNA network were obtained via the Enrichr tool (Chen et al., 2013; Kuleshov et al., 2016). In addition, BP and MF analyses were performed utilizing FuncPred (https://www.funcpred.com) (Perron et al., 2017) to get deeper insight into the DElncRNAs in ceRNA network.

Construction of protein-protein interaction network and hub genes identification

A PPI network was made by the online STRING database (https://string-db.org/) (Szklarczyk et al., 2019) to predict the interactive relevance of DEmRNAs encoding proteins. The Cytoscape software was used to visualize the PPI network and the hub genes analysis. Moreover, non-interacting genes were removed from the PPI network to simplify it. Subsequently, the top five genes of the PPI network were evaluated by Maximal Clique Centrality (MCC) method by use of CytoHubba (version 0.1) in Cytoscape (Chin et al., 2014). Compared to the other network scoring methods, MCC performed better (Chin et al., 2014).

Restoration of the subceRNA network

The DElncRNAs and miRNAs relevant to the primary structure of ceRNA were extracted and utilized for the network restoration of the ceRNA network hub genes.

Results

DEmRNAs and DElncRNAs identification

At first, background adjustment, normalization, and batch adjustment were performed. Following the normalization, box plots were illustrated to analyze the data distribution in the gene expression data (Appendix I). The correct adjustment was approved by the outcomes from distinct box plot arrays that show similar medians for expression levels.

Regarding to the chosen criteria (|log2FC | ≥ 1, and q < 0.01 for DEmRNAs and |log2FC | ≥ 0.585, and q < 0.01 for DElncRNAs), we detected 193 DEmRNAs and seven DElncRNAs between periplaque and normal-appearing white matter samples. The DEmRNAs volcano plot and DElncRNAs heat map are shown in Figures 1, 2, respectively. Furthermore, details of DEGs are summarized in Appendix II and III.

FIGURE 1

FIGURE 1. Differentially expressed mRNAs (DEmRNAs) volcano plot. A |(log2FC)| ≥ 1 and an p < 0.01 were utilized for screening DEmRNAs.

FIGURE 2

FIGURE 2. Differentially expressed long non-coding RNAs (DElncRNAs) Heatmap. Lowly expressed DElncRNAs are displayed in blue, while those that are highly expressed are displayed in red.

Predicting RNA interaction pairs

Six of the seven DElncRNAs might be targeted via the candidate miRNAs. Then, the interplay between candidate miRNAs (obtained from HMDD) and mRNAs was uncovered by miRTarBase. Eventually, 19 overlapping genes were obtained after comparing the candidate mRNAs and DEmRNAs.

Analysis of relationship between DElncRNAs and DEmRNAs, and lncRNA-associated ceRNA axes construction

At this level, to prove the assumption of ceRNA axes, i.e., lncRNA positively regulates (via miRNA interaction) mRNA’s expression; hence, we applied the Pearson correlation analysis among DEmRNAs and DElncRNAs (Figure 3). According to the interactions between DElncRNA-miRNA-DEmRNA and co-expression associations, we constructed the ceRNA regulatory axes to declare the related ceRNA axes with lncRNAs in periplaque regions in MS (Figure 4). Table 1 elucidates the ceRNA axes characteristics, including four DElncRNAs, 14 miRNAs, and 15 DEmRNAs, in total.

FIGURE 3

FIGURE 3. Pearson correlation analysis. Blue is used to represent positive correlations, while red is used to represent negative correlations. p greater than 0.05 are regarded as being insignificant (blank), and the intensity is related to correlation coefficients.

FIGURE 4

FIGURE 4. Long non-coding RNA (lncRNA)-associated competing endogenous RNA axes. Red represents upregulation of ceRNAs. LncRNAs, miRNAs, and mRNAs are represented by hexagon, round rectangle, and ellipse, respectively.

TABLE 1

TABLE 1. Details of competing endogenous RNA axes.

DEmRNAs and DElncRNAs enrichment analysis

Functional annotation obtained on whole DEmRNAs in the ceRNA network. Figures 5A,B indicated that target mRNAs were associated with “skeletal muscle tissue regeneration”, “cardiac conduction system development”, “positive regulation of nucleic acid-templated transcription”, “cardiac muscle tissue development”, “cellular response to reactive oxygen species”, “gap junction channel activity involved in cardiac conduction electrical coupling”, “calmodulin-dependent protein phosphatase activity”, “glutathione transmembrane transporter activity”, “tripeptide transmembrane transporter activity”, and “gap junction channel activity involved in cell communication by electrical coupling”. Moreover, the top enriched KEGG pathways were as follows: “MAPK signaling pathway”, “Kaposi sarcoma-associated herpes virus infection”, “Human immunodeficiency virus one infection”, “Lipid and atherosclerosis”, and “Amphetamine addiction” (Figure 5). The complete list of the enrichment analysis of DElncRNAs is included in Supplementary Files S1–S4.

FIGURE 5

FIGURE 5. Biological process (BP), molecular function (MF) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The top five enriched terms of (A) BP, (B) MF, and (C) KEGG are displayed along with the corresponding p-values in a bar graph. Those with significant p-values (i.e., <0.05) are matched with the colored bars. A p-value (shown with *) shows the significant adjusted p-value (<0.05).

Construction of PPI network and identification of hub genes

The PPI network was demonstrated in Figure 6A. The MCC method was applied to assess the top five genes in the PPI network (Figure 6B). The Hub genes with the strongest association are shown with nodes in red color. Besides, orange and yellow nodes illustrate the hub genes with moderate and poor connections, respectively. Fos Proto-Oncogene (FOS) and Gap Junction Protein Alpha 1 (GJA1) had the strongest correlation. The Neurotrophic Receptor Tyrosine Kinase 2 (NTRK2) and Catenin Delta 1 (CTNND1) had moderate relevance. In contrast, Sp3 Transcription Factor (SP3) demonstrated poor connection.

FIGURE 6

FIGURE 6. (A) The DEmRNAs interactions are shown as blue nodes in this PPI network. (B) The recognized hub genes in the PPI network. The yellow, orange, and red nodes, respectively, represent hub genes with weak, moderate, and strong connections.

Reconstruction of the subceRNA network

The DElncRNA-DEmiRNA-hub genes sub-network was reconstructed by use of the hub genes. There were three DElncRNAs (TUG1: Taurine upregulated 1, ASB16-AS1: ASB16 Antisense RNA 1, and LINC01094: Long Intergenic Non-Protein Coding RNA 1094), nine miRNAs (hsa-miR-145-5p, hsa-miR-200a-3p, hsa-miR-20a-5p, hsa-miR-22-3p, hsa-miR-23a-3p, hsa-miR-27a-3p, hsa-miR-29b-3p, hsa-miR-29c-3p, hsa-miR-34a-5p) and five hub genes (FOS, GJA1, NTRK2, CTNND1, SP3) overall (Figure 7).

FIGURE 7

FIGURE 7. The lncRNA-miRNA-hub gene subceRNA axes. Red represents upregulation of ceRNAs. LncRNAs, miRNAs, and mRNAs are represented by hexagon, round rectangle, and ellipse, respectively.

Discussion

Several studies have found that ceRNA regulatory axes and connected networks play an active role in a variety of developmental and pathological processes, including tumor development and a variety of brain-related illnesses (Ala, 2020; Moreno-García et al., 2020). On the basis of tissue-specific, cellular, and subcellular circumstances, the ceRNA expresses itself differently. In addition to mRNAs, a network may contain a variety of ceRNAs, such as lncRNAs, circRNAs, and pseudogenes. One of the primary RNA types, lncRNAs, is found in the ceRNA machinery and exerts a major impact on biological pathways in both healthy and diseased states (Cai and Wan, 2018). Currently, there is complete consensus that lncRNA expression varies depending on the tissue, cellular type, and developmental stage. In addition to subcellular dispersions, such a distinct tissue dependence is a clear sign that the expression of lncRNAs is tightly regulated (Gloss and Dinger, 2016). The theoretic ideas indicate that the ceRNA regulatory axes linked to lncRNAs can significantly contribute to MS pathogenicity. Based on these ceRNA axes and connections between DElncRNA and DEmRNA, we found the DElncRNA-miRNA-DEmRNA network in the current study. This network contained three important DElncRNAs (TUG1, ASB16-AS1, and LINC01094), nine important miRNAs (hsa-miR-145-5p, hsa-miR-200a-3p, hsa-miR-20a-5p, hsa-miR-22-3p, hsa-miR-23a-3p, hsa-miR-27a-3p, hsa-miR-29b-3p, hsa-miR-29c-3p, hsa-miR-34a-5p) and five hub genes (FOS, GJA1, NTRK2, CTNND1, SP3).

TUG1 is a 7.1-kb lncRNA whose function is connected to the normal development of the retina and nervous system. It was first found to be upregulated in response to taurine treatment of developing retinal cells (Young et al., 2005). In experimental autoimmune encephalomyelitis (EAE) mice and lipopolysaccharide-induced BV2 cells, Yue and colleagues discovered that down-regulation of TUG1 attenuates MS through the inhibition of inflammation by sponging miR-9-5p via targeting NF-κB1/p50 (Yue et al., 2019). Multiple lines of research have shown that patients with secondary progressive (Santoro et al., 2020) and relapsing-remitting (Santoro et al., 2016; Dastmalchi et al., 2018) MS exhibit increased TUG1 expression. For instance, reverse transcription quantitative PCR (RT-qPCR) analysis showed that TUG1 expression was significantly higher in peripheral serum (Santoro et al., 2016; Santoro et al., 2020) and whole blood (Dastmalchi et al., 2018) in patients with MS compared with controls. These results support our findings that TUG1 plays a role in the pathophysiology of MS. TUG1 was, however, downregulated in peripheral blood mononuclear cells (PBMCs) from relapsing-remitting MS patients, according to another expression study (Fenoglio et al., 2018). According to the authors, this contradictory outcome may be due to the various biological sources used for the analysis or a potentially active role for TUG1 in intercellular communication. Its increased levels in the bloodstream could result from exosomes, an enriched source of ncRNA, releasing it into the cellular environment (Fenoglio et al., 2018). ASB16-AS1 is localized to 17q21 and is approximately 2,275 bp (Li et al., 2021). The lncRNA ASB16-AS1 regulates cell proliferation, migration, invasion, and apoptosis through a ceRNA manner in a number of cancers, including glioma (Zhang et al., 2019) and cervical cancer, according to earlier research (Liu et al., 2020). A new long intergenic nonprotein coding RNA called LINC01094 is found on chromosome four and functions as a ceRNA in a variety of conditions, such as calcified aortic valve disease (Huang et al., 2022), clear cell renal cell carcinoma (Jiang et al., 2020), schizophrenia (Sabaie et al., 2021c), and glioma (Liu et al., 2022). Since our study is the first to identify a link between ASB16-AS1, LINC01094, and MS, the findings should be confirmed through further research. In agreement with our findings, glioma (Zhang et al., 2019) and schizophrenia (Sabaie et al., 2021c), two additional nervous system diseases, exhibit significantly higher expression levels of ASB16-AS1 and LINC01094, respectively.

We predicted that the sponging of nine key miRNAs (hsa-miR-145-5p, hsa-miR-200a-3p, hsa-miR-20a-5p, hsa-miR-22-3p, hsa-miR-23a-3p, hsa-miR-27a-3p, hsa-miR-29b-3p, hsa-miR-29c-3p, hsa-miR-34a-5p) by key lncRNAs might influence target genes. MiRNAs regulate target gene expression by binding to the non-transcript site of the targeted gene, thereby influencing signal transduction and biological pathways within a cell, which may have an impact on the beginning and progression of MS (Caputo et al., 2015). Based on the HMDD database, the majority of the mentioned miRNAs fall under the circulation biomarker class of literature evidence. However, their precise roles in MS are still unknown. Although our findings are supported by the dysregulation of these miRNAs in MS, the predicted ceRNA axes still need to be verified using molecular methods.

The enrichment analysis results indicated that DEmRNAs’ functions were focused on different BP and MF, including skeletal muscle characteristics, cardiac function, gap junction, protein phosphatase activity, and transmembrane transporter activity. Interestingly, the majority of the DEmRNAs in our study were found to be involved in the “MAPK signaling pathway”, “Kaposi sarcoma-associated herpesvirus infection”, “Human immunodeficiency virus one infection”, “Lipid and atherosclerosis”, and “Amphetamine addiction”, according to the results of our KEGG enrichment analysis. According to the data, the improper (over)activity of the mitogen-activated protein kinase pathway ERK (MAPK^ERK) causes microglial dysfunction. These findings concern both biochemistry and epigenetics, and they all point to the involvement of this pathway. Recent preclinical studies on neurodegeneration have already suggested that MAPK pathways, in particular MAPK^ERK, are involved in the process. This is crucial because it has been discovered that microglia with overactive MAPK disrupt local oligodendrocytes, which can result in locoregional demyelination, a characteristic of MS. This represents a novel idea in the pathophysiology of MS in contrast to the prevalent theory of autoimmunity. Smoking, hypovitaminosis D, and EBV infection are all known risk factors for MS that inhibit the MAPK^ERK negative feedback phosphatases that normally control MAPK^ERK activity. Due to inappropriate MAPK^ERK overactivity and subsequent neurodegeneration, these factors may be involved. A contributing factor in the pathophysiology of MS is MAPK^ERK overactivity in microglia, which may also explain why MS patients continue to experience neurodegeneration despite receiving optimal immunosuppressive or immunomodulatory therapy (ten Bosch et al., 2021). MS and Kaposi’s sarcoma, for many years, have been suspected to be of viral origin (Enbom, 2001; Tully et al., 2015; Walker and Brew, 2016). New lymphotropic herpesviruses have recently been linked to both of these diseases. Many potential oncogenes found in the viral genome are thought to have originated from the human genome during evolution. In conclusion, the full range of human diseases caused by these viruses is still not fully understood, and the fast-moving field of molecular biology will continue to provide fresh insights into how herpesviruses interact with their hosts (Enbom, 2001). HIV infection and MS are common and well-researched nosological entities. However, MS and HIV co-morbidity have only occasionally been documented in the medical literature. Recent epidemiological studies showed a decreased risk of MS in HIV-positive patients, indicating a negative association between MS and HIV. Increasing clinical evidence also points to a potential lower MS relapse rate in HIV patients. However, it is currently unknown whether this inverse correlation was caused by HIV infection, HIV treatment, or a combination of the two (Stefanou et al., 2019). When compared to the general population, MS patients were reported to have a higher mortality rate, which was attributed to a number of things, including the population’s higher prevalence of cardiovascular diseases (Manouchehrinia et al., 2016). Strong evidence links MS to an increased risk for stroke, myocardial infarction, and heart failure (Christiansen et al., 2010). The primary causes of vascular dysfunction in MS are not fully understood, but the presence of chronic inflammation and autoimmunity, oxidative stress, and the dysregulation of the cardiovascular autonomic in MS patients have all been proposed as potential causes. Atherosclerosis, arteriosclerosis, and endothelial dysfunction are all significantly influenced by inflammation. Additionally, it permits plaque rupture, which raises the danger of acute coronary syndrome (Crea and Libby, 2017). The amphetamines have central stimulant action and sympathomimetic properties. Clinical signs of toxicity include sudden death owing to psychosis, rhabdomyolysis, cardiac arrhythmia, stroke, and seizures, which are comparable to those of cocaine (Weis et al., 2011). Abuse of amphetamines can lead to neuropathologic changes such as cerebral infarcts and hemorrhages, which are frequently brought on by emboli connected to cardiac arrhythmias or myocardial infarction (Delaney and Estes, 1980). In addition, vasospasm or vasculitis resulting from the pharmacological effect of amphetamines may be contributing causes (Delaney and Estes, 1980; Weis et al., 2011). Weis et al. (Weis et al., 2011) reported a patient with amphetamine abuse-related acute demyelination that had not been previously described. The first report of amphetamine in the afflicted region of the brain strongly supports a link between amphetamine usage and demyelination; nevertheless, the pathogenetic processes of demyelination are yet unknown. In this study, BP and MF enrichment analyses were also performed on DElncRNAs in the ceRNA network. The findings suggested that the majority of DElncRNAs’ activities were focused on neuron recognition, synaptic function, ion channel activity, and transcriptional and post-transcriptional regulation.

We identified five hub genes in the ceRNA network, including FOS, GJA1, NTRK2, CTNND1, and SP3. FOS is a transcription factor linked to neuronal activity. Leucine zipper proteins from the FOS gene family can dimerize with JUN family proteins to form the transcription factor complex JUN (Shang et al., 2020). FOS may have essential regulatory roles in modulating gene expression in the pathophysiology of MS (Liu et al., 2013; Kotelnikova et al., 2015), as well as in the mechanism of action of IFNbeta1a and fingolimod (Howe et al., 2006; Kauffman et al., 2009; Yester et al., 2015; Anastasiadou and Knöll, 2016; Groves et al., 2018; Tran et al., 2020). In agreement with our findings, Yu et al. (1991) discovered a two-fold increase in FOS RNA in MS white matter relative to control tissue. The gap junctions (GJs), which form intercellular communication channels between two apposing cells or hemichannels with the extracellular environment, serve essential roles in maintaining small molecule homeostasis (Basu and Sarma, 2018). The GJs of the CNS are essential for maintaining myelin sheath and neuronal function. Connexin (Cx) proteins are the structural components of GJs (Basu and Sarma, 2018). The expression of GJA1 (Cx43) is essential for maintaining K+ buffering and nutritional homeostasis in oligodendrocytes, CNS myelin, and oligodendrocyte function, according to previous research (Basu and Sarma, 2018; Une et al., 2021). Furthermore, ablation of Cx43 in brain gray matter astroglia reduces the severity of EAE, an animal model of MS, by promoting an anti-inflammatory phenotype in astroglia and suppressing pro-inflammatory activation of spinal microglia partly through decreased cerebrospinal fluid pro-inflammatory cytokine/chemokine levels (Une et al., 2021). Therefore, brain astroglial Cx43 may be a potential MS therapeutic target (Une et al., 2021). GJA1 expression was dramatically elevated in the spinal cord periplaques of progressive MS patients and in EAE mice, according to previous research. (Lieury et al., 2014; Jin et al., 2017; Nataf et al., 2017). In addition, Nataf et al. (2017) found an astrocytosis-related co-expression module with the astrocyte gene GJA1 as its central hub. These findings are consistent with our result, indicating that the GJA1 gene is involved in the pathophysiology of MS. NTRK2 (also known as TrkB) is a membrane-bound receptor that belongs to the neurotrophic receptor kinase (NTRK) family. When neurotrophic proteins bind, members of the NTRK family and MAPK pathways are phosphorylated, resulting in cell differentiation via NTRK2 (Li et al., 2022). The Brain-Derived Neurotrophic Factor (BDNF) is one of the most important players in the biological pathways of brain plasticity. BDNF belongs to the neurotrophin family and binds to the NTRK2. The binding causes the MAPK pathway (Chao, 2003) to be activated, regulating synaptic plasticity and repair. On the other hand, fibrinogen (Fg)-containing plaques have been linked to memory loss in inflammatory neurodegenerative diseases such as Alzheimer’s, MS, stroke, and traumatic brain injury (Clark et al., 2018). Data indicate that Fg interacts with astrocytes, resulting in the upregulation of intercellular adhesion molecule 1 (ICAM-1) and TrkB and the phosphorylation of TrkB, and hence the activation of astrocytes (Clark et al., 2018). Since it is known that activated astrocytes express TrkB, which induces neuronal degeneration (Colombo et al., 2012), the interaction of Fg with astrocytes and the subsequent activation of TrkB may be a process involved in memory loss (Clark et al., 2018). Consistent with our findings, BDNF and its receptor truncated trkB tyrosine kinase receptor were detected in lesions of MS patients (Stadelmann et al., 2002) as well as EAE (Stadelmann et al., 2002; De Santi et al., 2009; Colombo et al., 2012). The CTNND1 gene encodes the cellular adhesion protein p120-catenin (p120), which is essential for myelinating Schwann cells, cell-cell contacts, and proper myelin sheath development. (Perrin-Tricaud et al., 2007). By altering Rho GTPase activity, CTNND1 is implicated in the regulation of cadherin-mediated adhesion and dynamic regulation of the actin cytoskeleton (Anastasiadis et al., 2000; Noren et al., 2000; Grosheva et al., 2001). CTNND1 was found among myelin basic protein (MBP) partners by Smirnova et al. (2021). Intrinsically aberrant MBP is one of the most important autoantigens in autoimmune neurodegeneration and MS (Benjamins and Morell, 1978; Garbay et al., 1988). These findings support our findings and show that the CTNND1 gene contributes to the pathophysiology of MS. Human transcription factor SP3, a member of the Sp1 family that binds GC/T box elements, has the ability to activate or suppress the expression of a number of immune system-related genes (Grekova et al., 2002). Multiple lines of evidence indicate a reduction in SP3 expression in PBMCs from patients with MS (Grekova et al., 1996; Grekova et al., 2000; Grekova et al., 2002; Lin et al., 2008). This discrepancy between these results and our findings may be explained as follows. The profile of gene expression in PBMCs does not reflect CNS components (Comabella and Martin, 2007; Dutta and Trapp, 2012). In addition, the gene expression profile derived from blood may potentially indicate the influence of variables unrelated to the illness process (Goertsches and Zettl, 2007).

It should be noted that a number of technical factors, including various methodologies, patient characteristics, a sample preparation, data analysis, and platforms, could have an impact on gene expression profiles. Additionally, small sample size may compromise statistical power. Of course, additional experimental studies and comparisons to reanalysis of modified microarray gene expression or other robust methodologies (e.g., WGCNA and ARACNE) are needed to confirm our findings.

Conclusion

Demyelination of the periplaque is likely to be the first stage of the lesion expansion process, making it a potential target for remyelinating therapies. Our research suggests potential ceRNA-mediated molecular mechanisms for this expansion involving ceRNA. Although the functions of the ceRNA axes need to be studied, this investigation provides potential research targets for examining molecular pathways that may be important for the pathogenesis of MS.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE52139).

Author contributions

HS, MR, and SA wrote the draft and revised it. MT and BH designed and supervised the study. MA, SK, GK, SH, and MS-B collected the data and performed the bio-informatic analysis. All the authors read and approved the submitted version.

Acknowledgments

We would like to thank the Clinical Research Development Unit of Tabriz Valiasr Hospital, Tabriz University of Medical Sciences, Tabriz, Iran for their assistance in this research.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1011350/full#supplementary-material

References

Ala, U. (2020). Competing endogenous RNAs, non-coding RNAs and diseases: An intertwined story. Cells 9 (7), 1574. doi:10.3390/cells9071574

PubMed Abstract | CrossRef Full Text | Google Scholar

Amato, M. P., Derfuss, T., Hemmer, B., Liblau, R., Montalban, X., Soelberg Sørensen, P., et al. (2018). Environmental modifiable risk factors for multiple sclerosis: Report from the 2016 ECTRIMS focused workshop. Mult. Scler. 24 (5), 590–603. doi:10.1177/1352458516686847

PubMed Abstract | CrossRef Full Text | Google Scholar

Anastasiadis, P. Z., Moon, S. Y., Thoreson, M. A., Mariner, D. J., Crawford, H. C., Zheng, Y., et al. (2000). Inhibition of RhoA by p120 catenin. Nat. Cell. Biol. 2 (9), 637–644. doi:10.1038/35023588

PubMed Abstract | CrossRef Full Text | Google Scholar

Anastasiadou, S., and Knöll, B. (2016). The multiple sclerosis drug fingolimod (FTY720) stimulates neuronal gene expression, axonal growth and regeneration. Exp. Neurol. 279, 243–260. doi:10.1016/j.expneurol.2016.03.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Basu, R., and Sarma, J. D. A. S. (2018). Connexin 43/47 channels are important for astrocyte/oligodendrocyte cross-talk in myelination and demyelination. J. Biosci. 43 (5), 1055–1068. doi:10.1007/s12038-018-9811-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamins, J. A., and Morell, P. (1978). Proteins of myelin and their metabolism. Neurochem. Res. 3 (2), 137–174. doi:10.1007/BF00964057

PubMed Abstract | CrossRef Full Text | Google Scholar

Bian, Z., Lei, W., Li, Q., Xue, W., Gao, Y., Zeng, Y., et al. (2020). Gm15575 functions as a ceRNA to up-regulate CCL7 expression through sponging miR-686 in Th17 cells. Mol. Immunol. 125, 32–42. doi:10.1016/j.molimm.2020.06.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Braschi, B., Denny, P., Gray, K., Jones, T., Seal, R., Tweedie, S., et al. (2019). Genenames.org: The HGNC and VGNC resources in 2019. Nucleic Acids Res. 47 (D1), D786-D792–D92. doi:10.1093/nar/gky930

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, Y., and Wan, J. (2018). Competing endogenous RNA regulations in neurodegenerative disorders: Current challenges and emerging insights. Front. Mol. Neurosci. 11, 370. doi:10.3389/fnmol.2018.00370

PubMed Abstract | CrossRef Full Text | Google Scholar

Caputo, V., Ciolfi, A., Macri, S., and Pizzuti, A. (2015). The emerging role of MicroRNA in schizophrenia. CNS Neurol. Disord. Drug Targets 14 (2), 208–221. doi:10.2174/1871527314666150116124253

PubMed Abstract | CrossRef Full Text | Google Scholar

Chao, M. V. (2003). Neurotrophins and their receptors: A convergence point for many signalling pathways. Nat. Rev. Neurosci. 4 (4), 299–309. doi:10.1038/nrn1078

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, E. Y., Tan, C. M., Kou, Y., Duan, Q., Wang, Z., Meirelles, G. V., et al. (2013). Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinforma. 14 (1), 128. doi:10.1186/1471-2105-14-128

PubMed Abstract | CrossRef Full Text | Google Scholar

Chin, C-H., Chen, S-H., Wu, H-H., Ho, C-W., Ko, M-T., and Lin, C-Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst. Biol. 8 (4), S11–S. doi:10.1186/1752-0509-8-S4-S11

PubMed Abstract | CrossRef Full Text | Google Scholar

Christiansen, C. F., Christensen, S., Farkas, D. K., Miret, M., Sørensen, H. T., and Pedersen, L. (2010). Risk of arterial cardiovascular diseases in patients with multiple sclerosis: A population-based cohort study. Neuroepidemiology 35 (4), 267–274. doi:10.1159/000320245

PubMed Abstract | CrossRef Full Text | Google Scholar

Clark, V. D., Layson, A., Charkviani, M., Muradashvili, N., and Lominadze, D. (2018). Hyperfibrinogenemia-mediated astrocyte activation. Brain Res. 1699, 158–165. doi:10.1016/j.brainres.2018.08.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Colombo, E., Cordiglieri, C., Melli, G., Newcombe, J., Krumbholz, M., Parada, L. F., et al. (2012). Stimulation of the neurotrophin receptor TrkB on astrocytes drives nitric oxide production and neurodegeneration. J. Exp. Med. 209 (3), 521–535. doi:10.1084/jem.20110698

PubMed Abstract | CrossRef Full Text | Google Scholar

Comabella, M., and Martin, R. (2007). Genomics in multiple sclerosis--current state and future directions. J. Neuroimmunol. 187 (1-2), 1–8. doi:10.1016/j.jneuroim.2007.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Crea, F., and Libby, P. (2017). Acute coronary syndromes: The way forward from mechanisms to precision treatment. Circulation 136 (12), 1155–1166. doi:10.1161/CIRCULATIONAHA.117.029870

PubMed Abstract | CrossRef Full Text | Google Scholar

Dal-Bianco, A., Grabner, G., Kronnerwetter, C., Weber, M., Höftberger, R., Berger, T., et al. (2017). Slow expansion of multiple sclerosis iron rim lesions: Pathology and 7 T magnetic resonance imaging. Acta Neuropathol. 133 (1), 25–42. doi:10.1007/s00401-016-1636-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Dastmalchi, R., Ghafouri-Fard, S., Omrani, M. D., Mazdeh, M., Sayad, A., and Taheri, M. (2018). Dysregulation of long non-coding RNA profile in peripheral blood of multiple sclerosis patients. Mult. Scler. Relat. Disord. 25, 219–226. doi:10.1016/j.msard.2018.07.044

PubMed Abstract | CrossRef Full Text | Google Scholar

De Santi, L., Annunziata, P., Sessa, E., and Bramanti, P. (2009). Brain-derived neurotrophic factor and TrkB receptor in experimental autoimmune encephalomyelitis and multiple sclerosis. J. Neurol. Sci. 287 (1-2), 17–26. doi:10.1016/j.jns.2009.08.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Delaney, P., and Estes, M. (1980). Intracranial hemorrhage with amphetamine abuse. Neurology 30 (10), 1125–1128. doi:10.1212/wnl.30.10.1125

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, Y., Li, T., Yan, X., Cui, M., Wang, C., Wang, S., et al. (2021). Identification of hub lncRNA ceRNAs in multiple sclerosis based on ceRNA mechanisms. Mol. Genet. Genomics. 296 (2), 423–435. doi:10.1007/s00438-020-01750-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, C., Liu, Y., Li, Y., Chen, H., Liu, X., Chen, X., et al. (2018). Sulfasalazine alters microglia phenotype by competing endogenous RNA effect of miR-136-5p and long non-coding RNA HOTAIR in cuprizone-induced demyelination. Biochem. Pharmacol. 155, 110–123. doi:10.1016/j.bcp.2018.06.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Dutta, R., and Trapp, B. D. (2012). Gene expression profiling in multiple sclerosis brain. Neurobiol. Dis. 45 (1), 108–114. doi:10.1016/j.nbd.2010.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Enbom, M. (2001). Multiple sclerosis and Kaposi's sarcoma--chronic diseases associated with new human herpesviruses? Scand. J. Infect. Dis. 33 (9), 648–658. doi:10.1080/00365540110027015

PubMed Abstract | CrossRef Full Text | Google Scholar

Fenoglio, C., Oldoni, E., Serpente, M., De Riz, M. A., Arcaro, M., D'Anca, M., et al. (2018). LncRNAs expression profile in peripheral blood mononuclear cells from multiple sclerosis patients. J. Neuroimmunol. 324, 129–135. doi:10.1016/j.jneuroim.2018.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Frischer, J. M., Weigand, S. D., Guo, Y., Kale, N., Parisi, J. E., Pirko, I., et al. (2015). Clinical and pathological insights into the dynamic nature of the white matter multiple sclerosis plaque. Ann. Neurol. 78 (5), 710–721. doi:10.1002/ana.24497

PubMed Abstract | CrossRef Full Text | Google Scholar

Garbay, B., Fournier, M., Sallafranque, M. L., Muller, S., Boiron, F., Heape, A., et al. (1988). Po, MBP, histone, and DNA levels in sciatic nerve. Postnatal accumulation studies in normal and trembler mice. Neurochem. Pathol. 8 (2), 91–107. doi:10.1007/BF03160138

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghafouri-Fard, S., and Taheri, M. (2020). A comprehensive review of non-coding RNAs functions in multiple sclerosis. Eur. J. Pharmacol. 879, 173127. doi:10.1016/j.ejphar.2020.173127

PubMed Abstract | CrossRef Full Text | Google Scholar

Gloss, B. S., and Dinger, M. E. (2016). The specificity of long noncoding RNA expression. Biochim. Biophys. Acta 1859 (1), 16–22. doi:10.1016/j.bbagrm.2015.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Goertsches, R., and Zettl, U. K. (2007). MS therapy research applying genome-wide RNA profiling of peripheral blood. Int. MS J. 14 (3), 98–107.

PubMed Abstract | Google Scholar

Grekova, M. C., Robinson, E. D., Faerber, M. A., Katz, P., McFarland, H. F., and Richert, J. R. (1996). Deficient expression in multiple sclerosis of the inhibitory transcription factor Sp3 in mononuclear blood cells. Ann. Neurol. 40 (1), 108–112. doi:10.1002/ana.410400117

PubMed Abstract | CrossRef Full Text | Google Scholar

Grekova, M. C., Salerno, K., Mikkilineni, R., and Richert, J. R. (2002). Sp3 expression in immune cells: A quantitative study. Lab. Investig. 82 (9), 1131–1138. doi:10.1097/01.lab.0000029149.38881.84

PubMed Abstract | CrossRef Full Text | Google Scholar

Grekova, M. C., Scherer, S. W., Trabb, J., and Richert, J. R. (2000). Localization of the human SP3 gene to chromosome 7p14-p15.2. The lack of expression in multiple sclerosis does not reflect abnormal gene organization. J. Neuroimmunol. 106 (1-2), 214–219. doi:10.1016/s0165-5728(99)00264-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Grosheva, I., Shtutman, M., Elbaum, M., and Bershadsky, A. D. (2001). p120 catenin affects cell motility via modulation of activity of Rho-family GTPases: a link between cell-cell contact formation and regulation of cell locomotion. J. Cell. Sci. 114 (4), 695–707. doi:10.1242/jcs.114.4.695

PubMed Abstract | CrossRef Full Text | Google Scholar

Groves, A., Kihara, Y., Jonnalagadda, D., Rivera, R., Kennedy, G., Mayford, M., et al. (2018). A functionally defined in vivo astrocyte population identified by c-fos activation in a mouse model of multiple sclerosis modulated by S1P signaling: Immediate-early astrocytes (ieAstrocytes). eNeuro 5, 1–13.

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, Y., He, M., Fu, Y., Zhao, C., Xiong, S., and Xu, X. (2022). Identification of novel key genes and pathways in multiple sclerosis based on weighted gene coexpression network analysis and long noncoding RNA-associated competing endogenous RNA network. Oxid. Med. Cell. Longev. 2022, 9328160. doi:10.1155/2022/9328160

PubMed Abstract | CrossRef Full Text | Google Scholar

Howe, C. L., Mayoral, S., and Rodriguez, M. (2006). Activated microglia stimulate transcriptional changes in primary oligodendrocytes via IL-1beta. Neurobiol. Dis. 23 (3), 731–739. doi:10.1016/j.nbd.2006.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, H. Y., Lin, Y. C., Li, J., Huang, K. Y., Shrestha, S., Hong, H. C., et al. (2020). miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 48 (D1), D148-D154–D54. doi:10.1093/nar/gkz896

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, K., Wu, L., Gao, Y., Li, Q., Wu, H., Liu, X., et al. (2022). Transcriptome sequencing data reveal LncRNA-miRNA-mRNA regulatory network in calcified aortic valve disease. Front. Cardiovasc. Med. 9, 886995. doi:10.3389/fcvm.2022.886995

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Z., Shi, J., Gao, Y., Cui, C., Zhang, S., Li, J., et al. (2019). HMDD v3.0: A database for experimentally supported human microRNA-disease associations. Nucleic Acids Res. 47 (D1), D1013-D1017–D7. doi:10.1093/nar/gky1010

PubMed Abstract | CrossRef Full Text | Google Scholar

Huber, W., Carey, V. J., Gentleman, R., Anders, S., Carlson, M., Carvalho, B. S., et al. (2015). Orchestrating high-throughput genomic analysis with Bioconductor. Nat. Methods 12 (2), 115–121. doi:10.1038/nmeth.3252

PubMed Abstract | CrossRef Full Text | Google Scholar

Irizarry, R. A., Hobbs, B., Collin, F., Beazer-Barclay, Y. D., Antonellis, K. J., Scherf, U., et al. (2003). Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostat. Oxf. Engl. 4 (2), 249–264. doi:10.1093/biostatistics/4.2.249

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Li, W., Yan, Y., Yao, X., Gu, W., and Zhang, H. (2020). LINC01094 triggers radio-resistance in clear cell renal cell carcinoma via miR-577/CHEK2/FOXM1 axis. Cancer Cell. Int. 20, 274. doi:10.1186/s12935-020-01306-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, Z., Ding, Y., Xue, R., Jia, Z., Huang, Z., Ding, Y., et al. (2017). Involvement of interstitial cells of Cajal in bladder dysfunction in mice with experimental autoimmune encephalomyelitis. Int. Urol. Nephrol. 49 (8), 1353–1359. doi:10.1007/s11255-017-1597-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Karagkouni, D., Paraskevopoulou, M. D., Tastsoglou, S., Skoufos, G., Karavangeli, A., Pierros, V., et al. (2020). DIANA-LncBase v3: Indexing experimentally supported miRNA targets on non-coding transcripts. Nucleic Acids Res. 48 (D1), D101-D110–D10. doi:10.1093/nar/gkz1036

PubMed Abstract | CrossRef Full Text | Google Scholar

Karimi, E., Azari, H., Tahmasebi, A., Nikpoor, A. R., Negahi, A. A., Sanadgol, N., et al. (2022). LncRNA-miRNA network analysis across the Th17 cell line reveals biomarker potency of lncRNA NEAT1 and KCNQ1OT1 in multiple sclerosis. J. Cell. Mol. Med. 26 (8), 2351–2362. doi:10.1111/jcmm.17256

PubMed Abstract | CrossRef Full Text | Google Scholar

Kauffman, M. A., Yankilevich, P., Barrero, P., Bello, R., Marangunich, L., Vidal, A., et al. (2009). Whole genome analysis of the action of interferon-beta. Int. J. Clin. Pharmacol. Ther. 47 (5), 328–357. doi:10.5414/cpp47328

PubMed Abstract | CrossRef Full Text | Google Scholar

Klistorner, S., Barnett, M. H., Yiannikas, C., Barton, J., Parratt, J., You, Y., et al. (2021). Expansion of chronic lesions is linked to disease progression in relapsing-remitting multiple sclerosis patients. Mult. Scler. 27 (10), 1533–1542. doi:10.1177/1352458520974357

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotelnikova, E., Bernardo-Faura, M., Silberberg, G., Kiani, N. A., Messinis, D., Melas, I. N., et al. (2015). Signaling networks in MS: A systems-based approach to developing new pharmacological therapies. Mult. Scler. 21 (2), 138–146. doi:10.1177/1352458514543339

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuleshov, M. V., Jones, M. R., Rouillard, A. D., Fernandez, N. F., Duan, Q., Wang, Z., et al. (2016). Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44 (W1), W90–W97. doi:10.1093/nar/gkw377

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Li, W., Zhang, L., He, J., Tang, L., Li, Z., et al. (2022). Comprehensive network analysis identified SIRT7, NTRK2, and CHI3L1 as new potential markers for intervertebral disc degeneration. J. Oncol. 2022, 4407541. doi:10.1155/2022/4407541

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Yin, B., Chen, M., Peng, J., Mu, X., Deng, Z., et al. (2021). Downregulation of the lncRNA ASB16-AS1 decreases LARP1 expression and promotes clear cell renal cell carcinoma progression via miR-185-5p/miR-214-3p. Front. Oncol. 10, 617105. doi:10.3389/fonc.2020.617105

PubMed Abstract | CrossRef Full Text | Google Scholar

Lieury, A., Chanal, M., Androdias, G., Reynolds, R., Cavagna, S., Giraudon, P., et al. (2014). Tissue remodeling in periplaque regions of multiple sclerosis spinal cord lesions. Glia 62 (10), 1645–1658. doi:10.1002/glia.22705

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, A. Y., Yang, Q. D., and Murong, S. X. (2008). Deficient mRNA expression of specific protein 3 gene in peripheral blood mononuclear cells from patients with multiple sclerosis. Zhonghua Yi Xue Yi Chuan Xue Za Zhi 25 (2), 187–189.

PubMed Abstract | Google Scholar

Liu, L., Xu, Q., Xiong, Y., Deng, H., and Zhou, J. (2022). LncRNA LINC01094 contributes to glioma progression by modulating miR-224-5p/CHSY1 axis. Hum. Cell. 35 (1), 214–225. doi:10.1007/s13577-021-00637-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, M., Hou, X., Zhang, P., Hao, Y., Yang, Y., Wu, X., et al. (2013). Microarray gene expression profiling analysis combined with bioinformatics in multiple sclerosis. Mol. Biol. Rep. 40 (5), 3731–3737. doi:10.1007/s11033-012-2449-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Zhuang, R., Feng, S., Bai, X., Jia, Z., Kapora, E., et al. (2020). Long non-coding RNA ASB16-AS1 enhances cell proliferation, migration and invasion via functioning as a ceRNA through miR-1305/Wnt/β-catenin axis in cervical cancer. Biomed. Pharmacother. 125, 109965. doi:10.1016/j.biopha.2020.109965

PubMed Abstract | CrossRef Full Text | Google Scholar

Manouchehrinia, A., Tanasescu, R., Tench, C. R., and Constantinescu, C. S. (2016). Mortality in multiple sclerosis: meta-analysis of standardised mortality ratios. J. Neurol. Neurosurg. Psychiatry 87 (3), 324–331. doi:10.1136/jnnp-2015-310361

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreno-García, L., López-Royo, T., Calvo, A. C., Toivonen, J. M., de la Torre, M., Moreno-Martínez, L., et al. (2020). Competing endogenous RNA networks as biomarkers in neurodegenerative diseases. Int. J. Mol. Sci. 21 (24), 9582. doi:10.3390/ijms21249582

CrossRef Full Text | Google Scholar

Nataf, S., Barritault, M., and Pays, L. (2017). A unique TGFB1-driven genomic program links astrocytosis, low-grade inflammation and partial demyelination in spinal cord periplaques from progressive multiple sclerosis patients. Int. J. Mol. Sci. 18 (10), 2097. doi:10.3390/ijms18102097

CrossRef Full Text | Google Scholar

Noren, N. K., Liu, B. P., Burridge, K., and Kreft, B. (2000). p120 catenin regulates the actin cytoskeleton via Rho family GTPases. J. Cell. Biol. 150 (3), 567–580. doi:10.1083/jcb.150.3.567

PubMed Abstract | CrossRef Full Text | Google Scholar

Perrin-Tricaud, C., Rutishauser, U., and Tricaud, N. (2007). P120 catenin is required for thickening of Schwann cell myelin. Mol. Cell. Neurosci. 35 (1), 120–129. doi:10.1016/j.mcn.2007.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Perron, U., Provero, P., and Molineris, I. (2017). In silico prediction of lncRNA function using tissue specific and evolutionary conserved expression. BMC Bioinforma. 18 (5), 144. doi:10.1186/s12859-017-1535-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Prineas, J. W., Kwon, E. E., Cho, E. S., Sharer, L. R., Barnett, M. H., Oleszak, E. L., et al. (2001). Immunopathology of secondary-progressive multiple sclerosis. Ann. Neurol. 50 (5), 646–657. doi:10.1002/ana.1255

PubMed Abstract | CrossRef Full Text | Google Scholar

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

Romero, L., and AgiMicroRna, P. (2020). Processing and differential expression analysis of agilent microRNA chips.

Google Scholar

Sabaie, H., Mazaheri Moghaddam, M., Mazaheri Moghaddam, M., Amirinejad, N., Asadi, M. R., Daneshmandpour, Y., et al. (2021). Long non-coding RNA-associated competing endogenous RNA axes in the olfactory epithelium in schizophrenia: A bioinformatics analysis. Sci. Rep. 11 (1), 24497. doi:10.1038/s41598-021-04326-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabaie, H., Moghaddam, M. M., Moghaddam, M. M., Ahangar, N. K., Asadi, M. R., Hussen, B. M., et al. (2021). Bioinformatics analysis of long non-coding RNA-associated competing endogenous RNA network in schizophrenia. Sci. Rep. 11 (1), 24413. doi:10.1038/s41598-021-03993-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Sabaie, H., Salkhordeh, Z., Asadi, M. R., Ghafouri-Fard, S., Amirinejad, N., Askarinejad Behzadi, M., et al. (2021). Long non-coding RNA- associated competing endogenous RNA axes in T-cells in multiple sclerosis. Front. Immunol. 12, 770679. doi:10.3389/fimmu.2021.770679

PubMed Abstract | CrossRef Full Text | Google Scholar

Salmena, L., Poliseno, L., Tay, Y., Kats, L., and Pandolfi, P. P. (2011). A ceRNA hypothesis: The rosetta stone of a hidden RNA language? Cell. 146 (3), 353–358. doi:10.1016/j.cell.2011.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Santoro, M., Nociti, V., Lucchini, M., De Fino, C., Losavio, F. A., and Mirabella, M. (2016). Expression profile of long non-coding RNAs in serum of patients with multiple sclerosis. J. Mol. Neurosci. 59 (1), 18–23. doi:10.1007/s12031-016-0741-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Santoro, M., Nociti, V., Lucchini, M., Loiodice, M., Centofanti, F., Botta, A., et al. (2020). A pilot study of lncRNAs expression profile in serum of progressive multiple sclerosis patients. Eur. Rev. Med. Pharmacol. Sci. 24 (6), 3267–3273. doi:10.26355/eurrev_202003_20694

PubMed Abstract | CrossRef Full Text | Google Scholar

Senousy, M. A., Shaker, O. G., Sayed, N. H., Fathy, N., and Kortam, M. A. (2020). LncRNA GAS5 and miR-137 polymorphisms and expression are associated with multiple sclerosis risk: Mechanistic insights and potential clinical impact. ACS Chem. Neurosci. 11 (11), 1651–1660. doi:10.1021/acschemneuro.0c00150

PubMed Abstract | CrossRef Full Text | Google Scholar

Shang, Z., Sun, W., Zhang, M., Xu, L., Jia, X., Zhang, R., et al. (2020). Identification of key genes associated with multiple sclerosis based on gene expression data from peripheral blood mononuclear cells. PeerJ 8, e8357. doi:10.7717/peerj.8357

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheng, W., Sheng, K., Zhao, Y., Li, H., Zhou, J., Yao, H., et al. (2015). Identifying the biomarkers of multiple sclerosis based on non-coding RNA signature. Eur. Rev. Med. Pharmacol. Sci. 19 (19), 3635–3642.

PubMed Abstract | Google Scholar

Smirnova, E. V., Rakitina, T. V., Ziganshin, R. H., Arapidi, G. P., Saratov, G. A., Kudriaeva, A. A., et al. (2021). Comprehensive atlas of the myelin basic protein interaction landscape. Biomolecules 11 (11), 1628. doi:10.3390/biom11111628

PubMed Abstract | CrossRef Full Text | Google Scholar

Stadelmann, C., Kerschensteiner, M., Misgeld, T., Brück, W., Hohlfeld, R., and Lassmann, H. (2002). BDNF and gp145trkB in multiple sclerosis brain lesions: Neuroprotective interactions between immune and neuronal cells? Brain 125 (1), 75–85. doi:10.1093/brain/awf015

PubMed Abstract | CrossRef Full Text | Google Scholar

Stefanou, M. I., Krumbholz, M., Ziemann, U., and Kowarik, M. C. (2019). Human immunodeficiency virus and multiple sclerosis: A review of the literature. Neurol. Res. Pract. 1, 24. doi:10.1186/s42466-019-0030-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47 (D1), D607-D613–d13. doi:10.1093/nar/gky1131

PubMed Abstract | CrossRef Full Text | Google Scholar

ten Bosch, G. J. A., Bolk, J., ‘t Hart, B. A., and Laman, J. D. (2021). Multiple sclerosis is linked to MAPKERK overactivity in microglia. J. Mol. Med. 99 (8), 1033–1042. doi:10.1007/s00109-021-02080-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, C., Heng, B., Teo, J. D., Humphrey, S. J., Qi, Y., Couttas, T. A., et al. (2020). Sphingosine 1-phosphate but not Fingolimod protects neurons against excitotoxic cell death by inducing neurotrophic gene expression in astrocytes. J. Neurochem. 153 (2), 173–188. doi:10.1111/jnc.14917

PubMed Abstract | CrossRef Full Text | Google Scholar

Tully, T., Barkley, A., and Silber, E. (2015). Kaposi sarcoma in a patient with relapsing-remitting multiple sclerosis receiving fingolimod. Neurology 84 (19), 1999–2001. doi:10.1212/WNL.0000000000001556

PubMed Abstract | CrossRef Full Text | Google Scholar

Une, H., Yamasaki, R., Nagata, S., Yamaguchi, H., Nakamuta, Y., Indiasari, U. C., et al. (2021). Brain gray matter astroglia-specific connexin 43 ablation attenuates spinal cord inflammatory demyelination. J. Neuroinflammation 18 (1), 126. doi:10.1186/s12974-021-02176-1

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Valk, P., and De Groot, C. J. (2000). Staging of multiple sclerosis (MS) lesions: Pathology of the time frame of MS. Neuropathol. Appl. Neurobiol. 26 (1), 2–10. doi:10.1046/j.1365-2990.2000.00217.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, S., and Brew, B. (2016). Kaposi sarcoma in a fingolimod-treated patient with multiple sclerosis. J. Clin. Neurosci. 31, 217–218. doi:10.1016/j.jocn.2016.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, T., Xu, S., Liu, L., Li, S., Zhang, H., Lu, X., et al. (2022). Integrated analysis of differentially expressed genes and a ceRNA network to identify hub lncRNAs and potential drugs for multiple sclerosis. Am. J. Transl. Res. 14 (2), 772–787.

PubMed Abstract | Google Scholar

Weis, S., Sachs, H., and Büttner, A. (2011). Acute demyelination in a person with amphetamine abuse. Case Rep. Pathol. 2011, 514613. doi:10.1155/2011/514613

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Wu, Y., Zhang, B., and Ni, B. (2018). Noncoding RNAs in multiple sclerosis. Clin. Epigenetics 10 (1), 149. doi:10.1186/s13148-018-0586-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Yester, J. W., Bryan, L., Waters, M. R., Mierzenski, B., Biswas, D. D., Gupta, A. S., et al. (2015). Sphingosine-1-phosphate inhibits IL-1-induced expression of C-C motif ligand 5 via c-Fos-dependent suppression of IFN-β amplification loop. Faseb J. 29 (12), 4853–4865. doi:10.1096/fj.15-275180

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, T. L., Matsuda, T., and Cepko, C. L. (2005). The noncoding RNA taurine upregulated gene 1 is required for differentiation of the murine retina. Curr. Biol. 15 (6), 501–512. doi:10.1016/j.cub.2005.02.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, A-M., Choi, Y. H., and Tu, M-J. (2020). RNA drugs and RNA targets for small molecules: Principles, progress, and challenges. Pharmacol. Rev. 72 (4), 862–898. doi:10.1124/pr.120.019554

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, J. S., Hayashi, T., Seboun, E., Sklar, R. M., Doolittle, T. H., and Hauser, S. L. (1991). Fos RNA accumulation in multiple sclerosis white matter tissue. J. Neurol. Sci. 103 (2), 209–215. doi:10.1016/0022-510x(91)90166-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, P., Jing, L., Zhao, X., Zhu, H., and Teng, J. (2019). Down-regulation of taurine-up-regulated gene 1 attenuates inflammation by sponging miR-9-5p via targeting NF-κB1/p50 in multiple sclerosis. Life Sci. 233, 116731. doi:10.1016/j.lfs.2019.116731

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Zhou, H., Liu, J., and Mao, J. (2019). Long noncoding RNA ASB16-AS1 promotes proliferation, migration, and invasion in glioma cells. Biomed. Res. Int. 2019, 5437531. doi:10.1155/2019/5437531

PubMed Abstract | CrossRef Full Text | Google Scholar