Meta-Analysis of Multiple Sclerosis Microarray Data Reveals Dysregulation in RNA Splicing Regulatory Genes

circRNAs as Epigenetic Regulators of Integrity in Blood-Brain Barrier Architecture: Mechanisms and Therapeutic Strategies in Multiple Sclerosis

Multiple sclerosis (MS) is a chronic inflammatory neurodegenerative disease leading to progressive demyelination and neuronal loss, with extensive neurological symptoms. As one of the most widespread neurodegenerative disorders, with an age onset of about 30 years, it turns out to be a socio-health and economic issue, thus necessitating therapeutic interventions currently unavailable. Loss of integrity in the blood-brain barrier (BBB) is one of the distinct MS hallmarks. Brain homeostasis is ensured by an endothelial cell-based monolayer at the interface between the central nervous system (CNS) and systemic bloodstream, acting as a selective barrier. MS results in enhanced barrier permeability, mainly due to the breakdown of tight (TJs) and adherens junctions (AJs) between endothelial cells. Specifically, proinflammatory mediator release causes failure in cytoplasmic exposure of junctions, resulting in compromised BBB integrity that enables blood cells to cross the barrier, establishing iron deposition and neuronal impairment. Cells with a compromised cytoskeletal protein network, fiber reorganization, and discontinuous junction structure can occur, resulting in BBB dysfunction. Recent investigations on spatial transcriptomics have proven circularRNAs (circRNAs) to be powerful multi-functional molecules able to epigenetically regulate transcription and structurally support proteins. In the present review, we provide an overview of the recent role ascribed to circRNAs in maintaining BBB integrity/permeability via cytoskeletal stability. Increased knowledge of the mechanisms responsible for impairment and circRNA's role in driving BBB damage and dysfunction might be helpful for the recognition of novel therapeutic targets to overcome BBB damage and unrestrained neurodegeneration.

Identifying inflammation-related targets of natural lactones using network pharmacology, molecular modeling and in vitro approaches

Natural lactones have been used in traditional and folklore medicine for centuries owing to their anti-inflammatory properties. The study uses a multifaceted approach to identify lead anti-inflammatory lactones from the SISTEMATX natural products database. The study analyzed the natural lactone database, revealing 18 lactones linked to inflammation targets. The primary targets were PTGES, PTGS1, COX-2, ALOX5 and IL1B. STX 12273 was the best hit, with the lowest binding energy and potential for inhibiting the COX-2 enzyme. The study suggested natural lactone, STX 12273, from the SISTEMATX database with anti-inflammatory potential and postulated its use for inflammation treatment or prevention.Communicated by Ramaswamy H. Sarma.

hnRNP A1 dysfunction alters RNA splicing and drives neurodegeneration in multiple sclerosis (MS)

Neurodegeneration is the primary driver of disease progression in multiple sclerosis (MS) resulting in permanent disability, creating an urgent need to discover its underlying mechanisms. Herein, we establish that dysfunction of the RNA binding protein heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) results in differential of binding to RNA targets causing alternative RNA splicing, which contributes to neurodegeneration in MS and its models. Using RNAseq of MS brains, we discovered differential expression and aberrant splicing of hnRNP A1 target RNAs involved in neuronal function and RNA homeostasis. We confirmed this in vivo in experimental autoimmune encephalomyelitis employing CLIPseq specific for hnRNP A1, where hnRNP A1 differentially binds and regulates RNA, including aberrantly spliced targets identified in human samples. Additionally, dysfunctional hnRNP A1 expression in neurons caused neurite loss and identical changes in splicing, corroborating hnRNP A1 dysfunction as a cause of neurodegeneration. Collectively, these data indicate hnRNP A1 dysfunction causes altered neuronal RNA splicing, resulting in neurodegeneration in MS.

TNF-α/STAT1/CXCL10 mutual inflammatory axis that contributes to the pathogenesis of experimental models of multiple sclerosis: A promising signaling pathway for targeted therapies

BACKGROUND

Identifying mutual neuroinflammatory axis in different experimental models of multiple sclerosis (MS) is essential to evaluate the de- and re-myelination processes and improve therapeutic interventions' reproducibility.

METHODS

The expression profile data set of EAE (GSE47900) and cuprizone (GSE100663) models were downloaded from the Gene Expression Omnibus database. The R package and GEO2R software processed these raw chip data. Gene Ontology (GO) functional analysis, KEGG pathway analysis, and protein-protein interaction network analysis were performed to investigate interactions between common differentially expressed genes (DEGs) in all models. Finally, the ELISA method assessed the protein level of highlighted mutual cytokines in serum.

RESULTS

Our data introduced 59 upregulated [CXCL10, CCL12, and GBP6 as most important] and 17 downregulated [Serpinb1a, Prr18, and Ugt8a as most important] mutual genes. The signal transducer and activator of transcription 1 (STAT1) and CXCL10 were the most crucial hub proteins among mutual upregulated genes. These mutual genes were found to be mainly involved in the TNF-α, TLRs, and complement cascade signaling, and animal models shared 26 mutual genes with MS individuals. Finally, significant upregulation of serum level of TNF-α/IL-1β/CXCL10 cytokines was confirmed in all models in a relatively similar pattern.

CONCLUSION

For the first time, our study revealed the common neuroinflammatory pathway in animal models of MS and introduced candidate hub genes for better evaluating the preclinical efficacy of pharmacological interventions and designing prospective targeted therapies.

The circular RNA landscape in multiple sclerosis: Disease-specific associated variants and exon methylation shape circular RNA expression profile

BACKGROUND

Circular RNAs (circRNAs) are a class of non-coding RNAs increasingly emerging as crucial actors in the pathogenesis of human diseases, including autoimmune and neurological disorders as multiple sclerosis (MS). Despite several efforts, the mechanisms regulating circRNAs expression are still largely unknown and the circRNA profile and regulation in MS-relevant cell models has not been completely investigated. In this work, we aimed at exploring the global landscape of circRNA expression in MS patients, also evaluating a possible correlation with their genetic and epigenetic background.

METHODS

We performed RNA-seq experiments on circRNA-enriched samples, derived from peripheral blood mononuclear cells (PBMCs) of 10 MS patients and 10 matched controls and performed differential circRNA expression. The genetic background was evaluated using array genotyping, and an expression quantitative trait loci (eQTL) analysis was carried out.

RESULTS

Expression analysis revealed 166 differentially expressed circRNAs in MS patients, 125 of which are downregulated. One of the top dysregulated circRNAs, hsa_circ_0007990, derives from the PGAP3 gene, encoding a protein relevant for the control of autoimmune responses. The downregulation of this circRNA was confirmed in two independent replication cohorts, suggesting its implementation as a possible RNA-based biomarker. The eQTL analysis evidenced a significant association between 89 MS-associated loci and the expression of at least one circRNA, suggesting that MS-associated variants could impact on disease pathogenesis by altering circRNA profiles. Finally, we found a significant correlation between exon methylation and circRNA expression levels, supporting the hypothesis that epigenetic features may play an important role in the definition of the cell circRNA pool.

CONCLUSION

We described the circRNA expression profile of PBMCs in MS patients, suggesting that MS-associated variants may tune the expression levels of circRNAs acting as "circ-QTLs", and proposing a role for exon-based DNA methylation in regulating circRNA expression.

Multiple sclerosis and computational biology (Review)

Multiple sclerosis (MS) is an autoimmune neurodegenerative disease whose prevalence has increased worldwide. The resultant symptoms may be debilitating and can substantially reduce the of patients. Computational biology, which involves the use of computational tools to answer biomedical questions, may provide the basis for novel healthcare approaches in the context of MS. The rapid accumulation of health data, and the ever-increasing computational power and evolving technology have helped to modernize and refine MS research. From the discovery of novel biomarkers to the optimization of treatment and a number of quality-of-life enhancements for patients, computational biology methods and tools are shaping the field of MS diagnosis, management and treatment. The final goal in such a complex disease would be personalized medicine, i.e., providing healthcare services that are tailored to the individual patient, in accordance to the particular biology of their disease and the environmental factors to which they are subjected. The present review article summarizes the current knowledge on MS, modern computational biology and the impact of modern computational approaches of MS.

Small noncoding RNA profiling across cellular and biofluid compartments and their implications for multiple sclerosis immunopathology

Dysregulation of microRNAs (miRNAs), a type of small noncoding RNAs (sncRNAs), has frequently been associated with multiple sclerosis (MS). However, most studies have focused on peripheral blood, and few investigated other classes of sncRNAs. To address this, we analyzed all classes of sncRNAs in matching peripheral blood mononuclear cells, plasma, cerebrospinal fluid (CSF) cells, and cell-free CSF from MS patients and controls. We demonstrate widespread alterations of small nuclear (snRNA)–derived RNAs, small nucleolar-derived RNAs (sdRNAs), transfer RNA–derived fragments, and miRNAs, particularly in CSF cells. The striking contrast between the periphery and central nervous system and between relapse and remission phases of disease highlights the importance of sncRNA-mediated mechanisms in MS, in particular alternative splicing and mRNA translation.

Multiple sclerosis (MS) is a chronic inflammatory demyelinating disease affecting the central nervous system (CNS). Small non-coding RNAs (sncRNAs) and, in particular, microRNAs (miRNAs) have frequently been associated with MS. Here, we performed a comprehensive analysis of all classes of sncRNAs in matching samples of peripheral blood mononuclear cells (PBMCs), plasma, cerebrospinal fluid (CSF) cells, and cell-free CSF from relapsing-remitting (RRMS, n = 12 in relapse and n = 11 in remission) patients, secondary progressive (SPMS, n = 6) MS patients, and noninflammatory and inflammatory neurological disease controls (NINDC, n = 11; INDC, n = 5). We show widespread changes in miRNAs and sncRNA-derived fragments of small nuclear, nucleolar, and transfer RNAs. In CSF cells, 133 out of 133 and 115 out of 117 differentially expressed sncRNAs were increased in RRMS relapse compared to remission and RRMS compared to NINDC, respectively. In contrast, 65 out of 67 differentially expressed PBMC sncRNAs were decreased in RRMS compared to NINDC. The striking contrast between the periphery and CNS suggests that sncRNA-mediated mechanisms, including alternative splicing, RNA degradation, and mRNA translation, regulate the transcriptome of pathogenic cells primarily in the CNS target organ.

[The role of mitochondria in pathological mechanisms of innate immunity in multiple]

The review discusses the role of mitochondria in multiple sclerosis (MS). Previously, damage to the mitochondria was regarded as a manifestation of secondary damage to axons and neurons, i.e. as a marker of neurodenegation. Recently, the role of mitochondria in the early stages of MS development, when they could participate in the activation of innate immunity and trigger activation of autoimmune responses of acquired immunity, has been increasingly discussed. The role of polymorphism mitochondrial DNA changes in MS is discussed.

Identification of common key genes and pathways between type 1 diabetes and multiple sclerosis using transcriptome and interactome analysis

PURPOSE

Type 1 diabetes (T1D) and multiple sclerosis (MS) are classified as T cell-mediated autoimmune diseases. Although convergent evidence proposed common genetic architecture for autoimmune diseases, it remains a challenge to identify them. This study aimed to determine common gene signature and pathways in T1D and MS via systems biology approach.

METHODS

Gene expression profiles of peripheral blood mononuclear cells (PBMCs) and pancreatic-β cells in T1D as well as PBMCs and cerebrospinal fluid (CSF) in MS were analyzed in our previous published data, and differential expressed genes were integrated with protein-protein interactions data to construct Query-Query PPI (QQPPI) networks. In this study, QQPPI networks were further analyzed to investigate more central genes, functional modules and complexes shared in T1D and MS progression. Lastly, the interaction of common genes with drugs was also explored.

RESULTS

Several cytokines such as IL-23A, IL-32, IL-34, and IL-37 tend to be differentially expressed in both diseases. In addition, PSMA1, MYC, SRPK1, YBX1, HNRNPM, NF-κB2, IKBKE, RAC1, FN1, ARRB2, ESR1, HSP90AB1, and PPP1CA were common high central genes in QQPPI networks corresponding to each disease. Proteasome, spliceosome, immune responses, apoptosis, cellular communication/signaling transduction mechanism, interaction with environment, and activity of intercellular mediators were shared biological processes in T1D and MS. Finally, azathioprine, melatonin, resveratrol, and geldanamycin identified as prioritized drugs for the treatment of patients with T1D and MS.

CONCLUSIONS

This study represented novel key genes and pathways shared between T1D and MS, which may facilitate the identification of potential therapeutic targets in these diseases.

Not only cancer: the long non-coding RNA MALAT1 affects the repertoire of alternatively spliced transcripts and circular RNAs in multiple sclerosis

Long non-coding RNAs (lncRNAs) are post-transcriptional and epigenetic regulators, whose implication in neurodegenerative and autoimmune diseases remains poorly understood. We analyzed publicly available microarray data sets to identify dysregulated lncRNAs in multiple sclerosis (MS), a neuroinflammatory autoimmune disease. We found a consistent upregulation in MS of the lncRNA MALAT1 (2.7-fold increase; meta-analysis, P = 1.3 × 10-8; 190 cases, 182 controls), known to regulate alternative splicing (AS). We confirmed MALAT1 upregulation in two independent MS cohorts (1.5-fold increase; P < 0.01; 59 cases, 50 controls). We hence performed MALAT1 overexpression/knockdown in cell lines, demonstrating that its modulation impacts on endogenous expression of splicing factors (HNRNPF and HNRNPH1) and on AS of MS-associated genes (IL7R and SP140). Minigene-based splicing assays upon MALAT1 modulation recapitulated IL7R and SP140 isoform unbalances observed in patients. RNA-sequencing of MALAT1-knockdown Jurkat cells further highlighted MALAT1 role in splicing (approximately 1100 significantly-modulated AS events) and revealed its contribution to backsplicing (approximately 50 differentially expressed circular RNAs). Our study proposes a possible novel role for MALAT1 dysregulation and the consequent AS alteration in MS pathogenesis, based on anomalous splicing/backsplicing profiles of MS-relevant genes.

"Bridging the Gap" Everything that Could Have Been Avoided If We Had Applied Gender Medicine, Pharmacogenetics and Personalized Medicine in the Gender-Omics and Sex-Omics Era

Gender medicine is the first step of personalized medicine and patient-centred care, an essential development to achieve the standard goal of a holistic approach to patients and diseases. By addressing the interrelation and integration of biological markers (i.e., sex) with indicators of psychological/cultural behaviour (i.e., gender), gender medicine represents the crucial assumption for achieving the personalized health-care required in the third millennium. However, ‘sex’ and ‘gender’ are often misused as synonyms, leading to frequent misunderstandings in those who are not deeply involved in the field. Overall, we have to face the evidence that biological, genetic, epigenetic, psycho-social, cultural, and environmental factors mutually interact in defining sex/gender differences, and at the same time in establishing potential unwanted sex/gender disparities. Prioritizing the role of sex/gender in physiological and pathological processes is crucial in terms of efficient prevention, clinical signs’ identification, prognosis definition, and therapy optimization. In this regard, the omics-approach has become a powerful tool to identify sex/gender-specific disease markers, with potential benefits also in terms of socio-psychological wellbeing for each individual, and cost-effectiveness for National Healthcare systems. “Being a male or being a female” is indeed important from a health point of view and it is no longer possible to avoid “sex and gender lens” when approaching patients. Accordingly, personalized healthcare must be based on evidence from targeted research studies aimed at understanding how sex and gender influence health across the entire life span. The rapid development of genetic tools in the molecular medicine approaches and their impact in healthcare is an example of highly specialized applications that have moved from specialists to primary care providers (e.g., pharmacogenetic and pharmacogenomic applications in routine medical practice). Gender medicine needs to follow the same path and become an established medical approach. To face the genetic, molecular and pharmacological bases of the existing sex/gender gap by means of omics approaches will pave the way to the discovery and identification of novel drug-targets/therapeutic protocols, personalized laboratory tests and diagnostic procedures (sex/gender-omics). In this scenario, the aim of the present review is not to simply resume the state-of-the-art in the field, rather an opportunity to gain insights into gender medicine, spanning from molecular up to social and psychological stances. The description and critical discussion of some key selected multidisciplinary topics considered as paradigmatic of sex/gender differences and sex/gender inequalities will allow to draft and design strategies useful to fill the existing gap and move forward.

Aberrant expression of alternative splicing variants in multiple sclerosis - A systematic review

OBJECTIVE

Alternative splicing is an important form of RNA processing that affects nearly all human genes. The differential expression of specific transcript and protein isoforms holds the potential of novel biomarkers for complex diseases. In this systematic review, we compiled the existing literature on aberrant alternative splicing events in multiple sclerosis (MS).

METHODS

A systematic literature search in the PubMed database was carried out and supplemented by screening the reference lists of the identified articles. We selected only MS-related original research studies which compared the levels of different isoforms of human protein-coding genes. A narrative synthesis of the research findings was conducted. Additionally, we performed a case-control analysis using high-density transcriptome microarray data to reevaluate the genes that were examined in the reviewed studies.

RESULTS

A total of 160 records were screened. Of those, 36 studies from the last two decades were included. Most commonly, peripheral blood samples were analyzed (32 studies), and PCR-based techniques were usually employed (27 studies) for measuring the expression of selected genes. Two studies used an exploratory genome-wide approach. Overall, 27 alternatively spliced genes were investigated. Nine of these genes appeared in at least two studies (CD40, CFLAR, FOXP3, IFNAR2, IL7R, MOG, PTPRC, SP140 and TNFRSF1A). The microarray data analysis confirmed differential alternative pre-mRNA splicing for 19 genes.

CONCLUSIONS

An altered RNA processing of genes mediating immune signaling pathways has been repeatedly implicated in MS. The analysis of individual exon-level expression patterns is stimulated by the advancement of transcriptome profiling technologies. In particular, the examination of genes encoded in MS-associated genetic regions may provide important insights into the pathogenesis of the disease and help to identify new biomarkers.

Interpreting Non-coding Genetic Variation in Multiple Sclerosis Genome-Wide Associated Regions

Multiple sclerosis (MS) is the most common neurological disorder in young adults. Despite extensive studies, only a fraction of MS heritability has been explained, with association studies focusing primarily on protein-coding genes, essentially for the difficulty of interpreting non-coding features. However, non-coding RNAs (ncRNAs) and functional elements, such as super-enhancers (SE), are crucial regulators of many pathways and cellular mechanisms, and they have been implicated in a growing number of diseases. In this work, we searched for possible enrichments in non-coding elements at MS genome-wide associated loci, with the aim to highlight their possible involvement in the susceptibility to the disease. We first reconstructed the linkage disequilibrium (LD) structure of the Italian population using data of 727,478 single-nucleotide polymorphisms (SNPs) from 1,668 healthy individuals. The genomic coordinates of the obtained LD blocks were intersected with those of the top hits identified in previously published MS genome-wide association studies (GWAS). By a bootstrapping approach, we hence demonstrated a striking enrichment of non-coding elements, especially of circular RNAs (circRNAs) mapping in the 73 LD blocks harboring MS-associated SNPs. In particular, we found a total of 482 circRNAs (annotated in publicly available databases) vs. a mean of 194 ± 65 in the random sets of LD blocks, using 1,000 iterations. As a proof of concept of a possible functional relevance of this observation, we experimentally verified that the expression levels of a circRNA derived from an MS-associated locus, i.e., hsa_circ_0043813 from the STAT3 gene, can be modulated by the three genotypes at the disease-associated SNP. Finally, by evaluating RNA-seq data of two cell lines, SH-SY5Y and Jurkat cells, representing tissues relevant for MS, we identified 18 (two novel) circRNAs derived from MS-associated genes. In conclusion, this work showed for the first time that MS-GWAS top hits map in LD blocks enriched in circRNAs, suggesting circRNAs as possible novel contributors to the disease pathogenesis.

Changes in expression profiles of internal jugular vein wall and plasma protein levels in multiple sclerosis

Background

Multiple sclerosis (MS) is an inflammatory, demyelinating and degenerative disorder of the central nervous system (CNS). Several observations support interactions between vascular and neurodegenerative mechanisms in multiple sclerosis (MS). To investigate the contribution of the extracranial venous compartment, we analysed expression profiles of internal jugular vein (IJV), which drains blood from CNS, and related plasma protein levels.

Methods

We studied a group of MS patients (n = 19), screened by echo-color Doppler and magnetic resonance venography, who underwent surgical reconstruction of IJV for chronic cerebrospinal venous insufficiency (CCSVI). Microarray-based transcriptome analysis was conducted on specimens of IJV wall from MS patients and from subjects undergoing carotid endarterectomy, as controls. Protein levels were determined by multiplex assay in: i) jugular and peripheral plasma from 17 MS/CCSVI patients; ii) peripheral plasma from 60 progressive MS patients, after repeated sampling and iii) healthy individuals.

Results

Of the differentially expressed genes (≥ 2 fold-change, multiple testing correction, P < 0.05), the immune-related CD86 (8.5 fold-change, P = 0.002) emerged among the up regulated genes (N = 409). Several genes encoding HOX transcription factors and histones potentially regulated by blood flow, were overexpressed. Smooth muscle contraction and cell adhesion processes emerged among down regulated genes (N = 515), including the neuronal cell adhesion L1CAM as top scorer (5 fold-change, P = 5 × 10^− 4).

Repeated measurements in jugular/peripheral plasma and overtime in peripheral plasma showed conserved individual plasma patterns for immune-inflammatory (CCL13, CCL18) and adhesion (NCAM1, VAP1, SELL) proteins, despite significant variations overtime (SELL P < 0.0001). Both age and MS disease phenotypes were determinants of VAP1 plasma levels.

Data supported cerebral related-mechanisms regulating ANGPT1 levels, which were remarkably lower in jugular plasma and correlated in repeated assays but not between jugular/peripheral compartments.

Conclusions

This study provides for the first time expression patterns of the IJV wall, suggesting signatures of altered vascular mRNA profiles in MS disease also independently from CCSVI. The combined transcriptome-protein analysis provides intriguing links between IJV wall transcript alteration and plasma protein expression, thus highlighting proteins of interest for MS pathophysiology.

Electronic supplementary material

The online version of this article (10.1186/s10020-018-0043-4) contains supplementary material, which is available to authorized users.

The Characterization of GSDMB Splicing and Backsplicing Profiles Identifies Novel Isoforms and a Circular RNA That Are Dysregulated in Multiple Sclerosis

Abnormalities in alternative splicing (AS) are emerging as recurrent features in autoimmune diseases (AIDs). In particular, a growing body of evidence suggests the existence of a pathogenic association between a generalized defect in splicing regulatory genes and multiple sclerosis (MS). Moreover, several studies have documented an unbalance in alternatively-spliced isoforms in MS patients possibly contributing to the disease etiology. In this work, using a combination of PCR-based techniques (reverse-transcription (RT)-PCR, fluorescent-competitive, real-time, and digital RT-PCR assays), we investigated the alternatively-spliced gene encoding Gasdermin B, GSDMB, which was repeatedly associated with susceptibility to asthma and AIDs. The in-depth characterization of GSDMB AS and backsplicing profiles led us to the identification of an exonic circular RNA (ecircRNA) as well as of novel GSDMB in-frame and out-of-frame isoforms. The non-productive splicing variants were shown to be downregulated by the nonsense-mediated mRNA decay (NMD) in human cell lines, suggesting that GSDMB levels are significantly modulated by NMD. Importantly, both AS isoforms and the identified ecircRNA were significantly dysregulated in peripheral blood mononuclear cells of relapsing-remitting MS patients compared to controls, further supporting the notion that aberrant RNA metabolism is a characteristic feature of the disease.

Systems Biomedicine of Rabies Delineates the Affected Signaling Pathways

The prototypical neurotropic virus, rabies, is a member of the Rhabdoviridae family that causes lethal encephalomyelitis. Although there have been a plethora of studies investigating the etiological mechanism of the rabies virus and many precautionary methods have been implemented to avert the disease outbreak over the last century, the disease has surprisingly no definite remedy at its late stages. The psychological symptoms and the underlying etiology, as well as the rare survival rate from rabies encephalitis, has still remained a mystery. We, therefore, undertook a systems biomedicine approach to identify the network of gene products implicated in rabies. This was done by meta-analyzing whole-transcriptome microarray datasets of the CNS infected by strain CVS-11, and integrating them with interactome data using computational and statistical methods. We first determined the differentially expressed genes (DEGs) in each study and horizontally integrated the results at the mRNA and microRNA levels separately. A total of 61 seed genes involved in signal propagation system were obtained by means of unifying mRNA and microRNA detected integrated DEGs. We then reconstructed a refined protein–protein interaction network (PPIN) of infected cells to elucidate the rabies-implicated signal transduction network (RISN). To validate our findings, we confirmed differential expression of randomly selected genes in the network using Real-time PCR. In conclusion, the identification of seed genes and their network neighborhood within the refined PPIN can be useful for demonstrating signaling pathways including interferon circumvent, toward proliferation and survival, and neuropathological clue, explaining the intricate underlying molecular neuropathology of rabies infection and thus rendered a molecular framework for predicting potential drug targets.