Although the causes of Multiple Sclerosis (MS) still remain largely unknown, multiple lines of evidence suggest that Epstein–Barr virus (EBV) infection may contribute to the development of MS. Here, we aimed to identify the potential contribution of EBV-encoded and host cellular miRNAs to MS pathogenesis. We identified differentially expressed host miRNAs in EBV infected B cells (LCLs) and putative host/EBV miRNA interactions with MS risk loci. We estimated the genotype effect of MS risk loci on the identified putative miRNA:mRNA interactions in silico. We found that the protective allele of MS risk SNP rs4808760 reduces the expression of hsa-mir-3188-3p. In addition, our analysis suggests that hsa-let-7b-5p may interact with ZC3HAV1 differently in LCLs compared to B cells. In vitro assays indicated that the protective allele of MS risk SNP rs10271373 increases ZC3HAV1 expression in LCLs, but not in B cells. The higher expression for the protective allele in LCLs is consistent with increased IFN response via ZC3HAV1 and so decreased immune evasion by EBV. Taken together, this provides evidence that EBV infection dysregulates the B cell miRNA machinery, including MS risk miRNAs, which may contribute to MS pathogenesis via interaction with MS risk genes either directly or indirectly.

Here we developed KARAJ, a fast and flexible Linux command-line tool to automate the end-to-end process of querying and downloading a wide range of genomic and transcriptomic sequence data types. The input to KARAJ is a list of PMCIDs or publication URLs or various types of accession numbers to automate four tasks as follows; firstly, it provides a summary list of accessible datasets generated by or used in these scientific articles, enabling users to select appropriate datasets; secondly, KARAJ calculates the size of files that users want to download and confirms the availability of adequate space on the local disk; thirdly, it generates a metadata table containing sample information and the experimental design of the corresponding study; and lastly, it enables users to download supplementary data tables attached to publications. Further, KARAJ provides a parallel downloading framework powered by Aspera connect which reduces the downloading time significantly.

Genome wide association studies have identified >200 susceptibility loci accounting for much of the heritability of Multiple Sclerosis (MS). Epstein Barr virus (EBV), a memory B cell tropic virus, has been identified as necessary but not sufficient for development of MS, with evidence for disease causation. The molecular and immunological basis for this has not been established. LCL proliferation is driven by signalling through the EBV produced cell surface protein LMP1, a homologue of the MS risk gene CD40. We show that the CD40 ligand, CD40L, potentially through competitive signalling with LMP1, reduces LCL proliferation (p<0.001). The MS risk variants of the LMP1 signalling inhibitor, TRAF3, had altered expression in B cells and LCLs. Both CD40 and TRAF3 risk SNPs are in binding sites for the EBV transcription factor EBNA2. We have investigated transcriptomes of B cells and EBV infected B cells at Latency III (LCLs) and identified 47 MS risk genes with altered expression, associated with the risk genotype. Overall these MS risk SNPs were overrepresented in target loci of the EBV transcription factor EBNA2 (p<10 ⁻¹⁶ ), in genes dysregulated between B and LCLs (p<10 ⁻⁵ ), and as targets for EBV miRNAs (p<10 ⁻⁴ ). The risk gene ZC3HAV1 is the putative target for multiple EBV miRNAs. It amplifies the interferon response, and was shown to have reduced expression in LCLs for the risk allele. These data indicate targeting EBV EBNA2, miRNAs, and MS risk genes on the LMP1/LMP2 pathways, and the pathways themselves, may be of therapeutic benefit in MS.