IMGT/GENE-DB: a comprehensive database for human and mouse immunoglobulin and T cell receptor genes

CLL stereotyped B-cell receptor immunoglobulin sequences are recurrent in the B-cell repertoire of healthy individuals: Apparent lack of central and early peripheral tolerance censoring

Introduction

The leukemic cells of patients with chronic lymphocytic leukemia (CLL) are often unique, expressing remarkably similar IGHV-IGHD-IGHJ gene rearrangements, “stereotyped BCRs”. The B-cell receptors (BCRs) on CLL cells are also distinctive in often deriving from autoreactive B lymphocytes, leading to the assumption of a defect in immune tolerance.

Results

Using bulk and single-cell immunoglobulin heavy and light chain variable domain sequencing, we enumerated CLL stereotype-like IGHV-IGHD-IGHJ sequences (CLL-SLS) in B cells from cord blood (CB) and adult peripheral blood (PBMC) and bone marrow (BM of healthy donors. CLL-SLS were found at similar frequencies among CB, BM, and PBMC, suggesting that age does not influence CLL-SLS levels. Moreover, the frequencies of CLL-SLS did not differ among B lymphocytes in the BM at early stages of development, and only re-circulating marginal zone B cells contained significantly higher CLL-SLS frequencies than other mature B-cell subpopulations. Although we identified CLL-SLS corresponding to most of the CLL major stereotyped subsets, CLL-SLS frequencies did not correlate with those found in patients. Interestingly, in CB samples, half of the CLL-SLS identified were attributed to two IGHV-mutated subsets. We also found satellite CLL-SLS among the same normal samples, and they were also enriched in naïve B cells but unexpectedly, these were ~10-fold higher than standard CLL-SLS. In general, IGHV-mutated CLL-SLS subsets were enriched among antigen-experienced B-cell subpopulations, and IGHV-unmutated CLL-SLS were found mostly in antigen-inexperienced B cells. Nevertheless, CLL-SLS with an IGHV-mutation status matching that of CLL clones varied among the normal B-cell subpopulations, suggesting that specific CLL-SLS could originate from distinct subpopulations of normal B cells. Lastly, using single-cell DNA sequencing, we identified paired IGH and IGL rearrangements in normal B lymphocytes resembling those of stereotyped BCRs in CLL, although some differed from those in patients based on IG isotype or somatic mutation.

Discussion

CLL-SLS are present in normal B-lymphocyte populations at all stages of development. Thus, despite their autoreactive profile they are not deleted by central tolerance mechanisms, possibly because the level of autoreactivity is not registered as dangerous by deletion mechanisms or because editing of L-chain variable genes occurred which our experimental approach could not identify.

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.

Identification of AL proteins from 10 λ-AL amyloidosis patients by mass spectrometry extracted from abdominal fat and heart tissue

BACKGROUND

Systemic AL amyloidosis arises from the misfolding of patient-specific immunoglobulin light chains (LCs). Potential drivers of LC amyloid formation are mutational changes and post-translational modifications (PTMs). However, little information is available on the exact primary structure of the AL proteins and their precursor LCs.

OBJECTIVE

We analyse the exact primary structure of AL proteins extracted from 10 λ AL amyloidosis patients and their corresponding precursor LCs.

MATERIALS AND METHODS

By cDNA sequencing of the precursor LC genes in combination with mass spectrometry of the AL proteins, the exact primary structure and PTMs were determined. This information was used to analyse their biochemical properties.

RESULTS

All AL proteins comprise the V_L and a small part of the C_L with a common C-terminal truncation region. While all AL proteins retain the conserved native disulphide bond of the V_L, we found no evidence for presence of other common PTMs. The analysis of the biochemical properties revealed that the isoelectric point of the V_L is significantly increased due to introduced mutations.

CONCLUSION

Our data imply that mutational changes influence the surface charge properties of the V_L and that common proteolytic processes are involved in the generation of the cleavage sites of AL proteins.

Type III interferon drives thymic B cell activation and regulatory T cell generation

The activation of thymic B cells is critical for their licensing as antigen presenting cells and resulting ability to mediate T cell central tolerance. The processes leading to licensing are still not fully understood. By comparing thymic B cells to activated Peyer's patch B cells at steady state, we found that thymic B cell activation starts during the neonatal period and is characterized by TCR/CD40-dependent activation, followed by immunoglobulin class switch recombination (CSR) without forming germinal centers. Transcriptional analysis also demonstrated a strong interferon signature, which was not apparent in the periphery. Thymic B cell activation and CSR were primarily dependent on type III IFN signaling, and loss of type III IFN receptor in thymic B cells resulted in reduced thymocyte regulatory T cell (Treg) development. Finally, from TCR deep sequencing, we estimate that licensed B cells induce development of a substantial fraction of the Treg cell repertoire. Together, these findings reveal the importance of steady-state type III IFN in generating licensed thymic B cells that induce T cell tolerance to activated B cells.

Reconstructing B cell lineage trees with minimum spanning tree and genotype abundances

B cell receptor (BCR) genes exposed to an antigen undergo somatic hypermutations and Darwinian antigen selection, generating a large BCR-antibody diversity. This process, known as B cell affinity maturation, increases antibody affinity, forming a specific B cell lineage that includes the unmutated ancestor and mutated variants. In a B cell lineage, cells with a higher antigen affinity will undergo clonal expansion, while those with a lower affinity will not proliferate and probably be eliminated. Therefore, cellular (genotype) abundance provides a valuable perspective on the ongoing evolutionary process. Phylogenetic tree inference is often used to reconstruct B cell lineage trees and represents the evolutionary dynamic of BCR affinity maturation. However, such methods should process B-cell population data derived from experimental sampling that might contain different cellular abundances. There are a few phylogenetic methods for tracing the evolutionary events occurring in B cell lineages; best-performing solutions are time-demanding and restricted to analysing a reduced number of sequences, while time-efficient methods do not consider cellular abundances. We propose ClonalTree, a low-complexity and accurate approach to construct B-cell lineage trees that incorporates genotype abundances into minimum spanning tree (MST) algorithms. Using both simulated and experimental data, we demonstrate that ClonalTree outperforms MST-based algorithms and achieves a comparable performance to a method that explores tree-generating space exhaustively. Furthermore, ClonalTree has a lower running time, being more convenient for building B-cell lineage trees from high-throughput BCR sequencing data, mainly in biomedical applications, where a lower computational time is appreciable. It is hundreds to thousands of times faster than exhaustive approaches, enabling the analysis of a large set of sequences within minutes or seconds and without loss of accuracy. The source code is freely available at github.com/julibinho/ClonalTree.

NSC243928 Treatment Induces Anti-Tumor Immune Response in Mouse Mammary Tumor Models

NSC243928 induces cell death in triple-negative breast cancer cells in a LY6K-dependent manner. NSC243928 has been reported as an anti-cancer agent in the NCI small molecule library. The molecular mechanism of NSC243928 as an anti-cancer agent in the treatment of tumor growth in the syngeneic mouse model has not been established. With the success of immunotherapies, novel anti-cancer drugs that may elicit an anti-tumor immune response are of high interest in the development of novel drugs to treat solid cancer. Thus, we focused on studying whether NSC243928 may elicit an anti-tumor immune response in the in vivo mammary tumor models of 4T1 and E0771. We observed that NSC243928 induced immunogenic cell death in 4T1 and E0771 cells. Furthermore, NSC243928 mounted an anti-tumor immune response by increasing immune cells such as patrolling monocytes, NKT cells, B1 cells, and decreasing PMN MDSCs in vivo. Further studies are required to understand the exact mechanism of NSC243928 action in inducing an anti-tumor immune response in vivo, which can be used to determine a molecular signature associated with NSC243928 efficacy. NSC243928 may be a good target for future immuno-oncology drug development for breast cancer.

Characterization of the immunoglobulin lambda chain locus from diverse populations reveals extensive genetic variation

Immunoglobulins (IGs), crucial components of the adaptive immune system, are encoded by three genomic loci. However, the complexity of the IG loci severely limits the effective use of short read sequencing, limiting our knowledge of population diversity in these loci. We leveraged existing long read whole-genome sequencing (WGS) data, fosmid technology, and IG targeted single-molecule, real-time (SMRT) long-read sequencing (IG-Cap) to create haplotype-resolved assemblies of the IG Lambda (IGL) locus from 6 ethnically diverse individuals. In addition, we generated 10 diploid assemblies of IGL from a diverse cohort of individuals utilizing IG-Cap. From these 16 individuals, we identified significant allelic diversity, including 36 novel IGLV alleles. In addition, we observed highly elevated single nucleotide variation (SNV) in IGLV genes relative to IGL intergenic and genomic background SNV density. By comparing SNV calls between our high quality assemblies and existing short read datasets from the same individuals, we show a high propensity for false-positives in the short read datasets. Finally, for the first time, we nucleotide-resolved common 5-10 Kb duplications in the IGLC region that contain functional IGLJ and IGLC genes. Together these data represent a significant advancement in our understanding of genetic variation and population diversity in the IGL locus.

PD-1highCXCR5-CD4+ peripheral helper T cells promote CXCR3+ plasmablasts in human acute viral infection

T cell-B cell interaction is the key immune response to protect the host from severe viral infection. However, how T cells support B cells to exert protective humoral immunity in humans is not well understood. Here, we use COVID-19 as a model of acute viral infections and analyze CD4+ T cell subsets associated with plasmablast expansion and clinical outcome. Peripheral helper T cells (Tph cells; denoted as PD-1highCXCR5-CD4+ T cells) are significantly increased, as are plasmablasts. Tph cells exhibit "B cell help" signatures and induce plasmablast differentiation in vitro. Interestingly, expanded plasmablasts show increased CXCR3 expression, which is positively correlated with higher frequency of activated Tph cells and better clinical outcome. Mechanistically, Tph cells help B cell differentiation and produce more interferon γ (IFNγ), which induces CXCR3 expression on plasmablasts. These results elucidate a role for Tph cells in regulating protective B cell response during acute viral infection.

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.

BTG1 mutation yields supercompetitive B cells primed for malignant transformation

Multicellular life requires altruistic cooperation between cells. The adaptive immune system is a notable exception, wherein germinal center B cells compete vigorously for limiting positive selection signals. Studying primary human lymphomas and developing new mouse models, we found that mutations affecting BTG1 disrupt a critical immune gatekeeper mechanism that strictly limits B cell fitness during antibody affinity maturation. This mechanism converted germinal center B cells into supercompetitors that rapidly outstrip their normal counterparts. This effect was conferred by a small shift in MYC protein induction kinetics but resulted in aggressive invasive lymphomas, which in humans are linked to dire clinical outcomes. Our findings reveal a delicate evolutionary trade-off between natural selection of B cells to provide immunity and potentially dangerous features that recall the more competitive nature of unicellular organisms.

Myeloma immunoglobulin rearrangement and translocation detection through targeted capture sequencing

Multiple myeloma is a plasma cell neoplasm characterized by clonal immunoglobulin V(D)J signatures and oncogenic immunoglobulin gene translocations. Additional subclonal genomic changes are acquired with myeloma progression and therapeutic selection. PCR-based methods to detect V(D)J rearrangements can have biases introduced by highly multiplexed reactions and primers undermined by somatic hypermutation, and are not readily extended to include mutation detection. Here, we report a hybrid-capture approach (CapIG-seq) targeting the 3' and 5' ends of the V and J segments of all immunoglobulin loci that enable the efficient detection of V(D)J rearrangements. We also included baits for oncogenic translocations and mutation detection. We demonstrate complete concordance with matched whole-genome sequencing and/or PCR clonotyping of 24 cell lines and report the clonal sequences for 41 uncharacterized cell lines. We also demonstrate the application to patient specimens, including 29 bone marrow and 39 cell-free DNA samples. CapIG-seq shows concordance between bone marrow and cfDNA blood samples (both contemporaneous and follow-up) with regard to the somatic variant, V(D)J, and translocation detection. CapIG-seq is a novel, efficient approach to examining genomic alterations in myeloma.

Antibody Sequence and Structure Analyses Using IMGT®: 30 Years of Immunoinformatics

IMGT^®, the international ImMunoGeneTics information system^®, http://www.imgt.org , the global reference in immunogenetics and immunoinformatics, was created in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS) to manage the huge diversity of the antigen receptors, immunoglobulins (IG) or antibodies, and T cell receptors (TR) of the adaptive immune responses. The founding of IMGT^® marked the advent of immunoinformatics, which emerged at the interface between immunogenetics and bioinformatics. IMGT^® standardized analysis of the IG, TR, and major histocompatibility (MH) genes and proteins bridges the gap between sequences and three-dimensional (3D) structures, for all jawed vertebrates from fish to humans. This is achieved through the IMGT Scientific chart rules, based on the IMGT-ONTOLOGY axioms, and primarily CLASSIFICATION (IMGT gene and allele nomenclature) and NUMEROTATION (IMGT unique numbering and IMGT Colliers de Perles). IMGT^® comprises seven databases (IMGT/LIGM-DB for nucleotide sequences, IMGT/GENE-DB for genes and alleles, etc.), 17 tools (IMGT/V-QUEST, IMGT/JunctionAnalysis, IMGT/HighV-QUEST for NGS, etc.), and more than 20,000 Web resources. In this chapter, the focus is on the tools for amino acid sequences per domain (IMGT/DomainGapAlign and IMGT/Collier-de-Perles), and on the databases for receptors (IMGT/2Dstructure-DB and IMGT/3D-structure-DB) described per receptor, chain, and domain and, for 3D, with contact analysis, paratope, and epitope. The IMGT/mAb-DB is the query interface for monoclonal antibodies (mAb), fusion proteins for immune applications (FPIA), composite proteins for clinical applications (CPCA), and related proteins of interest (RPI) with links to IMGT^® 2D and 3D databases and to the World Health Organization (WHO) International Nonproprietary Names (INN) program lists. The chapter includes the human IG allotypes and antibody engineered variants for effector properties used in the description of therapeutical mAb.

Characterization of heavy-chain antibody gene repertoires in Bactrian camels

Camelids are the only mammals that can produce functional heavy-chain antibodies (HCAbs). Although HCAbs were discovered over 30 years ago, the antibody gene repertoire of Bactrian camels remains largely underexplored. To characterize the diversity of variable genes of HCAbs (VHHs), germline and rearranged VHH repertoires are constructed. Phylogenetics analysis shows that all camelid VHH genes are derived from a common ancestor and the nucleotide diversity of VHHs is similar across all camelid species. While species-specific hallmark sites are identified, the non-canonical cysteines specific to VHHs are distinct in Bactrian camels and dromedaries compared with alpacas. Though low divergence at the germline repertoire between wild and domestic Bactrian camels, higher expression of VHHs is observed in some wild Bactrian camels than that of domestic ones. This study not only adds our understanding of VHH repertoire diversity across camelids, but also provides useful resources for HCAb engineering.

Adaptive immune responses to SARS-CoV-2 persist in the pharyngeal lymphoid tissue of children

Most studies of adaptive immunity to SARS-CoV-2 infection focus on peripheral blood, which may not fully reflect immune responses at the site of infection. Using samples from 110 children undergoing tonsillectomy and adenoidectomy during the COVID-19 pandemic, we identified 24 samples with evidence of previous SARS-CoV-2 infection, including neutralizing antibodies in serum and SARS-CoV-2-specific germinal center and memory B cells in the tonsils and adenoids. Single-cell B cell receptor (BCR) sequencing indicated virus-specific BCRs were class-switched and somatically hypermutated, with overlapping clones in the two tissues. Expanded T cell clonotypes were found in tonsils, adenoids and blood post-COVID-19, some with CDR3 sequences identical to previously reported SARS-CoV-2-reactive T cell receptors (TCRs). Pharyngeal tissues from COVID-19-convalescent children showed persistent expansion of germinal center and antiviral lymphocyte populations associated with interferon (IFN)-γ-type responses, particularly in the adenoids, and viral RNA in both tissues. Our results provide evidence for persistent tissue-specific immunity to SARS-CoV-2 in the upper respiratory tract of children after infection.

Ancestral reconstruction reveals catalytic inactivation of activation-induced cytidine deaminase concomitant with cold water adaption in the Gadiformes bony fish

BACKGROUND

Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID.

RESULTS

By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes.

CONCLUSIONS

Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function.

Transcriptome signatures preceding the induction of anti-stalk antibodies elicited after universal influenza vaccination

A phase 1 clinical trial to test the immunogenicity of a chimeric group 1 HA (cHA) universal influenza virus vaccine targeting the conserved stalk domain of the hemagglutinin of influenza viruses was carried out. Vaccination with adjuvanted-inactivated vaccines induced high anti-stalk antibody titers. We sought to identify gene expression signatures that correlate with such induction. Messenger-RNA sequencing in whole blood was performed on the peripheral blood of 53 vaccinees. We generated longitudinal data on the peripheral blood of 53 volunteers, at early (days 3 and 7) and late (28 days) time points after priming and boosting with cHAs. Differentially expressed gene analysis showed no differences between placebo and live-attenuated vaccine groups. However, an upregulation of genes involved in innate immune responses and type I interferon signaling was found at day 3 after vaccination with inactivated adjuvanted formulations. Cell type deconvolution analysis revealed a significant enrichment for monocyte markers and different subsets of dendritic cells as mediators for optimal B cell responses and significant increase of anti-stalk antibodies in sera. A significant upregulation of immunoglobulin-related genes was only observed after administration of adjuvanted vaccines (either as primer or booster) with specific induction of anti-stalk IGVH1-69. This approach informed of specific immune signatures that correlate with robust anti-stalk antibody responses, while also helping to understand the regulation of gene expression induced by cHA proteins under different vaccine regimens.

Potent monoclonal antibodies neutralize Omicron sublineages and other SARS-CoV-2 variants

The emergence and global spread of the SARS-CoV-2 Omicron variants, which carry an unprecedented number of mutations, raise serious concerns due to the reduced efficacy of current vaccines and resistance to therapeutic antibodies. Here, we report the generation and characterization of two potent human monoclonal antibodies, NA8 and NE12, against the receptor-binding domain of the SARS-CoV-2 spike protein. NA8 interacts with a highly conserved region and has a breadth of neutralization with picomolar potency against the Beta variant and the Omicron BA.1 and BA.2 sublineages and nanomolar potency against BA.2.12.1 and BA.4. Combination of NA8 and NE12 retains potent neutralizing activity against the major SARS-CoV-2 variants of concern. Cryo-EM analysis provides the structural basis for the broad and complementary neutralizing activity of these two antibodies. We confirm the in vivo protective and therapeutic efficacies of NA8 and NE12 in the hamster model. These results show that broad and potent human antibodies can overcome the continuous immune escape of evolving SARS-CoV-2 variants.

IgTreeZ, A Toolkit for Immunoglobulin Gene Lineage Tree-Based Analysis, Reveals CDR3s Are Crucial for Selection Analysis

Somatic hypermutation (SHM) is an important diversification mechanism that plays a part in the creation of immune memory. Immunoglobulin (Ig) variable region gene lineage trees were used over the last four decades to model SHM and the selection mechanisms operating on B cell clones. We hereby present IgTreeZ (Immunoglobulin Tree analyZer), a python-based tool that analyses many aspects of Ig gene lineage trees and their repertoires. Using simulations, we show that IgTreeZ can be reliably used for mutation and selection analyses. We used IgTreeZ on empirical data, found evidence for different mutation patterns in different B cell subpopulations, and gained insights into antigen-driven selection in corona virus disease 19 (COVID-19) patients. Most importantly, we show that including the CDR3 regions in selection analyses - which is only possible if these analyses are lineage tree-based - is crucial for obtaining correct results. Overall, we present a comprehensive lineage tree analysis tool that can reveal new biological insights into B cell repertoire dynamics.

IMGT ® Nomenclature of Engineered IGHG Variants Involved in Antibody Effector Properties and Formats

The constant region of the immunoglobulin (IG) or antibody heavy gamma chain is frequently engineered to modify the effector properties of the therapeutic monoclonal antibodies. These variants are classified in regards to their effects on effector functions, antibody-dependent cytotoxicity (ADCC), antibody-dependent phagocytosis (ADCP), complement-dependent cytotoxicity (CDC) enhancement or reduction, B cell inhibition by the coengagement of antigen and FcγR on the same cell, on half-life increase, and/or on structure such as prevention of IgG4 half-IG exchange, hexamerisation, knobs-into-holes and the heteropairing H-H of bispecific antibodies, absence of disulfide bridge inter H-L, absence of glycosylation site, and site-specific drug attachment engineered cysteine. The IMGT engineered variant identifier is comprised of the species and gene name (and eventually allele), the letter 'v' followed by a number (assigned chronologically), and for each concerned domain (e.g, CH1, h, CH2 and CH3), the novel AA (single letter abbreviation) and IMGT position according to the IMGT unique numbering for the C-domain and between parentheses, the Eu numbering. IMGT engineered variants are described with detailed amino acid changes, visualized in motifs based on the IMGT numbering bridging genes, sequences, and structures for higher order description.

Improving naive B cell isolation by absence of CD45RB glycosylation and CD27 expression in combination with BCR isotype

In past years ex vivo and in vivo experimental approaches involving human naive B cells have proven fundamental for elucidation of mechanisms promoting B cell differentiation in both health and disease. For such studies, it is paramount that isolation strategies yield a population of bona fide naive B cells, i.e., B cells that are phenotypically and functionally naive, clonally non-expanded, and have non-mutated BCR variable regions. In this study different combinations of common as well as recently identified B cell markers were compared to isolate naive B cells from human peripheral blood. High-throughput BCR sequencing was performed to analyze levels of somatic hypermutation and clonal expansion. Additionally, contamination from mature mutated B cells intrinsic to each cell-sorting strategy was evaluated and how this impacts the purity of obtained populations. Our results show that current naive B cell isolation strategies harbor contamination from non-naive B cells, and use of CD27-IgD⁺ is adequate but can be improved by including markers for CD45RB glycosylation and IgM. The finetuning of naive B cell classification provided herein will harmonize research lines using naive B cells, and will improve B cell profiling during health and disease, e.g. during diagnosis, treatment, and vaccination strategies.

Profiling genes encoding the adaptive immune receptor repertoire with gAIRR Suite

Adaptive immune receptor repertoire (AIRR) is encoded by T cell receptor (TR) and immunoglobulin (IG) genes. Profiling these germline genes encoding AIRR (abbreviated as gAIRR) is important in understanding adaptive immune responses but is challenging due to the high genetic complexity. Our gAIRR Suite comprises three modules. gAIRR-seq, a probe capture-based targeted sequencing pipeline, profiles gAIRR from individual DNA samples. gAIRR-call and gAIRR-annotate call alleles from gAIRR-seq reads and annotate whole-genome assemblies, respectively. We gAIRR-seqed TRV and TRJ of seven Genome in a Bottle (GIAB) DNA samples with 100% accuracy and discovered novel alleles. We also gAIRR-seqed and gAIRR-called the TR and IG genes of a subject from both the peripheral blood mononuclear cells (PBMC) and oral mucosal cells. The calling results from these two cell types have a high concordance (99% for all known gAIRR alleles). We gAIRR-annotated 36 genomes to unearth 325 novel TRV alleles and 29 novel TRJ alleles. We could further profile the flanking sequences, including the recombination signal sequence (RSS). We validated two structural variants for HG002 and uncovered substantial differences of gAIRR genes in references GRCh37 and GRCh38. gAIRR Suite serves as a resource to sequence, analyze, and validate germline TR and IG genes to study various immune-related phenotypes.

A multi-objective based clustering for inferring BCR clonal lineages from high-throughput B cell repertoire data

The adaptive B cell response is driven by the expansion, somatic hypermutation, and selection of B cell clonal lineages. A high number of clonal lineages in a B cell population indicates a highly diverse repertoire, while clonal size distribution and sequence diversity reflect antigen selective pressure. Identifying clonal lineages is fundamental to many repertoire studies, including repertoire comparisons, clonal tracking, and statistical analysis. Several methods have been developed to group sequences from high-throughput B cell repertoire data. Current methods use clustering algorithms to group clonally-related sequences based on their similarities or distances. Such approaches create groups by optimizing a single objective that typically minimizes intra-clonal distances. However, optimizing several objective functions can be advantageous and boost the algorithm convergence rate. Here we propose MobiLLe, a new method based on multi-objective clustering. Our approach requires V(D)J annotations to obtain the initial groups and iteratively applies two objective functions that optimize cohesion and separation within clonal lineages simultaneously. We show that our method greatly improves clonal lineage grouping on simulated benchmarks with varied mutation rates compared to other tools. When applied to experimental repertoires generated from high-throughput sequencing, its clustering results are comparable to the most performing tools and can reproduce the results of previous publications. The method based on multi-objective clustering can accurately identify clonally-related antibody sequences and presents the lowest running time among state-of-art tools. All these features constitute an attractive option for repertoire analysis, particularly in the clinical context. MobiLLe can potentially help unravel the mechanisms involved in developing and evolving B cell malignancies.

High-throughput sequencing can produce a large set of sequences and has profoundly changed our ability to study immune repertoires, particularly B cell receptor sequences. An important application is the analysis of the clonal lineage composition of B cell populations; it is the starting point of many immune repertoire studies, for instance, to differentiate between healthy individuals and those with lymphoid malignancies or other diseases. Several computational methods have been developed to identify clonal lineages from a set of B cell receptor sequences. Most of them apply clustering algorithms and optimize a single objective function that typically minimizes intra-clonal distances. However, optimizing several objective functions in parallel can benefit and increase the clustering performance and efficiency. We propose MobiLLe, the first multi-objective clonal lineage grouping method, which simultaneously optimizes two objective functions for minimizing intra-clonal diversity and maximizing inter-clonal differences. Our approach greatly improved clonal grouping on simulated benchmarks and performed comparably to the most powerful and recent methods on experimental samples. MobiLLe is computationally more efficient than existing tools and does not require any training process or hyper-parameter optimization. It can easily manage large-scale experimental repertoires, providing useful plots to help researchers detect clonally-related sequences in high-throughput B cell repertoire data.

A public antibody class recognizes an S2 epitope exposed on open conformations of SARS-CoV-2 spike

Delineating the origins and properties of antibodies elicited by SARS-CoV-2 infection and vaccination is critical for understanding their benefits and potential shortcomings. Therefore, we investigate the SARS-CoV-2 spike (S)-reactive B cell repertoire in unexposed individuals by flow cytometry and single-cell sequencing. We show that ∼82% of SARS-CoV-2 S-reactive B cells harbor a naive phenotype, which represents an unusually high fraction of total human naive B cells (∼0.1%). Approximately 10% of these naive S-reactive B cells share an IGHV1-69/IGKV3-11 B cell receptor pairing, an enrichment of 18-fold compared to the complete naive repertoire. Following SARS-CoV-2 infection, we report an average 37-fold enrichment of IGHV1-69/IGKV3-11 B cell receptor pairing in the S-reactive memory B cells compared to the unselected memory repertoire. This class of B cells targets a previously undefined non-neutralizing epitope on the S2 subunit that becomes exposed on S proteins used in approved vaccines when they transition away from the native pre-fusion state because of instability. These findings can help guide the improvement of SARS-CoV-2 vaccines.

Longitudinal analysis of anti-drug antibody development in multiple sclerosis patients treated with interferon beta-1a (Rebif™) using B cell receptor repertoire analysis

A significant proportion of multiple sclerosis (MS) patients treated with interferon beta-1a (Rebif™) develop anti-drug antibodies (ADA) with a negative impact on treatment efficacy. We hypothesized that high-throughput B-cell receptor (BCR) repertoire analysis could be used to predict and monitor ADA development. To study this we analyzed 228 peripheral blood samples from 68 longitudinally followed patients starting on interferon beta-1a. Our results show that whole blood BCR analysis does not reflect, and does not predict ADA development in MS patients treated with interferon beta-1a. We propose that BCR analysis of phenotypically selected cell subsets or tissues might be more informative.

Germinal center activity and B cell maturation are associated with protective antibody responses against Plasmodium pre-erythrocytic infection

Blocking Plasmodium, the causative agent of malaria, at the asymptomatic pre-erythrocytic stage would abrogate disease pathology and prevent transmission. However, the lack of well-defined features within vaccine-elicited antibody responses that correlate with protection represents a major roadblock to improving on current generation vaccines. We vaccinated mice (BALB/cJ and C57BL/6J) with Py circumsporozoite protein (CSP), the major surface antigen on the sporozoite, and evaluated vaccine-elicited humoral immunity and identified immunological factors associated with protection after mosquito bite challenge. Vaccination achieved 60% sterile protection and otherwise delayed blood stage patency in BALB/cJ mice. In contrast, all C57BL/6J mice were infected similar to controls. Protection was mediated by antibodies and could be passively transferred from immunized BALB/cJ mice into naïve C57BL/6J. Dissection of the underlying immunological features of protection revealed early deficits in antibody titers and polyclonal avidity in C57BL/6J mice. Additionally, PyCSP-vaccination in BALB/cJ induced a significantly higher proportion of antigen-specific B-cells and class-switched memory B-cell (MBCs) populations than in C57BL/6J mice. Strikingly, C57BL/6J mice also had markedly fewer CSP-specific germinal center experienced B cells and class-switched MBCs compared to BALB/cJ mice. Analysis of the IgG γ chain repertoires by next generation sequencing in PyCSP-specific memory B-cell repertoires also revealed higher somatic hypermutation rates in BALB/cJ mice than in C57BL/6J mice. These findings indicate that the development of protective antibody responses in BALB/cJ mice in response to vaccination with PyCSP was associated with increased germinal center activity and somatic mutation compared to C57BL/6J mice, highlighting the key role B cell maturation may have in the development of vaccine-elicited protective antibodies against CSP.

Identifying specific features of vaccine-elicited antibody responses that are associated with protection from malaria infection is a key step toward the development of a safe and effective vaccine. Here we compared antibody and B cell responses in two mouse strains that exhibited a differential ability to generate antibodies that protect from infection challenge. We found that protection was due to the presence of vaccine-elicited antibodies and could be transferred between strains, and that the ability of antibodies to neutralize the parasite was directly linked to the strength (affinity) with which it binds CSP. Thus, we sought to understand if there were differences in the two strains in the process of B cell maturation that leads to generation of high affinity, protective antibody responses after vaccination. Overall, our comparative analysis indicates that germinal center (GC) activity, a key process in B cell maturation, was significantly diminished in the non-protected strain. Further, we observed evidence of higher levels of somatic mutation, which is a result of germinal center activity, in protected mice. Thus, our results indicate that the ability to generate protective antibody responses was linked to enhanced B cell maturation in the protected strain, providing a key clue to the type of responses that should be generated by future vaccines.

AbLang: an antibody language model for completing antibody sequences

Motivation

General protein language models have been shown to summarize the semantics of protein sequences into representations that are useful for state-of-the-art predictive methods. However, for antibody specific problems, such as restoring residues lost due to sequencing errors, a model trained solely on antibodies may be more powerful. Antibodies are one of the few protein types where the volume of sequence data needed for such language models is available, e.g. in the Observed Antibody Space (OAS) database.

Results

Here, we introduce AbLang, a language model trained on the antibody sequences in the OAS database. We demonstrate the power of AbLang by using it to restore missing residues in antibody sequence data, a key issue with B-cell receptor repertoire sequencing, e.g. over 40% of OAS sequences are missing the first 15 amino acids. AbLang restores the missing residues of antibody sequences better than using IMGT germlines or the general protein language model ESM-1b. Further, AbLang does not require knowledge of the germline of the antibody and is seven times faster than ESM-1b.

Availability and implementation

AbLang is a python package available at https://github.com/oxpig/AbLang.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

IMGT/3Dstructure-DB: T-Cell Receptor TR Paratope and Peptide/Major Histocompatibility pMH Contact Sites and Epitope

T-cell receptors (TR), the antigen receptors of T cells, specifically recognize peptides presented by the major histocompatibility (MH) proteins, as peptide/MH (pMH), on the cell surface. The structure characterization of the trimolecular TR/pMH complexes is crucial to the fields of immunology, vaccination, and immunotherapy. IMGT/3Dstructure-DB is the three-dimensional (3-D) structure database of IMGT^®, the international ImMunoGenetics information system^®. By its creation, IMGT^® marks the advent of immunoinformatics, which emerged at the interface between immunogenetics and bioinformatics. The IMGT^® immunoglobulin (IG) and TR gene and allele nomenclature (CLASSIFICATION axiom) and the IMGT unique numbering and IMGT/Collier-de-Perles (NUMEROTATION axiom) are the two founding breakthroughs of immunoinformatics. IMGT-ONTOLOGY concepts and IMGT Scientific chart rules generated from these axioms allowed IMGT^® bridging genes, structures, and functions. IMGT/3Dstructure-DB contains 3-D structures of IG or antibodies, TR and MH proteins of the adaptive immune responses of jawed vertebrates (gnathostomata), IG or TR complexes with antigens (IG/Ag, TR/pMH), related proteins of the immune system of any species belonging to the IG and MH superfamilies, and fusion proteins for immune applications. The focus of this chapter is on the TR V domains and MH G domains and the contact analysis comparison in TR/pMH interactions. Standardized molecular characterization includes "IMGT pMH contact sites" for peptide and MH groove interactions and "IMGT paratopes and epitopes" for TR/pMH complexes. Data are available in the IMGT/3Dstructure database, at the IMGT Home page http://www.imgt.org .

Bulk Sequencing from mRNA with UMI for Evaluation of B-Cell Isotype and Clonal Evolution: A Method by the AIRR Community

During the course of an immune response to a virus such as influenza, B cells undergo activation, clonal expansion, isotype switching, and somatic hypermutation (SHM). Members of an antigen-experienced B-cell clone can have different sequence features including SHM in the immunoglobulin heavy-chain V (IGHV) gene and can use the same IGVH gene in combination with different constant regions or isotypes (e.g., IgM, IgG, IgA). To study these features of expanded clones in an immune response by AIRR-seq, we provide a bulk RNA-based sequencing experimental procedure with unique molecular identifiers (UMIs) and the accompanying bioinformatics analytical workflow.

Profiling of hMPV F-specific antibodies isolated from human memory B cells

Human metapneumovirus (hMPV) belongs to the Pneumoviridae family and is closely related to respiratory syncytial virus (RSV). The surface fusion (F) glycoprotein mediates viral fusion and is the primary target of neutralizing antibodies against hMPV. Here we report 113 hMPV-F specific monoclonal antibodies (mAbs) isolated from memory B cells of human donors. We characterize the antibodies' germline usage, epitopes, neutralization potencies, and binding specificities. We find that unlike RSV-F specific mAbs, antibody responses to hMPV F are less dominant against the apex of the antigen, and the majority of the potent neutralizing mAbs recognize epitopes on the side of hMPV F. Furthermore, neutralizing epitopes that differ from previously defined antigenic sites on RSV F are identified, and multiple binding modes of site V and II mAbs are discovered. Interestingly, mAbs that bind preferentially to the unprocessed prefusion F show poor neutralization potency. These results elucidate the immune recognition of hMPV infection and provide novel insights for future hMPV antibody and vaccine development.

Sensitive identification of neoantigens and cognate TCRs in human solid tumors

The identification of patient-specific tumor antigens is complicated by the low frequency of T cells specific for each tumor antigen. Here we describe NeoScreen, a method that enables the sensitive identification of rare tumor (neo)antigens and of cognate T cell receptors (TCRs) expressed by tumor-infiltrating lymphocytes. T cells transduced with tumor antigen-specific TCRs identified by NeoScreen mediate regression of established tumors in patient-derived xenograft mice.

Robust, persistent adaptive immune responses to SARS-CoV-2 in the oropharyngeal lymphoid tissue of children

SARS-CoV-2 infection triggers adaptive immune responses from both T and B cells. However, most studies focus on peripheral blood, which may not fully reflect immune responses in lymphoid tissues at the site of infection. To evaluate both local and systemic adaptive immune responses to SARS-CoV-2, we collected peripheral blood, tonsils, and adenoids from 110 children undergoing tonsillectomy/adenoidectomy during the COVID-19 pandemic and found 24 with evidence of prior SARS-CoV-2 infection, including detectable neutralizing antibodies against multiple viral variants. We identified SARS-CoV-2-specific germinal center (GC) and memory B cells; single cell BCR sequencing showed that these virus-specific B cells were class-switched and somatically hypermutated, with overlapping clones in the adenoids and tonsils. Oropharyngeal tissues from COVID-19-convalescent children showed persistent expansion of GC and anti-viral lymphocyte populations associated with an IFN-γ-type response, with particularly prominent changes in the adenoids, as well as evidence of persistent viral RNA in both tonsil and adenoid tissues of many participants. Our results show robust, tissue-specific adaptive immune responses to SARS-CoV-2 in the upper respiratory tract of children weeks to months after acute infection, providing evidence of persistent localized immunity to this respiratory virus.

Combining genotypes and T cell receptor distributions to infer genetic loci determining V(D)J recombination probabilities

Every T cell receptor (TCR) repertoire is shaped by a complex probabilistic tangle of genetically determined biases and immune exposures. T cells combine a random V(D)J recombination process with a selection process to generate highly diverse and functional TCRs. The extent to which an individual’s genetic background is associated with their resulting TCR repertoire diversity has yet to be fully explored. Using a previously published repertoire sequencing dataset paired with high-resolution genome-wide genotyping from a large human cohort, we infer specific genetic loci associated with V(D)J recombination probabilities using genome-wide association inference. We show that V(D)J gene usage profiles are associated with variation in the TCRB locus and, specifically for the functional TCR repertoire, variation in the major histocompatibility complex locus. Further, we identify specific variations in the genes encoding the Artemis protein and the TdT protein to be associated with biasing junctional nucleotide deletion and N-insertion, respectively. These results refine our understanding of genetically-determined TCR repertoire biases by confirming and extending previous studies on the genetic determinants of V(D)J gene usage and providing the first examples of trans genetic variants which are associated with modifying junctional diversity. Together, these insights lay the groundwork for further explorations into how immune responses vary between individuals.

IMGT® Biocuration and Analysis of the Rhesus Monkey IG Loci

The adaptive immune system, along with the innate immune system, are the two main biological processes that protect an organism from pathogens. The adaptive immune system is characterized by the specificity and extreme diversity of its antigen receptors. These antigen receptors are the immunoglobulins (IG) or antibodies of the B cells and the T cell receptors (TR) of the T cells. The IG are proteins that have a dual role in immunity: they recognize antigens and trigger elimination mechanisms, to rid the body of foreign cells. The synthesis of the immunoglobulin heavy and light chains requires gene rearrangements at the DNA level in the IGH, IGK, and IGL loci. The rhesus monkey (Macaca mulatta) is one of the most widely used nonhuman primate species in biomedical research. In this manuscript, we provide a thorough analysis of the three IG loci of the Mmul_10 assembly of rhesus monkey, integrating IMGT previously existing data. Detailed characterization of IG genes includes their localization and position in the loci, the determination of the allele functionality, and the description of the regulatory elements of their promoters as well as the sequences of the conventional recombination signals (RS). This complete annotation of the genomic IG loci of Mmul_10 assembly and the highly detailed IG gene characterization could be used as a model, in additional rhesus monkey assemblies, for the analysis of the IG allelic polymorphism and structural variation, which have been described in rhesus monkeys.

IMGT®Homo sapiens IG and TR Loci, Gene Order, CNV and Haplotypes: New Concepts as a Paradigm for Jawed Vertebrates Genome Assemblies

IMGT®, the international ImMunoGeneTics information system®, created in 1989, by Marie-Paule Lefranc (Université de Montpellier and CNRS), marked the advent of immunoinformatics, a new science which emerged at the interface between immunogenetics and bioinformatics for the study of the adaptive immune responses. IMGT® is based on a standardized nomenclature of the immunoglobulin (IG) and T cell receptor (TR) genes and alleles from fish to humans and on the IMGT unique numbering for the variable (V) and constant (C) domains of the immunoglobulin superfamily (IgSF) of vertebrates and invertebrates, and for the groove (G) domain of the major histocompatibility (MH) and MH superfamily (MhSF) proteins. IMGT® comprises 7 databases, 17 tools and more than 25,000 pages of web resources for sequences, genes and structures, based on the IMGT Scientific chart rules generated from the IMGT-ONTOLOGY axioms and concepts. IMGT® reference directories are used for the analysis of the NGS high-throughput expressed IG and TR repertoires (natural, synthetic and/or bioengineered) and for bridging sequences, two-dimensional (2D) and three-dimensional (3D) structures. This manuscript focuses on the IMGT®Homo sapiens IG and TR loci, gene order, copy number variation (CNV) and haplotypes new concepts, as a paradigm for jawed vertebrates genome assemblies.

The Organization of the Pig T-Cell Receptor γ (TRG) Locus Provides Insights into the Evolutionary Patterns of the TRG Genes across Cetartiodactyla

The domestic pig (Sus scrofa) is a species representative of the Suina, one of the four suborders within Cetartiodactyla. In this paper, we reported our analysis of the pig TRG locus in comparison with the loci of species representative of the Ruminantia, Tylopoda, and Cetacea suborders. The pig TRG genomic structure reiterates the peculiarity of the organization of Cetartiodactyla loci in TRGC “cassettes”, each containing the basic V-J-J-C unit. Eighteen genes arranged in four TRGC cassettes, form the pig TRG locus. All the functional TRG genes were expressed, and the TRGV genes preferentially rearrange with the TRGJ genes within their own cassette, which correlates the diversity of the γ-chain repertoire with the number of cassettes. Among them, the TRGC5, located at the 5′ end of the locus, is the only cassette that retains a marked homology with the corresponding TRGC cassettes of all the analyzed species. The preservation of the TRGC5 cassette for such a long evolutionary time presumes a highly specialized function of its genes, which could be essential for the survival of species. Therefore, the maintenance of this cassette in pigs confirms that it is the most evolutionarily ancient within Cetartiodactyla, and it has undergone a process of duplication to give rise to the other TRGC cassettes in the different artiodactyl species in a lineage-specific manner.

Extremely potent monoclonal antibodies neutralize Omicron and other SARS-CoV-2 variants

The ongoing coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has triggered a devastating global health, social and economic crisis. The RNA nature and broad circulation of this virus facilitate the accumulation of mutations, leading to the continuous emergence of variants of concern with increased transmissibility or pathogenicity ¹ . This poses a major challenge to the effectiveness of current vaccines and therapeutic antibodies ^{1, 2} . Thus, there is an urgent need for effective therapeutic and preventive measures with a broad spectrum of action, especially against variants with an unparalleled number of mutations such as the recently emerged Omicron variant, which is rapidly spreading across the globe ³ . Here, we used combinatorial antibody phage-display libraries from convalescent COVID-19 patients to generate monoclonal antibodies against the receptor-binding domain of the SARS-CoV-2 spike protein with ultrapotent neutralizing activity. One such antibody, NE12, neutralizes an early isolate, the WA-1 strain, as well as the Alpha and Delta variants with half-maximal inhibitory concentrations at picomolar level. A second antibody, NA8, has an unusual breadth of neutralization, with picomolar activity against both the Beta and Omicron variants. The prophylactic and therapeutic efficacy of NE12 and NA8 was confirmed in preclinical studies in the golden Syrian hamster model. Analysis by cryo-EM illustrated the structural basis for the neutralization properties of NE12 and NA8. Potent and broadly neutralizing antibodies against conserved regions of the SARS-CoV-2 spike protein may play a key role against future variants of concern that evade immune control.

Analysis of T-Cell Receptor Repertoire in Transplantation: Fingerprint of T Cell-mediated Alloresponse

T cells play a key role in determining allograft function by mediating allogeneic immune responses to cause rejection, and recent work pointed their role in mediating tolerance in transplantation. The unique T-cell receptor (TCR) expressed on the surface of each T cell determines the antigen specificity of the cell and can be the specific fingerprint for identifying and monitoring. Next-generation sequencing (NGS) techniques provide powerful tools for deep and high-throughput TCR profiling, and facilitate to depict the entire T cell repertoire profile and trace antigen-specific T cells in circulation and local tissues. Tailing T cell transcriptomes and TCR sequences at the single cell level provides a full landscape of alloreactive T-cell clones development and biofunction in alloresponse. Here, we review the recent advances in TCR sequencing techniques and computational tools, as well as the recent discovery in overall TCR profile and antigen-specific T cells tracking in transplantation. We further discuss the challenges and potential of using TCR sequencing-based assays to profile alloreactive TCR repertoire as the fingerprint for immune monitoring and prediction of rejection and tolerance.

Gut microbiota exaggerates triclosan-induced liver injury via gut-liver axis

Triclosan (TCS) is an antimicrobial ingredient that has been widely incorporated in consumer products. TCS can cause hepatic damage by disturbing lipid metabolism, which is often accompanied with gut microbiota dysbiosis. However, the effects of gut microbiota on the TCS-induced liver injury are still unknown. Therefore, we constructed a mouse model based on five-week-old male C57BL/6 mice to investigate the effects of dietary TCS exposure (40 ppm) on liver injury. We found that TCS treatment for 4 weeks dramatically disturbed gut microbiota homeostasis, resulting in overproduction of lipopolysaccharides (LPS) and deficiency of secondary bile acids such as deoxycholic acid (DCA) and lithocholic acid (LCA). In addition, TCS considerably increased intestinal permeability by reducing mucus excretion and expression of tight junction proteins (ZO-1, occludin and claudin 4), which facilitated translocation of LPS. The LPS accumulation in blood contributed to liver injury by triggering the inflammatory response via TLR4 pathway. In summary, this study provides novel insights into the underlying mechanisms of TCS-associated liver injury induced by gut microbiota via the gut-liver axis, and contributes to better interpretation of the health impact of the environmentally emerging contaminant TCS.

IMGT® Immunoinformatics Tools for Standardized V-DOMAIN Analysis

The variable domains (V-DOMAIN) of the antigen receptors, immunoglobulins (IG) or antibodies and T cell receptors (TR), which specifically recognize the antigens show a huge diversity in their sequences. This diversity results from the complex mechanisms involved in the synthesis of these domains at the DNA level (rearrangements of the variable (V), diversity (D), and joining (J) genes; N-diversity; and, for the IG, somatic hypermutations). The recognition of V, D, and J as "genes" and their entry in databases mark the creation of IMGT by Marie-Paule Lefranc, and the origin of immunoinformatics in 1989. For 30 years, IMGT^®, the international ImMunoGeneTics information system^® http://www.imgt.org , has implemented databases and developed tools for IG and TR immunoinformatics, based on the IMGT Scientific chart rules and IMGT-ONTOLOGY concepts and axioms, and more particularly, the princeps ones: IMGT genes and alleles (CLASSIFICATION axiom) and the IMGT unique numbering and IMGT Collier de Perles (NUMEROTATION axiom). This chapter describes the online tools for the characterization and annotation of the expressed V-DOMAIN sequences: (a) IMGT/V-QUEST analyzes in detail IG and TR rearranged nucleotide sequences, (b) IMGT/HighV-QUEST is its high throughput version, which includes a module for the identification of IMGT clonotypes and generates immunoprofiles of expressed V, D, and J genes and alleles, (c) IMGT/StatClonotype performs the pairwise comparison of IMGT/HighV-QUEST immunoprofiles, (d) IMGT/DomainGapAlign analyzes amino acid sequences and is frequently used in antibody engineering and humanization, and (e) IMGT/Collier-de-Perles provides two-dimensional (2D) graphical representations of V-DOMAIN, bridging the gap between sequences and 3D structures. These IMGT® tools are widely used in repertoire analyses of the adaptive immune responses in normal and pathological situations and in the design of engineered IG and TR for therapeutic applications.

OUP accepted manuscript

The study and manipulation of T cell receptors (TCRs) is central to multiple fields across basic and translational immunology research. Produced by V(D)J recombination, TCRs are often only recorded in the literature and data repositories as a combination of their V and J gene symbols, plus their hypervariable CDR3 amino acid sequence. However, numerous applications require full-length coding nucleotide sequences. Here we present Stitchr, a software tool developed to specifically address this limitation. Given minimal V/J/CDR3 information, Stitchr produces complete coding sequences representing a fully spliced TCR cDNA. Due to its modular design, Stitchr can be used for TCR engineering using either published germline or novel/modified variable and constant region sequences. Sequences produced by Stitchr were validated by synthesizing and transducing TCR sequences into Jurkat cells, recapitulating the expected antigen specificity of the parental TCR. Using a companion script, Thimble, we demonstrate that Stitchr can process a million TCRs in under ten minutes using a standard desktop personal computer. By systematizing the production and modification of TCR sequences, we propose that Stitchr will increase the speed, repeatability, and reproducibility of TCR research. Stitchr is available on GitHub.

The human antibody sequence space and structural design of the V, J regions, and CDRH3 with Rosetta

The human adaptive immune response enables the targeting of epitopes on pathogens with high specificity. Infection with a pathogen induces somatic hyper-mutation and B-cell selection processes that govern the shape and diversity of the antibody sequence landscape. To date, even the largest immunome repertoires of adaptive immune receptors acquired by next-generation sequencing cannot fully capture the vast antibody sequence space of a single individual, which is estimated to be at least 1012 potential sequences. Degeneracy of the genetic code means that the number of possible nucleotide triplets (64) is greater than the number of canonical amino acids (20), resulting in some amino acids being encoded by multiple triplets and different amino acids sharing the same nucleotide in 1 or 2 positions in the triplet. We hypothesize that the degeneracy of the genetic code can be used to statistically model an enlarged space of human antibody amino acid sequences, accommodating for the discrepancy between the observed and the hypothesized antibody sequence space. Facilitated by Bayesian statistics and immunome repertoire clustering, we calculated amino acid probabilities from single nucleotide frequencies to infer a human amino acid sequence space that is used to design human-like antibodies with Rosetta. We show that antibodies designed with our restraints are on average up to 16.6% more human-like in the V and J regions compared to the Rosetta designs produced without constraints. The human-likeness of the heavy-chain CDR3 region (CDRH3) could be increased for 8 of 27 antibodies compared to Rosetta designs with a similar number of mutations and could be successfully applied on Mus musculus antibodies to demonstrate humanization.

Adaptive Immune Receptor Repertoire (AIRR) Community Guide to TR and IG Gene Annotation

High-throughput sequencing of adaptive immune receptor repertoires (AIRR, i.e., IG and TR) has revolutionized the ability to carry out large-scale experiments to study the adaptive immune response. Since the method was first introduced in 2009, AIRR sequencing (AIRR-Seq) has been applied to survey the immune state of individuals, identify antigen-specific or immune-state-associated signatures of immune responses, study the development of the antibody immune response, and guide the development of vaccines and antibody therapies. Recent advancements in the technology include sequencing at the single-cell level and in parallel with gene expression, which allows the introduction of multi-omics approaches to understand in detail the adaptive immune response. Analyzing AIRR-seq data can prove challenging even with high-quality sequencing, in part due to the many steps involved and the need to parameterize each step. In this chapter, we outline key factors to consider when preprocessing raw AIRR-Seq data and annotating the genetic origins of the rearranged receptors. We also highlight a number of common difficulties with common AIRR-seq data processing and provide strategies to address them.

Clinical response to nivolumab in an INI1-deficient pediatric chordoma correlates with immunogenic recognition of brachyury

Poorly differentiated chordoma (PDC) is a recently recognized subtype of chordoma characterized by expression of the embryonic transcription factor, brachyury, and loss of INI1. PDC primarily affects children and is associated with a poor prognosis and limited treatment options. Here we describe the molecular and immune tumour microenvironment profiles of two paediatric PDCs produced using whole-genome, transcriptome and whole-genome bisulfite sequencing (WGBS) and multiplex immunohistochemistry. Our analyses revealed the presence of tumour-associated immune cells, including CD8+ T cells, and expression of the immune checkpoint protein, PD-L1, in both patient samples. Molecular profiling provided the rationale for immune checkpoint inhibitor (ICI) therapy, which resulted in a clinical and radiographic response. A dominant T cell receptor (TCR) clone specific for a brachyury peptide-MHC complex was identified from bulk RNA sequencing, suggesting that targeting of the brachyury tumour antigen by tumour-associated T cells may underlie this clinical response to ICI. Correlative analysis with rhabdoid tumours, another INI1-deficient paediatric malignancy, suggests that a subset of tumours may share common immune phenotypes, indicating the potential for a therapeutically targetable subgroup of challenging paediatric cancers.

IMGT® databases, related tools and web resources through three main axes of research and development

IMGT^®, the international ImMunoGeneTics information system^®, http://www.imgt.org/, is at the forefront of the immunogenetics and immunoinformatics fields with more than 30 years of experience. IMGT^® makes available databases and tools to the scientific community pertaining to the adaptive immune response, based on the IMGT-ONTOLOGY. We focus on the recent features of the IMGT^® databases, tools, reference directories and web resources, within the three main axes of IMGT^® research and development. Axis I consists in understanding the adaptive immune response, by deciphering the identification and characterization of the immunoglobulin (IG) and T cell receptor (TR) genes in jawed vertebrates. It is the starting point of the two other axes, namely the analysis and exploration of the expressed IG and TR repertoires based on comparison with IMGT reference directories in normal and pathological situations (Axis II) and the analysis of amino acid changes and functions of 2D and 3D structures of antibody and TR engineering (Axis III).

Individualized VDJ recombination predisposes the available Ig sequence space

The process of recombination between variable (V), diversity (D), and joining (J) immunoglobulin (Ig) gene segments determines an individual's naive Ig repertoire and, consequently, (auto)antigen recognition. VDJ recombination follows probabilistic rules that can be modeled statistically. So far, it remains unknown whether VDJ recombination rules differ between individuals. If these rules differed, identical (auto)antigen-specific Ig sequences would be generated with individual-specific probabilities, signifying that the available Ig sequence space is individual specific. We devised a sensitivity-tested distance measure that enables inter-individual comparison of VDJ recombination models. We discovered, accounting for several sources of noise as well as allelic variation in Ig sequencing data, that not only unrelated individuals but also human monozygotic twins and even inbred mice possess statistically distinguishable immunoglobulin recombination models. This suggests that, in addition to genetic, there is also nongenetic modulation of VDJ recombination. We demonstrate that population-wide individualized VDJ recombination can result in orders of magnitude of difference in the probability to generate (auto)antigen-specific Ig sequences. Our findings have implications for immune receptor-based individualized medicine approaches relevant to vaccination, infection, and autoimmunity.

Metabolite and thymocyte development defects in ADA-SCID mice receiving enzyme replacement therapy

Deficiency of adenosine deaminase (ADA, EC3.5.4.4), a housekeeping enzyme intrinsic to the purine salvage pathway, leads to severe combined immunodeficiency (SCID) both in humans and mice. Lack of ADA results in the intracellular accumulation of toxic metabolites which have effects on T cell development and function. While untreated ADA-SCID is a fatal disorder, there are different therapeutic options available to restore ADA activity and reconstitute a functioning immune system, including enzyme replacement therapy (ERT). Administration of ERT in the form of pegylated bovine ADA (PEG-ADA) has proved a life-saving though non-curative treatment for ADA-SCID patients. However, in many patients treated with PEG-ADA, there is suboptimal immune recovery with low T and B cell numbers. Here, we show reduced thymus cellularity in ADA-SCID mice despite weekly PEG-ADA treatment. This was associated with lack of effective adenosine (Ado) detoxification in the thymus. We also show that thymocyte development in ADA-deficient thymi is arrested at the DN3-to-DN4 stage transition with thymocytes undergoing dATP-induced apoptosis rather than defective TCRβ rearrangement or β-selection. Our studies demonstrate at a detailed level that exogenous once-a-week enzyme replacement does not fully correct intra-thymic metabolic or immunological abnormalities associated with ADA deficiency.

Human B cell lineages associated with germinal centers following influenza vaccination are measurably evolving

The poor efficacy of seasonal influenza virus vaccines is often attributed to pre-existing immunity interfering with the persistence and maturation of vaccine-induced B cell responses. We previously showed that a subset of vaccine-induced B cell lineages are recruited into germinal centers (GCs) following vaccination, suggesting that affinity maturation of these lineages against vaccine antigens can occur. However, it remains to be determined whether seasonal influenza vaccination stimulates additional evolution of vaccine-specific lineages, and previous work has found no significant increase in somatic hypermutation among influenza-binding lineages sampled from the blood following seasonal vaccination in humans. Here, we investigate this issue using a phylogenetic test of measurable immunoglobulin sequence evolution. We first validate this test through simulations and survey measurable evolution across multiple conditions. We find significant heterogeneity in measurable B cell evolution across conditions, with enrichment in primary response conditions such as HIV infection and early childhood development. We then show that measurable evolution following influenza vaccination is highly compartmentalized: while lineages in the blood are rarely measurably evolving following influenza vaccination, lineages containing GC B cells are frequently measurably evolving. Many of these lineages appear to derive from memory B cells. We conclude from these findings that seasonal influenza virus vaccination can stimulate additional evolution of responding B cell lineages, and imply that the poor efficacy of seasonal influenza vaccination is not due to a complete inhibition of vaccine-specific B cell evolution.

Germinal centre-driven maturation of B cell response to SARS-CoV-2 vaccination

Germinal centres (GC) are lymphoid structures where vaccine-responding B cells acquire affinity-enhancing somatic hypermutations (SHM), with surviving clones differentiating into memory B cells (MBCs) and long-lived bone marrow plasma cells (BMPCs) 1-4 . Induction of the latter is a hallmark of durable immunity after vaccination 5 . SARS-CoV-2 mRNA vaccination induces a robust GC response in humans 6-8 , but the maturation dynamics of GC B cells and propagation of their progeny throughout the B cell diaspora have not been elucidated. Here we show that anti-SARS-CoV-2 spike (S)-binding GC B cells were detectable in draining lymph nodes for at least six months in 10 out of 15 individuals who had received two doses of BNT162b2, a SARS-CoV-2 mRNA vaccine. Six months after vaccination, circulating S-binding MBCs were detected in all participants (n=42) and S-specific IgG-secreting BMPCs were detected in 9 out of 11 participants. Using a combined approach of single-cell RNA sequencing of responding blood and lymph node B cells from eight participants and expression of the corresponding monoclonal antibodies, we tracked the evolution of 1540 S-specific B cell clones. SHM accumulated along the B cell differentiation trajectory, with early blood plasmablasts showing the lowest frequencies, followed by MBCs and lymph node plasma cells whose SHM largely overlapped with GC B cells. By three months after vaccination, the frequency of SHM within GC B cells had doubled. Strikingly, S + BMPCs detected six months after vaccination accumulated the highest level of SHM, corresponding with significantly enhanced anti-S polyclonal antibody avidity in blood at that time point. This study documents the induction of affinity-matured BMPCs after two doses of SARS-CoV-2 mRNA vaccination in humans, providing a foundation for the sustained high efficacy observed with these vaccines.