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

B cell receptor repertoire abnormalities in autoimmune disease

B cells play a crucial role in the immune response and contribute to various autoimmune diseases. Recent studies have revealed abnormalities in the B cell receptor (BCR) repertoire of patients with autoimmune diseases, with distinct features observed among different diseases and B cell subsets. Classically, BCR repertoire was used as an identifier of distinct antigen-specific clonotypes, but the recent advancement of analyzing large-scale repertoire has enabled us to use it as a tool for characterizing cellular biology. In this review, we provide an overview of the BCR repertoire in autoimmune diseases incorporating insights from our latest research findings. In systemic lupus erythematosus (SLE), we observed a significant skew in the usage of VDJ genes, particularly in CD27+IgD+ unswitched memory B cells and plasmablasts. Notably, autoreactive clones within unswitched memory B cells were found to be increased and strongly associated with disease activity, underscoring the clinical significance of this subset. Similarly, various abnormalities in the BCR repertoire have been reported in other autoimmune diseases such as rheumatoid arthritis. Thus, BCR repertoire analysis holds potential for enhancing our understanding of the underlying mechanisms involved in autoimmune diseases. Moreover, it has the potential to predict treatment effects and identify therapeutic targets in autoimmune diseases.

nf-core/airrflow: an adaptive immune receptor repertoire analysis workflow employing the Immcantation framework

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. nf-core/airrflow is available free of charge, under the MIT license on GitHub (https://github.com/nf-core/airrflow). Detailed documentation and example results are available on the nf-core website at (https://nf-co.re/airrflow).

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

Lupus-associated innate receptors drive extrafollicular evolution of autoreactive B cells

In systemic lupus erythematosus, recent findings highlight the extrafollicular (EF) pathway as prominent origin of autoantibody-secreting cells (ASCs). CD21loCD11c+ B cells, associated with aging, infection, and autoimmunity, are contributors to autoreactive EF ASCs but have an obscure developmental trajectory. To study EF kinetics of autoreactive B cell in tissue, we adoptively transferred WT and gene knockout B cell populations into the 564Igi mice - an autoreactive host enriched with autoantigens and T cell help. Time-stamped analyses revealed TLR7 dependence in early escape of peripheral B cell tolerance and establishment of a pre-ASC division program. We propose CD21lo cells as precursors to EF ASCs due to their elevated TLR7 sensitivity and proliferative nature. Blocking receptor function reversed CD21 loss and reduced effector cell generation, portraying CD21 as a differentiation initiator and a possible target for autoreactive B cell suppression. Repertoire analysis further delineated proto-autoreactive B cell selection and receptor evolution toward self-reactivity. This work elucidates receptor and clonal dynamics in EF development of autoreactive B cells, and establishes modular, native systems to probe mechanisms of autoreactivity.

B cell phylogenetics in the single cell era

The widespread availability of single-cell RNA sequencing (scRNA-seq) has led to the development of new methods for understanding immune responses. Single-cell transcriptome data can now be paired with B cell receptor (BCR) sequences. However, RNA from BCRs cannot be analyzed like most other genes because BCRs are genetically diverse within individuals. In humans, BCRs are shaped through recombination followed by mutation and selection for antigen binding. As these processes co-occur with cell division, B cells can be studied using phylogenetic trees representing the mutations within a clone. B cell trees can link experimental timepoints, tissues, or cellular subtypes. Here, we review the current state and potential of how B cell phylogenetics can be combined with single-cell data to understand immune responses.

Systems immunology spanning tumors, lymph nodes, and periphery

The immune system defines a complex network of tissues and cell types that orchestrate responses across the body in a dynamic manner. The local and systemic interactions between immune and cancer cells contribute to disease progression. Lymphocytes are activated in lymph nodes, traffic through the periphery, and impact cancer progression through their interactions with tumor cells. As a result, therapeutic response and resistance are mediated across tissues, and a comprehensive understanding of lymphocyte dynamics requires a systems-level approach. In this review, we highlight experimental and computational methods that can leverage the study of leukocyte trafficking through an immunomics lens and reveal how adaptive immunity shapes cancer.

Spatial transcriptomics of B cell and T cell receptors reveals lymphocyte clonal dynamics

The spatial distribution of lymphocyte clones within tissues is critical to their development, selection, and expansion. We have developed spatial transcriptomics of variable, diversity, and joining (VDJ) sequences (Spatial VDJ), a method that maps B cell and T cell receptor sequences in human tissue sections. Spatial VDJ captures lymphocyte clones that match canonical B and T cell distributions and amplifies clonal sequences confirmed by orthogonal methods. We found spatial congruency between paired receptor chains, developed a computational framework to predict receptor pairs, and linked the expansion of distinct B cell clones to different tumor-associated gene expression programs. Spatial VDJ delineates B cell clonal diversity and lineage trajectories within their anatomical niche. Thus, Spatial VDJ captures lymphocyte spatial clonal architecture across tissues, providing a platform to harness clonal sequences for therapy.

Sustained Liver HBsAg Loss and Clonal T- and B-Cell Expansion upon Therapeutic DNA Vaccination Require Low HBsAg Levels

BACKGROUND

Suppression of HBV DNA, inhibition of HBV surface (HBsAg) production and therapeutic vaccination to reverse HBV-specific T-cell exhaustion in chronic HBV patients are likely required to achieve a functional cure. In the AAV-HBV mouse model, therapeutic vaccination can be effective in clearing HBV when HBsAg levels are low. Using a single-cell approach, we investigated the liver immune environment with different levels of HBsAg and sustained HBsAg loss through treatment with a GalNAc-HBV-siRNA followed by therapeutic vaccination.

METHODS

AAV-HBV-transduced C57BL/6 mice were treated with GalNAc-HBV-siRNA to lower HBsAg levels and then vaccinated using a DNA vaccine. We used single-cell RNA and V(D)J sequencing to understand liver immune microenvironment changes.

RESULTS

GalNAc-HBV-siRNA, followed by therapeutic vaccination, achieved sustained HBsAg loss in all mice. This was accompanied by CD4 follicular helper T-cell induction, polyclonal activation of CD8 T cells and clonal expansion of plasma cells that were responsible for antibody production.

CONCLUSIONS

This study provides novel insights into liver immune changes at the single-cell level, highlighting the correlation between induced reduction of HBsAg levels and clonal expansion of CD4, CD8 T cells and plasma cells in the liver upon HBV siRNA and subsequent therapeutic vaccination.

Enhanced SARS-CoV-2 humoral immunity following breakthrough infection builds upon the preexisting memory B cell pool

The human immune response must continuously adapt to newly emerging SARS-CoV-2 variants. To investigate how B cells respond to repeated SARS-CoV-2 antigen exposure by Wu01 booster vaccination and Omicron breakthrough infection, we performed a molecular longitudinal analysis of the memory B cell pool. We demonstrate that a subsequent breakthrough infection substantially increases the frequency of B cells encoding SARS-CoV-2-neutralizing antibodies. However, this is not primarily attributable to maturation, but to selection of preexisting B cell clones. Moreover, broadly reactive memory B cells arose early and even neutralized highly mutated variants like XBB.1.5 that the individuals had not encountered. Together, our data show that SARS-CoV-2 immunity is largely imprinted on Wu01 over the course of multiple antigen contacts but can respond to new variants through preexisting diversity.

Ultrasensitive allele inference from immune repertoire sequencing data with MiXCR

Allelic variability in the adaptive immune receptor loci, which harbor the gene segments that encode B cell and T cell receptors (BCR/TCR), has been shown to be of critical importance for immune responses to pathogens and vaccines. In recent years, B cell and T cell receptor repertoire sequencing (Rep-Seq) has become widespread in immunology research making it the most readily available source of information about allelic diversity in immunoglobulin (IG) and T cell receptor (TR) loci in different populations. Here we present a novel algorithm for extra-sensitive and specific variable (V) and joining (J) gene allele inference and genotyping allowing reconstruction of individual high-quality gene segment libraries. The approach can be applied for inferring allelic variants from peripheral blood lymphocyte BCR and TCR repertoire sequencing data, including hypermutated isotype-switched BCR sequences, thus allowing high-throughput genotyping and novel allele discovery from a wide variety of existing datasets. The developed algorithm is a part of the MiXCR software ( https://mixcr.com ) and can be incorporated into any pipeline utilizing upstream processing with MiXCR. We demonstrate the accuracy of this approach using Rep-Seq paired with long-read genomic sequencing data, comparing it to a widely used algorithm, TIgGER. We applied the algorithm to a large set of IG heavy chain (IGH) Rep-Seq data from 450 donors of ancestrally diverse population groups, and to the largest reported full-length TCR alpha and beta chain (TRA; TRB) Rep-Seq dataset, representing 134 individuals. This allowed us to assess the genetic diversity of genes within the IGH, TRA and TRB loci in different populations and demonstrate the connection between antibody repertoire gene usage and the number of allelic variants present in the population. Finally we established a database of allelic variants of V and J genes inferred from Rep-Seq data and their population frequencies with free public access at https://vdj.online .

Dynamic establishment and maintenance of the human intestinal B cell population and repertoire following transplantation

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 gut lymphocyte populations. Using polychromatic flow cytometry that includes HLA allele group-specific mAbs 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. We confirm the early presence of naïve donor B cells in the circulation and, for the first time, document the establishment of recipient B cell populations, including B resident memory cells, in the intestinal allograft mucosa. 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 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 healthy control adults. Collectively, our data demonstrate intestinal mucosal B cell repertoire establishment from a circulating pool, a process that continues for years without evidence of establishment of a stable mucosal B cell repertoire.

Single-cell analysis of memory B cells from top neutralizers reveals multiple sites of vulnerability within HCMV Trimer and Pentamer

Human cytomegalovirus (HCMV) can cause severe diseases in fetuses, newborns, and immunocompromised individuals. Currently, no vaccines are approved, and treatment options are limited. Here, we analyzed the human B cell response of four HCMV top neutralizers from a cohort of 9,000 individuals. By single-cell analyses of memory B cells targeting the pentameric and trimeric HCMV surface complexes, we identified vulnerable sites on the shared gH/gL subunits as well as complex-specific subunits UL128/130/131A and gO. Using high-resolution cryogenic electron microscopy, we revealed the structural basis of the neutralization mechanisms of antibodies targeting various binding sites. Moreover, we identified highly potent antibodies that neutralized a broad spectrum of HCMV strains, including primary clinical isolates, that outperform known antibodies used in clinical trials. Our study provides a deep understanding of the mechanisms of HCMV neutralization and identifies promising antibody candidates to prevent and treat HCMV infection.

Probabilities of developing HIV-1 bNAb sequence features in uninfected and chronically infected individuals

HIV-1 broadly neutralizing antibodies (bNAbs) are able to suppress viremia and prevent infection. Their induction by vaccination is therefore a major goal. However, in contrast to antibodies that neutralize other pathogens, HIV-1-specific bNAbs frequently carry uncommon molecular characteristics that might prevent their induction. Here, we perform unbiased sequence analyses of B cell receptor repertoires from 57 uninfected and 46 chronically HIV-1- or HCV-infected individuals and learn probabilistic models to predict the likelihood of bNAb development. We formally show that lower probabilities for bNAbs are predictive of higher HIV-1 neutralization activity. Moreover, ranking bNAbs by their probabilities allows to identify highly potent antibodies with superior generation probabilities as preferential targets for vaccination approaches. Importantly, we find equal bNAb probabilities across infected and uninfected individuals. This implies that chronic infection is not a prerequisite for the generation of bNAbs, fostering the hope that HIV-1 vaccines can induce bNAb development in uninfected people.

IL-4 downregulates BCL6 to promote memory B cell selection in germinal centers

Germinal center (GC)-derived memory B cells (MBCs) are critical for humoral immunity as they differentiate into protective antibody-secreting cells during re-infection. GC formation and cellular interactions within the GC have been studied in detail, yet the exact signals that allow for the selection and exit of MBCs are not understood. Here, we show that IL-4 signaling in GC B cells directly downregulates BCL6 via negative autoregulation to release cells from the GC program and promote MBC formation. This selection event requires additional survival cues and can therefore result in either GC exit or death. We demonstrate that both increasing IL-4 bioavailability or limiting IL-4 signaling disrupt MBC selection stringency. In this way, IL-4 control of BCL6 expression serves as a tunable switch within the GC to tightly regulate MBC selection and affinity maturation.

Recognition determinants of improved HIV-1 neutralization by a heavy chain matured pediatric antibody

The structural and characteristic features of HIV-1 broadly neutralizing antibodies (bnAbs) from chronically infected pediatric donors are currently unknown. Herein, we characterized a heavy chain matured HIV-1 bnAb 44m, identified from a pediatric elite-neutralizer. Interestingly, in comparison to its wild-type AIIMS-P01 bnAb, 44m exhibited moderately higher level of somatic hypermutations of 15.2%. The 44m neutralized 79% of HIV-1 heterologous viruses (n = 58) tested, with a geometric mean IC50 titer of 0.36 μg/mL. The cryo-EM structure of 44m Fab in complex with fully cleaved glycosylated native-like BG505.SOSIP.664.T332N gp140 envelope trimer at 4.4 Å resolution revealed that 44m targets the V3-glycan N332-supersite and GDIR motif to neutralize HIV-1 with improved potency and breadth, plausibly attributed by a matured heavy chain as compared to that of wild-type AIIMS-P01. This study further improves our understanding on pediatric HIV-1 bnAbs and structural basis of broad HIV-1 neutralization by 44m may be useful blueprint for vaccine design in future.

R2C2+UMI: Combining concatemeric consensus sequencing with unique molecular identifiers enables ultra-accurate sequencing of amplicons on Oxford Nanopore Technologies sequencers

The sequencing of PCR amplicons is a core application of high-throughput sequencing technology. Using unique molecular identifiers (UMIs), individual amplified molecules can be sequenced to very high accuracy on an Illumina sequencer. However, Illumina sequencers have limited read length and are therefore restricted to sequencing amplicons shorter than 600bp unless using inefficient synthetic long-read approaches. Native long-read sequencers from Pacific Biosciences and Oxford Nanopore Technologies can, using consensus read approaches, match or exceed Illumina quality while achieving much longer read lengths. Using a circularization-based concatemeric consensus sequencing approach (R2C2) paired with UMIs (R2C2+UMI) we show that we can sequence ~550nt antibody heavy-chain (IGH) and ~1500nt 16S amplicons at accuracies up to and exceeding Q50 (<1 error in 100,0000 sequenced bases), which exceeds accuracies of UMI-supported Illumina paired sequencing as well as synthetic long-read approaches.

Rational design of a booster vaccine against COVID-19 based on antigenic distance

Current COVID-19 vaccines are highly effective against symptomatic disease, but repeated booster doses using vaccines based on the ancestral strain offer limited additional protection against SARS-CoV-2 variants of concern (VOCs). To address this, we used antigenic distance to in silico select optimized booster vaccine seed strains effective against both current and future VOCs. Our model suggests that a SARS-CoV-1-based booster vaccine has the potential to cover a broader range of VOCs. Candidate vaccines including the spike protein from ancestral SARS-CoV-2, Delta, Omicron (BA.1), SARS-CoV-1, or MERS-CoV were experimentally evaluated in mice following two doses of the BNT162b2 vaccine. The SARS-CoV-1-based booster vaccine outperformed other candidates in terms of neutralizing antibody breadth and duration, as well as protective activity against Omicron (BA.2) challenge. This study suggests a unique strategy for selecting booster vaccines based on antigenic distance, which may be useful in designing future booster vaccines as new SARS-CoV-2 variants emerge.

SARS-CoV-2 infection and recovery in children: Distinct T cell responses in MIS-C compared to COVID-19

SARS-CoV-2 infection for most children results in mild or minimal symptoms, though in rare cases severe disease can develop, including a multisystem inflammatory syndrome (MIS-C) with myocarditis. Here, we present longitudinal profiling of immune responses during acute disease and following recovery in children who developed MIS-C, relative to children who experienced more typical symptoms of COVID-19. T cells in acute MIS-C exhibited transient signatures of activation, inflammation, and tissue residency which correlated with cardiac disease severity, while T cells in acute COVID-19 upregulated markers of follicular helper T cells for promoting antibody production. The resultant memory immune response in recovery showed increased frequencies of virus-specific memory T cells with pro-inflammatory functions in children with prior MIS-C compared to COVID-19 while both cohorts generated comparable antibody responses. Together our results reveal distinct effector and memory T cell responses in pediatric SARS-CoV-2 infection delineated by clinical syndrome, and a potential role for tissue-derived T cells in the immune pathology of systemic disease.

Germline-encoded specificities and the predictability of the B cell response

Antibodies result from the competition of B cell lineages evolving under selection for improved antigen recognition, a process known as affinity maturation. High-affinity antibodies to pathogens such as HIV, influenza, and SARS-CoV-2 are frequently reported to arise from B cells whose receptors, the precursors to antibodies, are encoded by particular immunoglobulin alleles. This raises the possibility that the presence of particular germline alleles in the B cell repertoire is a major determinant of the quality of the antibody response. Alternatively, initial differences in germline alleles' propensities to form high-affinity receptors might be overcome by chance events during affinity maturation. We first investigate these scenarios in simulations: when germline-encoded fitness differences are large relative to the rate and effect size variation of somatic mutations, the same germline alleles persistently dominate the response of different individuals. In contrast, if germline-encoded advantages can be easily overcome by subsequent mutations, allele usage becomes increasingly divergent over time, a pattern we then observe in mice experimentally infected with influenza virus. We investigated whether affinity maturation might nonetheless strongly select for particular amino acid motifs across diverse genetic backgrounds, but we found no evidence of convergence to similar CDR3 sequences or amino acid substitutions. These results suggest that although germline-encoded specificities can lead to similar immune responses between individuals, diverse evolutionary routes to high affinity limit the genetic predictability of responses to infection and vaccination.

Induction of bronchus-associated lymphoid tissue is an early life adaptation for promoting human B cell immunity

Infants and young children are more susceptible to common respiratory pathogens than adults but can fare better against novel pathogens like severe acute respiratory syndrome coronavirus 2. The mechanisms by which infants and young children mount effective immune responses to respiratory pathogens are unknown. Through investigation of lungs and lung-associated lymph nodes from infant and pediatric organ donors aged 0-13 years, we show that bronchus-associated lymphoid tissue (BALT), containing B cell follicles, CD4+ T cells and functionally active germinal centers, develop during infancy. BALT structures are prevalent around lung airways during the first 3 years of life, and their numbers decline through childhood coincident with the accumulation of memory T cells. Single-cell profiling and repertoire analysis reveals that early life lung B cells undergo differentiation, somatic hypermutation and immunoglobulin class switching and exhibit a more activated profile than lymph node B cells. Moreover, B cells in the lung and lung-associated lymph nodes generate biased antibody responses to multiple respiratory pathogens compared to circulating antibodies, which are mostly specific for vaccine antigens in the early years of life. Together, our findings provide evidence for BALT as an early life adaptation for mobilizing localized immune protection to the diverse respiratory challenges during this formative life stage.

T cell deletional tolerance restricts AQP4 but not MOG CNS autoimmunity

Aquaporin-4 (AQP4)-specific Th17 cells are thought to have a central role in neuromyelitis optica (NMO) pathogenesis. When modeling NMO, only AQP4-reactive Th17 cells from AQP4-deficient (AQP4-/-), but not wild-type (WT) mice, caused CNS autoimmunity in recipient WT mice, indicating that a tightly regulated mechanism normally ensures tolerance to AQP4. Here, we found that pathogenic AQP4 T cell epitopes bind MHC II with exceptionally high affinity. Examination of T cell receptor (TCR) α/β usage revealed that AQP4-specific T cells from AQP4-/- mice employed a distinct TCR repertoire and exhibited clonal expansion. Selective thymic AQP4 deficiency did not fully restore AQP4-reactive T cells, demonstrating that thymic negative selection alone did not account for AQP4-specific tolerance in WT mice. Indeed, AQP4-specific Th17 cells caused paralysis in recipient WT or B cell-deficient mice, which was followed by complete recovery that was associated with apoptosis of donor T cells. However, donor AQP4-reactive T cells survived and caused persistent paralysis in recipient mice deficient in both T and B cells or mice lacking T cells only. Thus, AQP4 CNS autoimmunity was limited by T cell-dependent deletion of AQP4-reactive T cells. In contrast, myelin oligodendrocyte glycoprotein (MOG)-specific T cells survived and caused sustained disease in WT mice. These findings underscore the importance of peripheral T cell deletional tolerance to AQP4, which may be relevant to understanding the balance of AQP4-reactive T cells in health and in NMO. T cell tolerance to AQP4, expressed in multiple tissues, is distinct from tolerance to MOG, an autoantigen restricted in its expression.

Genetic variation in the immunoglobulin heavy chain locus shapes the human antibody repertoire

Variation in the antibody response has been linked to differential outcomes in disease, and suboptimal vaccine and therapeutic responsiveness, the determinants of which have not been fully elucidated. Countering models that presume antibodies are generated largely by stochastic processes, we demonstrate that polymorphisms within the immunoglobulin heavy chain locus (IGH) impact the naive and antigen-experienced antibody repertoire, indicating that genetics predisposes individuals to mount qualitatively and quantitatively different antibody responses. We pair recently developed long-read genomic sequencing methods with antibody repertoire profiling to comprehensively resolve IGH genetic variation, including novel structural variants, single nucleotide variants, and genes and alleles. We show that IGH germline variants determine the presence and frequency of antibody genes in the expressed repertoire, including those enriched in functional elements linked to V(D)J recombination, and overlapping disease-associated variants. These results illuminate the power of leveraging IGH genetics to better understand the regulation, function, and dynamics of the antibody response in disease.

Vaccine-boosted CAR T crosstalk with host immunity to reject tumors with antigen heterogeneity

Chimeric antigen receptor (CAR) T cell therapy effectively treats human cancer, but the loss of the antigen recognized by the CAR poses a major obstacle. We found that in vivo vaccine boosting of CAR T cells triggers the engagement of the endogenous immune system to circumvent antigen-negative tumor escape. Vaccine-boosted CAR T promoted dendritic cell (DC) recruitment to tumors, increased tumor antigen uptake by DCs, and elicited the priming of endogenous anti-tumor T cells. This process was accompanied by shifts in CAR T metabolism toward oxidative phosphorylation (OXPHOS) and was critically dependent on CAR-T-derived IFN-γ. Antigen spreading (AS) induced by vaccine-boosted CAR T enabled a proportion of complete responses even when the initial tumor was 50% CAR antigen negative, and heterogeneous tumor control was further enhanced by the genetic amplification of CAR T IFN-γ expression. Thus, CAR-T-cell-derived IFN-γ plays a critical role in promoting AS, and vaccine boosting provides a clinically translatable strategy to drive such responses against solid tumors.

Dietary restriction mitigates the age-associated decline in mouse B cell receptor repertoire diversity

Aging impairs the capacity to respond to novel antigens, reducing immune protection against pathogens and vaccine efficacy. Dietary restriction (DR) extends life- and health span in diverse animals. However, little is known about the capacity of DR to combat the decline in immune function. Here, we study the changes in B cell receptor (BCR) repertoire during aging in DR and control mice. By sequencing the variable region of the BCR heavy chain in the spleen, we show that DR preserves diversity and attenuates the increase in clonal expansions throughout aging. Remarkably, mice starting DR in mid-life have repertoire diversity and clonal expansion rates indistinguishable from chronic DR mice. In contrast, in the intestine, these traits are unaffected by either age or DR. Reduced within-individual B cell repertoire diversity and increased clonal expansions are correlated with higher morbidity, suggesting a potential contribution of B cell repertoire dynamics to health during aging.

STAT3 signaling in B cells controls germinal center zone organization and recycling

Germinal centers (GCs), sites of antibody affinity maturation, are organized into dark (DZ) and light (LZ) zones. Here, we show a B cell-intrinsic role for signal transducer and activator of transcription 3 (STAT3) in GC DZ and LZ organization. Altered zonal organization of STAT3-deficient GCs dampens development of long-lived plasma cells (LL-PCs) but increases memory B cells (MBCs). In an abundant antigenic environment, achieved here by prime-boost immunization, STAT3 is not required for GC initiation, maintenance, or proliferation but is important for sustaining GC zonal organization by regulating GC B cell recycling. Th cell-derived signals drive STAT3 tyrosine 705 and serine 727 phosphorylation in LZ B cells, regulating their recycling into the DZ. RNA sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) analyses identified STAT3 regulated genes that are critical for LZ cell recycling and transiting through DZ proliferation and differentiation phases. Thus, STAT3 signaling in B cells controls GC zone organization and recycling, and GC egress of PCs, but negatively regulates MBC output.

Statistical inference reveals the role of length, GC content, and local sequence in V(D)J nucleotide trimming

To appropriately defend against a wide array of pathogens, humans somatically generate highly diverse repertoires of B cell and T cell receptors (BCRs and TCRs) through a random process called V(D)J recombination. Receptor diversity is achieved during this process through both the combinatorial assembly of V(D)J-genes and the junctional deletion and insertion of nucleotides. While the Artemis protein is often regarded as the main nuclease involved in V(D)J recombination, the exact mechanism of nucleotide trimming is not understood. Using a previously published TCRβ repertoire sequencing data set, we have designed a flexible probabilistic model of nucleotide trimming that allows us to explore various mechanistically interpretable sequence-level features. We show that local sequence context, length, and GC nucleotide content in both directions of the wider sequence, together, can most accurately predict the trimming probabilities of a given V-gene sequence. Because GC nucleotide content is predictive of sequence-breathing, this model provides quantitative statistical evidence regarding the extent to which double-stranded DNA may need to be able to breathe for trimming to occur. We also see evidence of a sequence motif that appears to get preferentially trimmed, independent of GC-content-related effects. Further, we find that the inferred coefficients from this model provide accurate prediction for V- and J-gene sequences from other adaptive immune receptor loci. These results refine our understanding of how the Artemis nuclease may function to trim nucleotides during V(D)J recombination and provide another step toward understanding how V(D)J recombination generates diverse receptors and supports a powerful, unique immune response in healthy humans.

Disrupted memory T cell expansion in HIV-exposed uninfected infants is preceded by premature skewing of T cell receptor clonality

While preventing vertical HIV transmission has been very successful, the increasing number of HIV-exposed uninfected infants (iHEU) experience an elevated risk to infections compared to HIV-unexposed and uninfected infants (iHUU). Immune developmental differences between iHEU and iHUU remains poorly understood and here we present a longitudinal multimodal analysis of infant immune ontogeny that highlights the impact of HIV/ARV exposure. Using mass cytometry, we show alterations and differences in the emergence of NK cell populations and T cell memory differentiation between iHEU and iHUU. Specific NK cells observed at birth were also predictive of acellular pertussis and rotavirus vaccine-induced IgG and IgA responses, respectively, at 3 and 9 months of life. T cell receptor Vβ clonotypic diversity was significantly and persistently lower in iHEU preceding the expansion of T cell memory. Our findings show that HIV/ARV exposure disrupts innate and adaptive immunity from birth which may underlie relative vulnerability to infections.

Abalign: a comprehensive multiple sequence alignment platform for B-cell receptor immune repertoires

The utilization of high-throughput sequencing (HTS) for B-cell receptor (BCR) immune repertoire analysis has become widespread in the fields of adaptive immunity and antibody drug development. However, the sheer volume of sequences generated by these experiments presents a challenge in data processing. Specifically, multiple sequence alignment (MSA), a critical aspect of BCR analysis, remains inadequate for handling massive BCR sequencing data and lacks the ability to provide immunoglobulin-specific information. To address this gap, we introduce Abalign, a standalone program specifically designed for ultrafast MSA of BCR/antibody sequences. Benchmark tests demonstrate that Abalign achieves comparable or even better accuracy than state-of-the-art MSA tools, and shows remarkable advantages in terms of speed and memory consumption, reducing the time required for high-throughput analysis from weeks to hours. In addition to its alignment capabilities, Abalign offers a broad range of BCR analysis features, including extracting BCRs, constructing lineage trees, assigning VJ genes, analyzing clonotypes, profiling mutations, and comparing BCR immune repertoires. With its user-friendly graphic interface, Abalign can be easily run on personal computers instead of computing clusters. Overall, Abalign is an easy-to-use and effective tool that enables researchers to analyze massive BCR/antibody sequences, leading to new discoveries in the field of immunoinformatics. The software is freely available at http://cao.labshare.cn/abalign/.

FLAIRR-Seq: A Method for Single-Molecule Resolution of Near Full-Length Antibody H Chain Repertoires

Current Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using short-read sequencing strategies resolve expressed Ab transcripts with limited resolution of the C region. In this article, we present the near-full-length AIRR-seq (FLAIRR-seq) method that uses targeted amplification by 5' RACE, combined with single-molecule, real-time sequencing to generate highly accurate (99.99%) human Ab H chain transcripts. FLAIRR-seq was benchmarked by comparing H chain V (IGHV), D (IGHD), and J (IGHJ) gene usage, complementarity-determining region 3 length, and somatic hypermutation to matched datasets generated with standard 5' RACE AIRR-seq using short-read sequencing and full-length isoform sequencing. Together, these data demonstrate robust FLAIRR-seq performance using RNA samples derived from PBMCs, purified B cells, and whole blood, which recapitulated results generated by commonly used methods, while additionally resolving H chain gene features not documented in IMGT at the time of submission. FLAIRR-seq data provide, for the first time, to our knowledge, simultaneous single-molecule characterization of IGHV, IGHD, IGHJ, and IGHC region genes and alleles, allele-resolved subisotype definition, and high-resolution identification of class switch recombination within a clonal lineage. In conjunction with genomic sequencing and genotyping of IGHC genes, FLAIRR-seq of the IgM and IgG repertoires from 10 individuals resulted in the identification of 32 unique IGHC alleles, 28 (87%) of which were previously uncharacterized. Together, these data demonstrate the capabilities of FLAIRR-seq to characterize IGHV, IGHD, IGHJ, and IGHC gene diversity for the most comprehensive view of bulk-expressed Ab repertoires to date.

SARS-CoV-2 Omicron boosting induces de novo B cell response in humans

The primary two-dose SARS-CoV-2 mRNA vaccine series are strongly immunogenic in humans, but the emergence of highly infectious variants necessitated additional doses and the development of vaccines aimed at the new variants1-4. SARS-CoV-2 booster immunizations in humans primarily recruit pre-existing memory B cells5-9. However, it remains unclear whether the additional doses induce germinal centre reactions whereby re-engaged B cells can further mature, and whether variant-derived vaccines can elicit responses to variant-specific epitopes. Here we show that boosting with an mRNA vaccine against the original monovalent SARS-CoV-2 mRNA vaccine or the bivalent B.1.351 and B.1.617.2 (Beta/Delta) mRNA vaccine induced robust spike-specific germinal centre B cell responses in humans. The germinal centre response persisted for at least eight weeks, leading to significantly more mutated antigen-specific bone marrow plasma cell and memory B cell compartments. Spike-binding monoclonal antibodies derived from memory B cells isolated from individuals boosted with either the original SARS-CoV-2 spike protein, bivalent Beta/Delta vaccine or a monovalent Omicron BA.1-based vaccine predominantly recognized the original SARS-CoV-2 spike protein. Nonetheless, using a more targeted sorting approach, we isolated monoclonal antibodies that recognized the BA.1 spike protein but not the original SARS-CoV-2 spike protein from individuals who received the mRNA-1273.529 booster; these antibodies were less mutated and recognized novel epitopes within the spike protein, suggesting that they originated from naive B cells. Thus, SARS-CoV-2 booster immunizations in humans induce robust germinal centre B cell responses and can generate de novo B cell responses targeting variant-specific epitopes.

Altered somatic hypermutation patterns in COVID-19 patients classifies disease severity

Introduction

The success of the human body in fighting SARS-CoV2 infection relies on lymphocytes and their antigen receptors. Identifying and characterizing clinically relevant receptors is of utmost importance.

Methods

We report here the application of a machine learning approach, utilizing B cell receptor repertoire sequencing data from severely and mildly infected individuals with SARS-CoV2 compared with uninfected controls.

Results

In contrast to previous studies, our approach successfully stratifies non-infected from infected individuals, as well as disease level of severity. The features that drive this classification are based on somatic hypermutation patterns, and point to alterations in the somatic hypermutation process in COVID-19 patients.

Discussion

These features may be used to build and adapt therapeutic strategies to COVID-19, in particular to quantitatively assess potential diagnostic and therapeutic antibodies. These results constitute a proof of concept for future epidemiological challenges.

Multiomics Empowers Predictive Pancreatic Cancer Immunotherapy

Advances in cancer immunotherapy, particularly immune checkpoint inhibitors, have dramatically improved the prognosis for patients with metastatic melanoma and other previously incurable cancers. However, patients with pancreatic ductal adenocarcinoma (PDAC) generally do not respond to these therapies. PDAC is exceptionally difficult to treat because of its often late stage at diagnosis, modest mutation burden, and notoriously complex and immunosuppressive tumor microenvironment. Simultaneously interrogating features of cancer, immune, and other cellular components of the PDAC tumor microenvironment is therefore crucial for identifying biomarkers of immunotherapeutic resistance and response. Notably, single-cell and multiomics technologies, along with the analytical tools for interpreting corresponding data, are facilitating discoveries of the systems-level cellular and molecular interactions contributing to the overall resistance of PDAC to immunotherapy. Thus, in this review, we will explore how multiomics and single-cell analyses provide the unprecedented opportunity to identify biomarkers of resistance and response to successfully sensitize PDAC to immunotherapy.

Characterisation of the immune repertoire of a humanised transgenic mouse through immunophenotyping and high-throughput sequencing

Immunoglobulin loci-transgenic animals are widely used in antibody discovery and increasingly in vaccine response modelling. In this study, we phenotypically characterised B-cell populations from the Intelliselect Transgenic mouse (Kymouse) demonstrating full B-cell development competence. Comparison of the naïve B-cell receptor (BCR) repertoires of Kymice BCRs, naïve human, and murine BCR repertoires revealed key differences in germline gene usage and junctional diversification. These differences result in Kymice having CDRH3 length and diversity intermediate between mice and humans. To compare the structural space explored by CDRH3s in each species' repertoire, we used computational structure prediction to show that Kymouse naïve BCR repertoires are more human-like than mouse-like in their predicted distribution of CDRH3 shape. Our combined sequence and structural analysis indicates that the naïve Kymouse BCR repertoire is diverse with key similarities to human repertoires, while immunophenotyping confirms that selected naïve B cells are able to go through complete development.

B cell M-CLL clones retain selection against replacement mutations in their immunoglobulin gene framework regions

Introduction

Chronic lymphocytic leukemia (CLL) is the most common adult leukemia, accounting for 30–40% of all adult leukemias. The dynamics of B-lymphocyte CLL clones with mutated immunoglobulin heavy chain variable region (IgHV) genes in their tumor (M-CLL) can be studied using mutational lineage trees.

Methods

Here, we used lineage tree-based analyses of somatic hypermutation (SHM) and selection in M-CLL clones, comparing the dominant (presumably malignant) clones of 15 CLL patients to their non-dominant (presumably normal) B cell clones, and to those of healthy control repertoires. This type of analysis, which was never previously published in CLL, yielded the following novel insights.

Results

CLL dominant clones undergo – or retain – more replacement mutations that alter amino acid properties such as charge or hydropathy. Although, as expected, CLL dominant clones undergo weaker selection for replacement mutations in the complementarity determining regions (CDRs) and against replacement mutations in the framework regions (FWRs) than non-dominant clones in the same patients or normal B cell clones in healthy controls, they surprisingly retain some of the latter selection in their FWRs. Finally, using machine learning, we show that even the non-dominant clones in CLL patients differ from healthy control clones in various features, most notably their expression of higher fractions of transition mutations.

Discussion

Overall, CLL seems to be characterized by significant loosening – but not a complete loss – of the selection forces operating on B cell clones, and possibly also by changes in SHM mechanisms.

Evidence of microbial activity in a uranium roll-front deposit: Unlocking their potential role as bioenhancers of the ore genesis

Uranium (U) roll-front deposits constitute a valuable source for an economical extraction by in situ recovery (ISR) mining. Such technology may induce changes in the subsurface microbiota, raising questions about the way their activities could build a functional ecosystem in such extreme environments (i.e.: oligotrophy and high SO₄ concentration and salinity). Additionally, more information is needed to dissipate the doubts about the microbial role in the genesis of such U orebodies. A U roll-front deposit hosted in an aquifer driven system (in Zoovch Ovoo, Mongolia), intended for mining by acid ISR, was previously explored and showed to be governed by a complex bacterial diversity, linked to the redox zonation and the geochemical conditions. Here for the first time, transcriptional activities of microorganisms living in such U ore deposits are determined and their metabolic capabilities allocated in the three redox-inherited compartments, naturally defined by the roll-front system. Several genes encoding for crucial metabolic pathways demonstrated a strong biological role controlling the subsurface cycling of many elements including nitrate, sulfate, metals and radionuclides (e.g.: uranium), through oxidation-reduction reactions. Interestingly, the discovered transcriptional behaviour gives important insights into the good microbial adaptation to the geochemical conditions and their active contribution to the stabilization of the U ore deposits. Overall, evidences on the importance of these microbial metabolic activities in the aquifer system are discussed that may clarify the doubts on the microbial role in the genesis of low-temperature U roll-front deposits, along the Zoovch Ovoo mine.

Heparan sulfate regulates IL-21 bioavailability and signal strength that control germinal center B cell selection and differentiation

In antibody responses, mutated germinal center B (B_GC) cells are positively selected for reentry or differentiation. As the products from GCs, memory B cells and antibody-secreting cells (ASCs) support high-affinity and long-lasting immunity. Positive selection of B_GC cells is controlled by signals received through the B cell receptor (BCR) and follicular helper T (T_FH) cell-derived signals, in particular costimulation through CD40. Here, we demonstrate that the T_FH cell effector cytokine interleukin-21 (IL-21) joins BCR and CD40 in supporting B_GC selection and reveal that strong IL-21 signaling prioritizes ASC differentiation in vivo. B_GC cells, compared with non-B_GC cells, show significantly reduced IL-21 binding and attenuated signaling, which is mediated by low cellular heparan sulfate (HS) sulfation. Mechanistically, N-deacetylase and N-sulfotransferase 1 (Ndst1)-mediated N-sulfation of HS in B cells promotes IL-21 binding and signal strength. Ndst1 is down-regulated in B_GC cells and up-regulated in ASC precursors, suggesting selective desensitization to IL-21 in B_GC cells. Thus, specialized biochemical regulation of IL-21 bioavailability and signal strength sets a balance between the stringency and efficiency of GC selection.

The temporal behavior of the murine T cell receptor repertoire following Immunotherapy

Immunotherapy is now an essential tool for cancer treatment, and the unique features of an individual's T cell receptor repertoire are known to play a key role in its effectiveness. The repertoire, famously vast due to a cascade of cellular mechanisms, can be quantified using repertoire sequencing. In this study, we sampled the repertoire over several time points following treatment with anti-CTLA-4, in a syngeniec mouse model for colorectal cancer, generating a longitudinal dataset of T cell clones and their abundance. The dynamics of the repertoire can be observed in response to treatment over a period of four weeks, as clonal expansion of specific clones ascends and descends. The data made available here can be used to determine treatment and predict its effect, while also providing a unique look at the behavior of the immune system over time.

Tissue-specific abundance of interferon-gamma drives regulatory T cells to restrain DC1-mediated priming of cytotoxic T cells against lung cancer

Local environmental factors influence CD8+ T cell priming in lymph nodes (LNs). Here, we sought to understand how factors unique to the tumor-draining mediastinal LN (mLN) impact CD8+ T cell responses toward lung cancer. Type 1 conventional dendritic cells (DC1s) showed a mLN-specific failure to induce robust cytotoxic T cells responses. Using regulatory T (Treg) cell depletion strategies, we found that Treg cells suppressed DC1s in a spatially coordinated manner within tissue-specific microniches within the mLN. Treg cell suppression required MHC II-dependent contact between DC1s and Treg cells. Elevated levels of IFN-γ drove differentiation Treg cells into Th1-like effector Treg cells in the mLN. In patients with cancer, Treg cell Th1 polarization, but not CD8+/Treg cell ratios, correlated with poor responses to checkpoint blockade immunotherapy. Thus, IFN-γ in the mLN skews Treg cells to be Th1-like effector Treg cells, driving their close interaction with DC1s and subsequent suppression of cytotoxic T cell responses.

Modeling and predicting the overlap of B- and T-cell receptor repertoires in healthy and SARS-CoV-2 infected individuals

Adaptive immunity's success relies on the extraordinary diversity of protein receptors on B and T cell membranes. Despite this diversity, the existence of public receptors shared by many individuals gives hope for developing population-wide vaccines and therapeutics. Using probabilistic modeling, we show many of these public receptors are shared by chance in healthy individuals. This predictable overlap is driven not only by biases in the random generation process of receptors, as previously reported, but also by their common functional selection. However, the model underestimates sharing between repertoires of individuals infected with SARS-CoV-2, suggesting strong specific antigen-driven convergent selection. We exploit this discrepancy to identify COVID-associated receptors, which we validate against datasets of receptors with known viral specificity. We study their properties in terms of sequence features and network organization, and use them to design an accurate diagnostic tool for predicting SARS-CoV-2 status from repertoire data.

GENTLE: a novel bioinformatics tool for generating features and building classifiers from T cell repertoire cancer data

BACKGROUND

In the global effort to discover biomarkers for cancer prognosis, prediction tools have become essential resources. TCR (T cell receptor) repertoires contain important features that differentiate healthy controls from cancer patients or differentiate outcomes for patients being treated with different drugs. Considering, tools that can easily and quickly generate and identify important features out of TCR repertoire data and build accurate classifiers to predict future outcomes are essential.

RESULTS

This paper introduces GENTLE (GENerator of T cell receptor repertoire features for machine LEarning): an open-source, user-friendly web-application tool that allows TCR repertoire researchers to discover important features; to create classifier models and evaluate them with metrics; and to quickly generate visualizations for data interpretations. We performed a case study with repertoires of TRegs (regulatory T cells) and TConvs (conventional T cells) from healthy controls versus patients with breast cancer. We showed that diversity features were able to distinguish between the groups. Moreover, the classifiers built with these features could correctly classify samples ('Healthy' or 'Breast Cancer')from the TRegs repertoire when trained with the TConvs repertoire, and from the TConvs repertoire when trained with the TRegs repertoire.

CONCLUSION

The paper walks through installing and using GENTLE and presents a case study and results to demonstrate the application's utility. GENTLE is geared towards any researcher working with TCR repertoire data and aims to discover predictive features from these data and build accurate classifiers. GENTLE is available on https://github.com/dhiego22/gentle and https://share.streamlit.io/dhiego22/gentle/main/gentle.py .

T-independent responses to polysaccharides in humans mobilize marginal zone B cells prediversified against gut bacterial antigens

Marginal zone (MZ) B cells are one of the main actors of T-independent (TI) responses in mice. To identify the B cell subset(s) involved in such responses in humans, we vaccinated healthy individuals with Pneumovax, a model TI vaccine. By high-throughput repertoire sequencing of plasma cells (PCs) isolated 7 days after vaccination and of different B cell subpopulations before and after vaccination, we show that the PC response mobilizes large clones systematically, including an immunoglobulin M component, whose diversification and amplification predated the pneumococcal vaccination. These clones could be mainly traced back to MZ B cells, together with clonally related IgA+ and, to a lesser extent, IgG+CD27+ B cells. Recombinant monoclonal antibodies isolated from large PC clones recognized a wide array of bacterial species from the gut flora, indicating that TI responses in humans largely mobilize MZ and switched B cells that most likely prediversified during mucosal immune responses against bacterial antigens and acquired pneumococcal cross-reactivity through somatic hypermutation.

The memory of pathogenic IgE is contained within CD23 + IgG1 + memory B cells poised to switch to IgE in food allergy

Food allergy is caused by allergen-specific IgE antibodies but little is known about the B cell memory of persistent IgE responses. Here we describe in human pediatric peanut allergy CD23 + IgG1 + memory B cells arising in type 2 responses that contain peanut specific clones and generate IgE cells on activation. These 'type2-marked' IgG1 + memory B cells differentially express IL-4/IL-13 regulated genes FCER2 / CD23, IL4R , and germline IGHE and carry highly mutated B cell receptors (BCRs). Further, high affinity memory B cells specific for the main peanut allergen Ara h 2 mapped to the population of 'type2-marked' IgG1 + memory B cells and included convergent BCRs across different individuals. Our findings indicate that CD23 + IgG1 + memory B cells transcribing germline IGHE are a unique memory population containing precursors of pathogenic IgE.

One-Sentence Summary

We describe a unique population of IgG + memory B cells poised to switch to IgE that contains high affinity allergen-specific clones in peanut allergy.

Distinct Genetically Determined Origins of Myd88/BCL2-Driven Aggressive Lymphoma Rationalize Targeted Therapeutic Intervention Strategies

Genomic profiling revealed the identity of at least 5 subtypes of diffuse large B-cell lymphoma (DLBCL), including the MCD/C5 cluster characterized by aberrations in MYD88, BCL2, PRDM1, and/or SPIB. We generated mouse models harboring B cell-specific Prdm1 or Spib aberrations on the background of oncogenic Myd88 and Bcl2 lesions. We deployed whole-exome sequencing, transcriptome, flow-cytometry, and mass cytometry analyses to demonstrate that Prdm1- or Spib-altered lymphomas display molecular features consistent with prememory B cells and light-zone B cells, whereas lymphomas lacking these alterations were enriched for late light-zone and plasmablast-associated gene sets. Consistent with the phenotypic evidence for increased B cell receptor signaling activity in Prdm1-altered lymphomas, we demonstrate that combined BTK/BCL2 inhibition displays therapeutic activity in mice and in five of six relapsed/refractory DLBCL patients. Moreover, Prdm1-altered lymphomas were immunogenic upon transplantation into immuno-competent hosts, displayed an actionable PD-L1 surface expression, and were sensitive to antimurine-CD19-CAR-T cell therapy, in vivo.

SIGNIFICANCE

Relapsed/refractory DLBCL remains a major medical challenge, and most of these patients succumb to their disease. Here, we generated mouse models, faithfully recapitulating the biology of MYD88-driven human DLBCL. These models revealed robust preclinical activity of combined BTK/BCL2 inhibition. We confirmed activity of this regimen in pretreated non-GCB-DLBCL patients. See related commentary by Leveille et al., p. 8. This article is highlighted in the In This Issue feature, p. 1.

Accurate estimation of molecular counts from amplicon sequence data with unique molecular identifiers

MOTIVATION

Amplicon sequencing is widely applied to explore heterogeneity and rare variants in genetic populations. Resolving true biological variants and quantifying their abundance is crucial for downstream analyses, but measured abundances are distorted by stochasticity and bias in amplification, plus errors during polymerase chain reaction (PCR) and sequencing. One solution attaches unique molecular identifiers (UMIs) to sample sequences before amplification. Counting UMIs instead of sequences provides unbiased estimates of abundance. While modern methods improve over naïve counting by UMI identity, most do not account for UMI reuse or collision, and they do not adequately model PCR and sequencing errors in the UMIs and sample sequences.

RESULTS

We introduce Deduplication and Abundance estimation with UMIs (DAUMI), a probabilistic framework to detect true biological amplicon sequences and accurately estimate their deduplicated abundance. DAUMI recognizes UMI collision, even on highly similar sequences, and detects and corrects most PCR and sequencing errors in the UMI and sampled sequences. DAUMI performs better on simulated and real data compared to other UMI-aware clustering methods.

AVAILABILITY AND IMPLEMENTATION

Source code is available at https://github.com/DormanLab/AmpliCI.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bioinformatic Analysis of Natively Paired VH:VL Antibody Repertoires for Antibody Discovery

Next-generation DNA sequencing (NGS) of human antibody repertoires has been extensively implemented to discover novel antibody drugs, to analyze B-cell developmental features, and to investigate antibody responses to infectious diseases and vaccination. Because the antibody repertoire encoded by human B cells is highly diverse, NGS analyses of antibody genes have provided a new window into understanding antibody responses for basic immunology, biopharmaceutical drug discovery, and immunotherapy. However, many antibody discovery protocols analyze the heavy and light chains separately due to the short-read nature of most NGS technologies, whereas paired heavy and light chain data are required for complete antibody characterization. Here, we describe a computational workflow to process millions of paired antibody heavy and light chain DNA sequence reads using the Illumina MiSeq 2x300 NGS platform. In this workflow, we describe raw NGS read processing and initial quality filtering, the annotation and assembly of antibody clonotypes relating to paired heavy and light chain antibody lineages, and the generation of complete heavy+light consensus sequences for the downstream cloning and expression of human antibody proteins.

Structural Modeling of Adaptive Immune Responses to Infection

Antibody and TCR modeling are becoming important as more and more sequence data becomes available to the public. One of the pressing questions now is how to use such data to understand adaptive immune responses to disease. Infectious disease is of particular interest because the antigens driving such responses are often known to some extent. Here, we describe tips for gathering data and cleaning it for use in downstream analysis. We present a method for high-throughput structural modeling of antibodies or TCRs using Repertoire Builder and its extensions. AbAdapt is an extension of Repertoire Builder for antibody-antigen docking from antibody and antigen sequences. ImmuneScape is a corresponding extension for TCR-pMHC 3D modeling. Together, these pipelines can help researchers to understand immune responses to infection from a structural point of view.

A somatic hypermutation-based machine learning model stratifies individuals with Crohn's disease and controls

Crohn's disease (CD) is a chronic relapsing-remitting inflammatory disorder of the gastrointestinal tract that is characterized by altered innate and adaptive immune function. Although massively parallel sequencing studies of the T cell receptor repertoire identified oligoclonal expansion of unique clones, much less is known about the B cell receptor (BCR) repertoire in CD. Here, we present a novel BCR repertoire sequencing data set from ileal biopsies from pediatric patients with CD and controls, and identify CD-specific somatic hypermutation (SHM) patterns, revealed by a machine learning (ML) algorithm trained on BCR repertoire sequences. Moreover, ML classification of a different data set from blood samples of adults with CD versus controls identified that V gene usage, clusters, or mutation frequencies yielded excellent results in classifying the disease (F1 > 90%). In summary, we show that an ML algorithm enables the classification of CD based on unique BCR repertoire features with high accuracy.

Immune repertoire profiling for disease pathobiology

Lymphocytes consist of highly heterogeneous populations, each expressing a specific cell surface receptor corresponding to a particular antigen. Lymphocytes are both the cause and regulator of various diseases, including autoimmune/allergic diseases, lifestyle diseases, neurodegenerative diseases, and cancers. Recently, immune repertoire sequencing has attracted much attention because it helps obtain global profiles of the immune receptor sequences of infiltrating T and B cells in specimens. Immune repertoire sequencing not only helps deepen our understanding of the molecular mechanisms of immune-related pathology but also assists in discovering novel therapeutic modalities for diseases, thereby shedding colorful light on otherwise tiny monotonous cells when observed under a microscope. In this review article, we introduce and detail the background and methodology of immune repertoire sequencing and summarize recent scientific achievements in association with human diseases. Future perspectives on this genetic technique in the field of histopathological research will also be discussed.

Cladribine treatment specifically affects peripheral blood memory B cell clones and clonal expansion in multiple sclerosis patients

Introduction

B cells are acknowledged as crucial players in the pathogenesis of multiple sclerosis (MS). Several disease modifying drugs including cladribine have been shown to exert differential effects on peripheral blood B cell subsets. However, little is known regarding functional changes within the peripheral B cell populations. In this study, we obtained a detailed picture of B cell repertoire changes under cladribine treatment on a combined immunoglobulin (Ig) transcriptome and proteome level.

Methods

We performed next-generation sequencing of Ig heavy chain (IGH) transcripts and Ig mass spectrometry in cladribine-treated patients with relapsing-remitting multiple sclerosis (n = 8) at baseline and after 6 and 12 months of treatment in order to generate Ig transcriptome and Ig peptide libraries. Ig peptides were overlapped with the corresponding IGH transcriptome in order to analyze B cell clones on a combined transcriptome and proteome level.

Results

The analysis of peripheral blood B cell percentages pointed towards a significant decrease of memory B cells and an increase of naive B cells following cladribine therapy. While basic IGH repertoire parameters (e.g. variable heavy chain family usage and Ig subclasses) were only slightly affected by cladribine treatment, a significantly decreased number of clones and significantly lower diversity in the memory subset was noticeable at 6 months following treatment which was sustained at 12 months. When looking at B-cell clones comprising sequences from the different time-points, clones spanning between all three time-points were significantly more frequent than clones including sequences from two time-points. Furthermore, Ig proteome analyses showed that Ig transcriptome specific peptides could mostly be equally aligned to all three time-points pointing towards a proportion of B-cell clones that are maintained during treatment.

Discussion

Our findings suggest that peripheral B cell related treatment effects of cladribine tablets might be exerted through a reduction of possibly disease relevant clones in the memory B cell subset without disrupting the overall clonal composition of B cells. Our results -at least partially- might explain the relatively mild side effects regarding infections and the sustained immune response after vaccinations during treatment. However, exact disease driving B cell subsets and their effects remain unknown and should be addressed in future studies.

T-cell repertoire diversity: friend or foe for protective antitumor response?

Profiling the T-Cell Receptor (TCR) repertoire is establishing as a potent approach to investigate autologous and treatment-induced antitumor immune response. Technical and computational breakthroughs, including high throughput next-generation sequencing (NGS) approaches and spatial transcriptomics, are providing unprecedented insight into the mechanisms underlying antitumor immunity. A precise spatiotemporal variation of T-cell repertoire, which dynamically mirrors the functional state of the evolving host-cancer interaction, allows the tracking of the T-cell populations at play, and may identify the key cells responsible for tumor eradication, the evaluation of minimal residual disease and the identification of biomarkers of response to immunotherapy. In this review we will discuss the relationship between global metrics characterizing the TCR repertoire such as T-cell clonality and diversity and the resultant functional responses. In particular, we will explore how specific TCR repertoires in cancer patients can be predictive of prognosis or response to therapy and in particular how a given TCR re-arrangement, following immunotherapy, can predict a specific clinical outcome. Finally, we will examine current improvements in terms of T-cell sequencing, discussing advantages and challenges of current methodologies.