Drugst.One - a plug-and-play solution for online systems medicine and network-based drug repurposing

In recent decades, the development of new drugs has become increasingly expensive and inefficient, and the molecular mechanisms of most pharmaceuticals remain poorly understood. In response, computational systems and network medicine tools have emerged to identify potential drug repurposing candidates. However, these tools often require complex installation and lack intuitive visual network mining capabilities. To tackle these challenges, we introduce Drugst.One, a platform that assists specialized computational medicine tools in becoming user-friendly, web-based utilities for drug repurposing. With just three lines of code, Drugst.One turns any systems biology software into an interactive web tool for modeling and analyzing complex protein-drug-disease networks. Demonstrating its broad adaptability, Drugst.One has been successfully integrated with 21 computational systems medicine tools. Available at https://drugst.one, Drugst.One has significant potential for streamlining the drug discovery process, allowing researchers to focus on essential aspects of pharmaceutical treatment research.

Single-Cell Multi-Omics Map of Cell Type-Specific Mechanistic Drivers of Multiple Sclerosis Lesions

BACKGROUND AND OBJECTIVES

In progressive multiple sclerosis (MS), compartmentalized inflammation plays a pivotal role in the complex pathology of tissue damage. The interplay between epigenetic regulation, transcriptional modifications, and location-specific alterations within white matter (WM) lesions at the single-cell level remains underexplored.

METHODS

We examined intracellular and intercellular pathways in the MS brain WM using a novel dataset obtained by integrated single-cell multi-omics techniques from 3 active lesions, 3 chronic active lesions, 3 remyelinating lesions, and 3 control WM of 6 patients with progressive MS and 3 non-neurologic controls. Single-nucleus RNA-seq and ATAC-seq were combined and additionally enriched with newly conducted spatial transcriptomics from 1 chronic active lesion. Functional gene modules were then validated in our previously published bulk tissue transcriptome data obtained from 73 WM lesions of patients with progressive MS and 25 WM of non-neurologic disease controls.

RESULTS

Our analysis uncovered an MS-specific oligodendrocyte genetic signature influenced by the KLF/SP gene family. This modulation has potential associations with the autocrine iron uptake signaling observed in transcripts of transferrin and its receptor LRP2. In addition, an inflammatory profile emerged within these oligodendrocytes. We observed unique cellular endophenotypes both at the periphery and within the chronic active lesion. These include a distinct metabolic astrocyte phenotype, the importance of FGF signaling among astrocytes and neurons, and a notable enrichment of mitochondrial genes at the lesion edge populated predominantly by astrocytes. Our study also identified B-cell coexpression networks indicating different functional B-cell subsets with differential location and specific tendencies toward certain lesion types.

DISCUSSION

The use of single-cell multi-omics has offered a detailed perspective into the cellular dynamics and interactions in MS. These nuanced findings might pave the way for deeper insights into lesion pathogenesis in progressive MS.

Expression of Bruton´s tyrosine kinase in different type of brain lesions of multiple sclerosis patients and during experimental demyelination

Background

Inhibition of Bruton's tyrosine kinase (BTK) is an emerging multiple sclerosis (MS) therapy. BTK inhibitors (BTKi) cross the blood-brain barrier and modulate B cells and microglia, major cellular players in active and chronic active lesions.

Objective

To assess potential lesional and cellular targets of BTKi, we examined BTK expression in different type of MS white matter (WM) lesions, in unmanipulated CNS resident cells, and in a degenerative MS model associated with microglia activation in vivo.

Methods

We examined BTK expression by next-generation RNA-sequencing in postmortem 25 control WM, 19 NAWM, 6 remyelinating, 18 active, 13 inactive and 17 chronic active lesions. Presence of B cells and microglia were examined by immunohistochemistry. CNS resident cells were isolated from the mouse brain by magnetic sorting. BTK expression was examined by quantitative PCR in isolated cells and dissected corpus callosum from mice treated with cuprizone (CPZ).

Results

BTK expression was significantly increased in active and chronic active lesions with upregulated complement receptors and Fcγ receptors. Active lesions contained high number of perivascular B cells, microglia, and macrophages. Chronic active lesions were characterized by microglia/macrophages in the rim. Microglia expressed BTK at high level (120-fold) in contrast to other CNS cell types (2-4-fold). BTK expression was increasing during CPZ treatment reaching significance after stopping CPZ.

Conclusion

Considering BTK expression in MS lesions and resident cells, BTKi may exert effect on B cells, microglia/macrophages in active lesions, and limit microglia activation in chronic active lesions, where tissue damage propagates.

Inference of differential key regulatory networks and mechanistic drug repurposing candidates from scRNA-seq data with SCANet

MOTIVATION

The reconstruction of small key regulatory networks that explain the differences in the development of cell (sub)types from single-cell RNA sequencing is a yet unresolved computational problem.

RESULTS

To this end, we have developed SCANet, an all-in-one package for single-cell profiling that covers the whole differential mechanotyping workflow, from inference of trait/cell-type-specific gene co-expression modules, driver gene detection, and transcriptional gene regulatory network reconstruction to mechanistic drug repurposing candidate prediction. To illustrate the power of SCANet, we examined data from two studies. First, we identify the drivers of the mechanotype of a cytokine storm associated with increased mortality in patients with acute respiratory illness. Secondly, we find 20 drugs for eight potential pharmacological targets in cellular driver mechanisms in the intestinal stem cells of obese mice.

AVAILABILITY AND IMPLEMENTATION

SCANet is a free, open-source, and user-friendly Python package that can be seamlessly integrated into single-cell-based systems medicine research and mechanistic drug discovery.

Whole blood miRNAs in relapsing MS patients treated with dimethyl fumarate in the phase 4 TREMEND trial

We investigated the impact of dimethyl fumarate (DMF), an oral therapy for relapsing multiple sclerosis (MS), on blood microRNA (miRNA) signatures and neurofilament light (NFL) levels. DMF normalized miR-660-5p and modulated various miRNAs associated with the NF-kB pathway. These alterations reached a peak 4-7 months after treatment. Notably, particular miRNAs correlated with high or low NFL levels, implying their potential role as markers of treatment efficacy. Our findings broaden the understanding of DMF's immunomodulatory effects and may aid in predicting treatment responses.

Drugst.One - A plug-and-play solution for online systems medicine and network-based drug repurposing

In recent decades, the development of new drugs has become increasingly expensive and inefficient, and the molecular mechanisms of most pharmaceuticals remain poorly understood. In response, computational systems and network medicine tools have emerged to identify potential drug repurposing candidates. However, these tools often require complex installation and lack intuitive visual network mining capabilities. To tackle these challenges, we introduce Drugst.One, a platform that assists specialized computational medicine tools in becoming user-friendly, web-based utilities for drug repurposing. With just three lines of code, Drugst.One turns any systems biology software into an interactive web tool for modeling and analyzing complex protein-drug-disease networks. Demonstrating its broad adaptability, Drugst.One has been successfully integrated with 21 computational systems medicine tools. Available at https://drugst.one, Drugst.One has significant potential for streamlining the drug discovery process, allowing researchers to focus on essential aspects of pharmaceutical treatment research.

De-novo reconstruction and identification of transcriptional gene regulatory network modules differentiating single-cell clusters

Single-cell RNA sequencing (scRNA-seq) technology provides an unprecedented opportunity to understand gene functions and interactions at single-cell resolution. While computational tools for scRNA-seq data analysis to decipher differential gene expression profiles and differential pathway expression exist, we still lack methods to learn differential regulatory disease mechanisms directly from the single-cell data. Here, we provide a new methodology, named DiNiro, to unravel such mechanisms de novo and report them as small, easily interpretable transcriptional regulatory network modules. We demonstrate that DiNiro is able to uncover novel, relevant, and deep mechanistic models that not just predict but explain differential cellular gene expression programs. DiNiro is available at https://exbio.wzw.tum.de/diniro/.

Hypothesis of a potential BrainBiota and its relation to CNS autoimmune inflammation

Infectious agents have been long considered to play a role in the pathogenesis of neurological diseases as part of the interaction between genetic susceptibility and the environment. The role of bacteria in CNS autoimmunity has also been highlighted by changes in the diversity of gut microbiota in patients with neurological diseases such as Parkinson's disease, Alzheimer disease and multiple sclerosis, emphasizing the role of the gut-brain axis. We discuss the hypothesis of a brain microbiota, the BrainBiota: bacteria living in symbiosis with brain cells. Existence of various bacteria in the human brain is suggested by morphological evidence, presence of bacterial proteins, metabolites, transcripts and mucosal-associated invariant T cells. Based on our data, we discuss the hypothesis that these bacteria are an integral part of brain development and immune tolerance as well as directly linked to the gut microbiome. We further suggest that changes of the BrainBiota during brain diseases may be the consequence or cause of the chronic inflammation similarly to the gut microbiota.

A Systematic Review of Tissue and Single Cell Transcriptome/Proteome Studies of the Brain in Multiple Sclerosis

Multiple sclerosis (MS) is an inflammatory demyelinating and degenerative disease of the central nervous system (CNS). Although inflammatory responses are efficiently treated, therapies for progression are scarce and suboptimal, and biomarkers to predict the disease course are insufficient. Cure or preventive measures for MS require knowledge of core pathological events at the site of the tissue damage. Novelties in systems biology have emerged and paved the way for a more fine-grained understanding of key pathological pathways within the CNS, but they have also raised questions still without answers. Here, we systemically review the power of tissue and single-cell/nucleus CNS omics and discuss major gaps of integration into the clinical practice. Systemic search identified 49 transcriptome and 11 proteome studies of the CNS from 1997 till October 2021. Pioneering molecular discoveries indicate that MS affects the whole brain and all resident cell types. Despite inconsistency of results, studies imply increase in transcripts/proteins of semaphorins, heat shock proteins, myelin proteins, apolipoproteins and HLAs. Different lesions are characterized by distinct astrocytic and microglial polarization, altered oligodendrogenesis, and changes in specific neuronal subtypes. In all white matter lesion types, CXCL12, SCD, CD163 are highly expressed, and STAT6- and TGFβ-signaling are increased. In the grey matter lesions, TNF-signaling seems to drive cell death, and especially CUX2-expressing neurons may be susceptible to neurodegeneration. The vast heterogeneity at both cellular and lesional levels may underlie the clinical heterogeneity of MS, and it may be more complex than the current disease phenotyping in the clinical practice. Systems biology has not solved the mystery of MS, but it has discovered multiple molecules and networks potentially contributing to the pathogenesis. However, these results are mostly descriptive; focused functional studies of the molecular changes may open up for a better interpretation. Guidelines for acceptable quality or awareness of results from low quality data, and standardized computational and biological pipelines may help to overcome limited tissue availability and the “snap shot” problem of omics. These may help in identifying core pathological events and point in directions for focus in clinical prevention.

Functional enrichment of alternative splicing events with NEASE reveals insights into tissue identity and diseases

Alternative splicing (AS) is an important aspect of gene regulation. Nevertheless, its role in molecular processes and pathobiology is far from understood. A roadblock is that tools for the functional analysis of AS-set events are lacking. To mitigate this, we developed NEASE, a tool integrating pathways with structural annotations of protein-protein interactions to functionally characterize AS events. We show in four application cases how NEASE can identify pathways contributing to tissue identity and cell type development, and how it highlights splicing-related biomarkers. With a unique view on AS, NEASE generates unique and meaningful biological insights complementary to classical pathways analysis.

CSF proteome in multiple sclerosis subtypes related to brain lesion transcriptomes

To identify markers in the CSF of multiple sclerosis (MS) subtypes, we used a two-step proteomic approach: (i) Discovery proteomics compared 169 pooled CSF from MS subtypes and inflammatory/degenerative CNS diseases (NMO spectrum and Alzheimer disease) and healthy controls. (ii) Next, 299 proteins selected by comprehensive statistics were quantified in 170 individual CSF samples. (iii) Genes of the identified proteins were also screened among transcripts in 73 MS brain lesions compared to 25 control brains. F-test based feature selection resulted in 8 proteins differentiating the MS subtypes, and secondary progressive (SP)MS was the most different also from controls. Genes of 7 out these 8 proteins were present in MS brain lesions: GOLM was significantly differentially expressed in active, chronic active, inactive and remyelinating lesions, FRZB in active and chronic active lesions, and SELENBP1 in inactive lesions. Volcano maps of normalized proteins in the different disease groups also indicated the highest amount of altered proteins in SPMS. Apolipoprotein C-I, apolipoprotein A-II, augurin, receptor-type tyrosine-protein phosphatase gamma, and trypsin-1 were upregulated in the CSF of MS subtypes compared to controls. This CSF profile and associated brain lesion spectrum highlight non-inflammatory mechanisms in differentiating CNS diseases and MS subtypes and the uniqueness of SPMS.

Unbiased examination of genome-wide human endogenous retrovirus transcripts in MS brain lesions

BACKGROUND

Human endogenous retrovirus (HERV) expression in multiple sclerosis (MS) brain lesions may contribute to chronic inflammation, but expression of genome-wide HERVs in different MS lesions is unknown.

OBJECTIVE

We examined the HERV expression landscape in different MS lesions compared to control brains.

METHODS

Transcripts from 71 MS brain samples and 25 control WM were obtained by next-generation RNA sequencing and mapped against HERV transcripts across the human genome. Differential expression of mapped HERV-W and HERV-H reads between MS lesion types and controls was analysed.

RESULTS

Out of 6.38 billion high-quality paired end reads, 174 million reads (2.73%) mapped to HERV transcripts. There was no difference in HERVs expression level between MS and control brains, but HERV-W transcripts were significantly reduced in chronic active lesions. Of the four HERV-W transcripts exclusively present in MS, ERV3633503 located on chromosome 7q21.13 close to the MS genetic risk locus had the highest number of reads. In the HERV-H family, 75% of transcripts located to nearby 7q21-22 were overrepresented in MS, and ERV3643914 was expressed more than 16 times in MS compared to control brains.

CONCLUSION

Novel HERV-W and HERV-H transcripts located at chromosome 7 regions were uniquely expressed in MS lesions, indicating their potential role in brain lesion evolution.

Multiple Sclerosis Atlas: A Molecular Map of Brain Lesion Stages in Progressive Multiple Sclerosis

Introduction: Multiple sclerosis (MS) is a chronic disorder of the central nervous system with an untreatable late progressive phase. Molecular maps of different stages of brain lesion evolution in patients with progressive multiple sclerosis (PMS) are missing but critical for understanding disease development and to identify novel targets to halt progression. Materials and Methods: The MS Atlas database comprises comprehensive high-quality transcriptomic profiles of 98 white matter (WM) brain samples of different lesion types (normal-appearing WM [NAWM], active, chronic active, inactive, remyelinating) from ten progressive MS patients and 25 WM areas from five non-neurological diseased cases. Results: We introduce the first MS brain lesion atlas (msatlas.dk), developed to address the current challenges of understanding mechanisms driving the fate on a lesion basis. The MS Atlas gives means for testing research hypotheses, validating biomarkers and drug targets. It comes with a user-friendly web interface, and it fosters bioinformatic methods for de novo network enrichment to extract mechanistic markers for specific lesion types and pathway-based lesion type comparison. We describe examples of how the MS Atlas can be used to extract systems medicine signatures and demonstrate the interface of MS Atlas. Conclusion: This compendium of mechanistic PMS WM lesion profiles is an invaluable resource to fuel future MS research and a new basis for treatment development.

Cladribine modifies functional properties of microglia

Cladribine (CdA), an oral prodrug approved for the treatment of relapsing multiple sclerosis, selectively depletes lymphocytes. CdA passes the blood–brain barrier, suggesting a potential effect on central nervous system (CNS) resident cells. We examined if CdA modifies the phenotype and function of naive and activated primary mouse microglia, when applied in the concentrations 0·1–1 μM that putatively overlap human cerebrospinal fluid (CSF) concentrations. Primary microglia cultures without stimulation or in the presence of proinflammatory lipopolysaccharide (LPS) or anti‐inflammatory interleukin (IL)‐4 were treated with different concentrations of CdA for 24 h. Viability was assessed by MTT [3‐(4,5‐dimethylthiazol‐2‐yl)‐2,5‐diphenyltetrazolium bromide] assay. Phagocytotic ability and morphology were examined by flow cytometry and random migration using IncuCyte Zoom and TrackMate. Change in gene expression was examined by quantitative polymerase chain reaction (qPCR) and protein secretion by Meso Scale Discovery. We found that LPS and IL‐4 up‐regulated deoxycytidine kinase (DCK) expression. Only activated microglia were affected by CdA, and this was unrelated to viability. CdA 0·1–1 μM significantly reduced granularity, phagocytotic ability and random migration of activated microglia. CdA 10 μM increased the IL‐4‐induced gene expression of arginase 1 (Arg1) and LPS‐induced expression of IL‐1β, tumor necrosis factor (TNF), inducible nitric oxide synthase (iNOS) and Arg1, but protein secretion remained unaffected. CdA 10 μM potentiated the increased expression of anti‐inflammatory TNF receptor 2 (TNF‐R2) but not TNF‐R1 induced by LPS. This suggests that microglia acquire a less activated phenotype when treated with 0·1–1 μM CdA that putatively overlaps human CSF concentrations. This may be related to the up‐regulated gene expression of DCK upon activation, and suggests a potential alternative mechanism of CdA with direct effect on CNS resident cells.

Absence of miRNA-146a Differentially Alters Microglia Function and Proteome

Background: MiR-146a is an important regulator of innate inflammatory responses and is also implicated in cell death and survival.

Methods: By sorting CNS resident cells, microglia were the main cellular source of miR-146a. Therefore, we investigated microglia function and phenotype in miR-146a knock-out (KO) mice, analyzed the proteome of KO and wild-type (WT) microglia by LC-MS/MS, and examined miR-146a expression in different brain lesions of patients with multiple sclerosis (MS).

Results: When stimulated with LPS or myelin in vitro, microglia from KO mice expressed higher levels of IL-1β, TNF, IL-6, IL-10, CCL3, and CCL2 compared to WT. Stimulation increased migration and phagocytosis of WT but not KO microglia. CD11c⁺ microglia were induced by cuprizone (CPZ) in the WT mice but less in the KO. The proteome of ex vivo microglia was not different in miR-146a KO compared to WT mice, but CPZ treatment induced differential and reduced protein responses in the KO: GOT1, COX5b, CRYL1, and cystatin-C were specifically changed in KO microglia. We explored discriminative features of microglia proteomes: sparse Partial Least Squares-Discriminant Analysis showed the best discrimination when control and CPZ-treated conditions were compared. Cluster of ten proteins separated WT and miR-146a KO microglia after CPZ: among them were sensomes allowing to perceive the environment, Atp1a3 that belongs to the signature of CD11c⁺ microglia, and proteins related to inflammatory responses (S100A9, Ppm1g). Finally, we examined the expression of miR-146a and its validated target genes in different brain lesions of MS patients. MiR-146 was upregulated in all lesion types, and the highest expression was in active lesions. Nineteen of 88 validated target genes were significantly changed in active lesions, while none were changed in NAWM.

Conclusion: Our data indicated that microglia is the major source of miR-146a in the CNS. The absence of miR-146a differentially affected microglia function and proteome, and miR-146a may play an important role in gene regulation of active MS lesions.

Proteomic changes during experimental de- and remyelination in the corpus callosum

BACKGROUND

In the cuprizone model of multiple sclerosis, de- and remyelination can be studied without major interference from the adaptive immune responses. Since previous proteomic studies did not focus on the corpus callosum, where cuprizone causes the most pronounced demyelination, we performed a bottom up proteomic analysis on this brain region.

METHODS

Eight week-old mice treated with 0.2% cuprizone, for 4 weeks and controls (C) were sacrificed after termination of the treatment (4wD), and 2 (2dR) or 14 (2wR) days later. Homogenates of dissected corpus callosum were analysed by quantitative proteomics. For data processing, clustering, gene ontology analysis, and regulatory network prediction, we used Perseus, PANTHER and Ingenuity Pathway Analysis softwares, respectively.

RESULTS

We identified 4886 unmodified, single- or multi phosphorylated and/or gycosylated (PTM) proteins. Out of them, 191 proteins were differentially regulated in at least one experimental group. We found 57 proteins specific for demyelination, 27 for early- and 57 for late remyelinationwhile 36 proteins were affected in two, and 23 proteins in all three groups. Phosphorylation represented 92% of the post translational modifications among differentially regulated modified (PTM) proteins with decreased level, while it was only 30% of the PTM proteins with increased level. Gene ontology analysis could not classify the demyelination specific proteins into any biological process category, while allocated the remyelination specific ones to nervous system development and myelination as the most specific subcategory. We also identified a protein network in experimental remyelination, and the gene orthologues of the network were differentially expressed in remyelinating multiple sclerosis brain lesions consistent with an early remyelination pattern.

CONCLUSION

Proteomic analysis seems more informative for remyelination than demyelination in the cuprizone model.

Molecular signature of different lesion types in the brain white matter of patients with progressive multiple sclerosis

To identify pathogenetic markers and potential drivers of different lesion types in the white matter (WM) of patients with progressive multiple sclerosis (PMS), we sequenced RNA from 73 different WM areas. Compared to 25 WM controls, 6713 out of 18,609 genes were significantly differentially expressed in MS tissues (FDR < 0.05). A computational systems medicine analysis was performed to describe the MS lesion endophenotypes. The cellular source of specific molecules was examined by RNAscope, immunohistochemistry, and immunofluorescence. To examine common lesion specific mechanisms, we performed de novo network enrichment based on shared differentially expressed genes (DEGs), and found TGFβ-R2 as a central hub. RNAscope revealed astrocytes as the cellular source of TGFβ-R2 in remyelinating lesions. Since lesion-specific unique DEGs were more common than shared signatures, we examined lesion-specific pathways and de novo networks enriched with unique DEGs. Such network analysis indicated classic inflammatory responses in active lesions; catabolic and heat shock protein responses in inactive lesions; neuronal/axonal specific processes in chronic active lesions. In remyelinating lesions, de novo analyses identified axonal transport responses and adaptive immune markers, which was also supported by the most heterogeneous immunoglobulin gene expression. The signature of the normal-appearing white matter (NAWM) was more similar to control WM than to lesions: only 465 DEGs differentiated NAWM from controls, and 16 were unique. The upregulated marker CD26/DPP4 was expressed by microglia in the NAWM but by mononuclear cells in active lesions, which may indicate a special subset of microglia before the lesion develops, but also emphasizes that omics related to MS lesions should be interpreted in the context of different lesions types. While chronic active lesions were the most distinct from control WM based on the highest number of unique DEGs (n = 2213), remyelinating lesions had the highest gene expression levels, and the most different molecular map from chronic active lesions. This may suggest that these two lesion types represent two ends of the spectrum of lesion evolution in PMS. The profound changes in chronic active lesions, the predominance of synaptic/neural/axonal signatures coupled with minor inflammation may indicate end-stage irreversible molecular events responsible for this less treatable phase.

Unique RNA signature of different lesion types in the brain white matter in progressive multiple sclerosis

The heterogeneity of multiple sclerosis is reflected by dynamic changes of different lesion types in the brain white matter (WM). To identify potential drivers of this process, we RNA-sequenced 73 WM areas from patients with progressive MS (PMS) and 25 control WM. Lesion endophenotypes were described by a computational systems medicine analysis combined with RNAscope, immunohistochemistry, and immunofluorescence. The signature of the normal-appearing WM (NAWM) was more similar to control WM than to lesions: one of the six upregulated genes in NAWM was CD26/DPP4 expressed by microglia. Chronic active lesions that become prominent in PMS had a signature that were different from all other lesion types, and were differentiated from them by two clusters of 62 differentially expressed genes (DEGs). An upcoming MS biomarker, CHI3L1 was among the top ten upregulated genes in chronic active lesions expressed by astrocytes in the rim. TGFβ-R2 was the central hub in a remyelination-related protein interaction network, and was expressed there by astrocytes. We used de novo networks enriched by unique DEGs to determine lesion-specific pathway regulation, i.e. cellular trafficking and activation in active lesions; healing and immune responses in remyelinating lesions characterized by the most heterogeneous immunoglobulin gene expression; coagulation and ion balance in inactive lesions; and metabolic changes in chronic active lesions. Because we found inverse differential regulation of particular genes among different lesion types, our data emphasize that omics related to MS lesions should be interpreted in the context of lesion pathology. Our data indicate that the impact of molecular pathways is substantially changing as different lesions develop. This was also reflected by the high number of unique DEGs that were more common than shared signatures. A special microglia subset characterized by CD26 may play a role in early lesion development, while astrocyte-derived TGFβ-R2 and TGFβ pathways may be drivers of repair in contrast to chronic tissue damage. The highly specific mechanistic signature of chronic active lesions indicates that as these lesions develop in PMS, the molecular changes are substantially skewed: the unique mitochondrial/metabolic changes and specific downregulation of molecules involved in tissue repair may reflect a stage of exhaustion.

Orthologous proteins of experimental de- and remyelination are differentially regulated in the CSF proteome of multiple sclerosis subtypes

OBJECTIVE

Here, we applied a multi-omics approach (i) to examine molecular pathways related to de- and remyelination in multiple sclerosis (MS) lesions; and (ii) to translate these findings to the CSF proteome in order to identify molecules that are differentially expressed among MS subtypes.

METHODS

To relate differentially expressed genes in MS lesions to de- and remyelination, we compared transcriptome of MS lesions to transcriptome of cuprizone (CPZ)-induced de- and remyelination. Protein products of the overlapping orthologous genes were measured within the CSF by quantitative proteomics, parallel reaction monitoring (PRM). Differentially regulated proteins were correlated with molecular markers of inflammation by using MesoScale multiplex immunoassay. Expression kinetics of differentially regulated orthologous genes and proteins were examined in the CPZ model.

RESULTS

In the demyelinated and remyelinated corpus callosum, we detected 1239 differentially expressed genes; 91 orthologues were also differentially expressed in MS lesions. Pathway analysis of these orthologues suggested that the TYROBP (DAP12)-TREM2 pathway, TNF-receptor 1, CYBA and the proteasome subunit PSMB9 were related to de- and remyelination. We designed 129 peptides representing 51 orthologous proteins, measured them by PRM in 97 individual CSF, and compared their levels between relapsing (n = 40) and progressive MS (n = 57). Four proteins were differentially regulated among relapsing and progressive MS: tyrosine protein kinase receptor UFO (UFO), TIMP-1, apolipoprotein C-II (APOC2), and beta-2-microglobulin (B2M). The orthologous genes/proteins in the mouse brain peaked during acute remyelination. UFO, TIMP-1 and B2M levels correlated inversely with inflammation in the CSF (IL-6, MCP-1/CCL2, TARC/CCL17). APOC2 showed positive correlation with IL-2, IL-16 and eotaxin-3/CCL26.

CONCLUSIONS

Pathology-based multi-omics identified four CSF markers that were differentially expressed in MS subtypes. Upregulated TIMP-1, UFO and B2M orthologues in relapsing MS were associated with reduced inflammation and reflected reparatory processes, in contrast to the upregulated orthologue APOC2 in progressive MS that reflected changes in lipid metabolism associated with increased inflammation.

Experimental Demyelination and Axonal Loss Are Reduced in MicroRNA-146a Deficient Mice

Background

The cuprizone (CPZ) model of multiple sclerosis (MS) was used to identify microRNAs (miRNAs) related to in vivo de- and remyelination. We further investigated the role of miR-146a in miR-146a-deficient (KO) mice: this miRNA is differentially expressed in MS lesions and promotes differentiation of oligodendrocyte precursor cells (OPCs) during remyelination, but its role has not been examined during demyelination.

Methods

MicroRNAs were examined by Agilent Mouse miRNA Microarray in the corpus callosum during CPZ-induced demyelination and remyelination. Demyelination, axonal loss, changes in number of oligodendrocytes, OPCs, and macrophages/microglia was compared by histology/immunohistochemistry between KO and WT mice. Differential expression of target genes and proteins of miR-146a was analyzed in the transcriptome (4 × 44K Agilent Whole Mouse Genome Microarray) and proteome (liquid chromatography tandem mass spectrometry) of CPZ-induced de- and remyelination in WT mice. Levels of proinflammatory molecules in the corpus callosum were compared in WT versus KO mice by Meso Scale Discovery multiplex protein analysis.

Results

miR-146a was increasingly upregulated during CPZ-induced de- and remyelination. The absence of miR-146a in KO mice protected against demyelination, axonal loss, body weight loss, and atrophy of thymus and spleen. The number of CNP⁺ oligodendrocytes was increased during demyelination in the miR-146a KO mice, while there was a trend of increased number of NG2⁺ OPCs in the WT mice. miR-146a target genes, SNAP25 and SMAD4, were downregulated in the proteome of demyelinating corpus callosum in WT mice. Higher levels of SNAP25 were measured by ELISA in the corpus callosum of miR-146a KO mice, but there was no difference between KO and WT mice during demyelination. Multiplex protein analysis of the corpus callosum lysate revealed upregulated TNF-RI, TNF-RII, and CCL2 in the WT mice in contrast to KO mice. The number of Mac3⁺ and Iba1⁺ macrophages/microglia was reduced in the demyelinating corpus callosum of the KO mice.

Conclusion

During demyelination, absence of miR-146a reduced inflammatory responses, demyelination, axonal loss, the number of infiltrating macrophages, and increased the number of myelinating oligodendrocytes. The number of OPCs was slightly higher in the WT mice during remyelination, indicating a complex role of miR-146a during in vivo de- and remyelination.

Immunolocalization of aquaporin-8 in rat kidney, gastrointestinal tract, testis, and airways

The purpose of this study was to determine the cellular and subcellular localization of aquaporin-8 (AQP8) in rat kidney and other organs by RT-PCR analyses and by immunoblotting and immunohistochemistry using peptide-derived rabbit antibodies to rat AQP8. RT-PCR and Southern blotting revealed the presence of AQP8 mRNA in all kidney zones. LLC-PK(1) cells transfected with a rat AQP8 construct exhibited strong labeling with the affinity-purified antibodies, whereas controls using cells transfected with the vector, but without the insert, were negative. The labeling was almost exclusively associated with intracellular vesicles. Immunoblotting of kidney membrane fractions revealed a predominant single band of 26-28 kDa. AQP8 immunoreactivity was mainly present in the cortex and outer stripe of the outer medulla. Sequential ultracentrifugation of rat kidney membrane revealed that AQP8 resides predominantly in intracellular vesicular fractions. Immunocytochemistry revealed modest labeling of proximal tubules and weak labeling of collecting ducts in cortex and medulla of rat kidney. The labeling was confined to cytoplasmic areas with no labeling of the brush border. Immunoblotting and RT-PCR/Southern blotting also revealed the presence of AQP8 protein and mRNA in rat liver, testis, epididymis, duodenum, jejunum, colon, and bronchi/trachea. Consistent with this, immunohistochemistry revealed AQP8 labeling in the hepatocytes and spermatogenic cells in testis and in the basal cells in ductus epididymis, trachea, and bronchial epithelia. Moreover, AQP8 labeling was observed in the myoepithelial cells in salivary, bronchial, and tracheal glands with no labeling of acini or ductal epithelial cells. AQP8 is also present in the surface epithelial cells in duodenum, jejunum, and colon. In conclusion, AQP8 is expressed at low levels in rat kidney proximal tubules and collecting ducts, and it is present in distinct cell types in liver, testis, epididymis, duodenum, jejunum, colon, trachea, and principal bronchi as well as in multiple glands, including salivary glands.

Detection of Duffy antigen in the plasma membranes and caveolae of vascular endothelial and epithelial cells of nonerythroid organs

The nonerythroid expression of the Duffy blood group protein (gp-Fy) was confined to certain cell types. Immunocytochemistry studies of the kidney showed gp-Fy in the endothelium of glomeruli, peritubular capillaries, vasa recta, and the principal cells (epithelial) of collecting ducts. Gp-Fy was also produced in the endothelial cells of large venules and epithelial cells (type-I) of pulmonary alveoli. In the thyroid, only the endothelial cells of capillaries produced gp-Fy. In the spleen, the endothelial cells of capillaries, high endothelial venule, and sinusoids produced abundant gp-Fy. Ultrastructural studies showed that apical and basolateral plasma membrane domains, including caveolae, had gp-Fy. Immunoblot analysis showed substantially less gp-Fy in nonerythroid cells than in erythrocytes. Moreover, the analyzed nonerythroid organs of Duffy-negative individuals did not produce more gp-Fy to compensate for the lack of this protein in their erythrocytes. The nucleotide sequence and the size of kidney mRNA from a Duffy-positive individual were the same as that of bone marrow. It is assumed, therefore, that nonerythroid Duffy protein is the product of the same gene as that of bone marrow. This notion is reinforced by the fact that nonerythroid and erythroid gp-Fy have the same antigenic domains.

Effects of microtubule disruption on endocytosis, membrane recycling and polarized distribution of Aquaporin-1 and gp330 in proximal tubule cells

Microtubules are critically involved in membrane traffic and maintenance of epithelial cell polarity. In this study we examined the effect of microtubule disruption by colchicine on 1) the subcellular organization of the apical endocytic apparatus, 2) apical endocytosis, and 3) subcellular distribution of gp330 and Aquaporin-1 water channels in renal proximal tubule cells. Rats were treated in vivo with colchicine for 5 h before fixation of the kidneys, and sections of proximal tubules were studied using electron microscopy, morphometry and immunocytochemistry. In proximal tubule cells from colchicine-treated animals, virtually no endocytic invaginations are present, indicating that colchicine blocks the formation of endocytic invaginations. Furthermore, no large endocytic vacuoles are present, also consistent with a block of endocytosis. This was confirmed by functional studies revealing a colchicine-induced arrest in apical endocytosis of peroxidase. There is a marked reduction in dense apical tubules (the exocytic vehicle for membrane recycling) but an extensive accumulation of small vesicles, suggesting a disruption in membrane recycling. This disruption may be primary or secondary to the block of endocytosis. Immunofluorescence and immunoelectron microscopy reveal extensive colchicine-induced changes in the subcellular distribution of gp330 and Aquaporin-1. The localization of gp330 is normally restricted to the apical endocytic apparatus. However, after colchicine treatment gp330 is localized to numerous vesicles distributed throughout the entire cytoplasm, as previously shown (Gutmann et al., Am. J. Physiol. 257, C397-C407 (1989). Also Aquaporin-1 water channels, which are normally almost exclusively present in the plasma membranes, are redistributed to small vesicles in addition to the plasma membrane localization after colchicine treatment. In summary, the present study suggests that cytoplasmic microtubules are critically involved in 1) formation and stabilization of endocytic invaginations, 2) formation of large endocytic vacuoles, 3) apical endocytosis, 4) maintaining the polarized distribution of vesicles in the apical part of the cell, and 5) maintaining the polarized organization of gp330 in the apical endocytic apparatus, and maintaining Aquaporin-1 water channels in the plasma membranes.