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.)

Human molecular genetics has brought incredible insights into the variants that confer risk for the development of tissue-specific autoimmune diseases, including type 1 diabetes. The hallmark cell-mediated immune destruction that is characteristic of type 1 diabetes is closely linked with risk conferred by the HLA class II gene locus, in combination with a broad array of additional candidate genes influencing islet-resident beta cells within the pancreas, as well as function, phenotype and trafficking of immune cells to tissues. In addition to the well-studied germline SNP variants, there are critical contributions conferred by T cell receptor (TCR) and B cell receptor (BCR) genes that undergo somatic recombination to yield the Adaptive Immune Receptor Repertoire (AIRR) responsible for autoimmunity in type 1 diabetes. We therefore created the T1D TCR/BCR Repository (The Type 1 Diabetes T Cell Receptor and B Cell Receptor Repository) to study these highly variable and dynamic gene rearrangements. In addition to processed TCR and BCR sequences, the T1D TCR/BCR Repository includes detailed metadata (e.g. participant demographics, disease-associated parameters and tissue type). We introduce the Type 1 Diabetes AIRR Consortium goals and outline methods to use and deposit data to this comprehensive repository. Our ultimate goal is to facilitate research community access to rich, carefully annotated immune AIRR datasets to enable new scientific inquiry and insight into the natural history and pathogenesis of type 1 diabetes., (© 2024. The Author(s).)

In the past decade, single-cell (sc) transcriptomics has overcome the limitations of bulk analysis by measuring gene expression in individual cells, not just a population average. This can identify diverse cell types and states within a sample with high resolution, even without prior purification. Various technologies exist, each with its own capture, barcoding, and library preparation methods. This chapter focuses on the analysis of normal and malignant mature B cells using the 10× Genomics 5' sc-gene expression in parallel with B cell immune repertoire profiling. By integrating the gene expression data from similar cells, the complete transcriptome for each population can be reconstructed, while the identification of the expressed immunoglobulin genes allows investigating clonotype evolution and the detection of tumor clones that share the same clonally rearranged B cell receptor sequence. Researchers are guided through both the experimental protocols and data analysis with a comprehensive, step-by-step walkthrough of how to use some of the more popular single-cell software tools., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

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.

Tumor-infiltrating B-lineage cells have become predictors of prognosis and immunotherapy responses in various cancers. However, limited knowledge about their infiltration and migration patterns has hindered the understanding of their anti-tumor functions. Here, we examined the immunoglobulin heavy chain (IGH) repertoires in 496 multi-regional tumor, 107 normal tissue, and 48 metastatic lymph node samples obtained from 107 patients with esophageal squamous cell carcinoma (ESCC). Our study revealed higher IgG-type B-lineage cells infiltration in tumors than in healthy tissue, which was associated with improved patient outcomes. Genes such as ACTN1, COL6A5, and pathways like focal adhesion, which shapes the physical structure of tumors, could affect B-lineage cell infiltration. Notably, the IGH sequence was used as an identity-tag to monitor B cell migration, and their infiltration schema within the tumor were depicted based on our multi-regional tumor specimens. This analysis revealed an escalation in B cell clones overlapped between metastatic lymph nodes and tumors. Therefore, the Lymph Node Activation Index was defined, which could predict the outcomes of patients with lymph node metastasis. This research introduces a novel framework for probing B cell infiltration and migration within the tumor microenvironment using large-scale transcriptome data, while simultaneously providing fresh perspectives on B cell immunology within ESCC., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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.)

It is a great challenge to isolate the broadly neutralizing antibodies (bnAbs) against foot-and-mouth disease virus (FMDV) due to its existence as seven distinct serotypes without cross-protection. Here, by vaccination of pig with FMDV serotype O and A whole virus antigens, we obtained 10 bnAbs against serotypes O, A and/or Asia1 by dissecting 216 common clonotypes of two serotype O and A specific porcine B-cell receptor (BCR) gene repertoires containing total 12720 B cell clones, indicating the induction of cross-serotype bnAbs after sequential vaccination with serotypes O and A antigens. The majority of porcine bnAbs (9/10) were derived from terminally differentiated B cells of different clonal lineages, which convergently targeted the conserved "RGDL" motif on structural protein VP1 of FMDV by mimicking receptor recognition to inhibit viral attachment to cells. Cryo-EM complex structures revealed that the other bnAb pOA-2 specifically targets a novel inter-pentamer antigen structure surrounding the viral three-fold axis, with a highly conserved determinant at residue 68 on VP2. This unique binding pattern enabled cross-serotype neutralization by destabilizing the viral particle. The evolutionary analysis of pOA-2 demonstrated its origin from an intermediate B-cell, emphasizing the crucial role of somatic hypermutations (SHMs) in balancing the breadth and potency of neutralization. However, excessive SHMs may deviate from the trajectory of broad neutralization. This study provides a strategy to uncover bnAbs against highly mutable pathogens and the cross-serotype antigenic structures to explore broadly protective FMDV vaccine. Author summary: The scarcity of broadly neutralizing antibodies (bnAbs) in the body significantly hinders their discovery and elucidation of the conserved antigenic structure of highly mutable pathogens. Foot-and-mouth disease virus (FMDV) severely infects pigs and exists as seven distinct serotypes without cross-protection. In this study, we present a strategy for rapid discovery of cross-serotype bnAbs by constructing of specific B-cell repertoires and exploring of common clonotypes. We originally isolated porcine bnAbs against FMDV serotypes O, A and/or Asia1, and proved for the first time that receptor binding region in VP1 can induce the production of cross-serotype bnAbs. The evolution of cross-serotype bnAbs reveals how somatic hypermutations balance the breadth and potency of neutralization. FMDV bnAbs employ dual neutralization mechanisms to neutralize viruses via receptor mimicry and destabilizing the viral particles. This study provides valuable guidelines that facilitate the isolation of bnAbs against the highly variable viruses and aid in designing a universal FMDV vaccine. [ABSTRACT FROM AUTHOR]

Single-cell RNA sequencing (scRNA-seq) can resolve transcriptional features from individual cells, but scRNA-seq techniques capable of resolving the variable regions of B cell receptors (BCRs) remain limited, especially from widely-used 3′-barcoded libraries. Here, we report a method that can recover paired, full-length variable region sequences of BCRs from 3′-barcoded scRNA-seq libraries. We first verify this method (B3E-seq) can produce accurate, full-length BCR sequences. We then apply this method to profile B cell responses elicited against the capsular polysaccharide of Streptococcus pneumoniae serotype 3 (ST3) by glycoconjugate vaccines in five infant rhesus macaques. We identify BCR features associated with specificity for the ST3 antigen which are present in multiple vaccinated monkeys, indicating a convergent response to vaccination. These results demonstrate the utility of our method to resolve key features of the B cell repertoire and profile antigen-specific responses elicited by vaccination. A method that recovers full-length, paired heavy- and light-chain variable regions of B cell receptor transcripts from 3'barcoded scRNA-seq libraries reveals a convergent response to pneumococcus vaccination in rhesus macaques. [ABSTRACT FROM AUTHOR]

Bats are reservoirs of many zoonotic viruses that are fatal in humans but do not cause disease in bats. Moreover, bats generate low neutralizing antibody titers in response to experimental viral infection, although more robust antibody responses have been observed in wild-caught bats during times of food stress. Here, we compared the antibody titers and B cell receptor (BCR) diversity of Jamaican fruit bats (Artibeus jamaicensis; JFBs) and BALB/c mice generated in response to T-dependent and T-independent antigens. We then manipulated the diet of JFBs and challenged them with H18N11 influenza A-like virus or a replication incompetent Nipah virus VSV (Nipah-riVSV). Under standard housing conditions, JFBs generated a lower avidity antibody response and possessed more BCR mRNA diversity compared to BALB/c mice. However, withholding protein from JFBs improved serum neutralization in response to Nipah-riVSV and improved serum antibody titers specific to H18 but reduced BCR mRNA diversity. Bats are known pathogen reservoirs and changes in their diet have been linked to spillover events. This study shows that bat antibodies have higher diversity but lower pathogen-binding capacity than those of mice, but dietary protein restriction enhances this capacity, potentially impacting spillover risk. [ABSTRACT FROM AUTHOR]

In recent years, rapid advances in research methods have made single cell analysis possible. Systemic sclerosis (SSc), a disease characterized by the triad of immune abnormalities, fibrosis, and vasculopathy, has also been the subject of various analyses. To summarize the results of single cell analysis in SSc accumulated to date and to deepen our understanding of SSc. Four databases were used to perform a database search on 23rd June 2023. Assessed Grading of Recommendations Assessment, Development and Evaluation certainty of evidence were performed according to PRISMA guidelines. The analysis was completed on July 2023. 17 studies with 358 SSc patients were included. Three studies used PBMCs, six used skin, nine used lung with SSc-interstitial lung diseases (ILDs), and one used lung with SSc-pulmonary arterial hypertension (PAH). The cells studied included immune cells such as T cells, natural killer cells, monocytes, macrophages, and dendritic cells, as well as endothelial cells, fibroblasts, keratinocytes, alveolar type I cells, basal epithelial cells, smooth muscle cells, mesothelial cells, etc. This systematic review revealed the results of single cell analysis, suggesting that PBMCs, skin, SSc-ILD, and SSc-PAH show activation and dysfunction of cells associated with immune-abnormalities, fibrosis, and vasculopathy, respectively. [ABSTRACT FROM AUTHOR]

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.

Background: An association has been observed between primary biliary cholangitis (PBC) and systemic rheumatic diseases (SRDs) in observational studies, however the exact causal link remains unclear. We aimed to evaluate the causal effects of PBC on SRDs through Mendelian randomization (MR) analysis. Methods: The genome-wide association study (GWAS) summary data were obtained from MRC IEU OpenGWAS and FinnGen databases. Independent genetic variants for PBC were selected as instrumental variables. Inverse variance weighted was used as the main approach to evaluate the causal effects of PBC on Sjögren syndrome (SS), rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), systemic sclerosis (SSc), mixed connective tissue disease (MCTD) and polymyositis (PM). Horizontal pleiotropy and heterogeneity were measured by MR‒Egger intercept test and Cochran's Q value, respectively. Results: PBC had causal effects on SS (OR = 1.177, P = 8.02e-09), RA (OR = 1.071, P = 9.80e-04), SLE (OR = 1.447, P = 1.04e-09), SSc (OR = 1.399, P = 2.52e-04), MCTD (OR = 1.306, P = 4.92e-14), and PM (OR = 1.416, P = 1.16e-04). Based on the MR‒Egger intercept tests, horizontal pleiotropy was absent (all P values > 0.05). The robustness of our results was further enhanced by the leave-one-out method. Conclusions: Our research has provided new insights into PBC and SRDs, indicating casual effects on various SRDs. [ABSTRACT FROM AUTHOR]

The leucine-rich repeat-based variable lymphocyte receptor B (VLRB) antibody system of jawless vertebrates is capable of generating an antibody repertoire equal to or exceeding the diversity of antibody repertoires of jawed vertebrates. Unlike immunoglobulin-based immune repertoires, the VLRB repertoire diversity is characterized by variable lengths of VLRB encoding transcripts, rendering conventional immunoreceptor repertoire sequencing approaches unsuitable for VLRB repertoire sequencing. Here we demonstrate that long-read single-molecule real-time (SMRT) sequencing (PacBio) approaches permit the efficient large-scale assessment of the VLRB repertoire. We present a computational pipeline for sequence data processing and provide the first repertoire-based analysis of VLRB protein characteristics including properties of its subunits and regions of diversity within each structural leucine-rich repeat subunit. Our study provides a template to explore changes in the VLRB repertoire during immune responses and to establish large scale VLRB repertoire databases for computational approaches aimed at isolating monoclonal VLRB reagents for biomedical research and clinical applications. [ABSTRACT FROM AUTHOR]

The upper airway is an important site of infection, but immune memory in the human upper airway is poorly understood, with implications for COVID-19 and many other human diseases1–4. Here we demonstrate that nasal and nasopharyngeal swabs can be used to obtain insights into these challenging problems, and define distinct immune cell populations, including antigen-specific memory B cells and T cells, in two adjacent anatomical sites in the upper airway. Upper airway immune cell populations seemed stable over time in healthy adults undergoing monthly swabs for more than 1 year, and prominent tissue resident memory T (TRM) cell and B (BRM) cell populations were defined. Unexpectedly, germinal centre cells were identified consistently in many nasopharyngeal swabs. In subjects with SARS-CoV-2 breakthrough infections, local virus-specific BRM cells, plasma cells and germinal centre B cells were identified, with evidence of local priming and an enrichment of IgA+ memory B cells in upper airway compartments compared with blood. Local plasma cell populations were identified with transcriptional profiles of longevity. Local virus-specific memory CD4+ TRM cells and CD8+ TRM cells were identified, with diverse additional virus-specific T cells. Age-dependent upper airway immunological shifts were observed. These findings provide new understanding of immune memory at a principal mucosal barrier tissue in humans.This study of immunological memory diversity in the human upper airway provides new understanding of immune memory at a major mucosal barrier tissue in humans. [ABSTRACT FROM AUTHOR]

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method. [ABSTRACT FROM AUTHOR]

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a valuable experimental tool to study the immune state in health and following immune challenges such as infectious diseases, (auto)immune diseases, and cancer. Several tools have been developed to reconstruct B cell and T cell receptor sequences from AIRR-seq data and infer B and T cell clonal relationships. However, currently available tools offer limited parallelization across samples, scalability or portability to high-performance computing infrastructures. To address this need, we developed nf-core/airrflow, an end-to-end bulk and single-cell AIRR-seq processing workflow which integrates the Immcantation Framework following BCR and TCR sequencing data analysis best practices. The Immcantation Framework is a comprehensive toolset, which allows the processing of bulk and single-cell AIRR-seq data from raw read processing to clonal inference. nf-core/airrflow is written in Nextflow and is part of the nf-core project, which collects community contributed and curated Nextflow workflows for a wide variety of analysis tasks. We assessed the performance of nf-core/airrflow on simulated sequencing data with sequencing errors and show example results with real datasets. To demonstrate the applicability of nf-core/airrflow to the high-throughput processing of large AIRR-seq datasets, we validated and extended previously reported findings of convergent antibody responses to SARS-CoV-2 by analyzing 97 COVID-19 infected individuals and 99 healthy controls, including a mixture of bulk and single-cell sequencing datasets. Using this dataset, we extended the convergence findings to 20 additional subjects, highlighting the applicability of nf-core/airrflow to validate findings in small in-house cohorts with reanalysis of large publicly available AIRR datasets. Author summary: We have created nf-core/airrflow, a workflow to help researchers study the immune system in healthy and disease states, such as infections, autoimmunity, and cancer. The adaptive immune system is responsible for the third line of defense responses, specific to each particular threat, after physical barriers have been compromised and the nonspecific innate immune response has failed to clear the danger. Two types of white blood cells are central players in the adaptive response, namely B cells and T cells. These cells have surface receptors that recognize suspicious elements (antigens). Learning what receptors bind to which antigens is of utmost interest to understand immune responses. Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is a technique that allows to determine the genetic sequence of these receptors. The amount of data generated in these experiments is large, and the analysis complex. nf-core/airrflow simplifies running a comprehensive analysis connecting tools from Immcantation, a specialized software project to analyze AIRR-seq data. The workflow can efficiently process large datasets on multiple computing platforms. We analyzed immune responses to COVID-19 in 97 infected individuals and 99 healthy people, and confirmed previous findings and provided new insights, demonstrating the workflow's applicability to reanalyzing large publicly available datasets. [ABSTRACT FROM AUTHOR]

Outbreaks of Ebolaviruses, such as Sudanvirus (SUDV) in Uganda in 2022, demonstrate that species other than the Zaire ebolavirus (EBOV), which is currently the sole virus represented in current licensed vaccines, remain a major threat to global health. There is a pressing need to develop effective pan-species vaccines and novel monoclonal antibody-based therapeutics for Ebolavirus disease. In response to recent outbreaks, the two dose, heterologous Ad26.ZEBOV/MVA-BN-Filo vaccine regimen was developed and was tested in a large phase II clinical trial (EBL2001) as part of the EBOVAC2 consortium. Here, we perform bulk sequencing of the variable heavy chain (VH) of B cell receptors (BCR) in forty participants from the EBL2001 trial in order to characterize the BCR repertoire in response to vaccination with Ad26.ZEBOV/MVA-BN-Filo. We develop a comprehensive database, EBOV-AbDab, of publicly available Ebolavirus-specific antibody sequences. We then use our database to predict the antigen-specific component of the vaccinee repertoires. Our results show striking convergence in VH germline gene usage across participants following the MVA-BN-Filo dose, and provide further evidence of the role of IGHV3--15 and IGHV3--13 antibodies in the B cell response to Ebolavirus glycoprotein. Furthermore, we found that previously described Ebola-specific mAb sequences present in EBOV-AbDab were sufficient to describe at least one of the ten most expanded BCR clonotypes in more than two thirds of our cohort of vaccinees following the boost, providing proof of principle for the utility of computational mining of immune repertoires. [ABSTRACT FROM AUTHOR]

Introduction: It is unknown how intestinal B cell populations and B cell receptor (BCR) repertoires are established and maintained over time in humans. Following intestinal transplantation (ITx), surveillance ileal mucosal biopsies provide a unique opportunity to map the dynamic establishment of recipient gut lymphocyte populations in immunosuppressed conditions. Methods: Using polychromatic flow cytometry that includes HLA allele groupspecific antibodies distinguishing donor from recipient cells along with high throughput BCR sequencing, we tracked the establishment of recipient B cell populations and BCR repertoire in the allograft mucosa of ITx recipients. Results: We confirm the early presence of naïve donor B cells in the circulation (donor age range: 1-14 years, median: 3 years) and, for the first time, document the establishment of recipient B cell populations, including B resident memory cells, in the intestinal allograft mucosa (recipient age range at the time of transplant: 1-44 years, median: 3 years). Recipient B cell repopulation of the allograft was most rapid in infant (<1 year old)-derived allografts and, unlike T cell repopulation, did not correlate with rejection rates. While recipient memory B cell populations were increased in graft mucosa compared to circulation, naïve recipient B cells remained detectable in the graft mucosa for years. Comparisons of peripheral and intra-mucosal B cell repertoires in the absence of rejection (recipient age range at the time of transplant: 1-9 years, median: 2 years) revealed increased BCR mutation rates and clonal expansion in graft mucosa compared to circulating B cells, but these parameters did not increase markedly after the first year post-transplant. Furthermore, clonal mixing between the allograft mucosa and the circulation was significantly greater in ITx recipients, even years after transplantation, than in deceased adult donors. In available pan-scope biopsies from pediatric recipients, we observed higher percentages of naïve recipient B cells in colon allograft compared to small bowel allograft and increased BCR overlap between native colon vs colon allograft compared to that between native colon vs ileum allograft in most cases, suggesting differential clonal distribution in large intestine vs small intestine. Discussion: Collectively, our data demonstrate intestinal mucosal B cell repertoire establishment from a circulating pool, a process that continues for years without evidence of stabilization of the mucosal B cell repertoire in pediatric ITx patients. [ABSTRACT FROM AUTHOR]

This manuscript describes the development of a resource module that is part of a learning platform named "NIGMS Sandbox for Cloud-based Learning" https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox at the beginning of this Supplement. This module delivers learning materials on RNA sequencing (RNAseq) data analysis in an interactive format that uses appropriate cloud resources for data access and analyses. Biomedical research is increasingly data-driven, and dependent upon data management and analysis methods that facilitate rigorous, robust, and reproducible research. Cloud-based computing resources provide opportunities to broaden the application of bioinformatics and data science in research. Two obstacles for researchers, particularly those at small institutions, are: (i) access to bioinformatics analysis environments tailored to their research; and (ii) training in how to use Cloud-based computing resources. We developed five reusable tutorials for bulk RNAseq data analysis to address these obstacles. Using Jupyter notebooks run on the Google Cloud Platform, the tutorials guide the user through a workflow featuring an RNAseq dataset from a study of prophage altered drug resistance in Mycobacterium chelonae. The first tutorial uses a subset of the data so users can learn analysis steps rapidly, and the second uses the entire dataset. Next, a tutorial demonstrates how to analyze the read count data to generate lists of differentially expressed genes using R/DESeq2. Additional tutorials generate read counts using the Snakemake workflow manager and Nextflow with Google Batch. All tutorials are open-source and can be used as templates for other analysis. [ABSTRACT FROM AUTHOR]