A gene regulatory network approach harmonizes genetic and epigenetic signals and reveals repurposable drug candidates for multiple sclerosis.

Multiple sclerosis (MS) is a complex dysimmune disorder of the central nervous system. Genome-wide association studies (GWAS) have identified 233 genetic variations associated with MS at the genome-wide significant level. Epigenetic studies have pinpointed differentially methylated CpG sites in MS patients. However, the interplay between genetic risk factors and epigenetic regulation remains elusive. Here, we employed a network model to integrate GWAS summary statistics of 14 802 MS cases and 26 703 controls with DNA methylation profiles from 140 MS cases and 139 controls and the human interactome. We identified differentially methylated genes by aggregating additive effects of differentially methylated CpG sites within promoter regions. We reconstructed a gene regulatory network (GRN) using literature-curated transcription factor knowledge. Colocalization of the MS GWAS and methylation quantitative trait loci (mQTL) was performed to assess the GRN. The resultant MS-associated GRN highlighted several single nucleotide polymorphisms with GWAS-mQTL colocalization: rs6032663, rs6065926 and rs2024568 of CD40 locus, rs9913597 of STAT3 locus, and rs887864 and rs741175 of CIITA locus. Moreover, synergistic mQTL and expression QTL signals were identified in CD40, suggesting gene expression alteration was likely induced by epigenetic changes. Web-based Cell-type Specific Enrichment Analysis of Genes (WebCSEA) indicated that the GRN was enriched in T follicular helper cells (P-value = 0.0016). Drug target enrichment analysis of annotations from the Therapeutic Target Database revealed the GRN was also enriched with drug target genes (P-value = 3.89 × 10-4), revealing repurposable candidates for MS treatment. These candidates included vorinostat (HDAC1 inhibitor) and sivelestat (ELANE inhibitor), which warrant further investigation.

An integrative study of genetic variants with brain tissue expression identifies viral etiology and potential drug targets of multiple sclerosis.

Multiple sclerosis (MS) is a neuroinflammatory disorder leading to chronic disability. Brain lesions in MS commonly arise in normal-appearing white matter (NAWM). Genome-wide association studies (GWAS) have identified genetic variants associated with MS. Transcriptome alterations have been observed in case-control studies of NAWM. We developed a Cross-Dataset Evaluation (CDE) function for our network-based tool, Edge-Weighted Dense Module Search of GWAS (EW_dmGWAS). We applied CDE to integrate publicly available MS GWAS summary statistics of 41,505 cases and controls with collectively 38 NAWM expression samples, using the human protein interactome as the reference network, to investigate biological underpinnings of MS etiology. We validated the resulting modules with colocalization of GWAS and expression quantitative trait loci (eQTL) signals, using GTEx Consortium expression data for MS-relevant tissues: 14 brain tissues and 4 immune-related tissues. Other network assessments included a drug target query and functional gene set enrichment analysis. CDE prioritized a MS NAWM network containing 55 unique genes. The gene list was enriched (p-value = 2.34 × 10-7) with GWAS-eQTL colocalized genes: CDK4, IFITM3, MAPK1, MAPK3, METTL12B and PIK3R2. The resultant network also included drug signatures of FDA-approved medications. Gene set enrichment analysis revealed the top functional term "intracellular transport of virus", among other viral pathways. We prioritize critical genes from the resultant network: CDK4, IFITM3, MAPK1, MAPK3, METTL12B and PIK3R2. Enriched drug signatures suggest potential drug targets and drug repositioning strategies for MS. Finally, we propose mechanisms of potential MS viral onset, based on prioritized gene set and functional enrichment analysis.

Dense module searching for gene networks associated with multiple sclerosis.

BACKGROUND: Multiple sclerosis (MS) is a complex disease in which the immune system attacks the central nervous system. The molecular mechanisms contributing to the etiology of MS remain poorly understood. Genome-wide association studies (GWAS) of MS have identified a small number of genetic loci significant at the genome level, but they are mainly non-coding variants. Network-assisted analysis may help better interpret the functional roles of the variants with association signals and potential translational medicine application. The Dense Module Searching of GWAS tool (dmGWAS version 2.4) developed in our team is applied to 2 MS GWAS datasets (GeneMSA and IMSGC GWAS) using the human protein interactome as the reference network. A dual evaluation strategy is used to generate results with reproducibility. RESULTS: Approximately 7500 significant network modules were identified for each independent GWAS dataset, and 20 significant modules were identified from the dual evaluation. The top modules included GRB2, HDAC1, JAK2, MAPK1, and STAT3 as central genes. Top module genes were enriched with functional terms such as "regulation of glial cell differentiation" (adjusted p-value = 2.58 × 10- 3), "T-cell costimulation" (adjusted p-value = 2.11 × 10- 6) and "virus receptor activity" (adjusted p-value = 1.67 × 10- 3). Interestingly, top gene networks included several MS FDA approved drug target genes HDAC1, IL2RA, KEAP1, and RELA, CONCLUSIONS: Our dmGWAS network analyses highlighted several genes (GRB2, HDAC1, IL2RA, JAK2, KEAP1, MAPK1, RELA and STAT3) in top modules that are promising to interpret GWAS signals and link to MS drug targets. The genes enriched with glial cell differentiation are important for understanding neurodegenerative processes in MS and for remyelination therapy investigation. Importantly, our identified genetic signals enriched in T cell costimulation and viral receptor activity supported the viral infection onset hypothesis for MS.