Jason A. Vander Heiden

B cells populating the multiple sclerosis brain mature in the draining cervical lymph nodes

In multiple sclerosis patients, B cells mature in the draining cervical lymph nodes before trafficking across the blood-brain barrier. B Cells Flip the Switch for MS B cells in multiple sclerosis (MS) patients may mature outside the central nervous system (CNS). Two complementary studies report that B cells found outside the CNS—in peripheral blood and draining cervical lymph nodes (CLNs)—share antigen specificity with intrathecal B cell repertoires. In patients with MS, immune cells attack the CNS; however, it remains unclear whether these cells mature in the CNS or traffic to the CNS as mature cells. Using paired tissues and high-throughput sequencing, Stern et al. found that clonally expanded B cells are found in both the CNS and CLNs but that founding members were more often found in the draining CLNs. Palanichamy et al. extend these findings by reporting a peripheral blood/CNS axis of mature B cells that have undergone class switch. These data support the therapeutic use of monoclonal antibodies that prevent lymphocytes from crossing the blood-brain barrier or induce peripheral B cell depletion in MS patients. Multiple sclerosis (MS) is an inflammatory disease of the central nervous system (CNS) characterized by autoimmune-mediated demyelination and neurodegeneration. The CNS of patients with MS harbors expanded clones of antigen-experienced B cells that reside in distinct compartments including the meninges, cerebrospinal fluid (CSF), and parenchyma. It is not understood whether this immune infiltrate initiates its development in the CNS or in peripheral tissues. B cells in the CSF can exchange with those in peripheral blood, implying that CNS B cells may have access to lymphoid tissue that may be the specific compartment(s) in which CNS-resident B cells encounter antigen and experience affinity maturation. Paired tissues were used to determine whether the B cells that populate the CNS mature in the draining cervical lymph nodes (CLNs). High-throughput sequencing of the antibody repertoire demonstrated that clonally expanded B cells were present in both compartments. Founding members of clones were more often found in the draining CLNs. More mature clonal members derived from these founders were observed in the draining CLNs and also in the CNS, including lesions. These data provide new evidence that B cells traffic freely across the tissue barrier, with the majority of B cell maturation occurring outside of the CNS in the secondary lymphoid tissue. Our study may aid in further defining the mechanisms of immunomodulatory therapies that either deplete circulating B cells or affect the intrathecal B cell compartment by inhibiting lymphocyte transmigration into the CNS.

Change-O: a toolkit for analyzing large-scale B cell immunoglobulin repertoire sequencing data

UNLABELLED Advances in high-throughput sequencing technologies now allow for large-scale characterization of B cell immunoglobulin (Ig) repertoires. The high germline and somatic diversity of the Ig repertoire presents challenges for biologically meaningful analysis, which requires specialized computational methods. We have developed a suite of utilities, Change-O, which provides tools for advanced analyses of large-scale Ig repertoire sequencing data. Change-O includes tools for determining the complete set of Ig variable region gene segment alleles carried by an individual (including novel alleles), partitioning of Ig sequences into clonal populations, creating lineage trees, inferring somatic hypermutation targeting models, measuring repertoire diversity, quantifying selection pressure, and calculating sequence chemical properties. All Change-O tools utilize a common data format, which enables the seamless integration of multiple analyses into a single workflow. AVAILABILITY AND IMPLEMENTATION Change-O is freely available for non-commercial use and may be downloaded from http://clip.med.yale.edu/changeo. CONTACT steven.kleinstein@yale.edu.

High-resolution antibody dynamics of vaccine-induced immune responses

Significance The immune system must constantly adapt to combat infections and other challenges. This is accomplished by continuously evolving the antibody repertoire, and by maintaining memory of prior challenges. By using next-generation DNA sequencing technology, we have examined the shear amount of antibody made by individuals during a flu vaccination trial. We demonstrate one of the first characterizations of the fast antibody dynamics through time in multiple individuals responding to an immune challenge. The adaptive immune system confers protection by generating a diverse repertoire of antibody receptors that are rapidly expanded and contracted in response to specific targets. Next-generation DNA sequencing now provides the opportunity to survey this complex and vast repertoire. In the present work, we describe a set of tools for the analysis of antibody repertoires and their application to elucidating the dynamics of the response to viral vaccination in human volunteers. By analyzing data from 38 separate blood samples across 2 y, we found that the use of the germ-line library of V and J segments is conserved between individuals over time. Surprisingly, there appeared to be no correlation between the use level of a particular VJ combination and degree of expansion. We found the antibody RNA repertoire in each volunteer to be highly dynamic, with each individual displaying qualitatively different response dynamics. By using combinatorial phage display, we screened selected VH genes paired with their corresponding VL library for affinity against the vaccine antigens. Altogether, this work presents an additional set of tools for profiling the human antibody repertoire and demonstrates characterization of the fast repertoire dynamics through time in multiple individuals responding to an immune challenge.

pRESTO: a toolkit for processing high-throughput sequencing raw reads of lymphocyte receptor repertoires

UNLABELLED Driven by dramatic technological improvements, large-scale characterization of lymphocyte receptor repertoires via high-throughput sequencing is now feasible. Although promising, the high germline and somatic diversity, especially of B-cell immunoglobulin repertoires, presents challenges for analysis requiring the development of specialized computational pipelines. We developed the REpertoire Sequencing TOolkit (pRESTO) for processing reads from high-throughput lymphocyte receptor studies. pRESTO processes raw sequences to produce error-corrected, sorted and annotated sequence sets, along with a wealth of metrics at each step. The toolkit supports multiplexed primer pools, single- or paired-end reads and emerging technologies that use single-molecule identifiers. pRESTO has been tested on data generated from Roche and Illumina platforms. It has a built-in capacity to parallelize the work between available processors and is able to efficiently process millions of sequences generated by typical high-throughput projects. AVAILABILITY AND IMPLEMENTATION pRESTO is freely available for academic use. The software package and detailed tutorials may be downloaded from http://clip.med.yale.edu/presto.

Models of Somatic Hypermutation Targeting and Substitution Based on Synonymous Mutations from High-Throughput Immunoglobulin Sequencing Data

Analyses of somatic hypermutation (SHM) patterns in B cell immunoglobulin (Ig) sequences contribute to our basic understanding of adaptive immunity, and have broad applications not only for understanding the immune response to pathogens, but also to determining the role of SHM in autoimmunity and B cell cancers. Although stochastic, SHM displays intrinsic biases that can confound statistical analysis, especially when combined with the particular codon usage and base composition in Ig sequences. Analysis of B cell clonal expansion, diversification, and selection processes thus critically depends on an accurate background model for SHM micro-sequence targeting (i.e., hot/cold-spots) and nucleotide substitution. Existing models are based on small numbers of sequences/mutations, in part because they depend on data from non-coding regions or non-functional sequences to remove the confounding influences of selection. Here, we combine high-throughput Ig sequencing with new computational analysis methods to produce improved models of SHM targeting and substitution that are based only on synonymous mutations, and are thus independent of selection. The resulting “S5F” models are based on 806,860 Synonymous mutations in 5-mer motifs from 1,145,182 Functional sequences and account for dependencies on the adjacent four nucleotides (two bases upstream and downstream of the mutation). The estimated profiles can explain almost half of the variance in observed mutation patterns, and clearly show that both mutation targeting and substitution are significantly influenced by neighboring bases. While mutability and substitution profiles were highly conserved across individuals, the variability across motifs was found to be much larger than previously estimated. The model and method source code are made available at http://clip.med.yale.edu/SHM

Salmonella Infection Drives Promiscuous B Cell Activation Followed by Extrafollicular Affinity Maturation

The B cell response to Salmonella typhimurium (STm) occurs massively at extrafollicular sites, without notable germinal centers (GCs). Little is known in terms of its specificity. To expand the knowledge of antigen targets, we screened plasmablast (PB)-derived monoclonal antibodies (mAbs) for Salmonella specificity, using ELISA, flow cytometry, and antigen microarray. Only a small fraction (0.5%-2%) of the response appeared to be Salmonella-specific. Yet, infection of mice with limited B cell receptor (BCR) repertoires impaired the response, suggesting that BCR specificity was important. We showed, using laser microdissection, that somatic hypermutation (SHM) occurred efficiently at extrafollicular sites leading to affinity maturation that in turn led to detectable STm Ag-binding. These results suggest a revised vision of how clonal selection and affinity maturation operate in response to Salmonella. Clonal selection initially is promiscuous, activating cells with virtually undetectable affinity, yet SHM and selection occur during the extrafollicular response yielding higher affinity, detectable antibodies.

Shared genetic basis for migraine and ischemic stroke A genome-wide analysis of common variants

Objective: To quantify genetic overlap between migraine and ischemic stroke (IS) with respect to common genetic variation. Methods: We applied 4 different approaches to large-scale meta-analyses of genome-wide data on migraine (23,285 cases and 95,425 controls) and IS (12,389 cases and 62,004 controls). First, we queried known genome-wide significant loci for both disorders, looking for potential overlap of signals. We then analyzed the overall shared genetic load using polygenic scores and estimated the genetic correlation between disease subtypes using data derived from these models. We further interrogated genomic regions of shared risk using analysis of covariance patterns between the 2 phenotypes using cross-phenotype spatial mapping. Results: We found substantial genetic overlap between migraine and IS using all 4 approaches. Migraine without aura (MO) showed much stronger overlap with IS and its subtypes than migraine with aura (MA). The strongest overlap existed between MO and large artery stroke (LAS; p = 6.4 × 10−28 for the LAS polygenic score in MO) and between MO and cardioembolic stroke (CE; p = 2.7 × 10−20 for the CE score in MO). Conclusions: Our findings indicate shared genetic susceptibility to migraine and IS, with a particularly strong overlap between MO and both LAS and CE pointing towards shared mechanisms. Our observations on MA are consistent with a limited role of common genetic variants in this subtype.

Individual heritable differences result in unique cell lymphocyte receptor repertoires of naive and antigen-experienced cells

The adaptive immune systems capability to protect the body requires a highly diverse lymphocyte antigen receptor repertoire. However, the influence of individual genetic and epigenetic differences on these repertoires is not typically measured. By leveraging the unique characteristics of B, CD4+ T and CD8+ T-lymphocyte subsets from monozygotic twins, we quantify the impact of heritable factors on both the V(D)J recombination process and on thymic selection. We show that the resulting biases in both V(D)J usage and N/P addition lengths, which are found in naïve and antigen experienced cells, contribute to significant variation in the CDR3 region. Moreover, we show that the relative usage of V and J gene segments is chromosomally biased, with ∼1.5 times as many rearrangements originating from a single chromosome. These data refine our understanding of the heritable mechanisms affecting the repertoire, and show that biases are evident on a chromosome-wide level.

Influence of seasonal exposure to grass pollen on local and peripheral blood IgE repertoires in patients with allergic rhinitis

Background Previous studies of immunoglobulin gene sequences in patients with allergic diseases using low-throughput Sanger sequencing have limited the analytic depth for characterization of IgE repertoires. Objectives We used a high-throughput, next-generation sequencing approach to characterize immunoglobulin heavy-chain gene (IGH) repertoires in patients with seasonal allergic rhinitis (AR) with the aim of better understanding the underlying disease mechanisms. Methods IGH sequences in matched peripheral blood and nasal biopsy specimens from nonallergic healthy control subjects (n = 3) and patients with grass pollen–related AR taken in season (n = 3) or out of season (n = 4) were amplified and pyrosequenced on the 454 GS FLX+ System. Results A total of 97,610 IGH (including 8,135 IgE) sequences were analyzed. Use of immunoglobulin heavy-chain variable region gene families 1 (IGHV1) and 5 (IGHV5) was higher in IgE clonotypic repertoires compared with other antibody classes independent of atopic status. IgE repertoires measured inside the grass pollen season were more diverse and more mutated (particularly in the biopsy specimens) and had more evidence of antigen-driven selection compared with those taken outside of the pollen season or from healthy control subjects. Clonal relatedness was observed for IgE between the blood and nasal biopsy specimens. Furthermore in patients with AR, but not healthy control subjects, we found clonal relatedness between IgE and IgG classes. Conclusion This is the first report that exploits next-generation sequencing to determine local and peripheral blood IGH repertoires in patients with respiratory allergic disease. We demonstrate that natural pollen exposure was associated with changes in IgE repertoires that were suggestive of ongoing germinal center reactions. Furthermore, these changes were more often apparent in nasal biopsy specimens compared with peripheral blood and in patients with AR compared with healthy control subjects.

Genetic analysis for a shared biological basis between migraine and coronary artery disease

Objective: To apply genetic analysis of genome-wide association data to study the extent and nature of a shared biological basis between migraine and coronary artery disease (CAD). Methods: Four separate methods for cross-phenotype genetic analysis were applied on data from 2 large-scale genome-wide association studies of migraine (19,981 cases, 56,667 controls) and CAD (21,076 cases, 63,014 controls). The first 2 methods quantified the extent of overlapping risk variants and assessed the load of CAD risk loci in migraineurs. Genomic regions of shared risk were then identified by analysis of covariance patterns between the 2 phenotypes and by querying known genome-wide significant loci. Results: We found a significant overlap of genetic risk loci for migraine and CAD. When stratified by migraine subtype, this was limited to migraine without aura, and the overlap was protective in that patients with migraine had a lower load of CAD risk alleles than controls. Genes indicated by 16 shared risk loci point to mechanisms with potential roles in migraine pathogenesis and CAD, including endothelial dysfunction (PHACTR1) and insulin homeostasis (GIP). Conclusions: The results suggest that shared biological processes contribute to risk of migraine and CAD, but surprisingly this commonality is restricted to migraine without aura and the impact is in opposite directions. Understanding the mechanisms underlying these processes and their opposite relationship to migraine and CAD may improve our understanding of both disorders.