Search

Background: Transforming growth factor β (TGF-β) is implicated as a key mediator of pathological fibrosis, but its pleiotropic activity in a range of homeostatic functions presents challenges to its safe and effective therapeutic targeting. There are three isoforms of TGF-β, TGF-β1, TGF-β2, and TGF-β3, which bind to a common receptor complex composed of TGF-βR1 and TGF-βR2 to induce similar intracellular signals in vitro. We have recently shown that the cellular expression patterns and activation thresholds of TGF-β2 and TGF-β3 are distinct from those of TGF-β1 and that selective short-term TGF-β2 and TGF-β3 inhibition can attenuate fibrosis in vivo without promoting excessive inflammation. Isoform-selective inhibition of TGF-β may therefore provide a therapeutic opportunity for patients with chronic fibrotic disorders., Methods: Transcriptomic profiling of skin biopsies from patients with systemic sclerosis (SSc) from multiple clinical trials was performed to evaluate the role of TGF-β3 in this disease. Antibody humanization, biochemical characterization, crystallization, and pre-clinical experiments were performed to further characterize an anti-TGF-β3 antibody., Findings: In the skin of patients with SSc, TGF-β3 expression is uniquely correlated with biomarkers of TGF-β signaling and disease severity. Crystallographic studies establish a structural basis for selective TGF-β3 inhibition with a potent and selective monoclonal antibody that attenuates fibrosis effectively in vivo at clinically translatable exposures. Toxicology studies suggest that, as opposed to pan-TGF-β inhibitors, this anti-TGF-β3 antibody has a favorable safety profile for chronic administration., Conclusion: We establish a rationale for targeting TGF-β3 in SSc with a favorable therapeutic index., Funding: This study was funded by Genentech, Inc., Competing Interests: Declaration of interests All authors are current or past employees of Genentech, a member of the Roche group, and may hold Roche stock or stock options. A clinical trial sponsored by Genentech, Inc., with the indicated therapeutic agent is currently ongoing (ClinicalTrials.gov: NCT05462522). A patent application titled “Isoform-selective anti-TGF-beta antibodies and methods of use,” relating to the subject matter of this manuscript, has been provisionally filed by Genentech, Inc. Currently, J.R.A. is an employee of Sonoma Biotherapeutics, Z.H. is an employee of Gilead Sciences, Inc., and P.L. is an employee of Synthekine., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Idiopathic pulmonary fibrosis (IPF) is a progressive fibrotic disease characterised by aberrant fibroblast/myofibroblast accumulation and excessive collagen matrix deposition in the alveolar areas of lungs. As the first approved IPF medication, pirfenidone (PFD) significantly decelerates lung function decline while its underlying anti-fibrotic mechanism remains elusive., Methods: We performed transcriptomic and immunofluorescence analyses of primary human IPF tissues., Results: We showed that myocardin-related transcription factor (MRTF) signalling is activated in myofibroblasts accumulated in IPF lungs. Furthermore, we showed that PFD inhibits MRTF activation in primary human lung fibroblasts at clinically achievable concentrations (half-maximal inhibitory concentration 50-150 µM, maximal inhibition >90%, maximal concentration of PFD in patients <100 µM). Mechanistically, PFD appears to exert its inhibitory effects by promoting the interaction between MRTF and actin indirectly. Finally, PFD-treated IPF lungs exhibit significantly less MRTF activation in fibroblast foci areas than naïve IPF lungs., Conclusions: Our results suggest MRTF signalling as a direct target for PFD and implicate that some of the anti-fibrotic effects of PFD may be due to MRTF inhibition in lung fibroblasts., Competing Interests: Conflict of interest: All authors except P.J. Wolters are or were employees of Genentech/Roche., (Copyright ©The authors 2023.)

Fibrosis is the major risk factor associated with morbidity and mortality in patients with non-alcoholic steatohepatitis (NASH)-driven chronic liver disease. Although numerous efforts have been made to identify the mediators of the initiation of liver fibrosis, the molecular underpinnings of fibrosis progression remain poorly understood, and therapies to arrest liver fibrosis progression are elusive. Here, we identify a pathway involving WNT1-inducible signaling pathway protein 1 (WISP1) and myocardin-related transcription factor (MRTF) as a central mechanism driving liver fibrosis progression through the integrin-dependent transcriptional reprogramming of myofibroblast cytoskeleton and motility. In mice, WISP1 deficiency protects against fibrosis progression, but not fibrosis onset. Moreover, the therapeutic administration of a novel antibody blocking WISP1 halted the progression of existing liver fibrosis in NASH models. These findings implicate the WISP1-MRTF axis as a crucial determinant of liver fibrosis progression and support targeting this pathway by antibody-based therapy for the treatment of NASH fibrosis., Competing Interests: Declaration of interests All authors except P.J.W. and A.E.F. are or were employees of Genentech, a member of the Roche group, and may hold Roche stock or stock options., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Partial reprogramming by expression of reprogramming factors (Oct4, Sox2, Klf4 and c-Myc) for short periods of time restores a youthful epigenetic signature to aging cells and extends the life span of a premature aging mouse model. However, the effects of longer-term partial reprogramming in physiologically aging wild-type mice are unknown. Here, we performed various long-term partial reprogramming regimens, including different onset timings, during physiological aging. Long-term partial reprogramming lead to rejuvenating effects in different tissues, such as the kidney and skin, and at the organismal level; duration of the treatment determined the extent of the beneficial effects. The rejuvenating effects were associated with a reversion of the epigenetic clock and metabolic and transcriptomic changes, including reduced expression of genes involved in the inflammation, senescence and stress response pathways. Overall, our observations indicate that partial reprogramming protocols can be designed to be safe and effective in preventing age-related physiological changes. We further conclude that longer-term partial reprogramming regimens are more effective in delaying aging phenotypes than short-term reprogramming., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Transforming growth factor-β (TGFβ) is a key driver of fibrogenesis. Three TGFβ isoforms (TGFβ1, TGFβ2, and TGFβ3) in mammals have distinct functions in embryonic development; however, the postnatal pathological roles and activation mechanisms of TGFβ2 and TGFβ3 have not been well characterized. Here, we show that the latent forms of TGFβ2 and TGFβ3 can be activated by integrin-independent mechanisms and have lower activation thresholds compared to TGFβ1. Unlike TGFB1 , TGFB2 and TGFB3 expression is increased in human lung and liver fibrotic tissues compared to healthy control tissues. Thus, TGFβ2 and TGFβ3 may play a pathological role in fibrosis. Inducible conditional knockout mice and anti-TGFβ isoform-selective antibodies demonstrated that TGFβ2 and TGFβ3 are independently involved in mouse fibrosis models in vivo, and selective TGFβ2 and TGFβ3 inhibition does not lead to the increased inflammation observed with pan-TGFβ isoform inhibition. A cocrystal structure of a TGFβ2-anti-TGFβ2/3 antibody complex reveals an allosteric isoform-selective inhibitory mechanism. Therefore, inhibiting TGFβ2 and/or TGFβ3 while sparing TGFβ1 may alleviate fibrosis without toxicity concerns associated with pan-TGFβ blockade., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Despite the development of effective therapies, a substantial proportion of asthmatics continue to have uncontrolled symptoms, airflow limitation, and exacerbations. Transient receptor potential cation channel member A1 (TRPA1) agonists are elevated in human asthmatic airways, and in rodents, TRPA1 is involved in the induction of airway inflammation and hyperreactivity. Here, the discovery and early clinical development of GDC-0334, a highly potent, selective, and orally bioavailable TRPA1 antagonist, is described. GDC-0334 inhibited TRPA1 function on airway smooth muscle and sensory neurons, decreasing edema, dermal blood flow (DBF), cough, and allergic airway inflammation in several preclinical species. In a healthy volunteer Phase 1 study, treatment with GDC-0334 reduced TRPA1 agonist-induced DBF, pain, and itch, demonstrating GDC-0334 target engagement in humans. These data provide therapeutic rationale for evaluating TRPA1 inhibition as a clinical therapy for asthma., (© 2021 Genentech.)

Progressive lung fibrosis is a major cause of mortality in systemic sclerosis (SSc) patients, but the underlying mechanisms remain unclear. We demonstrate that immune complexes (ICs) activate human monocytes to promote lung fibroblast migration partly via osteopontin (OPN) secretion, which is amplified by autocrine monocyte colony stimulating factor (MCSF) and interleukin-6 (IL-6) activity. Bulk and single-cell RNA sequencing demonstrate that elevated OPN expression in SSc lung tissue is enriched in macrophages, partially overlapping with CCL18 expression. Serum OPN is elevated in SSc patients with interstitial lung disease (ILD) and prognosticates future lung function deterioration in SSc cohorts. Serum OPN levels decrease following tocilizumab (monoclonal anti-IL-6 receptor) treatment, confirming the connection between IL-6 and OPN in SSc patients. Collectively, these data suggest a plausible link between autoantibodies and lung fibrosis progression, where circulating OPN serves as a systemic proxy for IC-driven profibrotic macrophage activity, highlighting its potential as a promising biomarker in SSc ILD., Competing Interests: All authors except P.J.W. and D.K. are current or past employees of Genentech, a member of the Roche group, and may hold Roche stock or stock options. D.K. is a consultant to Actelion, Acceleron, Arena, Bayer, Boehringer Ingelheim, Bristol-Myer Squibb, CSL Behring, Corbus, Cytori, GSK, Genentech/ Roche, Galapagos, Merck, Mitsubishi Tanabi, and UCB. He has received grants as part of investigator-initiated trials (to the University of Michigan) from Bayer, Bristol-Myer Squibb, and Pfizer and has stock options in Eicos Sciences, Inc. A.J. is a former employee of Genentech and is now employed by Gilead Sciences (Foster City, CA, USA)., (© 2020 The Author(s).)

The Adaptive Immune Receptor Repertoire (AIRR) Community is a research-driven group that is establishing a clear set of community-accepted data and metadata standards; standards-based reference implementation tools; and policies and practices for infrastructure to support the deposit, curation, storage, and use of high-throughput sequencing data from B-cell and T-cell receptor repertoires (AIRR-seq data). The AIRR Data Commons is a distributed system of data repositories that utilizes a common data model, a common query language, and common interoperability formats for storage, query, and downloading of AIRR-seq data. Here is described the principal technical standards for the AIRR Data Commons consisting of the AIRR Data Model for repertoires and rearrangements, the AIRR Data Commons (ADC) API for programmatic query of data repositories, a reference implementation for ADC API services, and tools for querying and validating data repositories that support the ADC API. AIRR-seq data repositories can become part of the AIRR Data Commons by implementing the data model and API. The AIRR Data Commons allows AIRR-seq data to be reused for novel analyses and empowers researchers to discover new biological insights about the adaptive immune system., (Copyright © 2020 Christley, Aguiar, Blanck, Breden, Bukhari, Busse, Jaglale, Harikrishnan, Laserson, Peters, Rocha, Schramm, Taylor, Vander Heiden, Zimonja, Watson, Corrie and Cowell.)

In order to produce effective antibodies, B cells undergo rapid somatic hypermutation (SHM) and selection for binding affinity to antigen via a process called affinity maturation. The similarities between this process and evolution by natural selection have led many groups to use phylogenetic methods to characterize the development of immunological memory, vaccination, and other processes that depend on affinity maturation. However, these applications are limited by the fact that most phylogenetic models are designed to be applied to individual lineages comprising genetically diverse sequences, while B cell repertoires often consist of hundreds to thousands of separate low-diversity lineages. Further, several features of affinity maturation violate important assumptions in standard phylogenetic models. Here, we introduce a hierarchical phylogenetic framework that integrates information from all lineages in a repertoire to more precisely estimate model parameters while simultaneously incorporating the unique features of SHM. We demonstrate the power of this repertoire-wide approach by characterizing previously undescribed phenomena in affinity maturation. First, we find evidence consistent with age-related changes in SHM hot-spot targeting. Second, we identify a consistent relationship between increased tree length and signs of increased negative selection, apparent in the repertoires of recently vaccinated subjects and those without any known recent infections or vaccinations. This suggests that B cell lineages shift toward negative selection over time as a general feature of affinity maturation. Our study provides a framework for undertaking repertoire-wide phylogenetic testing of SHM hypotheses and provides a means of characterizing dynamics of mutation and selection during affinity maturation., Competing Interests: The authors declare no competing interest., (Copyright © 2019 the Author(s). Published by PNAS.)

The adaptive immune system generates an incredible diversity of antigen receptors for B and T cells to keep dangerous pathogens at bay. The DNA sequences coding for these receptors arise by a complex recombination process followed by a series of productivity-based filters, as well as affinity maturation for B cells, giving considerable diversity to the circulating pool of receptor sequences. Although these datasets hold considerable promise for medical and public health applications, the complex structure of the resulting adaptive immune receptor repertoire sequencing (AIRR-seq) datasets makes analysis difficult. In this paper we introduce sumrep, an R package that efficiently performs a wide variety of repertoire summaries and comparisons, and show how sumrep can be used to perform model validation. We find that summaries vary in their ability to differentiate between datasets, although many are able to distinguish between covariates such as donor, timepoint, and cell type for BCR and TCR repertoires. We show that deletion and insertion lengths resulting from V(D)J recombination tend to be more discriminative characterizations of a repertoire than summaries that describe the amino acid composition of the CDR3 region. We also find that state-of-the-art generative models excel at recapitulating gene usage and recombination statistics in a given experimental repertoire, but struggle to capture many physiochemical properties of real repertoires., (Copyright © 2019 Olson, Moghimi, Schramm, Obraztsova, Ralph, Vander Heiden, Shugay, Shepherd, Lees and Matsen.)

Neuromyelitis optica spectrum disorders (NMOSD) constitute rare autoimmune disorders of the CNS that are primarily characterized by severe inflammation of the spinal cord and optic nerve. Approximately 75% of NMOSD patients harbour circulating pathogenic autoantibodies targeting the aquaporin-4 water channel (AQP4). The source of these autoantibodies remains unclear, but parallels between NMOSD and other autoantibody-mediated diseases posit compromised B cell tolerance checkpoints as common underlying and contributing factors. Using a well established assay, we assessed tolerance fidelity by creating recombinant antibodies from B cell populations directly downstream of each checkpoint and testing them for polyreactivity and autoreactivity. We examined a total of 863 recombinant antibodies. Those derived from three anti-AQP4-IgG seropositive NMOSD patients (n = 130) were compared to 733 antibodies from 15 healthy donors. We found significantly higher frequencies of poly- and autoreactive new emigrant/transitional and mature naïve B cells in NMOSD patients compared to healthy donors (P-values < 0.003), thereby identifying defects in both central and peripheral B cell tolerance checkpoints in these patients. We next explored whether pathogenic NMOSD anti-AQP4 autoantibodies can originate from the pool of poly- and autoreactive clones that populate the naïve B cell compartment of NMOSD patients. Six human anti-AQP4 autoantibodies that acquired somatic mutations were reverted back to their unmutated germline precursors, which were tested for both binding to AQP4 and poly- or autoreactivity. While the affinity of mature autoantibodies against AQP4 ranged from modest to strong (Kd 15.2-559 nM), none of the germline revertants displayed any detectable binding to AQP4, revealing that somatic hypermutation is required for the generation of anti-AQP4 autoantibodies. However, two (33.3%) germline autoantibody revertants were polyreactive and four (66.7%) were autoreactive, suggesting that pathogenic anti-AQP4 autoantibodies can originate from the pool of autoreactive naïve B cells, which develops as a consequence of impaired early B cell tolerance checkpoints in NMOSD patients., (© The Author(s) (2019). Published by Oxford University Press on behalf of the Guarantors of Brain. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

The adaptive immune receptor repertoire (AIRR) contains information on an individuals' immune past, present and potential in the form of the evolving sequences that encode the B cell receptor (BCR) repertoire. AIRR sequencing (AIRR-seq) studies rely on databases of known BCR germline variable (V), diversity (D), and joining (J) genes to detect somatic mutations in AIRR-seq data via comparison to the best-aligning database alleles. However, it has been shown that these databases are far from complete, leading to systematic misidentification of mutated positions in subsets of sample sequences. We previously presented TIgGER, a computational method to identify subject-specific V gene genotypes, including the presence of novel V gene alleles, directly from AIRR-seq data. However, the original algorithm was unable to detect alleles that differed by more than 5 single nucleotide polymorphisms (SNPs) from a database allele. Here we present and apply an improved version of the TIgGER algorithm which can detect alleles that differ by any number of SNPs from the nearest database allele, and can construct subject-specific genotypes with minimal prior information. TIgGER predictions are validated both computationally (using a leave-one-out strategy) and experimentally (using genomic sequencing), resulting in the addition of three new immunoglobulin heavy chain V (IGHV) gene alleles to the IMGT repertoire. Finally, we develop a Bayesian strategy to provide a confidence estimate associated with genotype calls. All together, these methods allow for much higher accuracy in germline allele assignment, an essential step in AIRR-seq studies.

Increased interest in the immune system's involvement in pathophysiological phenomena coupled with decreased DNA sequencing costs have led to an explosion of antibody and T cell receptor sequencing data collectively termed "adaptive immune receptor repertoire sequencing" (AIRR-seq or Rep-Seq). The AIRR Community has been actively working to standardize protocols, metadata, formats, APIs, and other guidelines to promote open and reproducible studies of the immune repertoire. In this paper, we describe the work of the AIRR Community's Data Representation Working Group to develop standardized data representations for storing and sharing annotated antibody and T cell receptor data. Our file format emphasizes ease-of-use, accessibility, scalability to large data sets, and a commitment to open and transparent science. It is composed of a tab-delimited format with a specific schema. Several popular repertoire analysis tools and data repositories already utilize this AIRR-seq data format. We hope that others will follow suit in the interest of promoting interoperable standards.

Epstein-Barr virus (EBV) ZEBRA protein activates the EBV lytic cycle. Cellular AP-1 proteins with alanine-to-serine [AP-1(A/S)] substitutions homologous to ZEBRA(S186) assume some functions of EBV ZEBRA. These AP-1(A/S) mutants bind methylated EBV DNA and activate expression of some EBV genes. Here, we compare expression of 67 viral genes induced by ZEBRA versus expression induced by AP-1(A/S) proteins. AP-1(A/S) activated 24 genes to high levels and 15 genes to intermediate levels; activation of 28 genes by AP-1(A/S) was severely impaired. We show that AP-1(A/S) proteins are defective at stimulating viral lytic DNA replication. The impairment of expression of many late genes compared to that of ZEBRA is likely due to the inability of AP-1(A/S) proteins to promote viral DNA replication. However, even in the absence of detectable viral DNA replication, AP-1(A/S) proteins stimulated expression of a subgroup of late genes that encode viral structural proteins and immune modulators. In response to ZEBRA, expression of this subgroup of late genes was inhibited by phosphonoacetic acid (PAA), which is a potent viral replication inhibitor. However, when the lytic cycle was activated by AP-1(A/S), PAA did not reduce expression of this subgroup of late genes. We also provide genetic evidence, using the BMRF1 knockout bacmid, that these genes are true late genes in response to ZEBRA. AP-1(A/S) binds to the promoter region of at least one of these late genes, BDLF3, encoding an immune modulator. IMPORTANCE Mutant c-Jun and c-Fos proteins selectively activate expression of EBV lytic genes, including a subgroup of viral late genes, in the absence of viral DNA replication. These findings indicate that newly synthesized viral DNA is not invariably required for viral late gene expression. While viral DNA replication may be obligatory for late gene expression driven by viral transcription factors, it does not limit the ability of cellular transcription factors to activate expression of some viral late genes. Our results show that expression of all late genes may not be strictly dependent on viral lytic DNA replication. The c-Fos A151S mutation has been identified in a human cancer. c-Fos A151S in combination with wild-type c-Jun activates the EBV lytic cycle. Our data provide proof of principle that mutant cellular transcription factors could cause aberrant regulation of viral lytic cycle gene expression and play important roles in EBV-associated diseases., (Copyright © 2018 American Society for Microbiology.)

High-throughput sequencing (HTS) of immunoglobulin (B-cell receptor, antibody) and T-cell receptor repertoires has increased dramatically since the technique was introduced in 2009 (1-3). This experimental approach explores the maturation of the adaptive immune system and its response to antigens, pathogens, and disease conditions in exquisite detail. It holds significant promise for diagnostic and therapy-guiding applications. New technology often spreads rapidly, sometimes more rapidly than the understanding of how to make the products of that technology reliable, reproducible, or usable by others. As complex technologies have developed, scientific communities have come together to adopt common standards, protocols, and policies for generating and sharing data sets, such as the MIAME protocols developed for microarray experiments. The Adaptive Immune Receptor Repertoire (AIRR) Community formed in 2015 to address similar issues for HTS data of immune repertoires. The purpose of this perspective is to provide an overview of the AIRR Community's founding principles and present the progress that the AIRR Community has made in developing standards of practice and data sharing protocols. Finally, and most important, we invite all interested parties to join this effort to facilitate sharing and use of these powerful data sets (join@airr-community.org).

CD4 + follicular regulatory T (T fr ) cells suppress B cell responses through modulation of follicular helper T (T fh ) cells and germinal center (GC) development. We found that T fr cells can also promote the GC response through provision of interleukin-10 (IL-10) after acute infection with lymphocytic choriomeningitis virus (LCMV). Sensing of IL-10 by B cells was necessary for optimal development of the GC response. GC B cells formed in the absence of T reg cell-derived IL-10 displayed an altered dark zone state and decreased expression of the transcription factor Forkhead box protein 1 (FOXO1). IL-10 promoted nuclear translocation of FOXO1 in activated B cells. These data indicate that T fr cells play a multifaceted role in the fine-tuning of the GC response and identify IL-10 as an important mediator by which T fr cells support the GC reaction., (Copyright © 2017 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Background: The B and T cells of the human adaptive immune system leverage a highly diverse repertoire of antigen-specific receptors to protect the human body from pathogens. The sequencing and analysis of immune repertoires is emerging as an important tool to understand immune responses, whether beneficial or harmful (in the case of autoimmunity). However, methods for studying these repertoires, and for directly comparing different immune repertoires, are lacking., Results: In this paper, we present a non-parametric method for directly comparing sequencing repertoires, with the goal of rigorously quantifying differences in V, D, and J gene segment utilization. This method, referred to as the Repertoire Dissimilarity Index (RDI), uses a bootstrapped subsampling approach to account for variance in sequencing depth, and, coupled with a data simulation approach, allows for direct quantification of the average variation between repertoires. We use the RDI method to recapitulate known differences in the formation of the CD4 + and CD8 + T cell repertoires, and further show that antigen-driven activation of naïve CD8 + T cells is more selective than in the CD4 + repertoire, resulting in a more specialized CD8 + memory repertoire., Conclusions: We prove that the RDI method is an accurate and versatile method for comparisons of immune repertoires. The RDI method has been implemented as an R package, and is available for download through Bitbucket.

Myasthenia gravis (MG) is a prototypical B cell-mediated autoimmune disease affecting 20-50 people per 100,000. The majority of patients fall into two clinically distinguishable types based on whether they produce autoantibodies targeting the acetylcholine receptor (AChR-MG) or muscle specific kinase (MuSK-MG). The autoantibodies are pathogenic, but whether their generation is associated with broader defects in the B cell repertoire is unknown. To address this question, we performed deep sequencing of the BCR repertoire of AChR-MG, MuSK-MG, and healthy subjects to generate ∼518,000 unique V H and V L sequences from sorted naive and memory B cell populations. AChR-MG and MuSK-MG subjects displayed distinct gene segment usage biases in both V H and V L sequences within the naive and memory compartments. The memory compartment of AChR-MG was further characterized by reduced positive selection of somatic mutations in the V H CDR and altered V H CDR3 physicochemical properties. The V L repertoire of MuSK-MG was specifically characterized by reduced V-J segment distance in recombined sequences, suggesting diminished V L receptor editing during B cell development. Our results identify large-scale abnormalities in both the naive and memory B cell repertoires. Particular abnormalities were unique to either AChR-MG or MuSK-MG, indicating that the repertoires reflect the distinct properties of the subtypes. These repertoire abnormalities are consistent with previously observed defects in B cell tolerance checkpoints in MG, thereby offering additional insight regarding the impact of tolerance defects on peripheral autoimmune repertoires. These collective findings point toward a deformed B cell repertoire as a fundamental component of MG., (Copyright © 2017 by The American Association of Immunologists, Inc.)

Analyses of somatic hypermutation (SHM) patterns in B cell Ig sequences have important basic science and clinical applications, but they are often confounded by the intrinsic biases of SHM targeting on specific DNA motifs (i.e., hot and cold spots). Modeling these biases has been hindered by the difficulty in identifying mutated Ig sequences in vivo in the absence of selection pressures, which skew the observed mutation patterns. To generate a large number of unselected mutations, we immunized B1-8 H chain transgenic mice with nitrophenyl to stimulate nitrophenyl-specific λ + germinal center B cells and sequenced the unexpressed κ L chains using next-generation methods. Most of these κ sequences had out-of-frame junctions and were presumably uninfluenced by selection. Despite being nonfunctionally rearranged, they were targeted by SHM and displayed a higher mutation frequency than functional sequences. We used 39,173 mutations to construct a quantitative SHM targeting model. The model showed targeting biases that were consistent with classic hot and cold spots, yet revealed additional highly mutable motifs. We observed comparable targeting for functional and nonfunctional sequences, suggesting similar biological processes operate at both loci. However, we observed species- and chain-specific targeting patterns, demonstrating the need for multiple SHM targeting models. Interestingly, the targeting of C/G bases and the frequency of transition mutations at C/G bases was higher in mice compared with humans, suggesting lower levels of DNA repair activity in mice. Our models of SHM targeting provide insights into the SHM process and support future analyses of mutation patterns., (Copyright © 2016 by The American Association of Immunologists, Inc.)

Background: The genes that produce antibodies and the immune receptors expressed on lymphocytes are not germline encoded; rather, they are somatically generated in each developing lymphocyte by a process called V(D)J recombination, which assembles specific, independent gene segments into mature composite genes. The full set of composite genes in an individual at a single point in time is referred to as the immune repertoire. V(D)J recombination is the distinguishing feature of adaptive immunity and enables effective immune responses against an essentially infinite array of antigens. Characterization of immune repertoires is critical in both basic research and clinical contexts. Recent technological advances in repertoire profiling via high-throughput sequencing have resulted in an explosion of research activity in the field. This has been accompanied by a proliferation of software tools for analysis of repertoire sequencing data. Despite the widespread use of immune repertoire profiling and analysis software, there is currently no standardized format for output files from V(D)J analysis. Researchers utilize software such as IgBLAST and IMGT/High V-QUEST to perform V(D)J analysis and infer the structure of germline rearrangements. However, each of these software tools produces results in a different file format, and can annotate the same result using different labels. These differences make it challenging for users to perform additional downstream analyses., Results: To help address this problem, we propose a standardized file format for representing V(D)J analysis results. The proposed format, VDJML, provides a common standardized format for different V(D)J analysis applications to facilitate downstream processing of the results in an application-agnostic manner. The VDJML file format specification is accompanied by a support library, written in C++ and Python, for reading and writing the VDJML file format., Conclusions: The VDJML suite will allow users to streamline their V(D)J analysis and facilitate the sharing of scientific knowledge within the community. The VDJML suite and documentation are available from https://vdjserver.org/vdjml/ . We welcome participation from the community in developing the file format standard, as well as code contributions.

Objective: Myasthenia gravis (MG) is an autoimmune condition in which neurotransmission is impaired by binding of autoantibodies to acetylcholine receptors (AChR) or, in a minority of patients, to muscle specific kinase (MuSK). There are differences in the dominant IgG subclass, pathogenic mechanisms, and treatment responses between the two MG subtypes (AChR or MuSK). The antibodies are thought to be T-cell dependent, but the mechanisms underlying their production are not well understood. One aspect not previously described is whether defects in central and peripheral tolerance checkpoints, which allow autoreactive B cells to accumulate in the naive repertoire, are found in both or either form of MG., Methods: An established set of assays that measure the frequency of both polyreactive and autoreactive B cell receptors (BCR) in naive populations was applied to specimens collected from patients with either AChR or MuSK MG and healthy controls. Radioimmuno- and cell-based assays were used to measure BCR binding to AChR and MuSK., Results: The frequency of polyreactive and autoreactive BCRs (n = 262) was higher in both AChR and MuSK MG patients than in healthy controls. None of the MG-derived BCRs bound AChR or MuSK., Interpretation: The results indicate that both these MG subtypes harbor defects in central and peripheral B cell tolerance checkpoints. Defective B cell tolerance may represent a fundamental contributor to autoimmunity in MG and is of particular importance when considering the durability of myasthenia gravis treatment strategies, particularly biologics that eliminate B cells.

The adaptive immune system's 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.

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., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

During the several-week course of an immune response, B cells undergo a process of clonal expansion, somatic hypermutation of the immunoglobulin (Ig) genes and affinity-dependent selection. Over a lifetime, each B cell may participate in multiple rounds of affinity maturation as part of different immune responses. These two time-scales for selection are apparent in the structure of B-cell lineage trees, which often contain a 'trunk' consisting of mutations that are shared across all members of a clone, and several branches that form a 'canopy' consisting of mutations that are shared by a subset of clone members. The influence of affinity maturation on the B-cell population can be inferred by analysing the pattern of somatic mutations in the Ig. While global analysis of mutation patterns has shown evidence of strong selection pressures shaping the B-cell population, the effect of different time-scales of selection and diversification has not yet been studied. Analysis of B cells from blood samples of three healthy individuals identifies a range of clone sizes with lineage trees that can contain long trunks and canopies indicating the significant diversity introduced by the affinity maturation process. We here show that observed mutation patterns in the framework regions (FWRs) are determined by an almost purely purifying selection on both short and long time-scales. By contrast, complementarity determining regions (CDRs) are affected by a combination of purifying and antigen-driven positive selection on the short term, which leads to a net positive selection in the long term. In both the FWRs and CDRs, long-term selection is strongly dependent on the heavy chain variable gene family., (© 2015 The Author(s) Published by the Royal Society. All rights reserved.)

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., (Copyright © 2015 Elsevier Inc. All rights reserved.)

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., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2014, American Association for the Advancement of Science.)

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., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

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.

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.