Transcriptional network analysis of peripheral blood leukocyte subsets in multiple sclerosis identifies a pathogenic role for a cytotoxicity-associated gene network in myeloid cells.

Multiple sclerosis (MS) is an autoimmune disease of the central nervous system affecting predominantly adults. It is a complex disease associated with both environmental and genetic risk factors. Although over 230 risk single-nucleotide polymorphisms have been associated with MS, all are common human variants. The mechanisms by which they increase the risk of MS, however, remain elusive. We hypothesized that a complex genetic phenotype such as MS could be driven by coordinated expression of genes controlled by transcriptional regulatory networks. We, therefore, constructed a gene coexpression network from microarray expression analyses of five purified peripheral blood leukocyte subsets of 76 patients with relapsing remitting MS and 104 healthy controls. These analyses identified a major network (or module) of expressed genes associated with MS that play key roles in cell-mediated cytotoxicity which was downregulated in monocytes of patients with MS. Manipulation of the module gene expression was achieved in vitro through small interfering RNA gene knockdown of identified drivers. In a mouse model, network gene knockdown modulated the autoimmune inflammatory MS model disease-experimental autoimmune encephalomyelitis. This research implicates a cytotoxicity-associated gene network in myeloid cells in the pathogenesis of MS.

Not all roads lead to the immune system: the genetic basis of multiple sclerosis severity.

Multiple sclerosis is a leading cause of neurological disability in adults. Heterogeneity in multiple sclerosis clinical presentation has posed a major challenge for identifying genetic variants associated with disease outcomes. To overcome this challenge, we used prospectively ascertained clinical outcomes data from the largest international multiple sclerosis registry, MSBase. We assembled a cohort of deeply phenotyped individuals of European ancestry with relapse-onset multiple sclerosis. We used unbiased genome-wide association study and machine learning approaches to assess the genetic contribution to longitudinally defined multiple sclerosis severity phenotypes in 1813 individuals. Our primary analyses did not identify any genetic variants of moderate to large effect sizes that met genome-wide significance thresholds. The strongest signal was associated with rs7289446 (ß = -0.4882, P = 2.73 × 10-7), intronic to SEZ6L on chromosome 22. However, we demonstrate that clinical outcomes in relapse-onset multiple sclerosis are associated with multiple genetic loci of small effect sizes. Using a machine learning approach incorporating over 62 000 variants together with clinical and demographic variables available at multiple sclerosis disease onset, we could predict severity with an area under the receiver operator curve of 0.84 (95% CI 0.79-0.88). Our machine learning algorithm achieved positive predictive value for outcome assignation of 80% and negative predictive value of 88%. This outperformed our machine learning algorithm that contained clinical and demographic variables alone (area under the receiver operator curve 0.54, 95% CI 0.48-0.60). Secondary, sex-stratified analyses identified two genetic loci that met genome-wide significance thresholds. One in females (rs10967273; ßfemale = 0.8289, P = 3.52 × 10-8), the other in males (rs698805; ßmale = -1.5395, P = 4.35 × 10-8), providing some evidence for sex dimorphism in multiple sclerosis severity. Tissue enrichment and pathway analyses identified an overrepresentation of genes expressed in CNS compartments generally, and specifically in the cerebellum (P = 0.023). These involved mitochondrial function, synaptic plasticity, oligodendroglial biology, cellular senescence, calcium and G-protein receptor signalling pathways. We further identified six variants with strong evidence for regulating clinical outcomes, the strongest signal again intronic to SEZ6L (adjusted hazard ratio 0.72, P = 4.85 × 10-4). Here we report a milestone in our progress towards understanding the clinical heterogeneity of multiple sclerosis outcomes, implicating functionally distinct mechanisms to multiple sclerosis risk. Importantly, we demonstrate that machine learning using common single nucleotide variant clusters, together with clinical variables readily available at diagnosis can improve prognostic capabilities at diagnosis, and with further validation has the potential to translate to meaningful clinical practice change.

Common and Low Frequency Variants in MERTK Are Independently Associated with Multiple Sclerosis Susceptibility with Discordant Association Dependent upon HLA-DRB1*15:01 Status.

Multiple Sclerosis (MS) is a chronic inflammatory demyelinating disease of the central nervous system. The risk of developing MS is strongly influenced by genetic predisposition, and over 100 loci have been established as associated with susceptibility. However, the biologically relevant variants underlying disease risk have not been defined for the vast majority of these loci, limiting the power of these genetic studies to define new avenues of research for the development of MS therapeutics. It is therefore crucial that candidate MS susceptibility loci are carefully investigated to identify the biological mechanism linking genetic polymorphism at a given gene to the increased chance of developing MS. MERTK has been established as an MS susceptibility gene and is part of a family of receptor tyrosine kinases known to be involved in the pathogenesis of demyelinating disease. In this study we have refined the association of MERTK with MS risk to independent signals from both common and low frequency variants. One of the associated variants was also found to be linked with increased expression of MERTK in monocytes and higher expression of MERTK was associated with either increased or decreased risk of developing MS, dependent upon HLA-DRB1*15:01 status. This discordant association potentially extended beyond MS susceptibility to alterations in disease course in established MS. This study provides clear evidence that distinct polymorphisms within MERTK are associated with MS susceptibility, one of which has the potential to alter MERTK transcription, which in turn can alter both susceptibility and disease course in MS patients.

Transcriptomics identifies blunted immunomodulatory effects of vitamin D in people with multiple sclerosis.

Vitamin D deficiency is a risk factor for developing multiple sclerosis (MS). However, the immune effects of vitamin D in people with MS are not well understood. We analyzed transcriptomic datasets generated by RNA sequencing of immune cell subsets (CD4+, CD8+ T cells, B cells, monocytes) from 33 healthy controls and 33 untreated MS cases. We utilized a traditional bioinformatic pipeline and weighted gene co-expression network analysis (WGCNA) to determine genes and pathways correlated with endogenous vitamin D. In controls, CD4+ and CD8+ T cells had 1079 and 1188 genes, respectively, whose expressions were correlated with plasma 25-hydroxyvitamin D level (P < 0.05). Functional enrichment analysis identified association with TNF-alpha and MAPK signaling. In CD4+ T cells of controls, vitamin D level was associated with expression levels of several genes proximal to multiple sclerosis risk loci (P = 0.01). Genes differentially associated with endogenous vitamin D by case-control status were enriched in TNF-alpha signaling via NF-κB. WGCNA suggested a blunted response to vitamin D in cases relative to controls. Collectively, our findings provide further evidence for the immune effects of vitamin D, and demonstrate a differential immune response to vitamin D in cases relative to controls, highlighting a possible mechanism contributing to MS pathophysiology.

Epigenetic differences at the HTR2A locus in progressive multiple sclerosis patients.

The pathology of progressive multiple sclerosis (MS) is poorly understood. We have previously assessed DNA methylation in the CD4+ T cells of relapsing-remitting (RR) MS patients compared to healthy controls and identified differentially methylated regions (DMRs) in HLA-DRB1 and RNF39. This study aimed to investigate the DNA methylation profiles of the CD4+ T cells of progressive MS patients. DNA methylation was measured in two separate case/control cohorts using the Illumina 450K/EPIC arrays and data was analysed with the Chip Analysis Methylation Pipeline (ChAMP). Single nucleotide polymorphisms (SNPs) were assessed using the Illumina Human OmniExpress24 arrays and analysed using PLINK. Expression was assessed using the Illumina HT12 array and analysed in R using a combination of Limma and Illuminaio. We identified three DMRs at HTR2A, SLC17A9 and HDAC4 that were consistent across both cohorts. The DMR at HTR2A is located within the bounds of a haplotype block; however, the DMR remained significant after accounting for SNPs in the region. No expression changes were detected in any DMRs. HTR2A is differentially methylated in progressive MS independent of genotype. This differential methylation is not evident in RRMS, making it a potential biomarker of progressive disease.

Multiple sclerosis risk variants regulate gene expression in innate and adaptive immune cells.

At least 200 single-nucleotide polymorphisms (SNPs) are associated with multiple sclerosis (MS) risk. A key function that could mediate SNP-encoded MS risk is their regulatory effects on gene expression. We performed microarrays using RNA extracted from purified immune cell types from 73 untreated MS cases and 97 healthy controls and then performed Cis expression quantitative trait loci mapping studies using additive linear models. We describe MS risk expression quantitative trait loci associations for 129 distinct genes. By extending these models to include an interaction term between genotype and phenotype, we identify MS risk SNPs with opposing effects on gene expression in cases compared with controls, namely, rs2256814 MYT1 in CD4 cells (q = 0.05) and rs12087340 RF00136 in monocyte cells (q = 0.04). The rs703842 SNP was also associated with a differential effect size on the expression of the METTL21B gene in CD8 cells of MS cases relative to controls (q = 0.03). Our study provides a detailed map of MS risk loci that function by regulating gene expression in cell types relevant to MS.

Epstein-Barr virus-encoded LMP1 regulates epithelial cell motility and invasion via the ERK-MAPK pathway.

The Epstein-Barr virus (EBV) latent membrane protein 1 (LMP1) is an oncogenic protein which has previously been shown to engage the NF-kappaB, stress-activated MAP kinase, phosphatidylinositol 3-kinase (PI 3-kinase), and extracellular-regulated kinase (ERK)-MAPK pathways. In this study, we demonstrate that LMP1 activates ERK-MAPK in epithelial cells via the canonical Raf-MEK-ERK-MAPK pathway but in a Ras-independent manner. In agreement with the results of a previous study (B. A. Mainou, D. N. Everly, Jr., and N. Raab-Traub, J. Virol. 81:9680-9692, 2007), we show that the ability of LMP1 to activate ERK-MAPK mapped to its CTAR1 domain, the TRAF binding domain previously implicated in PI 3-kinase activation. A role for ERK-MAPK in LMP1-induced epithelial cell motility was identified, as LMP1-expressing cells displayed increased rates of haptotactic migration compared to those of LMP1-negative cells. These data implicate the ERK-MAPK pathway in LMP1-induced effects associated with transformation, suggesting that this pathway may contribute to the oncogenicity of LMP1 through its ability to promote cell motility and to enhance the invasive properties of epithelial cells.

The EBV-Encoded Oncoprotein, LMP1, Induces an Epithelial-to-Mesenchymal Transition (EMT) via Its CTAR1 Domain through Integrin-Mediated ERK-MAPK Signalling.

The Epstein⁻Barr virus (EBV)-encoded latent membrane protein 1 (LMP1) oncogene can induce profound effects on epithelial growth and differentiation including many of the features of the epithelial-to-mesenchymal transition (EMT). To better characterise these effects, we used the well-defined Madin Darby Canine Kidney (MDCK) epithelial cell model and found that LMP1 expression in these cells induces EMT as defined by characteristic morphological changes accompanied by loss of E-cadherin, desmosomal cadherin and tight junction protein expression. The induction of the EMT phenotype required a functional CTAR1 domain of LMP1 and studies using pharmacological inhibitors revealed contributions from signalling pathways commonly induced by integrin⁻ligand interactions: extracellular signal-regulated kinases/mitogen-activated protein kinases (ERK-MAPK), PI3-Kinase and tyrosine kinases, but not transforming growth factor beta (TGFβ). More detailed analysis implicated the CTAR1-mediated induction of Slug and Twist in LMP1-induced EMT. A key role for β1 integrin signalling in LMP1-mediated ERK-MAPK and focal adhesion kianse (FAK) phosphorylation was observed, and β1 integrin activation was found to enhance LMP1-induced cell viability and survival. These findings support an important role for LMP1 in disease pathogenesis through transcriptional reprogramming that enhances tumour cell survival and leads to a more invasive, metastatic phenotype.

The Epstein-Barr virus encoded LMP1 oncoprotein modulates cell adhesion via regulation of activin A/TGFß and ß1 integrin signalling.

Approximately 20% of global cancer incidence is causally linked to an infectious agent. Epstein-Barr virus (EBV) accounts for around 1% of all virus-associated cancers and is associated with nasopharyngeal carcinoma (NPC). Latent membrane protein 1 (LMP1), the major oncoprotein encoded by EBV, behaves as a constitutively active tumour necrosis factor (TNF) receptor activating a variety of signalling pathways, including the three classic MAPKs (ERK-MAPK, p38 MAPK and JNK/SAPK). The present study identifies novel signalling properties for this integral membrane protein via the induction and secretion of activin A and TGFß1, which are both required for LMP1's ability to induce the expression of the extracellular matrix protein, fibronectin. However, it is evident that LMP1 is unable to activate the classic Smad-dependent TGFß signalling pathway, but rather elicits its effects through the non-Smad arm of TGFß signalling. In addition, there is a requirement for JNK/SAPK signalling in LMP1-mediated fibronectin induction. LMP1 also induces the expression and activation of the major fibronectin receptor, α5ß1 integrin, an effect that is accompanied by increased focal adhesion formation and turnover. Taken together, these findings support the putative role for LMP1 in the pathogenesis of NPC by contributing to the metastatic potential of epithelial cells.

Association of plasma levels of Protein S with disease severity in multiple sclerosis.

BACKGROUND: The TAM family of receptor tyrosine kinases (TYRO3, AXL and MERTK) play important roles in modulating innate immune responses and central demyelination. The TAM receptor ligand Protein S (PROS) has also been shown to modulate innate immune cell responses. OBJECTIVES: We assessed whether plasma levels of PROS are changed in multiple sclerosis (MS) patients and whether changes are associated with disease severity. METHODS: Plasma levels of total and free PROS were measured using enzyme-linked immunosorbent assay in a discovery cohort (MS: 65, control: 14) and an independent replication cohort (MS: 29, control: 29). The Multiple Sclerosis Severity Score (MSSS) was used to evaluate associations between plasma PROS levels and disease severity. RESULTS: We found plasma levels of total, but not free PROS, were decreased in MS patients compared with controls. In female MS patients, we observed decreases in total and free PROS levels compared with controls. In addition, we also observed higher MSSS in patients with very low levels of plasma free PROS. CONCLUSIONS: These data suggest PROS may represent a potential marker of disease severity in MS.

The MS Risk Allele of CD40 Is Associated with Reduced Cell-Membrane Bound Expression in Antigen Presenting Cells: Implications for Gene Function.

Human genetic and animal studies have implicated the costimulatory molecule CD40 in the development of multiple sclerosis (MS). We investigated the cell specific gene and protein expression variation controlled by the CD40 genetic variant(s) associated with MS, i.e. the T-allele at rs1883832. Previously we had shown that the risk allele is expressed at a lower level in whole blood, especially in people with MS. Here, we have defined the immune cell subsets responsible for genotype and disease effects on CD40 expression at the mRNA and protein level. In cell subsets in which CD40 is most highly expressed, B lymphocytes and dendritic cells, the MS-associated risk variant is associated with reduced CD40 cell-surface protein expression. In monocytes and dendritic cells, the risk allele additionally reduces the ratio of expression of full-length versus truncated CD40 mRNA, the latter encoding secreted CD40. We additionally show that MS patients, regardless of genotype, express significantly lower levels of CD40 cell-surface protein compared to unaffected controls in B lymphocytes. Thus, both genotype-dependent and independent down-regulation of cell-surface CD40 is a feature of MS. Lower expression of a co-stimulator of T cell activation, CD40, is therefore associated with increased MS risk despite the same CD40 variant being associated with reduced risk of other inflammatory autoimmune diseases. Our results highlight the complexity and likely individuality of autoimmune pathogenesis, and could be consistent with antiviral and/or immunoregulatory functions of CD40 playing an important role in protection from MS.

Axonally derived matrilin-2 induces proinflammatory responses that exacerbate autoimmune neuroinflammation.

In patients with multiple sclerosis (MS) and mice with experimental autoimmune encephalomyelitis (EAE), inflammatory axonal injury is a major determinant of disability; however, the drivers of this injury are incompletely understood. Here, we used the EAE model and determined that the extracellular matrix protein matrilin-2 (MATN2) is an endogenous neuronal molecule that is regulated in association with inflammatory axonal injury. Compared with WT mice, mice harboring a deletion of Matn2 exhibited reduced disease severity and axon damage following induction of EAE. Evaluation of neuron-macrophage cocultures revealed that exogenous MATN2 specifically signals through TLR4 and directly induces expression of proinflammatory genes in macrophages, promoting axonal damage. Moreover, the MATN2-induced proinflammatory response was attenuated greatly in macrophages from Myd88 KO mice. Examination of brain sections from patients with MS revealed that MATN2 is expressed in lesions but not in normal-appearing white matter. Together, our results indicate that MATN2 is a deleterious endogenous neuroaxonal injury response signal that activates innate immune cells and could contribute to early axonal damage in CNS inflammatory diseases like MS.