Transcriptomic meta-analysis of multiple sclerosis and its experimental models

Human neural cell type-specific extracellular vesicle proteome defines disease-related molecules associated with activated astrocytes in Alzheimer's disease brain

In neurodegenerative diseases, extracellular vesicles (EVs) transfer pathogenic molecules and are consequently involved in disease progression. We have investigated the proteomic profiles of EVs that were isolated from four different human-induced pluripotent stem cell-derived neural cell types (excitatory neurons, astrocytes, microglia-like cells, and oligodendrocyte-like cells). Novel cell type-specific EV protein markers were then identified for the excitatory neurons (ATP1A3, NCAM1), astrocytes (LRP1, ITGA6), microglia-like cells (ITGAM, LCP1), and oligodendrocyte-like cells (LAMP2, FTH1), as well as 16 pan-EV marker candidates, including integrins and annexins. To further demonstrate how cell-type-specific EVs may be involved in Alzheimer's disease (AD), we performed protein co-expression network analysis and conducted cell type assessments for the proteomes of brain-derived EVs from the control, mild cognitive impairment, and AD cases. A protein module enriched in astrocyte-specific EV markers was most significantly associated with the AD pathology and cognitive impairment, suggesting an important role in AD progression. The hub protein from this module, integrin-β1 (ITGB1), was found to be significantly elevated in astrocyte-specific EVs enriched from the total brain-derived AD EVs and associated with the brain β-amyloid and tau load in independent cohorts. Thus, our study provides a featured framework and rich resource for the future analyses of EV functions in neurodegenerative diseases in a cell type-specific manner.

Integrated analysis of human transcriptome data for Rett syndrome finds a network of involved genes

OBJECTIVES

Rett syndrome (RTT) is a rare disorder causing severe intellectual and physical disability. The cause is a mutation in the gene coding for the methyl-CpG binding protein 2 (MECP2), a multifunctional regulator protein. Purpose of the study was integration and investigation of multiple gene expression profiles in human cells with impaired MECP2 gene to obtain a robust, data-driven insight in molecular disease mechanisms.

METHODS

Information about changed gene expression was extracted from five previously published studies, integrated and the resulting differentially expressed genes were analysed using overrepresentation analysis of biological pathways and gene ontology, and network analysis.

RESULTS

We identified a set of genes, which are significantly changed not in all but several transcriptomics datasets and were not mentioned in the context of RTT before. We found that these genes are involved in several processes and molecular pathways known to be affected in RTT. Integrating transcription factors we identified a possible link how MECP2 regulates cytoskeleton organisation via MEF2C and CAPG.

CONCLUSIONS

Integrative analysis of omics data and prior knowledge databases is a powerful approach to identify links between mutation and phenotype especially in rare disease research where little data is available.

Revisiting the Pathoetiology of Multiple Sclerosis: Has the Tail Been Wagging the Mouse?

Multiple Sclerosis (MS) is traditionally considered an autoimmune-mediated demyelinating disease, the pathoetiology of which is unknown. However, the key question remains whether autoimmunity is the initiator of the disease (outside-in) or the consequence of a slow and as yet uncharacterized cytodegeneration (oligodendrocytosis), which leads to a subsequent immune response (inside-out). Experimental autoimmune encephalomyelitis has been used to model the later stages of MS during which the autoimmune involvement predominates. In contrast, the cuprizone (CPZ) model is used to model early stages of the disease during which oligodendrocytosis and demyelination predominate and are hypothesized to precede subsequent immune involvement in MS. Recent studies combining a boost, or protection, to the immune system with disruption of the blood brain barrier have shown CPZ-induced oligodendrocytosis with a subsequent immune response. In this Perspective, we review these recent advances and discuss the likelihood of an inside-out vs. an outside-in pathoetiology of MS.

A pathogenic and clonally expanded B cell transcriptome in active multiple sclerosis

B cells serve as a key weapon against infectious diseases. They also contribute to multiple autoimmune diseases, including multiple sclerosis (MS) where depletion of B cells is a highly effective therapy. We describe a comprehensive profile of central nervous system (CNS)-specific transcriptional B cell phenotypes in MS at single-cell resolution with paired immune repertoires. We reveal a polyclonal immunoglobulin M (IgM) and IgG1 cerebrospinal fluid B cell expansion polarized toward an inflammatory, memory and plasmablast/plasma cell phenotype, with differential up-regulation of specific proinflammatory pathways. We did not find evidence that CNS B cells harbor a neurotropic virus. These data support the targeting of activated resident B cells in the CNS as a potentially effective strategy for control of treatment-resistant chronic disease.

Central nervous system B cells have several potential roles in multiple sclerosis (MS): secretors of proinflammatory cytokines and chemokines, presenters of autoantigens to T cells, producers of pathogenic antibodies, and reservoirs for viruses that trigger demyelination. To interrogate these roles, single-cell RNA sequencing (scRNA-Seq) was performed on paired cerebrospinal fluid (CSF) and blood from subjects with relapsing-remitting MS (RRMS; n = 12), other neurologic diseases (ONDs; n = 1), and healthy controls (HCs; n = 3). Single-cell immunoglobulin sequencing (scIg-Seq) was performed on a subset of these subjects and additional RRMS (n = 4), clinically isolated syndrome (n = 2), and OND (n = 2) subjects. Further, paired CSF and blood B cell subsets (RRMS; n = 7) were isolated using fluorescence activated cell sorting for bulk RNA sequencing (RNA-Seq). Independent analyses across technologies demonstrated that nuclear factor kappa B (NF-κB) and cholesterol biosynthesis pathways were activated, and specific cytokine and chemokine receptors were up-regulated in CSF memory B cells. Further, SMAD/TGF-β1 signaling was down-regulated in CSF plasmablasts/plasma cells. Clonally expanded, somatically hypermutated IgM+ and IgG1+ CSF B cells were associated with inflammation, blood–brain barrier breakdown, and intrathecal Ig synthesis. While we identified memory B cells and plasmablast/plasma cells with highly similar Ig heavy-chain sequences across MS subjects, similarities were also identified with ONDs and HCs. No viral transcripts, including from Epstein–Barr virus, were detected. Our findings support the hypothesis that in MS, CSF B cells are driven to an inflammatory and clonally expanded memory and plasmablast/plasma cell phenotype.

Bassoon proteinopathy drives neurodegeneration in multiple sclerosis

Multiple sclerosis (MS) is characterized by inflammatory insults that drive neuroaxonal injury. However, knowledge about neuron-intrinsic responses to inflammation is limited. By leveraging neuron-specific messenger RNA profiling, we found that neuroinflammation leads to induction and toxic accumulation of the synaptic protein bassoon (Bsn) in the neuronal somata of mice and patients with MS. Neuronal overexpression of Bsn in flies resulted in reduction of lifespan, while genetic disruption of Bsn protected mice from inflammation-induced neuroaxonal injury. Notably, pharmacological proteasome activation boosted the clearance of accumulated Bsn and enhanced neuronal survival. Our study demonstrates that neuroinflammation initiates toxic protein accumulation in neuronal somata and advocates proteasome activation as a potential remedy.

Delayed Astrogliosis Associated with Reduced M1 Microglia Activation in Matrix Metalloproteinase 12 Knockout Mice during Theiler's Murine Encephalomyelitis

Theiler’s murine encephalomyelitis (TME) represents a versatile animal model for studying the pathogenesis of demyelinating diseases such as multiple sclerosis. Hallmarks of TME are demyelination, astrogliosis, as well as inflammation. Previous studies showed that matrix metalloproteinase 12 knockout (Mmp12^−/−) mice display an ameliorated clinical course associated with reduced demyelination. The present study aims to elucidate the impact of MMP12 deficiency in TME with special emphasis on astrogliosis, macrophages infiltrating the central nervous system (CNS), and the phenotype of microglia/macrophages (M1 or M2). SJL wild-type and Mmp12^−/− mice were infected with TME virus (TMEV) or vehicle (mock) and euthanized at 28 and 98 days post infection (dpi). Immunohistochemistry or immunofluorescence of cervical and thoracic spinal cord for detecting glial fibrillary acidic protein (GFAP), ionized calcium-binding adaptor molecule 1 (Iba1), chemokine receptor 2 (CCR2), CD107b, CD16/32, and arginase I was performed and quantitatively evaluated. Statistical analyses included the Kruskal–Wallis test followed by Mann–Whitney U post hoc tests. TMEV-infected Mmp12^−/− mice showed transiently reduced astrogliosis in association with a strong trend (p = 0.051) for a reduced density of activated/reactive microglia/macrophages compared with wild-type mice at 28 dpi. As astrocytes are an important source of cytokine production, including proinflammatory cytokines triggering or activating phagocytes, the origin of intralesional microglia/macrophages as well as their phenotype were determined. Only few phagocytes in wild-type and Mmp12^−/− mice expressed CCR2, indicating that the majority of phagocytes are represented by microglia. In parallel to the reduced density of activated/reactive microglia at 98 dpi, TMEV-infected Mmp12^−/− showed a trend (p = 0.073) for a reduced density of M1 (CD16/32- and CD107b-positive) microglia, while no difference regarding the density of M2 (arginase I- and CD107b-positive) cells was observed. However, a dominance of M1 cells was detected in the spinal cord of TMEV-infected mice at all time points. Reduced astrogliosis in Mmp12^−/− mice was associated with a reduced density of activated/reactive microglia and a trend for a reduced density of M1 cells. This indicates that MMP12 plays an important role in microglia activation, polarization, and migration as well as astrogliosis and microglia/astrocyte interaction.

Comparison of Reported Spinal Cord Lesions in Progressive Multiple Sclerosis with Theiler's Murine Encephalomyelitis Virus Induced Demyelinating Disease

BACKGROUND

Spinal cord (SC) lesions in Theiler's murine encephalomyelitis virus induced demyelinating disease (TMEV-IDD) resemble important features of brain lesions in progressive multiple sclerosis (MS) including inflammation, demyelination, and axonal damage. The aim of the present study was a comparison of SC lesions in MS and TMEV-IDD focusing on spatial and temporal distribution of demyelination, inflammation, SC atrophy (SCA), and axonal degeneration/loss in major descending motor pathways.

METHODS

TMEV and mock-infected mice were investigated clinically once a week. SC tissue was collected at 42, 98, 147, and 196 days post infection, and investigated using hematoxylin and eosin (HE) staining, immunohistochemistry targeting myelin basic protein (demyelination), Mac3 (microglia/macrophages), phosphorylated neurofilaments (axonal damage) and transmission electron microscopy.

RESULTS

Demyelination prevailed in SC white matter in TMEV-IDD, contrasting a predominant gray matter involvement in MS. TMEV-infected mice revealed a significant loss of axons similar to MS. Ultrastructural analysis in TMEV-IDD revealed denuded axons, degenerative myelin changes, axonal degeneration, as well as remyelination. SCA is a consistent finding in the SC of MS patients and was also detected at a late time point in TMEV-IDD.

CONCLUSION

This comparative study further indicates the suitability of TMEV-IDD as animal model also for the investigation of progressive SC lesions in MS.

Coming of Age for Autotaxin and Lysophosphatidate Signaling: Clinical Applications for Preventing, Detecting and Targeting Tumor-Promoting Inflammation

A quarter-century after the discovery of autotaxin in cell culture, the autotaxin-lysophosphatidate (LPA)-lipid phosphate phosphatase axis is now a promising clinical target for treating chronic inflammatory conditions, mitigating fibrosis progression, and improving the efficacy of existing cancer chemotherapies and radiotherapy. Nearly half of the literature on this axis has been published during the last five years. In cancer biology, LPA signaling is increasingly being recognized as a central mediator of the progression of chronic inflammation in the establishment of a tumor microenvironment which promotes cancer growth, immune evasion, metastasis, and treatment resistance. In this review, we will summarize recent advances made in understanding LPA signaling with respect to chronic inflammation and cancer. We will also provide perspectives on the applications of inhibitors of LPA signaling in preventing cancer initiation, as adjuncts extending the efficacy of current cancer treatments by blocking inflammation caused by either the cancer or the cancer therapy itself, and by disruption of the tumor microenvironment. Overall, LPA, a simple molecule that mediates a plethora of biological effects, can be targeted at its levels of production by autotaxin, LPA receptors or through LPA degradation by lipid phosphate phosphatases. Drugs for these applications will soon be entering clinical practice.

Canine dorsal root ganglia satellite glial cells represent an exceptional cell population with astrocytic and oligodendrocytic properties

Dogs can be used as a translational animal model to close the gap between basic discoveries in rodents and clinical trials in humans. The present study compared the species-specific properties of satellite glial cells (SGCs) of canine and murine dorsal root ganglia (DRG) in situ and in vitro using light microscopy, electron microscopy, and immunostainings. The in situ expression of CNPase, GFAP, and glutamine synthetase (GS) has also been investigated in simian SGCs. In situ, most canine SGCs (>80%) expressed the neural progenitor cell markers nestin and Sox2. CNPase and GFAP were found in most canine and simian but not murine SGCs. GS was detected in 94% of simian and 71% of murine SGCs, whereas only 44% of canine SGCs expressed GS. In vitro, most canine (>84%) and murine (>96%) SGCs expressed CNPase, whereas GFAP expression was differentially affected by culture conditions and varied between 10% and 40%. However, GFAP expression was induced by bone morphogenetic protein 4 in SGCs of both species. Interestingly, canine SGCs also stimulated neurite formation of DRG neurons. These findings indicate that SGCs represent an exceptional, intermediate glial cell population with phenotypical characteristics of oligodendrocytes and astrocytes and might possess intrinsic regenerative capabilities in vivo.

Morphologic, phenotypic, and transcriptomic characterization of classically and alternatively activated canine blood-derived macrophages in vitro

Macrophages are a heterogeneous cell population playing a pivotal role in tissue homeostasis and inflammation, and their phenotype strongly depends on the micromilieu. Despite its increasing importance as a translational animal model for human diseases, there is a considerable gap of knowledge with respect to macrophage polarization in dogs. The present study comprehensively investigated the morphologic, phenotypic, and transcriptomic characteristics of unstimulated (M0), M1- (GM-CSF, LPS, IFNγ-stimulated) and M2- (M-CSF, IL-4-stimulated)-polarized canine blood-derived macrophages in vitro. Scanning electron microscopy revealed distinct morphologies of polarized macrophages with formation of multinucleated cells in M2-macrophages, while immunofluorescence employing literature-based prototype-antibodies against CD16, CD32, iNOS, MHC class II (M1-markers), CD163, CD206, and arginase-1 (M2-markers) demonstrated that only CD206 was able to discriminate M2-macrophages from both other phenotypes, highlighting this molecule as a promising marker for canine M2-macrophages. Global microarray analysis revealed profound changes in the transcriptome of polarized canine macrophages. Functional analysis pointed out that M1-polarization was associated with biological processes such as "respiratory burst", whereas M2-polarization was associated with processes such as "mitosis". Literature-based marker gene selection revealed only minor overlaps in the gene sets of the dog compared to prototype markers of murine and human macrophages. Biomarker selection using supervised clustering suggested latexin (LXN) and membrane-spanning 4-domains, subfamily A, member 2 (MS4A2) to be the most powerful predicting biomarkers for canine M1- and M2-macrophages, respectively. Immunofluorescence for both markers demonstrated expression of both proteins by macrophages in vitro but failed to reveal differences between canine M1 and M2-macrophages. The present study provides a solid basis for future studies upon the role of macrophage polarization in spontaneous diseases of the dog, a species that has emerging importance for translational research.

Microarray-Based Gene Expression Analysis for Veterinary Pathologists: A Review

High-throughput, genome-wide transcriptome analysis is now commonly used in all fields of life science research and is on the cusp of medical and veterinary diagnostic application. Transcriptomic methods such as microarrays and next-generation sequencing generate enormous amounts of data. The pathogenetic expertise acquired from understanding of general pathology provides veterinary pathologists with a profound background, which is essential in translating transcriptomic data into meaningful biological knowledge, thereby leading to a better understanding of underlying disease mechanisms. The scientific literature concerning high-throughput data-mining techniques usually addresses mathematicians or computer scientists as the target audience. In contrast, the present review provides the reader with a clear and systematic basis from a veterinary pathologist's perspective. Therefore, the aims are (1) to introduce the reader to the necessary methodological background; (2) to introduce the sequential steps commonly performed in a microarray analysis including quality control, annotation, normalization, selection of differentially expressed genes, clustering, gene ontology and pathway analysis, analysis of manually selected genes, and biomarker discovery; and (3) to provide references to publically available and user-friendly software suites. In summary, the data analysis methods presented within this review will enable veterinary pathologists to analyze high-throughput transcriptome data obtained from their own experiments, supplemental data that accompany scientific publications, or public repositories in order to obtain a more in-depth insight into underlying disease mechanisms.

MicroRNA-142 regulates inflammation and T cell differentiation in an animal model of multiple sclerosis

Background

MicroRNAs have emerged as an important class of modulators of gene expression. These molecules influence protein synthesis through translational repression or degradation of mRNA transcripts. Herein, we investigated the potential role of miR-142a isoforms, miR-142a-3p and miR-142a-5p, in the context of autoimmune neuroinflammation.

Methods

The expression levels of two mature isoforms of miR-142 were measured in the brains of patients with multiple sclerosis (MS) and the CNS tissues from mice with experimental autoimmune encephalomyelitis (EAE), an animal model of MS. Expression analyses were also performed in mitogen and antigen-stimulated splenocytes, as well as macrophages and astrocytes using real-time RT-PCR. The role of the mature miRNAs was then investigated in T cell differentiation by transfection of CD4⁺ T cells, followed by flow cytometric analysis of intracellular cytokines. Luciferase assays using vectors containing the 3′UTR of predicted targets were performed to confirm the interaction of miRNA sequences with transcripts. Expression of targets were then analyzed in activated splenocytes and MS/EAE tissues.

Results

Expression of miR-142-5p was significantly increased in the frontal white matter from MS patients compared with white matter from non-MS controls. Likewise, expression levels of miR-142a-5p and miR-142a-3p showed significant upregulation in the spinal cords of EAE mice at days 15 and 25 post disease induction. Splenocytes stimulated with myelin oligodendrocyte glycoprotein (MOG) peptide or anti-CD3/anti-CD28 antibodies showed upregulation of miR-142a-5p and miR-142a-3p isoforms, whereas stimulated bone marrow-derived macrophages and primary astrocytes did not show any significant changes in miRNA expression levels. miR-142a-5p overexpression in activated lymphocytes shifted the pattern of T cell differentiation towards Th1 cells. Luciferase assays revealed SOCS1 and TGFBR1 as direct targets of miR-142a-5p and miR-142a-3p, respectively, and overexpression of miRNA mimic sequences suppressed the expression of these target transcripts in lymphocytes. SOCS1 levels were also diminished in MS white matter and EAE spinal cords.

Conclusions

Our findings suggest that increased expression of miR-142 isoforms might be involved in the pathogenesis of autoimmune neuroinflammation by influencing T cell differentiation, and this effect could be mediated by interaction of miR-142 isoforms with SOCS1 and TGFBR-1 transcripts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12974-017-0832-7) contains supplementary material, which is available to authorized users.

Persistent Morbillivirus Infection Leads to Altered Cortactin Distribution in Histiocytic Sarcoma Cells with Decreased Cellular Migration Capacity

Histiocytic sarcomas represent rare but fatal neoplasms in humans. Based on the absence of a commercially available human histiocytic sarcoma cell line the frequently affected dog displays a suitable translational model. Canine distemper virus, closely related to measles virus, is a highly promising candidate for oncolytic virotherapy. Therapeutic failures in patients are mostly associated with tumour invasion and metastasis often induced by misdirected cytoskeletal protein activities. Thus, the impact of persistent canine distemper virus infection on the cytoskeletal protein cortactin, which is frequently overexpressed in human cancers with poor prognosis, was investigated in vitro in a canine histiocytic sarcoma cell line (DH82). Though phagocytic activity, proliferation and apoptotic rate were unaltered, a significantly reduced migration activity compared to controls (6 hours and 1 day after seeding) accompanied by a decreased number of cortactin mRNA transcripts (1 day) was detected. Furthermore, persistently canine distemper virus infected DH82 cells showed a predominant diffuse intracytoplasmic cortactin distribution at 6 hours and 1 day compared to controls with a prominent membranous expression pattern (p ≤ 0.05). Summarized, persistent canine distemper virus infection induces reduced tumour cell migration associated with an altered intracellular cortactin distribution, indicating cytoskeletal changes as one of the major pathways of virus-associated inhibition of tumour spread.

Reduced angiogenic gene expression in morbillivirus-triggered oncolysis in a translational model for histiocytic sarcoma

Histiocytic sarcoma represents a rare malignant tumour with a short survival time, indicating the need of novel treatment strategies including oncolytic virotherapy. The underlying molecular mechanisms of viral oncolysis are largely unknown. As cancer in companion animals shares striking similarities with human counterparts, we chose a permanent canine histiocytic sarcoma cell line (DH82 cells) to identify global transcriptome changes following infection with canine distemper virus (CDV), a paramyxovirus closely related to human measles virus. Microarray analysis identified 3054 differentially expressed probe sets (DEPs), encoding for 892 up‐ and 869 down‐regulated unique canine genes, respectively, in DH82 cells persistently infected with the vaccine strain Onderstepoort of CDV (DH82‐Ond‐pi), compared to non‐infected DH82 cells. Up‐regulated genes were predominantly related to immune processes, as demonstrated by functional enrichment analysis. Moreover, there was substantial enrichment of genes characteristic for classically activated M1 and alternatively activated M2 macrophages in DH82‐Ond‐pi; however, significant polarization into either of both categories was lacking. ‘Angiogenesis’ was the dominant enriched functional term for the down‐regulated genes, highlighting decreased blood vessel generation as a potential mechanism of paramyxovirus‐induced oncolysis in DH82 cells. The anti‐angiogenic effect of infection was verified by immunohistochemistry, which revealed a lower blood vessel density in an in vivo mouse model, xenotransplanted with DH82‐Ond‐pi, compared to mice transplanted with non‐infected DH82 cells. Reduction in angiogenesis appears to be an important oncolytic mechanism of CDV in DH82 cells, suggesting that similar mechanisms might account for human histiocytic sarcoma and maybe other tumours in conjunction with measles virus.

Extracellular cues influencing oligodendrocyte differentiation and (re)myelination

There is an increasing number of neurologic disorders found to be associated with loss and/or dysfunction of the CNS myelin sheath, ranging from the classic demyelinating disease, multiple sclerosis, through CNS injury, to neuropsychiatric diseases. The disabling burden of these diseases has sparked a growing interest in gaining a better understanding of the molecular mechanisms regulating the differentiation of the myelinating cells of the CNS, oligodendrocytes (OLGs), and the process of (re)myelination. In this context, the importance of the extracellular milieu is becoming increasingly recognized. Under pathological conditions, changes in inhibitory as well as permissive/promotional cues are thought to lead to an overall extracellular environment that is obstructive for the regeneration of the myelin sheath. Given the general view that remyelination is, even though limited in human, a natural response to demyelination, targeting pathologically 'dysregulated' extracellular cues and their downstream pathways is regarded as a promising approach toward the enhancement of remyelination by endogenous (or if necessary transplanted) OLG progenitor cells. In this review, we will introduce the extracellular cues that have been implicated in the modulation of (re)myelination. These cues can be soluble, part of the extracellular matrix (ECM) or mediators of cell-cell interactions. Their inhibitory and permissive/promotional roles with regard to remyelination as well as their potential for therapeutic intervention will be discussed.

Characterisation of Transcriptional Changes in the Spinal Cord of the Progressive Experimental Autoimmune Encephalomyelitis Biozzi ABH Mouse Model by RNA Sequencing

Multiple sclerosis (MS) is a debilitating immune-mediated neurological disorder affecting young adults. MS is primarily relapsing-remitting, but neurodegeneration and disability accumulate from disease onset. The most commonly used mouse MS models exhibit a monophasic immune response with fast accumulation of neurological damage that does not allow the study of progressive neurodegeneration. The chronic relapsing and secondary progressive EAE (pEAE) Biozzi ABH mouse model of MS exhibits a reproducible relapsing-remitting disease course that slowly accumulates permanent neurological deficit and develops a post-relapsing progressive disease that permits the study of demyelination and neurodegeneration. RNA sequencing (RNAseq) was used to explore global gene expression in the pEAE Biozzi ABH mouse. Spinal cord tissue RNA from pEAE Biozzi ABH mice and healthy age-matched controls was sequenced. 2,072 genes were differentially expressed (q<0.05) from which 1,397 were significantly upregulated and 675 were significantly downregulated. This hypothesis-free investigation characterised the genomic changes that describe the pEAE mouse model. The differentially expressed genes revealed a persistent immunoreactant phenotype, combined with downregulation of the cholesterol biosynthesis superpathway and the LXR/RXR activation pathway. Genes differentially expressed include the myelination genes Slc17a7, Ugt8A and Opalin, the neuroprotective genes Sprr1A, Osm and Wisp2, as well as genes identified as MS risk factors, including RGs14 and Scap2. Novel genes with unestablished roles in EAE or MS were also identified. The identification of differentially expressed novel genes and genes involved in MS pathology, opens the door to their functional study in the pEAE mouse model which recapitulates some of the important clinical features of progressive MS.

Meta-Analysis of Gene Expression Patterns in Animal Models of Prenatal Alcohol Exposure Suggests Role for Protein Synthesis Inhibition and Chromatin Remodeling

BACKGROUND

Prenatal alcohol exposure (PAE) can result in an array of morphological, behavioral, and neurobiological deficits that can range in their severity. Despite extensive research in the field and a significant progress made, especially in understanding the range of possible malformations and neurobehavioral abnormalities, the molecular mechanisms of alcohol responses in development are still not well understood. There have been multiple transcriptomic studies looking at the changes in gene expression after PAE in animal models; however, there is a limited apparent consensus among the reported findings. In an effort to address this issue, we performed a comprehensive re-analysis and meta-analysis of all suitable, publically available expression data sets.

METHODS

We assembled 10 microarray data sets of gene expression after PAE in mouse and rat models consisting of samples from a total of 63 ethanol (EtOH)-exposed and 80 control animals. We re-analyzed each data set for differential expression and then used the results to perform meta-analyses considering all data sets together or grouping them by time or duration of exposure (pre- and postnatal, acute and chronic, respectively). We performed network and Gene Ontology enrichment analysis to further characterize the identified signatures.

RESULTS

For each subanalysis, we identified signatures of differential expressed genes that show support from multiple studies. Overall, the changes in gene expression were more extensive after acute EtOH treatment during prenatal development than in other models. Considering the analysis of all the data together, we identified a robust core signature of 104 genes down-regulated after PAE, with no up-regulated genes. Functional analysis reveals over representation of genes involved in protein synthesis, mRNA splicing, and chromatin organization.

CONCLUSIONS

Our meta-analysis shows that existing studies, despite superficial dissimilarity in findings, share features that allow us to identify a common core signature set of transcriptome changes in PAE. This is an important step to identifying the biological processes that underlie the etiology of fetal alcohol spectrum disorders.

Functional inflammatory profiles distinguish myelin-reactive T cells from patients with multiple sclerosis

Myelin-reactive T cells have been identified in patients with multiple sclerosis (MS) and healthy subjects with comparable frequencies, but the contribution of these autoreactive T cells to disease pathology remains unknown. A total of 13,324 T cell libraries generated from blood of 23 patients and 22 healthy controls were interrogated for reactivity to myelin antigens. Libraries derived from CCR6(+) myelin-reactive T cells from patients with MS exhibited significantly enhanced production of interferon-γ (IFN-γ), interleukin-17 (IL-17), and granulocyte-macrophage colony-stimulating factor (GM-CSF) compared to healthy controls. Single-cell clones isolated by major histocompatibility complex/peptide tetramers from CCR6(+) T cell libraries also secreted more proinflammatory cytokines, whereas clones isolated from controls secreted more IL-10. The transcriptomes of myelin-specific CCR6(+) T cells from patients with MS were distinct from those derived from healthy controls and, notably, were enriched in T helper cell 17 (TH17)-induced experimental autoimmune encephalitis gene signatures, and gene signatures derived from TH17 cells isolated other human autoimmune diseases. These data, although not causal, imply that functional differences between antigen-specific T cells from MS and healthy controls are fundamental to disease development and support the notion that IL-10 production from myelin-reactive T cells may act to limit disease progression or even pathogenesis.

The Autotaxin-Lysophosphatidic Acid Axis Modulates Histone Acetylation and Gene Expression during Oligodendrocyte Differentiation

During development, oligodendrocytes (OLGs), the myelinating cells of the CNS, undergo a stepwise progression during which OLG progenitors, specified from neural stem/progenitor cells, differentiate into fully mature myelinating OLGs. This progression along the OLG lineage is characterized by well synchronized changes in morphology and gene expression patterns. The latter have been found to be particularly critical during the early stages of the lineage, and they have been well described to be regulated by epigenetic mechanisms, especially by the activity of the histone deacetylases HDAC1 and HDAC2. The data presented here identify the extracellular factor autotaxin (ATX) as a novel upstream signal modulating HDAC1/2 activity and gene expression in cells of the OLG lineage. Using the zebrafish as an in vivo model system as well as rodent primary OLG cultures, this functional property of ATX was found to be mediated by its lysophospholipase D (lysoPLD) activity, which has been well characterized to generate the lipid signaling molecule lysophosphatidic acid (LPA). More specifically, the lysoPLD activity of ATX was found to modulate HDAC1/2 regulated gene expression during a time window coinciding with the transition from OLG progenitor to early differentiating OLG. In contrast, HDAC1/2 regulated gene expression during the transition from neural stem/progenitor to OLG progenitor appeared unaffected by ATX and its lysoPLD activity. Thus, together, our data suggest that an ATX-LPA-HDAC1/2 axis regulates OLG differentiation specifically during the transition from OLG progenitor to early differentiating OLG and via a molecular mechanism that is evolutionarily conserved from at least zebrafish to rodent.

SIGNIFICANCE STATEMENT

The formation of the axon insulating and supporting myelin sheath by differentiating oligodendrocytes (OLGs) in the CNS is considered an essential step during vertebrate development. In addition, loss and/or dysfunction of the myelin sheath has been associated with a variety of neurologic diseases in which repair is limited, despite the presence of progenitor cells with the potential to differentiate into myelinating OLGs. This study characterizes the autotaxin-lysophosphatidic acid signaling axis as a modulator of OLG differentiation in vivo in the developing zebrafish and in vitro in rodent OLGs in culture. These findings provide novel insight into the regulation of developmental myelination, and they are likely to lead to advancing studies related to the stimulation of myelin repair under pathologic conditions.

BubbleGUM: automatic extraction of phenotype molecular signatures and comprehensive visualization of multiple Gene Set Enrichment Analyses

Background

Recent advances in the analysis of high-throughput expression data have led to the development of tools that scaled-up their focus from single-gene to gene set level. For example, the popular Gene Set Enrichment Analysis (GSEA) algorithm can detect moderate but coordinated expression changes of groups of presumably related genes between pairs of experimental conditions. This considerably improves extraction of information from high-throughput gene expression data. However, although many gene sets covering a large panel of biological fields are available in public databases, the ability to generate home-made gene sets relevant to one’s biological question is crucial but remains a substantial challenge to most biologists lacking statistic or bioinformatic expertise. This is all the more the case when attempting to define a gene set specific of one condition compared to many other ones. Thus, there is a crucial need for an easy-to-use software for generation of relevant home-made gene sets from complex datasets, their use in GSEA, and the correction of the results when applied to multiple comparisons of many experimental conditions.

Result

We developed BubbleGUM (GSEA Unlimited Map), a tool that allows to automatically extract molecular signatures from transcriptomic data and perform exhaustive GSEA with multiple testing correction. One original feature of BubbleGUM notably resides in its capacity to integrate and compare numerous GSEA results into an easy-to-grasp graphical representation. We applied our method to generate transcriptomic fingerprints for murine cell types and to assess their enrichments in human cell types. This analysis allowed us to confirm homologies between mouse and human immunocytes.

Conclusions

BubbleGUM is an open-source software that allows to automatically generate molecular signatures out of complex expression datasets and to assess directly their enrichment by GSEA on independent datasets. Enrichments are displayed in a graphical output that helps interpreting the results. This innovative methodology has recently been used to answer important questions in functional genomics, such as the degree of similarities between microarray datasets from different laboratories or with different experimental models or clinical cohorts. BubbleGUM is executable through an intuitive interface so that both bioinformaticians and biologists can use it. It is available at http://www.ciml.univ-mrs.fr/applications/BubbleGUM/index.html.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2012-4) contains supplementary material, which is available to authorized users.

Neurodegeneration in Autoimmune Optic Neuritis Is Associated with Altered APP Cleavage in Neurons and Up-Regulation of p53

Multiple Sclerosis (MS) is a chronic autoimmune inflammatory disease of the central nervous system (CNS). Histopathological and radiological analysis revealed that neurodegeneration occurs early in the disease course. However, the pathological mechanisms involved in neurodegeneration are poorly understood. Myelin oligodendrocyte glycoprotein (MOG)-induced experimental autoimmune encephalomyelitis (EAE) in Brown Norway rats (BN-rats) is a well-established animal model, especially of the neurodegenerative aspects of MS. Previous studies in this animal model indicated that loss of retinal ganglion cells (RGCs), the neurons that form the axons of the optic nerve, occurs in the preclinical phase of the disease and is in part independent of overt histopathological changes of the optic nerve. Therefore, the aim of this study was to identify genes which are involved in neuronal cell loss at different disease stages of EAE. Furthermore, genes that are highly specific for autoimmune-driven neurodegeneration were compared to those regulated in RGCs after optic nerve axotomy at corresponding time points. Using laser capture micro dissection we isolated RNA from unfixed RGCs and performed global transcriptome analysis of retinal neurons. In total, we detected 582 genes sequentially expressed in the preclinical phase and 1150 genes in the clinical manifest EAE (P < 0.05, fold-induction >1.5). Furthermore, using ingenuity pathway analysis (IPA), we identified amyloid precursor protein (APP) as a potential upstream regulator of changes in gene expression in the preclinical EAE but neither in clinical EAE, nor at any time point after optic nerve transection. Therefore, the gene pathway analysis lead to the hypothesis that altered cleavage of APP in neurons in the preclinical phase of EAE leads to the enhanced production of APP intracellular domain (AICD), which in turn acts as a transcriptional regulator and thereby initiates an apoptotic signaling cascade via up-regulation of the target gene p53.

The immune response in the CNS in Theiler's virus induced demyelinating disease switches from an early adaptive response to a chronic innate-like response

Theiler's murine encephalomyelitis virus-induced demyelinating disease (TMEV-IDD) is an important model of the progressive disability caused by irreversible CNS tissue injury, and provides an example of how a CNS pathogen can cause inflammation, demyelination, and neuronal damage. We were interested in which molecules, especially inflammatory mediators, might be upregulated in the CNS throughout TMEV-IDD. We quantitated by a real-time RT-PCR multi-gene system the expression of a pathway-focused panel of genes at 30 and 165 days post infection, characterizing both the early inflammatory and the late neurodegenerative stages of TMEV-IDD. Also, we measured 32 cytokines/chemokines by multiplex Luminex analysis in CSF specimens from early and late TMEV-IDD as well as sham-treated mice. Results indicate that, in the later stage of TMEV-IDD, activation of the innate immune response is most prominent: TLRs, type I IFN response genes, and innate immunity-associated cytokines were highly expressed in late TMEV-IDD compared to sham (p ≤ 0.0001) and early TMEV-IDD (p < 0.05). Conversely, several molecular mediators of adaptive immune response were highly expressed in early TMEV-IDD (all p ≤ 0.001). Protein detection in the CSF was broadly concordant with mRNA abundance of the corresponding gene measured by real-time RT-PCR in the spinal cord, since several cytokines/chemokines were increased in the CSF of TMEV-IDD mice. Results show a clear shift from adaptive to innate immunity from early to late TMEV-IDD, indicating that adaptive and innate immune pathways are likely involved in the development and progression of the disease to different extents. CSF provides an optimal source of biomarkers of CNS neuroinflammation.

Transcriptional analysis of glial cell differentiation in the postnatal murine spinal cord

Postnatal murine spinal cord represents a good model system to study mammalian central nervous system myelination in vivo as a basis for further studies in demyelinating diseases. Transcriptional changes were analyzed in SJL/J mice on postnatal day 0, 14, 49 and 231 (P0, P14, P49, P231) employing Affymetrix GeneChip Mouse Genome 430 2.0 Arrays. Additionally, marker gene signatures for astrocyte and oligodendrocyte lineage-stages were defined to study their gene expression in more detail. In addition, immunohistochemistry was used to quantify the abundance of commonly used glial cell markers. 6092 differentially regulated genes (DEGs) were identified. The up-regulated DEGs at P14, P49 and P231 compared to P0 exhibited significantly enriched associations to gene ontology terms such as myelination and lipid metabolic transport and down-regulated DEGs to neurogenesis and axonogenesis. Expression values of marker gene signatures for neural stem cells, oligodendrocyte precursor cells, and developing astrocytes were constantly decreasing, whereas myelinating oligodendrocyte and mature astrocyte markers showed a steady increase. Molecular findings were substantiated by immunohistochemical observations. The transcriptional changes observed are an important reference for future analysis of degenerative and inflammatory conditions in the spinal cord.

Dynamic Changes of Microglia/Macrophage M1 and M2 Polarization in Theiler's Murine Encephalomyelitis

Microglia and macrophages play a central role for demyelination in Theiler's murine encephalomyelitis (TME) virus infection, a commonly used infectious model for chronic-progressive multiple sclerosis. In order to determine the dynamic changes of microglia/macrophage polarization in TME, the spinal cord of Swiss Jim Lambert (SJL) mice was investigated by gene expression profiling and immunofluorescence. Virus persistence and demyelinating leukomyelitis were confirmed by immunohistochemistry and histology. Electron microscopy revealed continuous myelin loss together with abortive myelin repair during the late chronic infection phase indicative of incomplete remyelination. A total of 59 genes out of 151 M1- and M2-related genes were differentially expressed in TME virus-infected mice over the study period. The onset of virus-induced demyelination was associated with a dominating M1 polarization, while mounting M2 polarization of macrophages/microglia together with sustained prominent M1-related gene expression was present during the chronic-progressive phase. Molecular results were confirmed by immunofluorescence, showing an increased spinal cord accumulation of CD16/32(+) M1-, arginase-1(+) M2- and Ym1(+) M2-type cells associated with progressive demyelination. The present study provides a comprehensive database of M1-/M2-related gene expression involved in the initiation and progression of demyelination supporting the hypothesis that perpetuating interaction between virus and macrophages/microglia induces a vicious circle with persistent inflammation and impaired myelin repair in TME.