MSPrecise: A molecular diagnostic test for multiple sclerosis using next generation sequencing

The clinical applications of immunosequencing

Technological advances in high-throughput sequencing have opened the door for the interrogation of adaptive immune responses at unprecedented scale. It is now possible to determine the sequences of antibodies or T-cell receptors produced by individual B and T cells in a sample. This capability, termed immunosequencing, has transformed the study of both infectious and non-infectious diseases by allowing the tracking of dynamic changes in B and T cell clonal populations over time. This has improved our understanding of the pathology of cancers, autoimmune diseases, and infectious diseases. However, to date there has been only limited clinical adoption of the technology. Advances over the last decade and on the horizon that reduce costs and improve interpretability could enable widespread clinical use. Many clinical applications have been proposed and, while most are still undergoing research and development, some methods relying on immunosequencing data have been implemented, the most widespread of which is the detection of measurable residual disease. Here, we review the diagnostic, prognostic, and therapeutic applications of immunosequencing for both infectious and non-infectious diseases.

Clonally expanded B cells in multiple sclerosis bind EBV EBNA1 and GlialCAM

Multiple sclerosis (MS) is a heterogenous autoimmune disease in which autoreactive lymphocytes attack the myelin sheath of the central nervous system. B lymphocytes in the cerebrospinal fluid (CSF) of patients with MS contribute to inflammation and secrete oligoclonal immunoglobulins^1,2. Epstein-Barr virus (EBV) infection has been epidemiologically linked to MS, but its pathological role remains unclear³. Here we demonstrate high-affinity molecular mimicry between the EBV transcription factor EBV nuclear antigen 1 (EBNA1) and the central nervous system protein glial cell adhesion molecule (GlialCAM) and provide structural and in vivo functional evidence for its relevance. A cross-reactive CSF-derived antibody was initially identified by single-cell sequencing of the paired-chain B cell repertoire of MS blood and CSF, followed by protein microarray-based testing of recombinantly expressed CSF-derived antibodies against MS-associated viruses. Sequence analysis, affinity measurements and the crystal structure of the EBNA1-peptide epitope in complex with the autoreactive Fab fragment enabled tracking of the development of the naive EBNA1-restricted antibody to a mature EBNA1-GlialCAM cross-reactive antibody. Molecular mimicry is facilitated by a post-translational modification of GlialCAM. EBNA1 immunization exacerbates disease in a mouse model of MS, and anti-EBNA1 and anti-GlialCAM antibodies are prevalent in patients with MS. Our results provide a mechanistic link for the association between MS and EBV and could guide the development of new MS therapies.

LncRNA MIAT enhances systemic lupus erythematosus by upregulating CFHR5 expression via miR-222 degradation

Systemic lupus erythematosus (SLE), a complex polygenic autoimmune disease, is associated with increased complement activation. Complement factor H related protein 5 (CFHR5) may contribute to dysfunctional complement activation, thus predisposing to SLE. The expression levels of anti-dsDNA, C3 and CFHR5 in blood samples from 50 SLE patients and 50 healthy individuals were evaluated, and also their expression levels were measured in an MRL/lpr mouse model and control MRL/MPJ mice. The results showed that CFHR5 expression increased in SLE patients together with the increase of anti-dsDNA in comparison with the healthy control. Furthermore, CFHR5 expression was inversely correlated with C3, down-regulation of which was associated with worse SLE. Previous studies indicated that long noncoding RNA (lncRNA) regulates mRNA synthesis via microRNA (miRNA) inhibition. The present bioinformatics analysis revealed that the target miRNA (miR-222) was combined with both lncRNA MIAT and mRNA CFHR5. H&E staining of the kidney tissues of the MRL/lpr mice revealed that lncRNA MIAT, as a competitive inhibitor of miR-222, enhanced SLE by upregulating CFHR5 expression through the degradation of miR-222 in vivo. Thus, our study revealed for the first time the role of lncRNA MIAT in regulating CFHR5 expression in SLE in vivo and provided new insights into the role of lncRNA in regulation and complement function of SLE pathogenesis.

Evaluating Infectious, Neoplastic, Immunological, and Degenerative Diseases of the Central Nervous System with Cerebrospinal Fluid-Based Next-Generation Sequencing

Cerebrospinal fluid (CSF) is a clear and paucicellular fluid that circulates within the ventricular system and the subarachnoid space of the central nervous system (CNS), and diverse CNS disorders can impact its composition, volume, and flow. As conventional CSF testing suffers from suboptimal sensitivity, this review aimed to evaluate the role of next-generation sequencing (NGS) in the work-up of infectious, neoplastic, neuroimmunological, and neurodegenerative CNS diseases. Metagenomic NGS showed improved sensitivity—compared to traditional methods—to detect bacterial, viral, parasitic, and fungal infections, while the overall performance was maximized in some studies when all diagnostic modalities were used. In patients with primary CNS cancer, NGS findings in the CSF were largely concordant with the molecular signatures derived from tissue-based molecular analysis; of interest, additional mutations were identified in the CSF in some glioma studies, reflecting intratumoral heterogeneity. In patients with metastasis to the CNS, NGS facilitated diagnosis, prognosis, therapeutic management, and monitoring, exhibiting higher sensitivity than neuroimaging, cytology, and plasma-based molecular analysis. Although evidence is still rudimentary, NGS could enhance the diagnosis and pathogenetic understanding of multiple sclerosis in addition to Alzheimer and Parkinson disease. To conclude, NGS has shown potential to aid the research, facilitate the diagnostic approach, and improve the management outcomes of all the aforementioned CNS diseases. However, to establish its role in clinical practice, the clinical validity and utility of each NGS protocol should be determined. Lastly, as most evidence has been derived from small and retrospective studies, results from randomized control trials could be of significant value.

Single-Cell High-Throughput Technologies in Cerebrospinal Fluid Research and Diagnostics

High-throughput single-cell technologies have recently emerged as essential tools in biomedical research with great potential for clinical pathology when studying liquid and solid biopsies. We provide an update on current single-cell methods in cerebrospinal fluid research and diagnostics, focusing on high-throughput cell-type specific proteomic and genomic technologies. Proteomic methods comprising flow cytometry and mass cytometry as well as genomic approaches including immune cell repertoire and single-cell transcriptomic studies are critically reviewed and future directions discussed.

Persistence of intrathecal oligoclonal B cells and IgG in multiple sclerosis

In multiple sclerosis (MS), B cells are trafficking across the blood-brain barrier, but it is not known how this relates to the synthesis of oligoclonal IgG. We used quantitative mass spectrometry of oligoclonal bands and high-throughput sequencing of immunoglobulin heavy-chain variable transcripts to study the longitudinal B cell response in the cerebrospinal fluid (CSF) and blood of two MS patients. Twenty of 22 (91%) and 25 of 28 (89%) of oligoclonal band peptides persisted in samples collected 18 months apart, in spite of a dynamic exchange across the blood-CSF barrier of B lineage cells connecting to oligoclonal IgG.

VDJServer: A Cloud-Based Analysis Portal and Data Commons for Immune Repertoire Sequences and Rearrangements

Background

Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation.

Results

VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity.

Conclusion

VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.

VDJPipe: a pipelined tool for pre-processing immune repertoire sequencing data

Background

Pre-processing of high-throughput sequencing data for immune repertoire profiling is essential to insure high quality input for downstream analysis. VDJPipe is a flexible, high-performance tool that can perform multiple pre-processing tasks with just a single pass over the data files.

Results

Processing tasks provided by VDJPipe include base composition statistics calculation, read quality statistics calculation, quality filtering, homopolymer filtering, length and nucleotide filtering, paired-read merging, barcode demultiplexing, 5′ and 3′ PCR primer matching, and duplicate reads collapsing. VDJPipe utilizes a pipeline approach whereby multiple processing steps are performed in a sequential workflow, with the output of each step passed as input to the next step automatically. The workflow is flexible enough to handle the complex barcoding schemes used in many immunosequencing experiments. Because VDJPipe is designed for computational efficiency, we evaluated this by comparing execution times with those of pRESTO, a widely-used pre-processing tool for immune repertoire sequencing data. We found that VDJPipe requires <10% of the run time required by pRESTO.

Conclusions

VDJPipe is a high-performance tool that is optimized for pre-processing large immune repertoire sequencing data sets.

Statistical classifiers for diagnosing disease from immune repertoires: a case study using multiple sclerosis

Background

Deep sequencing of lymphocyte receptor repertoires has made it possible to comprehensively profile the clonal composition of lymphocyte populations. This opens the door for novel approaches to diagnose and prognosticate diseases with a driving immune component by identifying repertoire sequence patterns associated with clinical phenotypes. Indeed, recent studies support the feasibility of this, demonstrating an association between repertoire-level summary statistics (e.g., diversity) and patient outcomes for several diseases. In our own prior work, we have shown that six codons in VH4-containing genes in B cells from the cerebrospinal fluid of patients with relapsing remitting multiple sclerosis (RRMS) have higher replacement mutation frequencies than observed in healthy controls or patients with other neurological diseases. However, prior methods to date have been limited to focusing on repertoire-level summary statistics, ignoring the vast amounts of information in the millions of individual immune receptors comprising a repertoire. We have developed a novel method that addresses this limitation by using innovative approaches for accommodating the extraordinary sequence diversity of immune receptors and widely used machine learning approaches. We applied our method to RRMS, an autoimmune disease that is notoriously difficult to diagnose.

Results

We use the biochemical features encoded by the complementarity determining region 3 of each B cell receptor heavy chain in every patient repertoire as input to a detector function, which is fit to give the correct diagnosis for each patient using maximum likelihood optimization methods. The resulting statistical classifier assigns patients to one of two diagnosis categories, RRMS or other neurological disease, with 87% accuracy by leave-one-out cross-validation on training data (N = 23) and 72% accuracy on unused data from a separate study (N = 102).

Conclusions

Our method is the first to apply statistical learning to immune repertoires to aid disease diagnosis, learning repertoire-level labels from the set of individual immune repertoire sequences. This method produced a repertoire-based statistical classifier for diagnosing RRMS that provides a high degree of diagnostic capability, rivaling the accuracy of diagnosis by a clinical expert. Additionally, this method points to a diagnostic biochemical motif in the antibodies of RRMS patients, which may offer insight into the disease process.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1814-6) contains supplementary material, which is available to authorized users.

Experimental Autoimmune Encephalomyelitis (EAE)-Induced Elevated Expression of the E1 Isoform of Methyl CpG Binding Protein 2 (MeCP2E1): Implications in Multiple Sclerosis (MS)-Induced Neurological Disability and Associated Myelin Damage

Multiple sclerosis (MS) is a chronic neurological disease characterized by the destruction of central nervous system (CNS) myelin. At present, there is no cure for MS due to the inability to repair damaged myelin. Although the neurotrophin brain derived neurotrophic factor (BDNF) has a beneficial role in myelin repair, these effects may be hampered by the over-expression of a transcriptional repressor isoform of methyl CpG binding protein 2 (MeCP2) called MeCP2E1. We hypothesize that following experimental autoimmune encephalomyelitis (EAE)-induced myelin damage, the immune system induction of the pathogenic MeCP2E1 isoform hampers the myelin repair process by repressing BDNF expression. Using an EAE model of MS, we identify the temporal gene and protein expression changes of MeCP2E1, MeCP2E2 and BDNF. The expression changes of these key biological targets were then correlated with the temporal changes in neurological disability scores (NDS) over the entire disease course. Our results indicate that MeCP2E1 mRNA levels are elevated in EAE animals relative to naïve control (NC) and active control (AC) animals during all time points of disease progression. Our results suggest that the EAE-induced elevations in MeCP2E1 expression contribute to the repressed BDNF production in the spinal cord (SC). The sub-optimal levels of BDNF result in sustained NDS and associated myelin damage throughout the entire disease course. Conversely, we observed no significant differences in the expression patterns displayed for the MeCP2E2 isoform amongst our experimental groups. However, our results demonstrate that baseline protein expression ratios between the MeCP2E1 versus MeCP2E2 isoforms in the SC are higher than those identified within the dorsal root ganglia (DRG). Thus, the DRG represents a more conducive environment than that of the SC for BDNF production and transport to the CNS to assist in myelin repair. Henceforth, the sub-optimal BDNF levels we report in the SC may arise from the elevated MeCP2E1 vs. MeCP2E2 ratio in the SC that creates a more hostile environment thereby preventing local BDNF production. At the level of transcript, we demonstrate that EAE-induces the pathological enhanced expression of MeCP2E1 that contributes to enhanced NDS during the entire disease course. Thus, the pathological induction of the MeCP2E1 isoform contributes to the disruption of the normal homeostatic signaling equilibrium network that exists between cytokines, neurotrophins and chemokines that regulate the myelin repair process by repressing BDNF. Our research suggests that the elevated ratio of MeCP2E1 relative to MeCP2E2 may be a useful diagnostic marker that clinicians can utilize to determine the degree of neurological disability with associated myelin damage. The elevated MeCP2E1 vs. MeCP2E2 ratios (E1/E2) in the SC prevent BDNF from reaching optimal levels required for myelin repair. Thus, the lower E1/E2 ratios in the DRG, allow the DRG to serve as a weak secondary compensatory mechanism for enhanced production and delivery of BDNF to the SC to try to assist in myelin repair.

Viruses and Multiple Sclerosis: From Mechanisms and Pathways to Translational Research Opportunities

Viruses are directly or indirectly implicated in multiple sclerosis (MS). Here, we review the evidence on the virus-related pathophysiology of MS, introduce common experimental models, and explore the ways in which viruses cause demyelination. By emphasizing knowledge gaps, we highlight future research directions for effective MS diagnostics and therapies: (i) identifying biomarkers for at-risk individuals, (ii) searching for direct evidence of specific causative viruses, (iii) establishing the contribution of host genetic factors and viruses, and (iv) investigating the contribution of immune regulation at extra-CNS sites. Research in these areas is likely to be facilitated by the application of high-throughput technologies, the development of systems-based bioinformatic approaches, careful selection of experimental models, and the acquisition of high-quality clinical material for tissue-based research.

Peripheral VH4+ plasmablasts demonstrate autoreactive B cell expansion toward brain antigens in early multiple sclerosis patients

Plasmablasts are a highly differentiated, antibody secreting B cell subset whose prevalence correlates with disease activity in Multiple Sclerosis (MS). For most patients experiencing partial transverse myelitis (PTM), plasmablasts are elevated in the blood at the first clinical presentation of disease (known as a clinically isolated syndrome or CIS). In this study we found that many of these peripheral plasmablasts are autoreactive and recognize primarily gray matter targets in brain tissue. These plasmablasts express antibodies that over-utilize immunoglobulin heavy chain V-region subgroup 4 (VH4) genes, and the highly mutated VH4+ plasmablast antibodies recognize intracellular antigens of neurons and astrocytes. Most of the autoreactive, highly mutated VH4+ plasmablast antibodies recognize only a portion of cortical neurons, indicating that the response may be specific to neuronal subgroups or layers. Furthermore, CIS-PTM patients with this plasmablast response also exhibit modest reactivity toward neuroantigens in the plasma IgG antibody pool. Taken together, these data indicate that expanded VH4+ peripheral plasmablasts in early MS patients recognize brain gray matter antigens. Peripheral plasmablasts may be participating in the autoimmune response associated with MS, and provide an interesting avenue for investigating the expansion of autoreactive B cells at the time of the first documented clinical event.

High-throughput sequencing of immune repertoires in multiple sclerosis

T cells and B cells are crucial in the initiation and maintenance of multiple sclerosis (MS), and the activation of these cells is believed to be mediated through specific recognition of antigens by the T‐ and B‐cell receptors. The antigen receptors are highly polymorphic due to recombination (T‐ and B‐cell receptors) and mutation (B‐cell receptors) of the encoding genes, which can therefore be used as fingerprints to track individual T‐ and B‐cell clones. Such studies can shed light on mechanisms driving the immune responses and provide new insights into the pathogenesis. Here, we summarize studies that have explored the T‐ and B‐cell receptor repertoires using earlier methodological approaches, and we focus on how high‐throughput sequencing has provided new knowledge by surveying the immune repertoires in MS in even greater detail and with unprecedented depth.

Practical guidelines for B-cell receptor repertoire sequencing analysis

High-throughput sequencing of B-cell immunoglobulin repertoires is increasingly being applied to gain insights into the adaptive immune response in healthy individuals and in those with a wide range of diseases. Recent applications include the study of autoimmunity, infection, allergy, cancer and aging. As sequencing technologies continue to improve, these repertoire sequencing experiments are producing ever larger datasets, with tens- to hundreds-of-millions of sequences. These data require specialized bioinformatics pipelines to be analyzed effectively. Numerous methods and tools have been developed to handle different steps of the analysis, and integrated software suites have recently been made available. However, the field has yet to converge on a standard pipeline for data processing and analysis. Common file formats for data sharing are also lacking. Here we provide a set of practical guidelines for B-cell receptor repertoire sequencing analysis, starting from raw sequencing reads and proceeding through pre-processing, determination of population structure, and analysis of repertoire properties. These include methods for unique molecular identifiers and sequencing error correction, V(D)J assignment and detection of novel alleles, clonal assignment, lineage tree construction, somatic hypermutation modeling, selection analysis, and analysis of stereotyped or convergent responses. The guidelines presented here highlight the major steps involved in the analysis of B-cell repertoire sequencing data, along with recommendations on how to avoid common pitfalls.

A Distinct Class of Antibodies May Be an Indicator of Gray Matter Autoimmunity in Early and Established Relapsing Remitting Multiple Sclerosis Patients

*These authors contributed equally to the work in this manuscript.We have previously identified a distinct class of antibodies expressed by B cells in the cerebrospinal fluid (CSF) of early and established relapsing remitting multiple sclerosis (RRMS) patients that is not observed in healthy donors. These antibodies contain a unique pattern of mutations in six codons along VH4 antibody genes that we termed the antibody gene signature (AGS). In fact, patients who have such B cells in their CSF are identified as either having RRMS or developing RRMS in the future. As mutations in antibody genes increase antibody affinity for particular antigens, the goal for this study was to investigate whether AGS(+) antibodies bind to brain tissue antigens. Single B cells were isolated from the CSF of 10 patients with early or established RRMS. We chose 32 of these B cells that expressed antibodies enriched for the AGS for further study. We generated monoclonal full-length recombinant human antibodies (rhAbs) and used both immunological assays and immunohistochemistry to investigate the capacity of these AGS(+) rhAbs to bind brain tissue antigens. AGS(+) rhAbs did not recognize myelin tracts in the corpus callosum. Instead, AGS(+) rhAbs recognized neuronal nuclei and/or astrocytes, which are prevalent in the cortical gray matter. This pattern was unique to the AGS(+) antibodies from early and established RRMS patients, as AGS(+) antibodies from an early neuromyelitis optica patient did not display the same reactivity. Prevalence of CSF-derived B cells expressing AGS(+) antibodies that bind to these cell types may be an indicator of gray matter-directed autoimmunity in early and established RRMS patients.