Next Article in Journal
Development of Novel Residual-Dense-Attention (RDA) U-Net Network Architecture for Hepatocellular Carcinoma Segmentation
Previous Article in Journal
Detecting Endotracheal Tube and Carina on Portable Supine Chest Radiographs Using One-Stage Detector with a Coarse-to-Fine Attention
Previous Article in Special Issue
MiRNAs in Lung Cancer: Diagnostic, Prognostic, and Therapeutic Potential
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated Approaches to Identify miRNA Biomarkers Associated with Cognitive Dysfunction in Multiple Sclerosis Using Text Mining, Gene Expression, Pathways, and GWAS

1
R&D Division, Eriks-Precision Components India Pvt Ltd., Mohali 160055, Punjab, India
2
Regenerative Biology, Morgridge Institute for Research, Madison, WI 53715, USA
*
Author to whom correspondence should be addressed.
Diagnostics 2022, 12(8), 1914; https://doi.org/10.3390/diagnostics12081914
Submission received: 5 July 2022 / Revised: 3 August 2022 / Accepted: 3 August 2022 / Published: 8 August 2022
(This article belongs to the Special Issue MicroRNA in Diagnosis and Prognosis)

Abstract

:
Multiple sclerosis (MS), a chronic autoimmune disorder, affects the central nervous system of many young adults. More than half of MS patients develop cognition problems. Although several genomic and transcriptomic studies are currently reported in MS cognitive impairment, a comprehensive repository dealing with all the experimental data is still underdeveloped. In this study, we combined text mining, gene regulation, pathway analysis, and genome-wide association studies (GWAS) to identify miRNA biomarkers to explore the cognitive dysfunction in MS, and to understand the genomic etiology of the disease. We first identified the dysregulated miRNAs associated with MS and cognitive dysfunction using PubTator (text mining), HMDD (experimental associations), miR2Disease, and PhenomiR database (differentially expressed miRNAs). Our results suggest that miRNAs such as hsa-mir-148b-3p, hsa-mir-7b-5p, and hsa-mir-7a-5p are commonly associated with MS and cognitive dysfunction. Next, we retrieved GWAS signals from GWAS Catalog, and analyzed the enrichment analysis of association signals in genes/miRNAs and their association networks. Then, we identified susceptible genetic loci, rs17119 (chromosome 6; p = 1 × 10−10), rs1843938 (chromosome 7; p = 1 × 10−10), and rs11637611 (chromosome 15; p = 1.00 × 10−15), associated with significant genetic risk. Lastly, we conducted a pathway analysis for the susceptible genetic variants and identified novel risk pathways. The ECM receptor signaling pathway (p = 3.98 × 10−8) and PI3K/Akt signaling pathway (p = 5.98 × 10−5) were found to be associated with differentially expressed miRNA biomarkers.

1. Introduction

Multiple sclerosis (MS) is a long-lasting autoimmune disease of the central nervous system (brain and spinal cord) that leads to demyelination and neurodegeneration [1]. Nearly 2.5 million patients are suffering with MS worldwide, and individuals in the age group of 20 to 40 are the most affected [2]. MS is termed as a medical mystery because the cause of the disease remains elusive. A common assumption is that MS is caused by an abnormal immune response to ecological factors in individuals with a genetic predisposition [3]. A neurological manifestation of MS comprises motor and visual deficits, limb weakness, sensory problems, gait ataxia, uncontrolled bladder, and cognitive decline [4,5]. In addition to the known motor and sensory problems, the effect of cognitive dysfunction is very common in more than half of all the individuals with MS [6,7].
Apart from MS, cognitive dysfunction is a common condition in patients with neurodegenerative diseases. It causes focal impairment to brain structures mediating mild to moderate cognition, and it affects personality, behavior, and decision-making abilities [8,9]. Cognitive dysfunction in MS patients specifically includes declined information processing, memory, concentration, reasoning, and verbal fluency [10,11]. The pathological characteristics of MS include multifocal lesions, inflammation, loss of oligodendrocyte, various ranges of axonal injury, gliosis, and temporary dysfunction to permanent loss [12].
In our study, we aimed to identify biomarkers in MS cognitive impairment analysis using miRNA pathway associations involved in MS cognition. Biomarkers play an important role in the primary assessment and evaluation of a disease, which helps in early detection and treatment. They act as the indicators of changes in normal biological processes and pathological processes, helping to extend our knowledge about the underlying disease conditions [13]. Biomarkers identify neurological disease at an early stage and help to determine the disease condition. If the disease worsens or improves in the patient, the concentration of the biomarker might increase or decrease accordingly; hence, biomarkers act as a noninvasive method of detection for the patients [14].
Gaetani et al. investigated cerebrospinal fluid biomarkers and determined different neuropsychological patterns of cognitive deficits in MS [15]. The cerebral hypoperfusion in MS is associated with chronic hypoxia, focal lesion formation, diffuse axonal degeneration, cognitive dysfunction, and fatigue. MS patients with pathological vascular abnormalities may have reduced quality of life and worsened cerebral perfusion, which may contribute to greater disease progression. Patients with cerebrovascular diseases have poor cognitive performance in MS. These findings suggest that restoring the cerebral blood flow may act as a new therapeutic target in MS [16,17].
The identification of disease-specific biomarkers provides an opportunity to the patients to start the treatment early enough to delay the disease progression [11]. Examples of biomarkers in the progression of MS are (i) MxA expression defining the pathological effect of interferon-β (IFN-β) in vivo [18], (ii) TNF-related apoptosis-inducing ligand (TRAIL), a potential biomarker of IFN-β therapeutic efficacy [19], (iii) complement regulator factor H, a serum biomarker predicted to be effective in disease progression [20], and (iv) the genes CHI3L1 and NF-L, recently identified as biomarkers for mild cognitive impairment in early stages of MS at genomic levels [21]. MicroRNAs (miRNAs) play an important role in various cellular functions and disease-causing mechanisms [22,23]. miRNAs are known to play a critical role in neuron development and maintenance of the central nervous system, including neural differentiation, synaptic plasticity, and cognitive functions [24,25].
An MS cognition study performed by Liguori et al. demonstrated the expression of miRNA in association with cognitive impairment in pediatric MS and reported the associated genes, BST1, NTNG2, SPTB, and STAB1 [26]. Due to the complexity of the disease, various clinical phenotypes are observed in patients at different stages of MS, which led the researchers to analyze various biological and pathological pathways responding to the therapy at both genomic and transcriptomic levels [27,28]. Despite the occurrence of cognitive symptoms in MS, the role of miRNAs and their target mRNAs associated with cognitive decline and neuropsychiatric problems in MS disease is less reported. A study by Scapoli et al. highlighted the influence of vascular components in MS. The study combined the transcriptome analysis of the internal jugular vein (IJV) walls from MS patients and whole-exome sequencing (WES) [29]. Studies on the expression patterns of the IJV wall in MS through combined transcriptome–protein analysis highlighted the proteins of interest for MS pathophysiology [30].
GWAS and other mapping studies have identified 250 loci for autoimmune diseases, some of which are common to two or more diseases. Most of the variants are in the noncoding regions, especially in miRNA target sites [31]. miRNAs act by imperfect base paring with target messenger RNA (mRNA), leading to the negative regulation of mRNA. Single-nucleotide polymorphisms (SNPs) in mRNA sites complementary to miRNA are referred to as miRSNPs [32]. miRSNPs have functional consequences for autoimmune diseases including MS [31]. Statistical approaches and meta-analysis of neurological diseases A recent study showed that miRSNPs have an expression quantitative trait locus (eQTL) effect on genes reported in GWAS [33], and the number of miRSNPs associated with human diseases keep increasing [34,35]; A recent study integrated GWAS with pathways to report novel risk pathways and genes [36]. We believe that the existing analysis on MS cognitive dysfunction can reveal the pathways and genes with common risks that are associated with the cognitive dysfunction in MS. This study can clarify the disease mechanism and its genetic etiology.
This paper extracted the biomarkers for various cognitive function changes in MS patients reported by high-throughput experiments and other experimental techniques from PubMed abstracts using PubTator [37], a text-mining tool. Cognitive problems occurring due to alterations in normal functions of memory, attention and concentration, information processing, executive function, and verbal fluency that are associated with cognitive dysfunction, cognitive impairment, cognitive decline, etc. are reported here. MS cognition and potential disease associations were identified for a wide range of clinical use. While many MS cognition-associated biomarkers have been reported in recent years, this is the first attempt to present an integrated approach to study all the genes, miRNAs, and associated pathways.

2. Materials and Methods

2.1. Data Collection

We extracted the research articles from PubMed (https://www.ncbi.nlm.nih.gov/pubmed, accessed on 30 June 2019) using the search terms multiple sclerosis and cognition: “multiple sclerosis” [MeSH terms] OR (“multiple” [all fields] AND “sclerosis” [all fields]) OR “multiple sclerosis” [all fields]) AND (“cognitive” [all fields] OR “cognition” [all fields]) AND “humans” [MeSH terms]. We retrieved 4102 publications (June 2019) and exported the PubMed IDs of retrieved articles to a text file for further processing. The retrieved articles were used as input for text-mining analysis.

2.2. Gene Concept Finding Using PubTator

PubTator [37] supports five kinds of biological concepts: gene, disease, chemical, species, and mutation. We utilized the gene concept finding, a module within PubTator, for recognizing miRNA and gene mentions in PubMed articles. (https://www.ncbi.nlm.nih.gov/research/bionlp/APIs/usage/, accessed on 30 June 2019) [38]. The gene concept finding module contains information on genes (and gene IDs) reported in each PubMed article. Text-mining results from PubTator were retrieved using libwww-perl (LWP), a Perl script. The retrieved information was manually verified, and a nonredundant set of data was mapped to Entrez gene name. For each gene, the official gene symbol approved by the HUGO gene nomenclature committee (HGNC) [39] was retrieved from the NCBI database [40].

2.3. Predicting Gene–miRNA Targets

We used DIANA-TarBase v8.0, a database which provides high-quality manually curated experimentally validated miRNA–gene interactions [41], to determine the experimentally validated miRNA–gene interactions. The database includes information on the methodology, experimental conditions such as cell/tissue type, and study design. The results are enhanced with detailed metadata. The genes associated with MS cognition in human were used to determine the corresponding miRNA targets and the associated mRNAs. A single miRNA can bind to hundreds of target mRNAs, or a single mRNA can be targeted by multiple miRNAs. These miRNA–mRNA associations play an important role in the regulation of genome [42,43,44].

2.4. Identification and Analysis of Genetic Risk Variants

GWAS Catalog [45] was used to identify the risk variants associated with MS and cognitive dysfunction. For MS, we used the keyword “multiple sclerosis”. For cognitive dysfunction, we used two keywords, “cognitive impairment” (CI) and “cognitive decline” (CD). A list of 458 associated genetic variants were retrieved for MS, 150 associated genetic variants were retrieved for CI, and 33 associated genetic variants were retrieved for CD. We selected 77 genetic variants for MS and 25 genetic variants for CI with a p-value less than 1 × 10−5. For MS, the genetic variants were associated with MS, MS pleiotropy, MS neuromyelitis optica, and MS cholesterol measurements. For CI, the genetic variants were associated with cognitive decline rate in late mild cognitive impairment. For CD, no genetic variants with a p-value less than 1 × 10−5 were retrieved.

2.5. Function Analysis

HaploReg tool (v4.1) [46] was used to annotate the genetic variants, 77 for MS and 25 for CI, and to identify both coding and noncoding mutations. PhenoScanner (v2) [47] database, a collection of human genotype–phenotype associations, was used to perform eQTL analysis on the selected genetic variants. In addition, we used miRDSNP (mirdsnp.ccr.buffalo.edu), a database of disease-associated SNPs and microRNA target sites on 3’UTRs of human genes and MSDD (MiRNA SNP Disease Database, http://bio-bigdata.hrbmu.edu.cn/msdd/, accessed on 15 July 2019) [48] to determine miRNA.

2.6. Pathway Analysis

DIANA-mirPath (mirPath v.3) [49] is a webserver for miRNA pathway analysis that utilizes KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways for analysis. This webserver utilizes experimentally validated miRNA interactions derived from DIANA-Tarbase database [41]. In this study, we used a combined list of MS cognition genes by using overrepresentation statistical analysis. This enrichment analysis identified the pathways significantly enriched with genes targeted by at least one of the selected miRNAs. The network was visualized using Cytoscape Software 3.7.2 [50]. This study used the cytoHubba plugin in Cytoscape to determine the potential association of pathways to genes associated with MS cognitive dysfunction [51]. The key genes responsible for MS cognitive dysfunction were identified by the cytoHubba plugin and miRNA–gene network.

2.7. Functional Annotation

We performed functional annotation in two phases, functional annotation of genes from text mining and other databases, and identification of risk pathways from the list of risk variants from GWAS. We used DAVID, a Database for Annotation, Visualization, and Integrated Discovery (https://david.ncifcrf.gov/, accessed on 1 August 2019), for identifying the functional annotation of MS cognition genes. The input to DAVID was a list of official gene symbols with a p-value cutoff of 0.05. We identified key biological processes (BP), molecular functions (MF), and cellular components (CC) using DAVID [52,53]. For this analysis, p < 0.05 was considered to be statistically significant, and the calculated false discovery rate (FDR) score was taken for the study.

3. Results

3.1. Collection of MS Cognition Genomic and Transcriptomic Data from Literature

PubTator identified 16 genes and 20 miRNA mentions in 4102 PubMed articles retrieved for MS and cognition. A total of 318 gene–miRNA associations were retrieved using the gene concept identification module of PubTator. All these associations were manually curated and verified to get an accurate list of genes/miRNAs. A nonredundant list of genes/miRNA and PubMed IDs were collected using an in-house Perl program. The genes associated with MS cognition were further mapped to their official gene symbols approved by HGNC, and a set of 16 genes were found to be involved in the cognition. A set of 20 nonredundant miRNAs associated with cognitive dysfunction in MS were identified and mapped to miRNA IDs from miRBase [54]. The identified list of genes and miRNAs with PubMed reference is shown in Supplementary Table S1 (gene-pubmed).

3.2. Identification of miRNA–mRNA Targets

Using DIANA-TarBase v8.0, experimentally validated miRNA–gene interactions were identified for each of the gene entity from our list of 16 genes. A total of 248 miRNA–mRNA associations were found along with the experimental evidence and related PubMed resource. This miRNA–mRNA pairs contained 83 nonredundant miRNAs, and these were visualized using Cytoscape software. The list of miRNA–mRNA interactions is shown in Supplementary Table S1 (gene-miRNA).

3.3. Functional Analysis of Genetic Variants

Functional analysis of genetic variants using Haploreg tool identified 63 eQTLS from MS and one eQTL associated with CI. A total of 57 coding and 19 noncoding mutations were found in MS, while 14 coding and 11 noncoding mutations were found in CI. Among the noncoding mutations, 16 were novel transcripts and three were long noncoding RNAs in MS, while eight were novel transcripts and three were long noncoding RNAs in CI. Using PhenoScanner, we identified that 77 genetic variants were eQTLs in MS, and intronic and intergenic locations were identified. The detailed descriptions are shown in Supplementary Table S2 for MS and Supplementary Table S3 for CI. Allelic variations and more detailed information are also provided.

3.4. KEGG Pathway Analysis

miRNAs involved in MS cognition dysfunctions were subjected to pathway analysis using DIANA-mirPath (mirPath v.3). KEGG pathway analysis utilizes overrepresentation statistical analysis and determines the experimentally validated results with a p-value. The network was visualized using Cytoscape 3.2, and the cytoHubba plugin was used to identify the hub genes associated with the reported pathways [55]. The miRNA–pathway associations are shown in Supplementary Figure S1. The complete list of 185 miRNA pathway interactions is given in Supplementary Table S4 (miR-KEGG).
Our results show the most significant pathways identified using the highest degree of distribution of the nodes. hsa-mir-155-5p has a higher number of pathway interactions; hence, this miRNA entity may serve as the main target for MS cognitive dysfunction therapy and treatment. hsa-mir-182-5p, hsa-mir-320a, hsa-mir-148b-3p, and has-mir-301a-5p also serve as the biomarkers in MS cognitive impairment analysis. Further analysis of these miRNA pathway interactions determined the most significant pathways, as shown in Figure 1. The ECM receptor pathway is highly involved in the miRNA pathway network.

3.5. Identification of KEGG Pathways Based on Genetic Variants from GWAS

Using these 76 mapped genes of MS, we identified 100 pathways associated with MS. Using 26 mapped genes of these genetic variants in CI, we identified 88 pathways associated with CI. These pathways were determined using KEGG pathway mapper https://www.genome.jp/kegg/tool/map_pathway1.html/, accessed on 15 August 2019). From these results, we determined that the genes mapped by the genetic variants are also involved in the ECM signaling pathway, PIK3/Akt signaling pathway, estrogen signaling pathway, and TGF-beta signaling pathway. In addition, we performed the functional annotation of the mapped genes of MS and CI.

3.6. DAVID Functional Annotation

All the target genes in the list were given as input to the DAVID tool to identify the related KEGG pathways. Gene Ontology (GO) analysis results revealed that the target genes were expressively enriched in biological processes (BP), including regulation of neurological system process, defense response, regulation of neuron apoptosis, regulation of synaptic transmission, regulation of signaling, neuron development, etc., as shown in Table 1. Molecular functions (MF) in GO included cytokine activity, growth factor activity, neurotrophin receptor binding, interferon-alpha/beta receptor binding, tumor necrosis factor receptor superfamily binding, and structural constituent of cytoskeleton, as shown in Table 2. Cellular components (CC) in GO involved the extracellular region, extracellular region part, and extracellular space functions, as shown in Table 3. Data with the complete list of reported GO functional annotation are shown in Supplementary Table S4. The genes involved in each of GO term functional enrichment analysis reflect the enriched genes involved in various biological pathways, as shown in Figure 2.
The functional annotation of 63 mapped genes from MS genetic variants and 26 mapped genes from CI genetic variants was determined using the DAVID tool. The results of KEGG enrichment analysis and GO pathway analysis are shown in Supplementary Table S5 for MS and Supplementary Table S6 for CI.

3.7. Identification of Key Risk Genes in miRNA–Gene Network

The miRNA–mRNA interactions involving 16 genes and 83 miRNAs were studied to determine the most significant genes that may act as risk genes in the complete MS cognition network. Among these interactions, brain-derived neurotrophic factor (BDNF) was the most significant gene, followed by interferon beta (IFNB), neurofilament L (NEFL), neurotropin (NTF3), and interleukin 10 (IL10). These genes were found to be the top five nodes ranked by cytoHubba, as shown in Figure 3. Other important genes included nerve growth factor (NGF) and chitinase-3-like protein 1 (CH3L1).

4. Discussion

We presented an integrated bioinformatics approach to analyze and report the most significant genes that could act as biomarkers in the treatment of MS cognitive impairment. Our approach used text mining, GWAS signals, and pathway analysis to identify the biomarkers associated with MS. Although several genomic and transcriptomic studies are currently reported in MS cognitive impairment, a comprehensive repository dealing with all the experimental data is still underdeveloped. Automated approaches like our approach will be useful to identify the biomarkers associated with MS cognitive impairment in an efficient way.
Cognitive impairment is known to be involved in the mTOR signaling pathway [56,57], PI3K/Akt/mTOR signaling pathway [58], cholesterol biosynthesis pathway [59,60], cholinergic pathway [61,62,63], visual pathway [64], etc. Management of cognitive dysfunction and decline in MS patients is more important to improve the lifestyle of people with MS. Several clinical trials reported the use of simvastatin in MS (MS-STAT trial), affecting brain atrophy, clinical, and cognitive outcome measures. These studies implicate the involvement of metabolites in the cholesterol synthesis pathway.
The development of computational models might elucidate the causal architecture underlying the treatment effects in clinical trials of progressive MS. Physical activity is associated with the motor function and mobility of healthy people, and MS patients have poor physical activity. MS and aging predict the motor function in people, but their cognitive function is not modified. Hence, there is huge need to evaluate and treat cognitive dysfunction in these patients who suffer from MS.
Studies suggest that visual pathway measures are known to bring improved outcomes in MS cognition during neuroprotection trials [65,66]. MS cognition is associated with the cholinergic pathway, but studies suggest that MS does not involve the reduction of cholinergic neurons as in other diseases like Alzheimer’s [67]. Acetyl choline esterase (AChE) inhibitors are used to treat cognitively impaired MS patients, but their therapeutic effects are poor and unsatisfactory [67,68].
CI genetic variants determined from GWAS signals reveal the presence of eQTLs that could regulate the expression of genes with p < 0.01, e.g., mapping gene PARP6 with the variant rs11637611 and risk allele C. This study was earlier reported by Hu et al. (rs11637611 with SNPs at chr15q23 and locus p = 1.07 × 10−15 [69]. Risk genes associated with CI are involved in risk pathways such as the TGF-beta signaling pathway, ECM receptor interaction, PI3k/Akt signaling pathway, mTOR signaling pathway, and estrogen signaling pathway. These results are in correlation with gene/miRNA functional analysis, thus confirming that these risk genes are highly associated with MS cognitive dysfunction.
Risk variant rs703842 located on chromosome 12 was found to be in association with CYP27B1 with a p-value < 1.00 × 10−5. However, this association was earlier reported by Jiang et al. in the association between genetic polymorphism rs703842 in CYP27B1 and MS [70]. Although MS genes were associated with cognition genes, the genes reported by GWAS variants in MS and CI were not correlated. When functional annotation was performed on these genes, we found that their pathways were correlated (see Supplementary Table S7).
In addition to the identification of risk genes from GWAS, we identified miRSNPS from the MSDD database [47] and found that hsa-miR-223 and hsa-miR-23a were associated with genetic variants such as rs1044165 and rs3745453, respectively [70]. These associations were reported in the literature [70,71]. However, our GWAS analysis did not validate these results. Variants such as rs17119 and rs1843938, located on chromosomes 6 and 7, were mapped with novel transcripts with a p-value of 1.00 × 10−10. These variants were highly correlated. Their associations were earlier validated by GWAS studies performed by Disanto et al. [72].

5. Conclusions

Our study showed that BDNF, IFNB, IL10, NEFL, and CHI3L1 could be the potential core genes that play an important role in cognitive dysfunction and impairment in MS. GWAS signals identified the risk variants and alleles associated with MS and cognition, as well as their association with the mapped genes and risk pathways. To further confirm these studies, we utilized miRNAs and their target genes to determine the risk pathways associated with miRNAs. The results from both miRNA functional annotation and GWAS risk gene functional annotation were correlated to determine the key functional genes and risk pathways. In the future, we will evaluate our findings using gene expression datasets in MS patient samples who are cognitively impaired or have cognitive decline.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/diagnostics12081914/s1: Table S1. List of genes and miRNAs associated with MS cognitive dysfunction; Table S2. Genetic variants associated with MS; Table S3. Genetic variants associated with cognitive impairment (CI); Table S4. miRNA pathway interactions from KEGG for MS cognitive dysfunction; Table S5. KEGG enrichment analysis and GO pathway analysis for MS; Table S6. KEGG enrichment analysis and GO pathway analysis for CI; Table S7. List of KEGG correlated pathways in MS and CI. Figure S1. miRNA-pathway associations involved in MS-cognition (Orange nodes represents miRNA and the green nodes represent KEGG pathways).

Author Contributions

Conceptualization, A.P. and K.R.; methodology, A.P.; software, A.P.; validation, A.P. and K.R.; formal analysis, A.P. and K.R.; investigation, A.P.; resources, A.P.; data curation, A.P.; writing—original draft preparation, A.P.; writing—review and editing, A.P. and K.R.; visualization, A.P.; supervision, K.R.; project administration, A.P and K.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Patient consent was not required because the study used publicly available data only.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Karussis, D. The diagnosis of multiple sclerosis and the various related demyelinating syndromes: A critical review. J. Autoimmun. 2014, 48, 134–142. [Google Scholar] [CrossRef] [PubMed]
  2. Koch-Henriksen, N.; Sørensen, P.S. The changing demographic pattern of multiple sclerosis epidemiology. Lancet Neurol. 2010, 9, 520–532. [Google Scholar] [CrossRef]
  3. Skarke, C. Dissecting the Mystery Around Multiple Sclerosis. Sci. Transl. Med. 2012, 4, 162ec216. [Google Scholar] [CrossRef]
  4. de Sa, J.C.C.; Airas, L.; Bartholome, E.; Grigoriadis, N.; Mattle, H.; Oreja-Guevara, C.; O’Riordan, J.; Sellebjerg, F.; Stankoff, B.; Vass, K. Symptomatic therapy in multiple sclerosis: A review for a multimodal approach in clinical practice. Ther. Adv. Neurol. Disord. 2011, 4, 139–168. [Google Scholar] [CrossRef] [Green Version]
  5. Gil-González, I.; Martín-Rodríguez, A.; Conrad, R.; Pérez-San-Gregorio, M.Á. Quality of life in adults with multiple sclerosis: A systematic review. BMJ Open 2020, 10, e041249. [Google Scholar] [CrossRef]
  6. Benedict, R.H.; Zivadinov, R. Risk factors for and management of cognitive dysfunction in multiple sclerosis. Nat. Rev. Neurol. 2011, 7, 332–342. [Google Scholar] [CrossRef]
  7. Rocca, M.A.; Amato, M.P.; De Stefano, N.; Enzinger, C.; Geurts, J.J.; Penner, I.-K.; Rovira, A.; Sumowski, J.F.; Valsasina, P.; Filippi, M. Clinical and imaging assessment of cognitive dysfunction in multiple sclerosis. Lancet Neurol. 2015, 14, 302–317. [Google Scholar] [CrossRef]
  8. Shany-Ur, T.; Rankin, K.P. Personality and social cognition in neurodegenerative disease. Curr. Opin. Neurol. 2011, 24, 550–555. [Google Scholar] [CrossRef] [Green Version]
  9. Neary, D.; Snowden, J. Fronto-temporal dementia: Nosology, neuropsychology, and neuropathology. Brain Cogn. 1996, 31, 176–187. [Google Scholar] [CrossRef] [Green Version]
  10. Migliore, S.; Ghazaryan, A.; Simonelli, I.; Pasqualetti, P.; Squitieri, F.; Curcio, G.; Landi, D.; Palmieri, M.; Moffa, F.; Filippi, M. Cognitive impairment in relapsing-remitting multiple sclerosis patients with very mild clinical disability. Behav. Neurol. 2017, 2017, 7404289. [Google Scholar] [CrossRef]
  11. Harris, V.K.; Sadiq, S.A. Biomarkers of therapeutic response in multiple sclerosis: Current status. Mol. Diagn. Ther. 2014, 18, 605–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, H.; Zhang, P.; Chen, W.; Feng, D.; Jia, Y.; Xie, L.-X. Evidence for serum miR-15a and miR-16 levels as biomarkers that distinguish sepsis from systemic inflammatory response syndrome in human subjects. Clin. Chem. Lab. Med. 2012, 50, 1423–1428. [Google Scholar] [CrossRef] [PubMed]
  13. Sahab, Z.J.; Semaan, S.M.; Sang, Q.-X.A. Methodology and applications of disease biomarker identification in human serum. Biomark. Insights 2007, 2, 21–43. [Google Scholar] [CrossRef] [PubMed]
  14. Comabella, M.; Montalban, X. Body fluid biomarkers in multiple sclerosis. Lancet Neurol. 2014, 13, 113–126. [Google Scholar] [CrossRef]
  15. Gaetani, L.; Salvadori, N.; Chipi, E.; Gentili, L.; Borrelli, A.; Parnetti, L.; Di Filippo, M. Cognitive impairment in multiple sclerosis: Lessons from cerebrospinal fluid biomarkers. Neural Regen. Res. 2021, 16, 36. [Google Scholar] [PubMed]
  16. Jakimovski, D.; Bergsland, N.; Dwyer, M.G.; Traversone, J.; Hagemeier, J.; Fuchs, T.A.; Ramasamy, D.P.; Weinstock-Guttman, B.; Benedict, R.H.; Zivadinov, R. Cortical and deep gray matter perfusion associations with physical and cognitive performance in multiple sclerosis patients. Front. Neurol. 2020, 11, 700. [Google Scholar] [CrossRef]
  17. Jakimovski, D.S.; Gandhi, I.; Paunkoski, N.; Bergsland, J.; Hagemeier, D.P.; Ramasamy, D.; Hojnacki, D.; Kolb, C.; Benedict, R.H.B.; Weinstock-Guttman, B.; et al. Hypertension and heart disease are associated with development of brain atrophy in multiple sclerosis: A 5-year longitudinal study. Eur. J. Neurol. 2019, 26, 87-e8. [Google Scholar] [CrossRef] [Green Version]
  18. Wandinger, K.P.; Stürzebecher, C.S.; Bielekova, B.; Detore, G.; Rosenwald, A.; Staudt, L.M.; McFarland, H.F.; Martin, R. Complex immunomodulatory effects of interferon-β in multiple sclerosis include the upregulation of T helper 1-associated marker genes. Ann. Neurol. Off. J. Am. Neurol. Assoc. Child Neurol. Soc. 2001, 50, 349–357. [Google Scholar] [CrossRef]
  19. Wandinger, K.-P.; Lünemann, J.D.; Wengert, O.; Bellmann-Strobl, J.; Aktas, O.; Weber, A.; Grundström, E.; Ehrlich, S.; Wernecke, K.-D.; Volk, H.-D. TNF-related apoptosis inducing ligand (TRAIL) as a potential response marker for interferon-beta treatment in multiple sclerosis. Lancet 2003, 361, 2036–2043. [Google Scholar] [CrossRef]
  20. Ingram, G.; Hakobyan, S.; Hirst, C.L.; Harris, C.L.; Pickersgill, T.P.; Cossburn, M.D.; Loveless, S.; Robertson, N.P.; Morgan, B.P. Complement regulator factor H as a serum biomarker of multiple sclerosis disease state. Brain 2010, 133, 1602–1611. [Google Scholar] [CrossRef] [Green Version]
  21. Quintana, E.; Coll, C.; Salavedra-Pont, J.; Muñoz-San Martín, M.; Robles-Cedeño, R.; Tomàs-Roig, J.; Buxó, M.; Matute-Blanch, C.; Villar, L.; Montalban, X. Cognitive impairment in early stages of multiple sclerosis is associated with high cerebrospinal fluid levels of chitinase 3-like 1 and neurofilament light chain. Eur. J. Neurol. 2018, 25, 1189–1191. [Google Scholar] [CrossRef] [PubMed]
  22. Cortez, M.A.; Calin, G.A. MicroRNA identification in plasma and serum: A new tool to diagnose and monitor diseases. Expert Opin. Biol. Ther. 2009, 9, 703–711. [Google Scholar] [CrossRef] [PubMed]
  23. Li, Y.; Zhang, Z.; Liu, F.; Vongsangnak, W.; Jing, Q.; Shen, B. Performance comparison and evaluation of software tools for microRNA deep-sequencing data analysis. Nucleic Acids Res. 2012, 40, 4298–4305. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Bhalala, O.G.; Srikanth, M.; Kessler, J.A. The emerging roles of microRNAs in CNS injuries. Nat. Rev. Neurol. 2013, 9, 328–339. [Google Scholar] [CrossRef] [PubMed]
  25. Moon, J.-m.; Xu, L.; Giffard, R.G. Inhibition of microRNA-181 reduces forebrain ischemia-induced neuronal loss. J. Cereb. Blood Flow Metab. 2013, 33, 1976–1982. [Google Scholar] [CrossRef] [Green Version]
  26. Liguori, M.; Nuzziello, N.; Simone, M.; Amoroso, N.; Viterbo, R.G.; Tangaro, S.; Consiglio, A.; Giordano, P.; Bellotti, R.; Trojano, M. Association between miRNAs expression and cognitive performances of Pediatric Multiple Sclerosis patients: A pilot study. Brain Behav. 2019, 9, e01199. [Google Scholar] [CrossRef] [Green Version]
  27. Comi, G.; Filippi, M.; Barkhof, F.; Durelli, L.; Edan, G.; Fernández, O.; Hartung, H.-P.; Seeldrayers, P.; Sørensen, P.S.; Rovaris, M. Effect of early interferon treatment on conversion to definite multiple sclerosis: A randomised study. Lancet 2001, 357, 1576–1582. [Google Scholar] [CrossRef]
  28. Siva, A. The spectrum of multiple sclerosis and treatment decisions. Clin. Neurol. Neurosurg. 2006, 108, 333–338. [Google Scholar] [CrossRef] [PubMed]
  29. Scapoli, C.; Ziliotto, N.; Lunghi, B.; Menegatti, E.; Salvi, F.; Zamboni, P.; Baroni, M.; Mascoli, F.; Bernardi, F.; Marchetti, G. Combination of genomic and transcriptomic approaches highlights vascular and circadian clock components in multiple sclerosis. Int. J. Mol. Sci. 2021, 23, 310. [Google Scholar] [CrossRef]
  30. Marchetti, G.; Ziliotto, N.; Meneghetti, S.; Baroni, M.; Lunghi, B.; Menegatti, E.; Pedriali, M.; Salvi, F.; Bartolomei, I.; Straudi, S. Changes in expression profiles of internal jugular vein wall and plasma protein levels in multiple sclerosis. Mol. Med. 2018, 24, 1–11. [Google Scholar] [CrossRef]
  31. Ricaño-Ponce, I.; Wijmenga, C. Mapping of immune-mediated disease genes. Annu. Rev. Genom. Hum. Genet. 2013, 14, 325–353. [Google Scholar] [CrossRef] [PubMed]
  32. Liu, C.; Zhang, F.; Li, T.; Lu, M.; Wang, L.; Yue, W.; Zhang, D. MirSNP, a database of polymorphisms altering miRNA target sites, identifies miRNA-related SNPs in GWAS SNPs and eQTLs. BMC Genom. 2012, 13, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Klein, K.; Tremmel, R.; Winter, S.; Fehr, S.; Battke, F.; Scheurenbrand, T.; Schaeffeler, E.; Biskup, S.; Schwab, M.; Zanger, U.M. A new panel-based next-generation sequencing method for ADME genes reveals novel associations of common and rare variants with expression in a human liver cohort. Front. Genet. 2019, 10, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Stegeman, S.; Amankwah, E.; Klein, K.; O’Mara, T.A.; Kim, D.; Lin, H.-Y.; Permuth-Wey, J.; Sellers, T.A.; Srinivasan, S.; Eeles, R. A large-scale analysis of genetic variants within putative miRNA binding sites in prostate cancer. Cancer Discov. 2015, 5, 368–379. [Google Scholar] [CrossRef] [Green Version]
  35. Cipolla, G.A.; Park, J.K.; de Oliveira, L.A.; Lobo-Alves, S.C.; de Almeida, R.C.; Farias, T.D.; Lemos, D.d.S.; Malheiros, D.; Lavker, R.M.; Petzl-Erler, M.L. A 3′ UTR polymorphism marks differential KLRG1 mRNA levels through disruption of a miR-584-5p binding site and associates with pemphigus foliaceus susceptibility. Biochim. Biophys. Acta (BBA) Gene Regul. Mech. 2016, 1859, 1306–1313. [Google Scholar] [CrossRef] [Green Version]
  36. Gao, F.; Yao, Y.; Zhang, Y.; Tian, J. Integrating genome-wide association studies with pathway analysis and gene expression analysis highlights novel osteoarthritis risk pathways and genes. Front. Genet. 2019, 10, 827. [Google Scholar] [CrossRef] [Green Version]
  37. Wei, C.-H.; Kao, H.-Y.; Lu, Z. PubTator: A web-based text mining tool for assisting biocuration. Nucleic Acids Res. 2013, 41, W518–W522. [Google Scholar] [CrossRef]
  38. Wei, C.-H.; Leaman, R.; Lu, Z. Beyond accuracy: Creating interoperable and scalable text-mining web services. Bioinformatics 2016, 32, 1907–1910. [Google Scholar] [CrossRef] [Green Version]
  39. Povey, S.; Lovering, R.; Bruford, E.; Wright, M.; Lush, M.; Wain, H. The HUGO gene nomenclature committee (HGNC). Human Genet. 2001, 109, 678–680. [Google Scholar] [CrossRef]
  40. Wheeler, D.L.; Barrett, T.; Benson, D.A.; Bryant, S.H.; Canese, K.; Chetvernin, V.; Church, D.M.; DiCuccio, M.; Edgar, R.; Federhen, S. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2007, 36, D13–D21. [Google Scholar] [CrossRef] [Green Version]
  41. Karagkouni, D.; Paraskevopoulou, M.D.; Chatzopoulos, S.; Vlachos, I.S.; Tastsoglou, S.; Kanellos, I.; Papadimitriou, D.; Kavakiotis, I.; Maniou, S.; Skoufos, G. DIANA-TarBase v8: A decade-long collection of experimentally supported miRNA–gene interactions. Nucleic Acids Res. 2018, 46, D239–D245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Bartel, D.P. MicroRNAs: Target recognition and regulatory functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Lewis, B.P.; Burge, C.B.; Bartel, D.P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005, 120, 15–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Prabahar, A.; Natarajan, J. Prediction of microRNAs involved in immune system diseases through network based features. J. Biomed. Inform. 2017, 65, 34–45. [Google Scholar] [CrossRef] [PubMed]
  45. Buniello, A.; MacArthur, J.A.L.; Cerezo, M.; Harris, L.W.; Hayhurst, J.; Malangone, C.; McMahon, A.; Morales, J.; Mountjoy, E.; Sollis, E. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019, 47, D1005–D1012. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Ward, L.D.; Kellis, M. HaploReg: A resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012, 40, D930–D934. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Staley, J.R.; Blackshaw, J.; Kamat, M.A.; Ellis, S.; Surendran, P.; Sun, B.B.; Paul, D.S.; Freitag, D.; Burgess, S.; Danesh, J. PhenoScanner: A database of human genotype–phenotype associations. Bioinformatics 2016, 32, 3207–3209. [Google Scholar] [CrossRef] [Green Version]
  48. Yue, M.; Zhou, D.; Zhi, H.; Wang, P.; Zhang, Y.; Gao, Y.; Guo, M.; Li, X.; Wang, Y.; Zhang, Y. MSDD: A manually curated database of experimentally supported associations among miRNAs, SNPs and human diseases. Nucleic Acids Res. 2018, 46, D181–D185. [Google Scholar] [CrossRef]
  49. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-miRPath v3. 0: Deciphering microRNA function with experimental support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef]
  50. Cline, M.S.; Smoot, M.; Cerami, E.; Kuchinsky, A.; Landys, N.; Workman, C.; Christmas, R.; Avila-Campilo, I.; Creech, M.; Gross, B. Integration of biological networks and gene expression data using Cytoscape. Nat. Protoc. 2007, 2, 2366–2382. [Google Scholar] [CrossRef] [Green Version]
  51. Chin, C.-H.; Chen, S.-H.; Wu, H.-H.; Ho, C.-W.; Ko, M.-T.; Lin, C.-Y. cytoHubba: Identifying hub objects and sub-networks from complex interactome. BMC Syst. Biol. 2014, 8, 1–7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  53. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Kozomara, A.; Griffiths-Jones, S. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42, D68–D73. [Google Scholar] [CrossRef] [Green Version]
  55. Ni, W.-J.; Leng, X.-M. Dynamic miRNA–mRNA paradigms: New faces of miRNAs. Biochem. Biophys. Rep. 2015, 4, 337–341. [Google Scholar] [CrossRef] [Green Version]
  56. Yang, L.; Zhang, Y.; Yan, Z.; Tian, F. The role of mTOR signaling pathway on cognitive functions in cerebral ischemia-reperfusion. Exp. Ther. Med. 2017, 14, 2839–2844. [Google Scholar] [CrossRef] [Green Version]
  57. Cho, C.; Michailidis, V.; Martin, L.J. Revealing brain mechanisms of mTOR-mediated translational regulation: Implications for chronic pain. Neurobiol. Pain 2018, 4, 27–34. [Google Scholar] [CrossRef] [PubMed]
  58. Mammana, S.; Bramanti, P.; Mazzon, E.; Cavalli, E.; Basile, M.S.; Fagone, P.; Petralia, M.C.; McCubrey, J.A.; Nicoletti, F.; Mangano, K. Preclinical evaluation of the PI3K/Akt/mTOR pathway in animal models of multiple sclerosis. Oncotarget 2018, 9, 8263. [Google Scholar] [CrossRef] [Green Version]
  59. Uher, T.; Fellows, K.; Horakova, D.; Zivadinov, R.; Vaneckova, M.; Sobisek, L.; Tyblova, M.; Seidl, Z.; Krasensky, J.; Bergsland, N. Serum lipid profile changes predict neurodegeneration in interferon-β1a-treated multiple sclerosis patients. J. Lipid Res. 2017, 58, 403–411. [Google Scholar] [CrossRef] [Green Version]
  60. Eshaghi, A.; Kievit, R.A.; Prados, F.; Sudre, C.H.; Nicholas, J.; Cardoso, M.J.; Chan, D.; Nicholas, R.; Ourselin, S.; Greenwood, J. Applying causal models to explore the mechanism of action of simvastatin in progressive multiple sclerosis. Proc. Natl. Acad. Sci. USA 2019, 116, 11020–11027. [Google Scholar] [CrossRef] [Green Version]
  61. Jiang, Y.-Y.; Shao, S.; Zhang, Y.; Zheng, J.; Chen, X.; Cui, S.; Liu, F.-Y.; Wan, Y.; Yi, M. Neural pathways in medial septal cholinergic modulation of chronic pain: Distinct contribution of the anterior cingulate cortex and ventral hippocampus. Pain 2018, 159, 1550–1561. [Google Scholar] [CrossRef] [PubMed]
  62. Cotter, J.; Muhlert, N.; Talwar, A.; Granger, K. Examining the effectiveness of acetylcholinesterase inhibitors and stimulant-based medications for cognitive dysfunction in multiple sclerosis: A systematic review and meta-analysis. Neurosci. Biobehav. Rev. 2018, 86, 99–107. [Google Scholar] [CrossRef] [PubMed]
  63. Di Filippo, M.; Portaccio, E.; Mancini, A.; Calabresi, P. Multiple sclerosis and cognition: Synaptic failure and network dysfunction. Nat. Rev. Neurosci. 2018, 19, 599–609. [Google Scholar] [CrossRef]
  64. Nguyen, J.; Rothman, A.; Fitzgerald, K.; Whetstone, A.; Syc-Mazurek, S.; Aquino, J.; Balcer, L.J.; Frohman, E.M.; Frohman, T.C.; Crainiceanu, C. Visual pathway measures are associated with neuropsychological function in multiple sclerosis. Curr. Eye Res. 2018, 43, 941–948. [Google Scholar] [CrossRef]
  65. Motl, R.W.; Sandroff, B.M. Exercise as a countermeasure to declining central nervous system function in multiple sclerosis. Clin. Ther. 2018, 40, 16–25. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Sandroff, B.M.; Motl, R.W.; Kam, J.P.; Pula, J.H. Accelerometer measured physical activity and the integrity of the anterior visual pathway in multiple sclerosis. Mult. Scler. Relat. Disord. 2014, 3, 117–122. [Google Scholar] [CrossRef] [PubMed]
  67. Schliebs, R.; Arendt, T. The cholinergic system in aging and neuronal degeneration. Behav. Brain Res. 2011, 221, 555–563. [Google Scholar] [CrossRef]
  68. Kimura, Y.; Sato, N.; Ota, M.; Maikusa, N.; Maekawa, T.; Sone, D.; Enokizono, M.; Sugiyama, A.; Imabayashi, E.; Matsuda, H. A structural MRI study of cholinergic pathways and cognition in multiple sclerosis. eNeurologicalSci 2017, 8, 11–16. [Google Scholar] [CrossRef]
  69. Hu, X.; Pickering, E.; Hall, S.; Naik, S.; Liu, Y.; Soares, H.; Katz, E.; Paciga, S.; Liu, W.; Aisen, P. Genome-wide association study identifies multiple novel loci associated with disease progression in subjects with mild cognitive impairment. Transl. Psychiatry 2011, 1, e54. [Google Scholar] [CrossRef] [Green Version]
  70. Jiang, T.; Li, L.; Wang, Y.; Zhao, C.; Yang, J.; Ma, D.; Guan, Y.; Zhao, D.; Bao, Y.; Wang, Y. The association between genetic polymorphism rs703842 in CYP27B1 and multiple sclerosis: A meta-analysis. Medicine 2016, 95, e3612. [Google Scholar] [CrossRef]
  71. Ridolfi, E.; Fenoglio, C.; Cantoni, C.; Calvi, A.; De Riz, M.; Pietroboni, A.; Villa, C.; Serpente, M.; Bonsi, R.; Vercellino, M. Expression and genetic analysis of microRNAs involved in multiple sclerosis. Int. J. Mol. Sci. 2013, 14, 4375–4384. [Google Scholar] [CrossRef] [PubMed]
  72. Disanto, G.; Dobson, R.; Pakpoor, J.; Elangovan, R.I.; Adiutori, R.; Kuhle, J.; Giovannoni, G. The refinement of genetic predictors of multiple sclerosis. PLoS ONE 2014, 9, e96578. [Google Scholar] [CrossRef] [PubMed]
Figure 1. miRNA pathway interactions of most significant pathways.
Figure 1. miRNA pathway interactions of most significant pathways.
Diagnostics 12 01914 g001
Figure 2. The X and Y-axes represent pathways (GO Terms) and the −log10 p-value of the corresponding GO term, respectively. GO terms are represented as follows: BP (biological process), CC (cellular component), and MF (molecular function).
Figure 2. The X and Y-axes represent pathways (GO Terms) and the −log10 p-value of the corresponding GO term, respectively. GO terms are represented as follows: BP (biological process), CC (cellular component), and MF (molecular function).
Diagnostics 12 01914 g002
Figure 3. Top five ranked gene nodes and their interactions from gene–miRNA network.
Figure 3. Top five ranked gene nodes and their interactions from gene–miRNA network.
Diagnostics 12 01914 g003
Table 1. GO analysis biological process (BP) of genes associated with cognitive dysfunction in MS.
Table 1. GO analysis biological process (BP) of genes associated with cognitive dysfunction in MS.
GO IDDescriptionCountp-ValueGenes
GO:0031644Regulation of neurological system process64.75 × 10−7BDNF, TNF, NTF4, NTF3, IL10, NGF
GO:0006952Defense response81.81 × 10−6IL17A, IFNA1, TNF, STAB1, IFNB1, IL17F, IL10, NGF
GO:0043523Regulation of neuron apoptosis52.36 × 10−6BDNF, TNF, NTF3, NEFL, NGF
GO:0050804Regulation of synaptic transmission51.22 × 10−5BDNF, TNF, NTF4, NTF3, NGF
GO:0009611Response to wounding71.30 × 10−5IL17A, TNF, STAB1, IL17F, NEFL, IL10, NGF
GO:0044057Regulation of system process61.50 × 10−5BDNF, TNF, NTF4, NTF3, IL10, NGF
GO:0051969Regulation of transmission of nerve impulse51.66 × 10−5BDNF, TNF, NTF4, NTF3, NGF
GO:0006954Inflammatory response61.91 × 10−5IL17A, TNF, STAB1, IL17F, IL10, NGF
GO:0009617Response to bacterium54.85 × 10−5TNF, STAB1, IFNB1, IL10, NGF
GO:0051384Response to glucocorticoid stimulus47.98 × 10−5TNF, NEFL, IL10, NGF
GO:0031960Response to corticosteroid stimulus41.03 × 10−4TNF, NEFL, IL10, NGF
GO:0042981Regulation of apoptosis71.36 × 10−4BDNF, TNF, NTF3, IFNB1, NEFL, IL10, NGF
GO:0043067Regulation of programmed cell death71.44 × 10−4BDNF, TNF, NTF3, IFNB1, NEFL, IL10, NGF
GO:0031175Neuron projection development51.45 × 10−4BDNF, NTF3, NTNG2, NEFL, NGF
GO:0010941Regulation of cell death71.47 × 10−4BDNF, TNF, NTF3, IFNB1, NEFL, IL10, NGF
GO:0051094Positive regulation of developmental process51.99 × 10−4BDNF, TNF, NTF3, NEFL, NGF
GO:0042742Defense response to bacterium42.34 × 10−4TNF, STAB1, IFNB1, IL10
GO:0048666neuron development54.25 × 10−4BDNF, NTF3, NTNG2, NEFL, NGF
GO:0043066Negative regulation of apoptosis55.01 × 10−4BDNF, TNF, NEFL, IL10, NGF
GO:0043069Negative regulation of programmed cell death55.28 × 10−4BDNF, TNF, NEFL, IL10, NGF
GO:0060548Negative regulation of cell death55.33 × 10−4BDNF, TNF, NEFL, IL10, NGF
GO:0030030Cell projection organization55.80 × 10−4BDNF, NTF3, NTNG2, NEFL, NGF
Table 2. GO analysis molecular function (MF) of genes associated with cognitive dysfunction in MS.
Table 2. GO analysis molecular function (MF) of genes associated with cognitive dysfunction in MS.
GO IDDescriptionCountp ValueGenes
GO:0005125Cytokine activity61.30 × 10−6IL17A, IFNA1, TNF, IFNB1, IL17F, IL10
GO:0008083Growth factor activity52.07 × 10−5BDNF, NTF4, NTF3, IL10, NGF
GO:0005165Neurotrophin receptor binding20.002156NTF3, NGF
GO:0005132Interferon-alpha/beta receptor binding20.009666IFNA1, IFNB1
GO:0032813Tumor necrosis factor receptor superfamily binding20.032931TNF, NGF
GO:0005200Structural constituent of cytoskeleton20.076944NEFL, SPTB
Table 3. GO analysis cell component (CC) of genes associated with cognitive dysfunction in MS.
Table 3. GO analysis cell component (CC) of genes associated with cognitive dysfunction in MS.
GO IDDescriptionCountp ValueGenes
GO:0005576Extracellular region121.04 × 10−6IL17A, IFNA1, BDNF, TNF, NTF4, NTF3, IFNB1, IL17F, CHI3L1, NTNG2, IL10, NGF
GO:0044421Extracellular region part93.93 × 10−6IL17A, IFNA1, TNF, IFNB1, IL17F, CHI3L1, NTNG2, IL10, NGF
GO:0005615Extracellular space85.43 × 10−6IL17A, IFNA1, TNF, IFNB1, IL17F, CHI3L1, IL10, NGF
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Prabahar, A.; Raja, K. Integrated Approaches to Identify miRNA Biomarkers Associated with Cognitive Dysfunction in Multiple Sclerosis Using Text Mining, Gene Expression, Pathways, and GWAS. Diagnostics 2022, 12, 1914. https://doi.org/10.3390/diagnostics12081914

AMA Style

Prabahar A, Raja K. Integrated Approaches to Identify miRNA Biomarkers Associated with Cognitive Dysfunction in Multiple Sclerosis Using Text Mining, Gene Expression, Pathways, and GWAS. Diagnostics. 2022; 12(8):1914. https://doi.org/10.3390/diagnostics12081914

Chicago/Turabian Style

Prabahar, Archana, and Kalpana Raja. 2022. "Integrated Approaches to Identify miRNA Biomarkers Associated with Cognitive Dysfunction in Multiple Sclerosis Using Text Mining, Gene Expression, Pathways, and GWAS" Diagnostics 12, no. 8: 1914. https://doi.org/10.3390/diagnostics12081914

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

Article Metrics

Back to TopTop