IMGT®, the international ImMunoGeneTics information system® 25 years on

Susceptibility to Plasmodium falciparum Malaria: Influence of Combined Polymorphisms of IgG3 Gm Allotypes and Fc Gamma Receptors IIA, IIIA, and IIIB

The binding of immunoglobulin (Ig) to Fc gamma receptors (FcgR) at the immune cell surface is an important step to initiate immunological defense against malaria. However, polymorphisms in receptors and/or constant regions of the IgG heavy chains may modulate this binding. Here, we investigated whether polymorphisms located in FcgR and constant regions of the heavy chain of IgG are associated with susceptibility to P. falciparum malaria. For this purpose, a clinical and parasitological follow-up on malaria was conducted among 656 infants in southern Benin. G3m allotypes (from total IgG3) were determined by a serological method of hemagglutination inhibition. FcgRIIA 131R/H and FcgRIIIA 176F/V genotypes were determined using the TaqMan method and FcgRIIIB NA1/NA2 genotypes were assessed by polymerase chain reaction using allele-specific primers. Association analyses between the number of malaria infections during the follow-up and polymorphisms in IgG G3m allotypes and FcgR were studied independently by zero inflated binomial negative regression. The influence of combinations of G3m allotypes and FcgRIIA/FcgRIIIA/FcgRIIIB polymorphisms on the number of P. falciparum infections, and their potential interaction with environmental exposure to malaria was assessed by using the generalized multifactor dimensionality reduction (GMDR) method. Results showed that individual carriage of G3m24 single allotype and of G3m5,6,10,11,13,14,24 phenotype was independently associated with a high risk of malaria infection. A risk effect for G3m6 was observed only under high environmental exposure. FcgRIIIA 176VV single genotype and combined carriage of FcgRIIA 131RH/FcgRIIIA 176VV/FcgRIIIB NA1NA2, FcgRIIA 131HH/FcgRIIIA 176FF/FcgRIIIB NA1NA1, FcgRIIA 131HH/FcgRIIIA 176VV/FcgRIIIB NA2NA2 and FcgRIIA 131HH/FcgRIIIA 176VV/FcgRIIIB NA1NA2 genotypes were related to a high number of malaria infections. The risk was accentuated for FcgRIIIA 176VV when considering the influence of environmental exposure to malaria. Finally, the GMDR analysis including environmental exposure showed strengthened associations with a malaria risk when FcgRIIA/FcgRIIIA/FcgRIIIB genotypes were combined to G3m5,6,11,24 and G3m5,6,10,11,13,15,24 phenotypes or G3m10 and G3m13 single allotypes. Our results highlight the relevance of studying IgG heavy chain and FcgR polymorphisms, independently as well as in combination, in relation to the individual susceptibility to P. falciparum infection. The intensity of individual exposure to mosquito bites was demonstrated to impact the relationships found.

The adaptive immune receptor repertoire community as a model for FAIR stewardship of big immunology data

Systems biology involves network-oriented, computational approaches to modeling biological systems through analysis of big biological data. To contribute maximally to scientific progress, big biological data should be FAIR: findable, accessible, interoperable, and reusable. Here, we describe high-throughput sequencing data that characterize the vast diversity of B- and T-cell clones comprising the adaptive immune receptor repertoire (AIRR-seq data) and its contribution to our understanding of COVID-19 (coronavirus disease 19). We describe the accomplishments of the AIRR community, a grass-roots network of interdisciplinary laboratory scientists, bioinformaticians, and policy wonks, in creating and publishing standards, software and repositories for AIRR-seq data based on the FAIR principles.

SARS-CoV-2 neutralizing antibody structures inform therapeutic strategies

The coronavirus disease 2019 (COVID-19) pandemic presents an urgent health crisis. Human neutralizing antibodies that target the host ACE2 receptor-binding domain (RBD) of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) spike protein1-5 show promise therapeutically and are being evaluated clinically6-8. Here, to identify the structural correlates of SARS-CoV-2 neutralization, we solved eight new structures of distinct COVID-19 human neutralizing antibodies5 in complex with the SARS-CoV-2 spike trimer or RBD. Structural comparisons allowed us to classify the antibodies into categories: (1) neutralizing antibodies encoded by the VH3-53 gene segment with short CDRH3 loops that block ACE2 and bind only to 'up' RBDs; (2) ACE2-blocking neutralizing antibodies that bind both up and 'down' RBDs and can contact adjacent RBDs; (3) neutralizing antibodies that bind outside the ACE2 site and recognize both up and down RBDs; and (4) previously described antibodies that do not block ACE2 and bind only to up RBDs9. Class 2 contained four neutralizing antibodies with epitopes that bridged RBDs, including a VH3-53 antibody that used a long CDRH3 with a hydrophobic tip to bridge between adjacent down RBDs, thereby locking the spike into a closed conformation. Epitope and paratope mapping revealed few interactions with host-derived N-glycans and minor contributions of antibody somatic hypermutations to epitope contacts. Affinity measurements and mapping of naturally occurring and in vitro-selected spike mutants in 3D provided insight into the potential for SARS-CoV-2 to escape from antibodies elicited during infection or delivered therapeutically. These classifications and structural analyses provide rules for assigning current and future human RBD-targeting antibodies into classes, evaluating avidity effects and suggesting combinations for clinical use, and provide insight into immune responses against SARS-CoV-2.

The respiratory syncytial virus fusion protein-specific B cell receptor repertoire reshaped by post-fusion subunit vaccination

Respiratory syncytial virus (RSV) is the major cause of acute lower respiratory illness in children of less than 5 years of age which usually results in hospitalization or even in death. Vaccine development is hampered in consequence of a failed vaccine trial with fatalities in the 1960s. Even though research has been more focused on the RSV fusion protein in its pre-fusion conformation, maternal vaccination with post-fusion protein (post F) was considered as a promising vaccine strategy for passive immunization of babies, because post F preserves very potent neutralizing epitopes. We extensively analyzed post F-binding B cell receptor (BCR) repertoires of three vaccinees who received a post F-subunit vaccine in the context of a first-in-human, Phase 1, randomized, observer-blind, placebo-controlled clinical trial (ClinicalTrials.gov Identifier: NCT02298179). In order to compare the vaccine-induced BCR repertoires with BCR repertoires induced by natural infection, we also analyzed pre F- and post F-binding BCRs isolated from a healthy blood donor with relatively high F-binding memory B cell (MBC) frequencies. Analysis of the vaccine-induced repertoires revealed that preferentially V_H4-encoded BCRs were expanded in response to vaccination. Estimation of antigen-driven selection further demonstrated that expanded BCRs accumulated positively selected replacement mutations which substantiated the hypothesis that post F-vaccination induces diversification of V_H4-encoded BCRs in germinal centers. Comparison of the vaccine-induced BCR repertoires with clonally related pre and post F-binding BCRs of the healthy blood donor suggested that the vaccine expanded pre/post F cross-reactive MBCs. Interestingly, several vaccine-induced BCRs shared stereotypic VDJ gene junctions with known neutralizing Abs. Once expressed for functional characterization, the selected monoclonal Abs demonstrated the predicted neutralization activities in plaque reduction neutralization assays indicating that the post F-vaccine induced expansion of neutralizing BCRs.

Dataset of antibody variable region sequence features inferred from a respiratory syncytial virus fusion protein-specific B cell receptor repertoire induced by natural infection of a healthy adult

Respiratory syncytial virus (RSV) is the primary cause for acute lower respiratory syndrome in children younger than 5 years. Research on B cell repertoires and antibodies binding the RSV fusion protein (RSV F) is of major interest in the development of potential vaccine candidates and therapies. B cell receptors (BCRs) which have higher affinities for a specific antigen are preferentially selected for B cell clonal expansion in germinal center reactions. Consequently, antigen-specific BCR repertoires share common features, as for instance preferential variable gene usage, variable region mutation levels or lengths of the heavy chain complementarity-determining region 3. Since RSV repeatedly infects every person throughout life, memory B cells (MBC) expressing RSV F-binding BCRs circulate in the blood of healthy adults. This dataset of BCR variable region sequence features was derived from single cell-sorted RSV F-directed MBCs of a healthy adult blood donor [1]. The dataset was produced with publicly available data analysis software programs and scripts, which facilitates integration or comparison with antibody sequence repertoire data of different individuals derived with the same or comparable data analysis approaches and tools.

Gut γδ T cells as guardians, disruptors, and instigators of cancer

Colorectal cancer is the third most common cancer worldwide with nearly 2 million cases per year. Immune cells and inflammation are a critical component of colorectal cancer progression, and they are used as reliable prognostic indicators of patient outcome. With the growing appreciation for immunology in colorectal cancer, interest is growing on the role γδ T cells have to play, as they represent one of the most prominent immune cell populations in gut tissue. This group of cells consists of both resident populations-γδ intraepithelial lymphocytes (γδ IELs)-and transient populations that each has unique functions. The homeostatic role of these γδ T cell subsets is to maintain barrier integrity and prevent microorganisms from breaching the mucosal layer, which is accomplished through crosstalk with enterocytes and other immune cells. Recent years have seen a surge in discoveries regarding the regulation of γδ IELs in the intestine and the colon with particular new insights into the butyrophilin family. In this review, we discuss the development, specialities, and functions of γδ T cell subsets during cancer progression. We discuss how these cells may be used to predict patient outcome, as well as how to exploit their behavior for cancer immunotherapy.

TRIP - T cell receptor/immunoglobulin profiler

Background

Antigen receptors are characterized by an extreme diversity of specificities, which poses major computational and analytical challenges, particularly in the era of high-throughput immunoprofiling by next generation sequencing (NGS). The T cell Receptor/Immunoglobulin Profiler (TRIP) tool offers the opportunity for an in-depth analysis based on the processing of the output files of the IMGT/HighV-Quest tool, a standard in NGS immunoprofiling, through a number of interoperable modules. These provide detailed information about antigen receptor gene rearrangements, including variable (V), diversity (D) and joining (J) gene usage, CDR3 amino acid and nucleotide composition and clonality of both T cell receptors (TR) and B cell receptor immunoglobulins (BcR IG), and characteristics of the somatic hypermutation within the BcR IG genes. TRIP is a web application implemented in R shiny.

Results

Two sets of experiments have been performed in order to evaluate the efficiency and performance of the TRIP tool. The first used a number of synthetic datasets, ranging from 250k to 1M sequences, and established the linear response time of the tool (about 6 h for 1M sequences processed through the entire BcR IG data pipeline). The reproducibility of the tool was tested comparing the results produced by the main TRIP workflow with the results from a previous pipeline used on the Galaxy platform. As expected, no significant differences were noted between the two tools; although the preselection process seems to be stricter within the TRIP pipeline, about 0.1% more rearrangements were filtered out, with no impact on the final results.

Conclusions

TRIP is a software framework that provides analytical services on antigen receptor gene sequence data. It is accurate and contains functions for data wrangling, cleaning, analysis and visualization, enabling the user to build a pipeline tailored to their needs. TRIP is publicly available at https://bio.tools/TRIP_-_T-cell_Receptor_Immunoglobulin_Profiler.

Computational approaches to therapeutic antibody design: established methods and emerging trends

Antibodies are proteins that recognize the molecular surfaces of potentially noxious molecules to mount an adaptive immune response or, in the case of autoimmune diseases, molecules that are part of healthy cells and tissues. Due to their binding versatility, antibodies are currently the largest class of biotherapeutics, with five monoclonal antibodies ranked in the top 10 blockbuster drugs. Computational advances in protein modelling and design can have a tangible impact on antibody-based therapeutic development. Antibody-specific computational protocols currently benefit from an increasing volume of data provided by next generation sequencing and application to related drug modalities based on traditional antibodies, such as nanobodies. Here we present a structured overview of available databases, methods and emerging trends in computational antibody analysis and contextualize them towards the engineering of candidate antibody therapeutics.

The expansion of the TRB and TRG genes in domestic goats (Capra hircus) is characteristic of the ruminant species

BACKGROUND

Goats (Capra hircus), one of the first domesticated species, are economically important for milk and meat production, and their broad geographical distribution reflects their successful adaptation to diverse environmental conditions. Despite the relevance of this species, the genetic research on the goat traits is limited compared to other domestic species. Thanks to the latest goat reference genomic sequence (ARS1), which is considered to be one of the most continuous assemblies in livestock, we deduced the genomic structure of the T cell receptor beta (TRB) and gamma (TRG) loci in this ruminant species.

RESULTS

Our analyses revealed that although the organization of the goat TRB locus is broadly similar to that of the other artiodactyl species, with three in-tandem D-J-C clusters located at the 3' end, a complex and extensive series of duplications have occurred in the V genes at the 5' end, leading to a marked expansion in the number of the TRBV genes. This phenomenon appears to be a feature of the ruminant lineage since similar gene expansions have also occurred in sheep and cattle. Likewise, the general organization of the goat TRG genes is typical of ruminant species studied so far, with two paralogous TRG loci, TRG1 and TRG2, located in two distinct and distant positions on the same chromosome as result of a split in the ancestral locus. Each TRG locus consists of reiterated V-J-J-C cassettes, with the goat TRG2 containing an additional cassette relative to the corresponding sheep and cattle loci.

CONCLUSIONS

Taken together, these findings demonstrate that strong evolutionary pressures in the ruminant lineage have selected for the development of enlarged sets of TRB and TRG genes that contribute to a diverse T cell receptor repertoire. However, differences observed among the goat, sheep and cattle TRB and TRG genes indicate that distinct evolutionary histories, with independent expansions and/or contractions, have also affected each ruminant species.

Activated PI3Kδ breaches multiple B cell tolerance checkpoints and causes autoantibody production

In patients, gain-of-function (GOF) mutations in PIK3CD break tolerance, causing highly penetrant secretion of autoreactive IgM. Mouse models reveal that Pik3cd GOF subverts the response to self-antigen, preventing the induction of anergy and instead stimulating plasmablast and GC formation.

Antibody-mediated autoimmune diseases are a major health burden. However, our understanding of how self-reactive B cells escape self-tolerance checkpoints to secrete pathogenic autoantibodies remains incomplete. Here, we demonstrate that patients with monogenic immune dysregulation caused by gain-of-function mutations in PIK3CD, encoding the p110δ catalytic subunit of phosphoinositide 3-kinase (PI3K), have highly penetrant secretion of autoreactive IgM antibodies. In mice with the corresponding heterozygous Pik3cd activating mutation, self-reactive B cells exhibit a cell-autonomous subversion of their response to self-antigen: instead of becoming tolerized and repressed from secreting autoantibody, Pik3cd gain-of-function B cells are activated by self-antigen to form plasmablasts that secrete high titers of germline-encoded IgM autoantibody and hypermutating germinal center B cells. However, within the germinal center, peripheral tolerance was still enforced, and there was selection against B cells with high affinity for self-antigen. These data show that the strength of PI3K signaling is a key regulator of pregerminal center B cell self-tolerance and thus represents a druggable pathway to treat antibody-mediated autoimmunity.

Progress, Challenges, and Surprises in Annotating the Human Genome

Our understanding of the human genome has continuously expanded since its draft publication in 2001. Over the years, novel assays have allowed us to progressively overlay layers of knowledge above the raw sequence of A's, T's, G's, and C's. The reference human genome sequence is now a complex knowledge base maintained under the shared stewardship of multiple specialist communities. Its complexity stems from the fact that it is simultaneously a template for transcription, a record of evolution, a vehicle for genetics, and a functional molecule. In short, the human genome serves as a frame of reference at the intersection of a diversity of scientific fields. In recent years, the progressive fall in sequencing costs has given increasing importance to the quality of the human reference genome, as hundreds of thousands of individuals are being sequenced yearly, often for clinical applications. Also, novel sequencing-based assays shed light on novel functions of the genome, especially with respect to gene expression regulation. Keeping the human genome annotation up to date and accurate is therefore an ongoing partnership between reference annotation projects and the greater community worldwide.

Immunoglobulins or Antibodies: IMGT® Bridging Genes, Structures and Functions

IMGT®, the international ImMunoGeneTics® information system founded in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS), marked the advent of immunoinformatics, a new science at the interface between immunogenetics and bioinformatics. For the first time, the immunoglobulin (IG) or antibody and T cell receptor (TR) genes were officially recognized as 'genes' as well as were conventional genes. This major breakthrough has allowed the entry, in genomic databases, of the IG and TR variable (V), diversity (D) and joining (J) genes and alleles of Homo sapiens and of other jawed vertebrate species, based on the CLASSIFICATION axiom. The second major breakthrough has been the IMGT unique numbering and the IMGT Collier de Perles for the V and constant (C) domains of the IG and TR and other proteins of the IG superfamily (IgSF), based on the NUMEROTATION axiom. IMGT-ONTOLOGY axioms and concepts bridge genes, sequences, structures and functions, between biological and computational spheres in the IMGT® system (Web resources, databases and tools). They provide the IMGT Scientific chart rules to identify, to describe and to analyse the IG complex molecular data, the huge diversity of repertoires, the genetic (alleles, allotypes, CNV) polymorphisms, the IG dual function (paratope/epitope, effector properties), the antibody humanization and engineering.

Structural classification of neutralizing antibodies against the SARS-CoV-2 spike receptor-binding domain suggests vaccine and therapeutic strategies

The COVID-19 pandemic presents an urgent health crisis. Human neutralizing antibodies (hNAbs) that target the host ACE2 receptor-binding domain (RBD) of the SARS-CoV-2 spike1-5 show therapeutic promise and are being evaluated clincally6-8. To determine structural correlates of SARS-CoV-2 neutralization, we solved 8 new structures of distinct COVID-19 hNAbs5 in complex with SARS-CoV-2 spike trimer or RBD. Structural comparisons allowed classification into categories: (1) VH3-53 hNAbs with short CDRH3s that block ACE2 and bind only to "up" RBDs, (2) ACE2-blocking hNAbs that bind both "up" and "down" RBDs and can contact adjacent RBDs, (3) hNAbs that bind outside the ACE2 site and recognize "up" and "down" RBDs, and (4) Previously-described antibodies that do not block ACE2 and bind only "up" RBDs9. Class 2 comprised four hNAbs whose epitopes bridged RBDs, including a VH3-53 hNAb that used a long CDRH3 with a hydrophobic tip to bridge between adjacent "down" RBDs, thereby locking spike into a closed conformation. Epitope/paratope mapping revealed few interactions with host-derived N-glycans and minor contributions of antibody somatic hypermutations to epitope contacts. Affinity measurements and mapping of naturally-occurring and in vitro-selected spike mutants in 3D provided insight into the potential for SARS-CoV-2 escape from antibodies elicited during infection or delivered therapeutically. These classifications and structural analyses provide rules for assigning current and future human RBD-targeting antibodies into classes, evaluating avidity effects, suggesting combinations for clinical use, and providing insight into immune responses against SARS-CoV-2.

Structures of Human Antibodies Bound to SARS-CoV-2 Spike Reveal Common Epitopes and Recurrent Features of Antibodies

Neutralizing antibody responses to coronaviruses mainly target the receptor-binding domain (RBD) of the trimeric spike. Here, we characterized polyclonal immunoglobulin Gs (IgGs) and Fabs from COVID-19 convalescent individuals for recognition of coronavirus spikes. Plasma IgGs differed in their focus on RBD epitopes, recognition of alpha- and beta-coronaviruses, and contributions of avidity to increased binding/neutralization of IgGs over Fabs. Using electron microscopy, we examined specificities of polyclonal plasma Fabs, revealing recognition of both S1A and RBD epitopes on SARS-CoV-2 spike. Moreover, a 3.4 Å cryo-electron microscopy (cryo-EM) structure of a neutralizing monoclonal Fab-spike complex revealed an epitope that blocks ACE2 receptor binding. Modeling based on these structures suggested different potentials for inter-spike crosslinking by IgGs on viruses, and characterized IgGs would not be affected by identified SARS-CoV-2 spike mutations. Overall, our studies structurally define a recurrent anti-SARS-CoV-2 antibody class derived from VH3-53/VH3-66 and similarity to a SARS-CoV VH3-30 antibody, providing criteria for evaluating vaccine-elicited antibodies.

De novo prediction of cancer-associated T cell receptors for noninvasive cancer detection

The adaptive immune system recognizes tumor antigens at an early stage to eradicate cancer cells. This process is accompanied by systemic proliferation of the tumor antigen-specific T lymphocytes. While detection of asymptomatic early-stage cancers is challenging due to small tumor size and limited somatic alterations, tracking peripheral T cell repertoire changes may provide an attractive solution to cancer diagnosis. Here, we developed a deep learning method called DeepCAT to enable de novo prediction of cancer-associated T cell receptors (TCRs). We validated DeepCAT using cancer-specific or non-cancer TCRs obtained from multiple major histocompatibility complex I (MHC-I) multimer-sorting experiments and demonstrated its prediction power for TCRs specific to cancer antigens. We blindly applied DeepCAT to distinguish over 250 patients with cancer from over 600 healthy individuals using blood TCR sequences and observed high prediction accuracy, with area under the curve (AUC) ≥ 0.95 for multiple early-stage cancers. This work sets the stage for using the peripheral blood TCR repertoire for noninvasive cancer detection.

Beyond Allotypes: The Influence of Allelic Diversity in Antibody Constant Domains

Polymorphic diversity in antibody constant domains has long been defined by allotypic motifs that cross react with the sera of other individuals. Improvements in sequencing technologies have led to the discovery of a large number of new allelic sequences that underlie this diversity. Many of the point mutations lie outside traditional allotypic motifs suggesting they do not elicit immunogenic responses. As antibodies play an important role in immune defense and biotechnology, understanding how this newly resolved diversity influences the function of antibodies is important. This review investigates the current known diversity of antibody alleles at a protein level for each antibody isotype as well as the kappa and lambda light chains. We focus on evidence emerging for how these mutations perturb antibody interactions with antigens and Fc receptors that are critical for function, as well as the influence this might have on the use of antibodies as therapeutics and reagents.

A minority of T cells recognizing tumor-associated antigens presented in self-HLA can provoke antitumor reactivity

Tumor-associated antigens (TAAs) are monomorphic self-antigens that are proposed as targets for immunotherapeutic approaches to treat malignancies. We investigated whether T cells with sufficient avidity to recognize naturally overexpressed self-antigens in the context of self-HLA can be found in the T-cell repertoire of healthy donors. Minor histocompatibility antigen (MiHA)-specific T cells were used as a model, as the influence of thymic selection on the T-cell repertoire directed against MiHA can be studied in both self (MiHApos donors) and non-self (MiHAneg donors) backgrounds. T-cell clones directed against the HLA*02:01-restricted MiHA HA-1H were isolated from HA-1Hneg/HLA-A*02:01pos and HA-1Hpos/HLA-A*02:01pos donors. Of the 16 unique HA-1H–specific T-cell clones, five T-cell clones derived from HA-1Hneg/HLA-A*02:01pos donors and one T-cell clone derived from an HA-1Hpos/HLA-A*02:01pos donor showed reactivity against HA-1Hpos target cells. In addition, in total, 663 T-cell clones (containing at least 91 unique clones expressing different T-cell receptors) directed against HLA*02:01-restricted peptides of TAA WT1-RMF, RHAMM-ILS, proteinase-3-VLQ, PRAME-VLD, and NY-eso-1-SLL were isolated from HLA-A*02:01pos donors. Only 3 PRAME-VLD–specific and one NY-eso-1-SLL–specific T-cell clone provoked interferon-γ production and/or cytolysis upon stimulation with HLA-A*02:01pos malignant cell lines (but not primary malignant samples) naturally overexpressing the TAA. These results show that self-HLA–restricted T cells specific for self-antigens such as MiHA in MiHApos donors and TAAs are present in peripheral blood of healthy individuals. However, clinical efficacy would require highly effective in vivo priming by peptide vaccination in the presence of proper adjuvants or in vitro expansion of the low numbers of self-antigen–specific T cells of sufficient avidity to recognize endogenously processed antigen.

Specificity of the T Cell Response to Protein Biopharmaceuticals

The anti-drug antibody (ADA) response is an undesired humoral response raised against protein biopharmaceuticals (BPs) which can dramatically disturb their therapeutic properties. One particularity of the ADA response resides in the nature of the immunogens, which are usually human(ized) proteins and are therefore expected to be tolerated. CD4 T cells initiate, maintain and regulate the ADA response and are therefore key players of this immune response. Over the last decade, advances have been made in characterizing the T cell responses developed by patients treated with BPs. Epitope specificity and phenotypes of BP-specific T cells have been reported and highlight the effector and regulatory roles of T cells in the ADA response. BP-specific T cell responses are assessed in healthy subjects to anticipate the immunogenicity of BP prior to their testing in clinical trials. Immunogenicity prediction, also called preclinical immunogenicity assessment, aims at identifying immunogenic BPs and immunogenic BP sequences before any BP injection in humans. All of the approaches that have been developed to date rely on the detection of BP-specific T cells in donors who have never been exposed to BPs. The number of BP-specific T cells circulating in the blood of these donors is therefore limited. T cell assays using cells collected from healthy donors might reveal the weak tolerance induced by BPs, whose endogenous form is expressed at a low level. These BPs have a complete human sequence, but the level of their endogenous form appears insufficient to promote the negative selection of autoreactive T cell clones. Multiple T cell epitopes have also been identified in therapeutic antibodies and some other BPs. The pattern of identified T cell epitopes differs across the antibodies, notwithstanding their humanized, human or chimeric nature. However, in all antibodies, the non-germline amino acid sequences mainly found in the CDRs appear to be the main driver of immunogenicity, provided they can be presented by HLA class II molecules. Considering the fact that the BP field is expanding to include new formats and gene and cell therapies, we face new challenges in understanding and mastering the immunogenicity of new biological products.

PyIR: a scalable wrapper for processing billions of immunoglobulin and T cell receptor sequences using IgBLAST

BACKGROUND

Recent advances in DNA sequencing technologies have enabled significant leaps in capacity to generate large volumes of DNA sequence data, which has spurred a rapid growth in the use of bioinformatics as a means of interrogating antibody variable gene repertoires. Common tools used for annotation of antibody sequences are often limited in functionality, modularity and usability.

RESULTS

We have developed PyIR, a Python wrapper and library for IgBLAST, which offers a minimal setup CLI and API, FASTQ support, file chunking for large sequence files, JSON and Python dictionary output, and built-in sequence filtering.

CONCLUSIONS

PyIR offers improved processing speed over multithreaded IgBLAST (version 1.14) when spawning more than 16 processes on a single computer system. Its customizable filtering and data encapsulation allow it to be adapted to a wide range of computing environments. The API allows for IgBLAST to be used in customized bioinformatics workflows.

A single mutation increases heavy-chain heterodimer assembly of bispecific antibodies by inducing structural disorder in one homodimer species

We previously reported efficient heavy-chain assembly of heterodimeric bispecific antibodies by exchanging the interdomain protein interface of the human IgG1 CH3 dimer with the protein interface of the constant α and β domains of the human T-cell receptor, a technology known as bispecific engagement by antibodies based on the T-cell receptor (BEAT). Efficient heterodimerization in mammalian cell transient transfections was observed, but levels were influenced by the nature of the binding arms, particularly in the Fab-scFv-Fc format. In this study, we report a single amino acid change that significantly and consistently improved the heterodimerization rate of this format (≥95%) by inducing partial disorder in one homodimer species without affecting the heterodimer. Correct folding and assembly of the heterodimer were confirmed by the high-resolution (1.88-1.98 Å) crystal structure presented here. Thermal stability and 1-anilinonaphthalene-8-sulfonic acid-binding experiments, comparing original BEAT, mutated BEAT, and "knobs-into-holes" interfaces, suggested a cooperative assembly process of heavy chains in heterodimers. The observed gain in stability of the interfaces could be classified in the following rank order: mutated BEAT > original BEAT > knobs-into-holes. We therefore propose that the superior cooperativity found in BEAT interfaces is the key driver of their greater performance. Furthermore, we show how the mutated BEAT interface can be exploited for the routine preparation of drug candidates, with minimal risk of homodimer contamination using a single Protein A chromatography step.

Characterization of Human Mucosal-associated Invariant T (MAIT) Cells

Mucosal-associated invariant T (MAIT) cells are a subset of unconventional T cells restricted by the major histocompatibility complex (MHC) class I-like molecule MHC-related protein 1 (MR1). MAIT cells are found throughout the body, especially in human blood and liver. Unlike conventional T cells, which are stimulated by peptide antigens presented by MHC molecules, MAIT cells recognize metabolite antigens derived from an intermediate in the microbial biosynthesis of riboflavin. MAIT cells mediate protective immunity to infections by riboflavin-producing microbes via the production of cytokines and cytotoxicity. The discovery of stimulating MAIT cell antigens allowed for the development of an analytical tool, the MR1 tetramer, that binds specifically to the MAIT T cell receptor (TCR) and is becoming the gold standard for identification of MAIT cells by flow cytometry. This article describes protocols to characterize the phenotype of human MAIT cells in blood and tissues by flow cytometry using fluorescently labeled human MR1 tetramers alongside antibodies specific for MAIT cell markers. © 2019 by John Wiley & Sons, Inc. The main protocols include: Basic Protocol 1: Determining the frequency and steady-state surface phenotype of human MAIT cells Basic Protocol 2: Determining the activation phenotype of human MAIT cells in blood Basic Protocol 3: Characterizing MAIT cell TCRs using TCR-positive reporter cell lines Alternate protocols are provided for determining the absolute number, transcription factor phenotype, and TCR usage of human MAIT cells; and determining activation phenotype by staining for intracellular markers, measuring secreted cytokines, and measuring fluorescent dye dilution due to proliferation. Additional methods are provided for determining the capacity of MAIT cells to produce cytokine independently of antigen using plate-bound or bead-immobilized CD3/CD28 stimulation; and determining the MR1-Ag dependence of MAIT cell activation using MR1-blocking antibody or competitive inhibition. For TCR-positive reporter cell lines, methods are also provided for evaluating the MAIT TCR-mediated MR1-Ag response, determining the capacity of the reporter lines to produce cytokine independently of antigen, determining the MR1-Ag dependence of the reporter lines, and evaluating the MR1-Ag response of the reporter lines using IL-2 secretion. Support Protocols describe the preparation of PBMCs from human blood, the preparation of single-cell suspensions from tissue, the isolation of MAIT cells by FACS and MACS, cloning MAIT TCRα and β chain genes and MR1 genes for transduction, generating stably and transiently transfected cells lines, generating a stable MR1 knockout antigen-presenting cell line, and generating monocyte-derived dendritic cells.

AncesTree: An interactive immunoglobulin lineage tree visualizer

High-throughput sequencing of human immunoglobulin genes allows analysis of antibody repertoires and the reconstruction of clonal lineage evolution. The study of antibodies (Abs) affinity maturation is of specific interest to understand the generation of Abs with high affinity or broadly neutralizing activities. Moreover, phylogenic analysis enables the identification of the key somatic mutations required to achieve optimal antigen binding. The Immcantation framework provides a start-to-finish set of analytical methods for high-throughput adaptive immune receptor repertoire sequencing (AIRR-Seq; Rep-Seq) data. Furthermore, Immcantation's Change-O package has developed IgPhyML, an algorithm designed to build specifically immunoglobulin (Ig) phylogenic trees. Meanwhile Phylip, an algorithm that has been originally developed for applications in ecology and macroevolution, can also be used for the phylogenic reconstruction of antibodies maturation pathway. To complement Ig lineages made by IgPhyML or Dnaml (Phylip), we developed AncesTree, a graphic user interface (GUI) that aims to give researchers the opportunity to interactively explore antibodies clonal evolution. AncesTree displays interactive immunoglobulins phylogenic tree, Ig related mutations and sequence alignments using additional information coming from specialized antibody tools. The GUI is a Java standalone application allowing interaction with Ig tree that can run under Windows, Linux and Mac OS.

Landscape of variable domain of heavy-chain-only antibody repertoire from alpaca

Heavy-chain-only antibodies (HCAbs), which are devoid of light chains, have been found naturally occurring in various species including camelids and cartilaginous fish. Because of their high thermostability, refoldability and capacity for cell permeation, the variable regions of the heavy chain of HCAbs (VHHs) have been widely used in diagnosis, bio-imaging, food safety and therapeutics. Most immunogenetic and functional studies of HCAbs are based on case studies or a limited number of low-throughput sequencing data. A complete picture derived from more abundant high-throughput sequencing (HTS) data can help us gain deeper insights. We cloned and sequenced the full-length coding region of VHHs in Alpaca (Vicugna pacos) via HTS in this study. A new pipeline was developed to conduct an in-depth analysis of the HCAb repertoires. Various critical features, including the length distribution of complementarity-determining region 3 (CDR3), V(D)J usage, VJ pairing, germline-specific mutation rate and germline-specific scoring profiles (GSSPs), were systematically characterized. The quantitative data show that V(D)J usage and VHH recombination are highly biased. Interestingly, we found that the average CDR3 length of classical VHHs is longer than that of non-classical ones, whereas the mutation rates are similar in both kinds of VHHs. Finally, GSSPs were built to quantitatively describe and compare sequences that originate from each VJ pair. Overall, this study presents a comprehensive landscape of the HCAb repertoire, which can provide useful guidance for the modeling of somatic hypermutation and the design of novel functional VHHs or VHH repertoires via evolutionary profiles.

Profiling immunoglobulin repertoires across multiple human tissues using RNA sequencing

Profiling immunoglobulin (Ig) receptor repertoires with specialized assays can be cost-ineffective and time-consuming. Here we report ImReP, a computational method for rapid and accurate profiling of the Ig repertoire, including the complementary-determining region 3 (CDR3), using regular RNA sequencing data such as those from 8,555 samples across 53 tissues types from 544 individuals in the Genotype-Tissue Expression (GTEx v6) project. Using ImReP and GTEx v6 data, we generate a collection of 3.6 million Ig sequences, termed the atlas of immunoglobulin repertoires (TAIR), across a broad range of tissue types that often do not have reported Ig repertoires information. Moreover, the flow of Ig clonotypes and inter-tissue repertoire similarities across immune-related tissues are also evaluated. In summary, TAIR is one of the largest collections of CDR3 sequences and tissue types, and should serve as an important resource for studying immunological diseases.

Information on immune receptor repertoire provides important insights on disease progression and therapy development, but can be expensive and time-consuming to obtain. Here the authors report ImReP, a computational method that can extract detailed immune repertoire information from existing tissue-specific RNA sequencing data.

Sequencing B cell receptors from ferrets (Mustela putorius furo)

The domestic ferret (Mustela putorius furo) provides a critical animal model to study human respiratory diseases. However immunological insights are restricted due to a lack of ferret-specific reagents and limited genetic information about ferret B and T cell receptors. Here, variable, diversity and joining genes within the ferret kappa, lambda and heavy chain immunoglobulin loci were annotated using available genomic information. A multiplex PCR approach was derived that facilitated the recovery of paired heavy and light chain immunoglobulin sequences from single sorted ferret B cells, allowing validation of predicted germline gene sequences and the identification of putative novel germlines. Eukaryotic expression vectors were developed that enabled the generation of recombinant ferret monoclonal antibodies. This work advances the ferret as an informative immunological model for viral diseases by allowing the in-depth interrogation of antibody-based immunity.

Structures of human antibodies bound to SARS-CoV-2 spike reveal common epitopes and recurrent features of antibodies

Neutralizing antibody responses to coronaviruses focus on the trimeric spike, with most against the receptor-binding domain (RBD). Here we characterized polyclonal IgGs and Fabs from COVID-19 convalescent individuals for recognition of coronavirus spikes. Plasma IgGs differed in their degree of focus on RBD epitopes, recognition of SARS-CoV, MERS-CoV, and mild coronaviruses, and how avidity effects contributed to increased binding/neutralization of IgGs over Fabs. Electron microscopy reconstructions of polyclonal plasma Fab-spike complexes showed recognition of both S1A and RBD epitopes. A 3.4Å cryo-EM structure of a neutralizing monoclonal Fab-S complex revealed an epitope that blocks ACE2 receptor-binding on "up" RBDs. Modeling suggested that IgGs targeting these sites have different potentials for inter-spike crosslinking on viruses and would not be greatly affected by identified SARS-CoV-2 spike mutations. These studies structurally define a recurrent anti-SARS-CoV-2 antibody class derived from VH3-53/VH3-66 and similarity to a SARS-CoV VH3-30 antibody, providing criteria for evaluating vaccine-elicited antibodies.

The ABCD database: a repository for chemically defined antibodies

The ABCD (for AntiBodies Chemically Defined) database is a repository of sequenced antibodies, integrating curated information about the antibody and its antigen with cross-links to standardized databases of chemical and protein entities. It is freely available to the academic community, accessible through the ExPASy server (https://web.expasy.org/abcd/). The ABCD database aims at helping to improve reproducibility in academic research by providing a unique, unambiguous identifier associated to each antibody sequence. It also allows to determine rapidly if a sequenced antibody is available for a given antigen.

Evaluation of Somatic Hypermutation Status in Chronic Lymphocytic Leukemia (CLL) in the Era of Next Generation Sequencing

Somatic hypermutation (SHM) status provides an important prognostic indicator for chronic lymphocytic leukemia (CLL), a very common type of mature B-cell leukemia. Owing to the adverse prognosis associated with an unmutated immunoglobulin heavy chain variable (IGHV) status, SHM testing is performed as a standard of care in CLL. Conventionally, SHM testing has been performed using labor intensive and primarily analog Sanger sequencing method following PCR amplification of the clonal immunoglobulin heavy chain gene rearrangements in CLL cells. In comparison, recent availability of next generation sequencing (NGS) allows more versatile detection and direct identification of clonal immunoglobulin gene rearrangements in neoplastic B-cell populations. The ability to identify specific clonal IGHV signature(s) in both baseline (diagnostic) and post-treatment settings enables unique clinical applications of NGS such as determination of SHM status, minimal residual disease (MRD) monitoring, clonal heterogeneity and B cell receptor IG stereotypy. We provide a review of current practices and recommendations for SHM determination using NGS including examples of difficult cases.

OGRDB: a reference database of inferred immune receptor genes

High-throughput sequencing of the adaptive immune receptor repertoire (AIRR-seq) is providing unprecedented insights into the immune response to disease and into the development of immune disorders. The accurate interpretation of AIRR-seq data depends on the existence of comprehensive germline gene reference sets. Current sets are known to be incomplete and unrepresentative of the degree of polymorphism and diversity in human and animal populations. A key issue is the complexity of the genomic regions in which they lie, which, because of the presence of multiple repeats, insertions and deletions, have not proved tractable with short-read whole genome sequencing. Recently, tools and methods for inferring such gene sequences from AIRR-seq datasets have become available, and a community approach has been developed for the expert review and publication of such inferences. Here, we present OGRDB, the Open Germline Receptor Database (https://ogrdb.airr-community.org), a public resource for the submission, review and publication of previously unknown receptor germline sequences together with supporting evidence.

Modeling and mitigation of high-concentration antibody viscosity through structure-based computer-aided protein design

For an antibody to be a successful therapeutic many competing factors require optimization, including binding affinity, biophysical characteristics, and immunogenicity risk. Additional constraints may arise from the need to formulate antibodies at high concentrations (>150 mg/ml) to enable subcutaneous dosing with reasonable volume (ideally <1.0 mL). Unfortunately, antibodies at high concentrations may exhibit high viscosities that place impractical constraints (such as multiple injections or large needle diameters) on delivery and impede efficient manufacturing. Here we describe the optimization of an anti-PDGF-BB antibody to reduce viscosity, enabling an increase in the formulated concentration from 80 mg/ml to greater than 160 mg/ml, while maintaining the binding affinity. We performed two rounds of structure guided rational design to optimize the surface electrostatic properties. Analysis of this set demonstrated that a net-positive charge change, and disruption of negative charge patches were associated with decreased viscosity, but the effect was greatly dependent on the local surface environment. Our work here provides a comprehensive study exploring a wide sampling of charge-changes in the Fv and CDR regions along with targeting multiple negative charge patches. In total, we generated viscosity measurements for 40 unique antibody variants with full sequence information which provides a significantly larger and more complete dataset than has previously been reported.

IMGT® Biocuration and Comparative Study of the T Cell Receptor Beta Locus of Veterinary Species Based on Homo sapiens TRB

IMGT®, the international ImMunoGeneTics information system® is the global reference in immunogenetics and immunoinformatics. By its creation in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS), IMGT® marked the advent of immunoinformatics, which emerged at the interface between immunogenetics and bioinformatics. IMGT® is specialized in the immunoglobulins (IG) or antibodies, T cell receptors (TR), major histocompatibility (MH), and proteins of the IgSF and MhSF superfamilies. T cell receptors are divided into two groups, αβ and γδ TR, which express distinct TR containing either α and β, or γ and δ chains, respectively. The TRβ locus (TRB) was recently described and annotated by IMGT® biocurators for several veterinary species, i.e., cat (Felis catus), dog (Canis lupus familiaris), ferret (Mustela putorius furo), pig (Sus scrofa), rabbit (Oryctolagus cuniculus), rhesus monkey (Macaca mulatta), and sheep (Ovis aries). The aim of the present study is to compare the genes of the TRB locus among these different veterinary species based on Homo sapiens. The results reveal that there are similarities but also differences including the number of genes by subgroup which may demonstrate duplications and/or deletions during evolution.

Murine gammaherpesvirus infection is skewed toward Igλ+ B cells expressing a specific heavy chain V-segment

One of the defining characteristics of the B cell receptor (BCR) is the extensive diversity in the repertoire of immunoglobulin genes that make up the BCR, resulting in broad range of specificity. Gammaherpesviruses are B lymphotropic viruses that establish life-long infection in B cells, and although the B cell receptor plays a central role in B cell biology, very little is known about the immunoglobulin repertoire of gammaherpesvirus infected cells. To begin to characterize the Ig genes expressed by murine gammaherpesvirus 68 (MHV68) infected cells, we utilized single cell sorting to sequence and clone the Ig variable regions of infected germinal center (GC) B cells and plasma cells. We show that MHV68 infection is biased towards cells that express the Igλ light chain along with a single heavy chain variable gene, IGHV10-1*01. This population arises through clonal expansion but is not viral antigen specific. Furthermore, we show that class-switching in MHV68 infected cells differs from that of uninfected cells. Fewer infected GC B cells are class-switched compared to uninfected GC B cells, while more infected plasma cells are class-switched compared to uninfected plasma cells. Additionally, although they are germinal center derived, the majority of class switched plasma cells display no somatic hypermutation regardless of infection status. Taken together, these data indicate that selection of infected B cells with a specific BCR, as well as virus mediated manipulation of class switching and somatic hypermutation, are critical aspects in establishing life-long gammaherpesvirus infection.

ASAP-SML: An antibody sequence analysis pipeline using statistical testing and machine learning

Antibodies are capable of potently and specifically binding individual antigens and, in some cases, disrupting their functions. The key challenge in generating antibody-based inhibitors is the lack of fundamental information relating sequences of antibodies to their unique properties as inhibitors. We develop a pipeline, Antibody Sequence Analysis Pipeline using Statistical testing and Machine Learning (ASAP-SML), to identify features that distinguish one set of antibody sequences from antibody sequences in a reference set. The pipeline extracts feature fingerprints from sequences. The fingerprints represent germline, CDR canonical structure, isoelectric point and frequent positional motifs. Machine learning and statistical significance testing techniques are applied to antibody sequences and extracted feature fingerprints to identify distinguishing feature values and combinations thereof. To demonstrate how it works, we applied the pipeline on sets of antibody sequences known to bind or inhibit the activities of matrix metalloproteinases (MMPs), a family of zinc-dependent enzymes that promote cancer progression and undesired inflammation under pathological conditions, against reference datasets that do not bind or inhibit MMPs. ASAP-SML identifies features and combinations of feature values found in the MMP-targeting sets that are distinct from those in the reference sets.

VSV-Displayed HIV-1 Envelope Identifies Broadly Neutralizing Antibodies Class-Switched to IgG and IgA

The HIV-1 envelope (Env) undergoes conformational changes during infection. Broadly neutralizing antibodies (bNAbs) are typically isolated by using soluble Env trimers, which do not capture all Env states. To address these limitations, we devised a vesicular stomatitis virus (VSV)-based probe to display membrane-embedded Env trimers and isolated five bNAbs from two chronically infected donors, M4008 and M1214. Donor B cell receptor (BCR) repertoires identified two bNAb lineages, M4008_N1 and M1214_N1, that class-switched to immunoglobulin G (IgG) and IgA. Variants of these bNAbs reconstituted as IgA demonstrated broadly neutralizing activity, and the IgA fraction of M1214 plasma conferred neutralization. M4008_N1 epitope mapping revealed a glycan-independent V3 epitope conferring tier 2 virus neutralization. A 4.86-Å-resolution cryogenic electron microscopy (cryo-EM) structure of M1214_N1 complexed with CH505 SOSIP revealed another elongated epitope, the V2V5 corridor, extending from V2 to V5. Overall, the VSV_ENV probe identified bNAb lineages with neutralizing IgG and IgA members targeting distinct sites of HIV-1 Env vulnerability.

•

VSV-displayed HIV-1 envelope trimers identified five HIV-1 bNAbs

•

BCR repertoires identified two bNAb lineages class-switched to both IgG and IgA

•

The V3 crown-targeting bNAb M4008_N1 conferred tier 2 virus neutralization

•

Cryo-EM structure of bNAb M1214_N1 with CH505 SOSIP defined a V2V5 corridor epitope

Development of a RACE-based RNA-Seq approach to characterize the T-cell receptor repertoire of porcine γδ T cells

Recent data suggest that porcine γδ T cells exhibit a similar degree of functional plasticity as human and murine γδ T cells. Due to the high frequency of TCR-γδ⁺ cells in blood and secondary lymphatic organs, the pig is an attractive model to study these cells, especially their combined features of the innate and the adaptive immune system. Using a 5' RACE-like approach, we translated a human/murine NGS library preparation strategy to capture full-length V-(D)-J TRG and TRD clonotypes in swine. After oligo(dT) primed conversion of input RNA, the cDNA population was enriched for full-length V(D)J TCR transcripts with porcine-specific primers including Illumina adaptor sequences as overhangs for Illumina MiSeq analysis. After quality control and processing by FastQC and ea-utils, porcine TRG and TRD sequences were mapped against the human IMGT reference directory. Porcine blood-derived CD2⁺ and CD2‾ TCR-γδ⁺ cells exhibited two distinct clonotypes Vγ11JγP1 (74.6%) and Vγ10JγP1 (57.7%), respectively. Despite the high TCR-δ diversity among CD2⁺ cells (39 clonotypes), both subsets shared the same abundant Vδ1DδxJδ4 clonotype at approximately identically frequencies (CD2⁺: 31.2%; CD2‾: 37.0%). The flexible nature of this approach will facilitate the assessment of organ-specific phenotypes of γδ T cell subsets alongside with their respective TCR diversity at single cell resolution.

Genotyping and Copy Number Analysis of Immunoglobin Heavy Chain Variable Genes Using Long Reads

One of the remaining challenges to describing an individual's genetic variation lies in the highly heterogeneous and complex genomic regions that impede the use of classical reference-guided mapping and assembly approaches. Once such region is the Immunoglobulin heavy chain locus (IGH), which is critical for the development of antibodies and the adaptive immune system. We describe ImmunoTyper, the first PacBio-based genotyping and copy number calling tool specifically designed for IGH V genes (IGHV). We demonstrate that ImmunoTyper's multi-stage clustering and combinatorial optimization approach represents the most comprehensive IGHV genotyping approach published to date, through validation using gold-standard IGH reference sequence. This preliminary work establishes the feasibility of fine-grained genotype and copy number analysis using error-prone long reads in complex multi-gene loci and opens the door for in-depth investigation into IGHV heterogeneity using accessible and increasingly common whole-genome sequence.

Chromosome-level assembly of wild Bactrian camel genome reveals organization of immune gene loci

Camelids are characterized by their unique adaptive immune system that exhibits the generation of homodimeric heavy-chain immunoglobulins, somatic hypermutation of T-cell receptors, and low genetic diversity of major histocompatibility complex (MHC) genes. However, short-read assemblies are typically highly fragmented in these gene loci owing to their repetitive and polymorphic nature. Here, we constructed a chromosome-level assembly of wild Bactrian camel genome based on high-coverage long-read sequencing and chromatin interaction mapping. The assembly with a contig N50 of 5.37 Mb and a scaffold N50 of 76.03 Mb, represents the most contiguous camelid genome to date. The genomic organization of immunoglobulin heavy-chain locus was similar between the wild Bactrian camel and alpaca, and genes encoding for conventional and heavy-chain antibodies were intermixed. The organizations of two immunoglobulin light-chain loci and four T cell receptor loci were also fully deciphered using the new assembly. Additionally, the complete classical MHC region was resolved into a single contig. The high-quality assembly presented here provides an essential reference for future investigations examining the camelid immune system.

Salmonella Typhimurium biofilm disruption by a human antibody that binds a pan-amyloid epitope on curli

Bacterial biofilms, especially those associated with implanted medical devices, are difficult to eradicate. Curli amyloid fibers are important components of the biofilms formed by the Enterobacteriaceae family. Here, we show that a human monoclonal antibody with pan-amyloid-binding activity (mAb 3H3) can disrupt biofilms formed by Salmonella enterica serovar Typhimurium in vitro and in vivo. The antibody disrupts the biofilm structure, enhancing biofilm eradication by antibiotics and immune cells. In mice, 3H3 injections allow antibiotic-mediated clearance of catheter-associated S. Typhimurium biofilms. Thus, monoclonal antibodies that bind a pan-amyloid epitope have potential to prevent or eradicate bacterial biofilms.

Curli amyloid fibers are important components of bacterial biofilms formed by E. coli and Salmonella. Here, Tursi et al. show that a human monoclonal antibody with pan-amyloid binding activity can disrupt biofilms formed by Salmonella Typhimurium in vitro and in vivo.

Enterovirus 71 Infection Shapes Host T Cell Receptor Repertoire and Presumably Expands VP1-Specific TCRβ CDR3 Cluster

Enterovirus 71 (EV71) has become an important public health problem in the Asia-Pacific region in the past decades. EV71 infection might cause neurological and psychiatric complications and even death. Although an EV71 vaccine has been currently approved, there is no effective therapy for treating EV71-infected patients. Virus infections have been reported to shape host T cell receptor (TCR) repertoire. Therefore, understanding of host TCR repertoire in EV71 infection could better the knowledge in viral pathogenesis and further benefit the anti-viral therapy development. In this study, we used a mouse-adapted EV71 (mEV71) model to observe changes of host TCR repertoire in an EV71-infected central nervous system. Neonate mice were infected with mEV71 and mouse brainstem TCRβ repertoires were explored. Here, we reported that mEV71 infection impacted host brainstem TCRβ repertoire, where mEV71 infection skewed TCRβ diversity, changed VJ combination usages, and further expanded specific TCRβ CDR3 clones. Using bioinformatics analysis and ligand-binding prediction, we speculated the expanded TCRβ CDR3 clone harboring CASSLGANSDYTF sequence was capable of binding cleaved EV71 VP1 peptides in concert with major histocompatibility complex (MHC) molecules. We observed that mEV71 infection shaped host TCRβ repertoire and presumably expanded VP1-specific TCRβ CDR3 in mEV71-infected mouse brainstem that integrated EV71 pathogenesis in central nervous system.

Versatile Bispidine-Based Bifunctional Chelators for 64 CuII -Labelling of Biomolecules

Bifunctional chelators as parts of modular metal-based radiopharmaceuticals are responsible for stable complexation of the radiometal ion and for covalent linkage between the complex and the targeting vector. To avoid loss of complex stability, the bioconjugation strategy should not interfere with the radiometal chelation by occupying coordinating groups. The C9 position of the very stable CuII chelator 3,7-diazabicyclo[3.3.1]nonane (bispidine) is virtually predestined to introduce functional groups for facile bioconjugation as this functionalisation does not disturb the metal binding centre. We describe the preparation and characterisation of a set of novel bispidine derivatives equipped with suitable functional groups for diverse bioconjugation reactions, including common amine coupling strategies (bispidine-isothiocyanate) and the Cu-free strain-promoted alkyne-azide cycloaddition. We demonstrate their functionality and versatility in an exemplary way by conjugation to an antibody-based biomolecule and validate the obtained conjugate in vitro and in vivo.

A large peptidome dataset improves HLA class I epitope prediction across most of the human population

Prediction of HLA epitopes is important for the development of cancer immunotherapies and vaccines. However, current prediction algorithms have limited predictive power, in part because they were not trained on high-quality epitope datasets covering a broad range of HLA alleles. To enable prediction of endogenous HLA class I-associated peptides across a large fraction of the human population, we used mass spectrometry to profile >185,000 peptides eluted from 95 HLA-A, -B, -C and -G mono-allelic cell lines. We identified canonical peptide motifs per HLA allele, unique and shared binding submotifs across alleles and distinct motifs associated with different peptide lengths. By integrating these data with transcript abundance and peptide processing, we developed HLAthena, providing allele-and-length-specific and pan-allele-pan-length prediction models for endogenous peptide presentation. These models predicted endogenous HLA class I-associated ligands with 1.5-fold improvement in positive predictive value compared with existing tools and correctly identified >75% of HLA-bound peptides that were observed experimentally in 11 patient-derived tumor cell lines.

Aberrant B cell repertoire selection associated with HIV neutralizing antibody breadth

A goal of HIV vaccine development is to elicit antibodies with neutralizing breadth. Broadly neutralizing antibodies (bNAbs) to HIV often have unusual sequences with long heavy-chain complementarity-determining region loops, high somatic mutation rates and polyreactivity. A subset of HIV-infected individuals develops such antibodies, but it is unclear whether this reflects systematic differences in their antibody repertoires or is a consequence of rare stochastic events involving individual clones. We sequenced antibody heavy-chain repertoires in a large cohort of HIV-infected individuals with bNAb responses or no neutralization breadth and uninfected controls, identifying consistent features of bNAb repertoires, encompassing thousands of B cell clones per individual, with correlated T cell phenotypes. These repertoire features were not observed during chronic cytomegalovirus infection in an independent cohort. Our data indicate that the development of numerous B cell lineages with antibody features associated with autoreactivity may be a key aspect in the development of HIV neutralizing antibody breadth.

High-Throughput Prediction of MHC Class I and II Neoantigens with MHCnuggets

Computational prediction of binding between neoantigen peptides and major histocompatibility complex (MHC) proteins can be used to predict patient response to cancer immunotherapy. Current neoantigen predictors focus on in silico estimation of MHC binding affinity and are limited by low predictive value for actual peptide presentation, inadequate support for rare MHC alleles, and poor scalability to high-throughput data sets. To address these limitations, we developed MHCnuggets, a deep neural network method that predicts peptide-MHC binding. MHCnuggets can predict binding for common or rare alleles of MHC class I or II with a single neural network architecture. Using a long short-term memory network (LSTM), MHCnuggets accepts peptides of variable length and is faster than other methods. When compared with methods that integrate binding affinity and MHC-bound peptide (HLAp) data from mass spectrometry, MHCnuggets yields a 4-fold increase in positive predictive value on independent HLAp data. We applied MHCnuggets to 26 cancer types in The Cancer Genome Atlas, processing 26.3 million allele-peptide comparisons in under 2.3 hours, yielding 101,326 unique predicted immunogenic missense mutations (IMM). Predicted IMM hotspots occurred in 38 genes, including 24 driver genes. Predicted IMM load was significantly associated with increased immune cell infiltration (P < 2 × 10^-16), including CD8⁺ T cells. Only 0.16% of predicted IMMs were observed in more than 2 patients, with 61.7% of these derived from driver mutations. Thus, we describe a method for neoantigen prediction and its performance characteristics and demonstrate its utility in data sets representing multiple human cancers.

A Proposed New Nomenclature for the Immunoglobulin Genes of Mus musculus

Mammalian immunoglobulin (IG) genes are found in complex loci that contain hundreds of highly similar pseudogenes, functional genes and repetitive elements, which has made their investigation particularly challenging. High-throughput sequencing has provided new avenues for the investigation of these loci, and has recently been applied to study the IG genes of important inbred mouse strains, revealing unexpected differences between their IG loci. This demonstrated that the structural differences are of such magnitude that they call into question the merits of the current mouse IG gene nomenclatures. Three nomenclatures for the mouse IG heavy chain locus (Igh) are presently in use, and they are all positional nomenclatures using the C57BL/6 genome reference sequence as their template. The continued use of these nomenclatures requires that genes of other inbred strains be confidently identified as allelic variants of C57BL/6 genes, but this is clearly impossible. The unusual breeding histories of inbred mouse strains mean that, regardless of the genetics of wild mice, no single ancestral origin for the IG loci exists for laboratory mice. Here we present a general discussion of the challenges this presents for any IG nomenclature. Furthermore, we describe principles that could be followed in the formulation of a solution to these challenges. Finally, we propose a non-positional nomenclature that accords with the guidelines of the International Mouse Nomenclature Committee, and outline strategies that can be adopted to meet the nomenclature challenges if three systems are to give way to a new one.

Investigation of Antigen-Specific T-Cell Receptor Clusters in Human Cancers

PURPOSE

Cancer antigen-specific T cells are key components in antitumor immune response, yet their identification in the tumor microenvironment remains challenging, as most cancer antigens are unknown. Recent advance in immunology suggests that similar T-cell receptor (TCR) sequences can be clustered to infer shared antigen specificity. This study aims to identify antigen-specific TCRs from the tumor genomics sequencing data.

EXPERIMENTAL DESIGN

We used the TRUST (Tcr Repertoire Utilities for Solid Tissue) algorithm to assemble the TCR hypervariable CDR3 regions from 9,700 bulk tumor RNA-sequencing (RNA-seq) samples, and developed a computational method, iSMART, to group similar TCRs into antigen-specific clusters. Integrative analysis on the TCR clusters with multi-omics datasets was performed to profile cancer-associated T cells and to uncover novel cancer antigens.

RESULTS

Clustered TCRs are associated with signatures of T-cell activation after antigen encounter. We further elucidated the phenotypes of clustered T cells using single-cell RNA-seq data, which revealed a novel subset of tissue-resident memory T-cell population with elevated metabolic status. An exciting application of the TCR clusters is to identify novel cancer antigens, exemplified by our identification of a candidate cancer/testis gene, HSFX1, through integrated analysis of HLA alleles and genomics data. The target was further validated using vaccination of humanized HLA-A*02:01 mice and ELISpot assay. Finally, we showed that clustered tumor-infiltrating TCRs can differentiate patients with early-stage cancer from healthy donors, using blood TCR repertoire sequencing data, suggesting potential applications in noninvasive cancer detection.

CONCLUSIONS

Our analysis on the antigen-specific TCR clusters provides a unique resource for alternative antigen discovery and biomarker identification for cancer immunotherapies.

PTMProphet: Fast and Accurate Mass Modification Localization for the Trans-Proteomic Pipeline

Spectral matching sequence database search engines commonly used on mass spectrometry-based proteomics experiments excel at identifying peptide sequence ions, and in addition, possible sequence ions carrying post-translational modifications (PTMs), but most do not provide confidence metrics for the exact localization of those PTMs when several possible sites are available. Localization is absolutely required for downstream molecular cell biology analysis of PTM function in vitro and in vivo. Therefore, we developed PTMProphet, a free and open-source software tool integrated into the Trans-Proteomic Pipeline, which reanalyzes identified spectra from any search engine for which pepXML output is available to provide localization confidence to enable appropriate further characterization of biologic events. Localization of any type of mass modification (e.g., phosphorylation) is supported. PTMProphet applies Bayesian mixture models to compute probabilities for each site/peptide spectrum match where a PTM has been identified. These probabilities can be combined to compute a global false localization rate at any threshold to guide downstream analysis. We describe the PTMProphet tool, its underlying algorithms, and demonstrate its performance on ground-truth synthetic peptide reference data sets, one previously published small data set, one new larger data set, and also on a previously published phosphoenriched data set where the correct sites of modification are unknown. Data have been deposited to ProteomeXchange with identifier PXD013210.

A generalized HIV vaccine design strategy for priming of broadly neutralizing antibody responses

Vaccine induction of broadly neutralizing antibodies (bnAbs) to HIV remains a major challenge. Germline-targeting immunogens hold promise for initiating the induction of certain bnAb classes; yet for most bnAbs, a strong dependence on antibody heavy chain complementarity-determining region 3 (HCDR3) is a major barrier. Exploiting ultradeep human antibody sequencing data, we identified a diverse set of potential antibody precursors for a bnAb with dominant HCDR3 contacts. We then developed HIV envelope trimer-based immunogens that primed responses from rare bnAb-precursor B cells in a mouse model and bound a range of potential bnAb-precursor human naïve B cells in ex vivo screens. Our repertoire-guided germline-targeting approach provides a framework for priming the induction of many HIV bnAbs and could be applied to most HCDR3-dominant antibodies from other pathogens.