IMGT unique numbering for immunoglobulin and T cell receptor constant domains and Ig superfamily C-like domains

For antibody sequence generative modeling, mixture models may be all you need

MOTIVATION

Antibody therapeutic candidates must exhibit not only tight binding to their target but also good developability properties, especially low risk of immunogenicity.

RESULTS

In this work, we fit a simple generative model, SAM, to sixty million human heavy and seventy million human light chains. We show that the probability of a sequence calculated by the model distinguishes human sequences from other species with the same or better accuracy on a variety of benchmark datasets containing >400 million sequences than any other model in the literature, outperforming large language models (LLMs) by large margins. SAM can humanize sequences, generate new sequences, and score sequences for humanness. It is both fast and fully interpretable. Our results highlight the importance of using simple models as baselines for protein engineering tasks. We additionally introduce a new tool for numbering antibody sequences which is orders of magnitude faster than existing tools in the literature.

AVAILABILITY AND IMPLEMENTATION

All tools developed in this study are available at https://github.com/Wang-lab-UCSD/AntPack.

Prediction of polyspecificity from antibody sequence data by machine learning

Antibodies are generated with great diversity in nature resulting in a set of molecules, each optimized to bind a specific target. Taking advantage of their diversity and specificity, antibodies make up for a large part of recently developed biologic drugs. For therapeutic use antibodies need to fulfill several criteria to be safe and efficient. Polyspecific antibodies can bind structurally unrelated molecules in addition to their main target, which can lead to side effects and decreased efficacy in a therapeutic setting, for example via reduction of effective drug levels. Therefore, we created a neural-network-based model to predict polyspecificity of antibodies using the heavy chain variable region sequence as input. We devised a strategy for enriching antibodies from an immunization campaign either for antigen-specific or polyspecific binding properties, followed by generation of a large sequencing data set for training and cross-validation of the model. We identified important physico-chemical features influencing polyspecificity by investigating the behaviour of this model. This work is a machine-learning-based approach to polyspecificity prediction and, besides increasing our understanding of polyspecificity, it might contribute to therapeutic antibody development.

De novo antibody discovery in human blood from full-length single B cell transcriptomics and matching haplotyped-resolved germline assemblies

Immunoglobulin (IGH, IGK, IGL) loci in the human genome are highly polymorphic regions that encode the building blocks of the light and heavy chain IG proteins that dimerize to form antibodies. The processes of V(D)J recombination and somatic hypermutation in B cells are responsible for creating an enormous reservoir of highly specific antibodies capable of binding a vast array of possible antigens. However, the antibody repertoire is fundamentally limited by the set of variable (V), diversity (D), and joining (J) alleles present in the germline IG loci. To better understand how the germline IG haplotypes contribute to the expressed antibody repertoire, we combined genome sequencing of the germline IG loci with single-cell transcriptome sequencing of B cells from the same donor. Sequencing and assembly of the germline IG loci captured the IGH locus in a single fully-phased contig where the maternal and paternal contributions to the germline V, D, and J repertoire can be fully resolved. The B cells were collected following a measles, mumps, and rubella (MMR) vaccination, resulting in a population of cells that were activated in response to this specific immune challenge. Single-cell, full-length transcriptome sequencing of these B cells resulted in whole transcriptome characterization of each cell, as well as highly-accurate consensus sequences for the somatically rearranged and hypermutated light and heavy chain IG transcripts. A subset of antibodies synthesized based on their consensus heavy and light chain transcript sequences demonstrated binding to measles antigens and neutralization of measles live virus.

Identification of TRDV-TRAJ V domains in human and mouse T-cell receptor repertoires

Here, we describe the identification of two T-cell receptors (TRs) containing TRDV genes in their TRA chains, the first one in human and the second one in mouse. First, using 5'RACE on a mixed lymphocyte-tumor cell culture (MLTC), we identified TRDV1 5'-untranslated region (UTR) and complete coding sequence rearranged productively to TRAJ24. Single-cell TR RNA sequencing (RNA-seq) of the MLTC, conducted to identify additional clonotypes, revealed that the analysis software detected the hybrid TRDV-TRAJ TRA (TRA) chain but excluded it from the final results. In a separate project, we performed TR sequencing of tumor-infiltrating lymphocytes (TILs) in a murine tumor model. Here, the predominant clonotype contained a TRA chain with a TRDV2-2-TRAJ49 rearrangement. Again, the hybrid TRA chain was not reported in the final results. Transfection of both TR cDNAs resulted in cell surface localization of TR together with CD3, suggesting a productive protein in both cases. Tumor recognition of the Homo sapiens (Homsap) TRDV1-containing TR could be demonstrated by IFN Gamma ELISA ELISpot kit, whereas the Mus musculus (Musmus) TR did not recognize a tumor-derived cell line. To determine whether the TRDV-containing TRA chains we detected were rare events or whether TRDV genes are commonly incorporated into TRA chains, we queried the NCBI Sequence Read Archive for TR single-cell RNA-seq data and analyzed 21 human and 23 murine datasets. We found that especially Homsap TRDV1, Musmus TRDV1, and to some extent Musmus TRDV2-2 are more commonly incorporated into TRA chains than several TRAV genes, making those TRDV genes a relevant contribution to TRA diversity. TRDV-containing TRA chains are currently excluded from the final results of V-(D)-J dataset analyses with the CellRanger software. We provide a work-around to avoid exclusion of those hybrid TRA chains from the final analysis results.

Genomic and comparative analysis of the T cell receptor gamma locus in two Equus species

The genus Equus is the only extant genus of the Equidae family, which belongs to Perissodactyla, an order of mammals characterized by an odd number of toes (odd-toes ungulates). Taking advantage of the latest release of the genome assembly, we studied, for the first time in two organisms belonging to the Equus genus, the horse (Equus caballus) and the donkey (Equus asinus), the T cell receptor gamma (TRG) locus encoding the gamma chain of the γδ T cell receptor. Forty-five Variable (TRGV) genes belonging to the seven IMGT-NC validated mammalian TRGV subgroups, 25 Joining (TRGJ) and 17 Constant (TRGC) genes organized in 17 V-J-(J)-C cassettes, in tandem on about 1100 Kb, characterize the horse TRG locus, making the horse TRG locus the one with the greatest extension and with a significantly higher number of genes than the orthologous loci of the other mammalian species. A clonotype analysis of an RNA-seq transcriptomic dataset derived from spleen of an adult healthy horse, using the complete set of the horse TRGJ germline gene sequences as a probe, revealed that, in addition to the most prominent V-J rearrangements within each cassette, there is a relevant proportion of trans-cassette V-J recombination, whereby the same TRGV genes can recombine with different TRGJ genes spliced to the corresponding TRGC genes. This recombinant event strongly contributes to the diversity of the γ chain repertoire. In the donkey TRG locus, 34 TRGV, 21 TRGJ and 14 TRGC genes distributed in 14 V-J-(J)-C cassettes were found in a region of approximately 860 kb. Although the donkey's TRG is smaller than that of the horse, in Equus genus, this is still the second largest locus so far found in any mammalian species. Finally, the comparative analysis highlighted differences in size and gene content between the horse and donkey TRG loci, despite belonging to the same genus, indicating a good level of diversification within Equus. These data is in agreement with the evolutionary idea of the existence of a Equus recent common ancestor in rapid evolution, for which a mutation rate between horses and donkeys is more comparable to that between species belonging to different genera rather than to species of the same genus.

Evolution of T cell receptor beta loci in salmonids

T-cell mediated immunity relies on a vast array of antigen specific T cell receptors (TR). Characterizing the structure of TR loci is essential to study the diversity and composition of T cell responses in vertebrate species. The lack of good-quality genome assemblies, and the difficulty to perform a reliably mapping of multiple highly similar TR sequences, have hindered the study of these loci in non-model organisms. High-quality genome assemblies are now available for the two main genera of Salmonids, Salmo and Oncorhynchus. We present here a full description and annotation of the TRB loci located on chromosomes 19 and 25 of rainbow trout (Oncorhynchus mykiss). To get insight about variations of the structure and composition of TRB locus across salmonids, we compared rainbow trout TRB loci with other salmonid species and confirmed that the basic structure of salmonid TRB locus is a double set of two TRBV-D-J-C loci in opposite orientation on two different chromosomes. Our data shed light on the evolution of TRB loci in Salmonids after their whole genome duplication (WGD). We established a coherent nomenclature of salmonid TRB loci based on comprehensive annotation. Our work provides a fundamental basis for monitoring salmonid T cell responses by TRB repertoire sequencing.

KA-Search, a method for rapid and exhaustive sequence identity search of known antibodies

Antibodies with similar amino acid sequences, especially across their complementarity-determining regions, often share properties. Finding that an antibody of interest has a similar sequence to naturally expressed antibodies in healthy or diseased repertoires is a powerful approach for the prediction of antibody properties, such as immunogenicity or antigen specificity. However, as the number of available antibody sequences is now in the billions and continuing to grow, repertoire mining for similar sequences has become increasingly computationally expensive. Existing approaches are limited by either being low-throughput, non-exhaustive, not antibody specific, or only searching against entire chain sequences. Therefore, there is a need for a specialized tool, optimized for a rapid and exhaustive search of any antibody region against all known antibodies, to better utilize the full breadth of available repertoire sequences. We introduce Known Antibody Search (KA-Search), a tool that allows for the rapid search of billions of antibody variable domains by amino acid sequence identity across either the variable domain, the complementarity-determining regions, or a user defined antibody region. We show KA-Search in operation on the [Formula: see text]2.4 billion antibody sequences available in the OAS database. KA-Search can be used to find the most similar sequences from OAS within 30 minutes and a representative subset of 10 million sequences in less than 9 seconds. We give examples of how KA-Search can be used to obtain new insights about an antibody of interest. KA-Search is freely available at https://github.com/oxpig/kasearch .

Evolution of the Immunoglobulin Isotypes-Variations of Biophysical Properties among Animal Classes

The adaptive immune system arose around 500 million years ago in jawed fish, and, since then, it has mediated the immune defense against pathogens in all vertebrates. Antibodies play a central role in the immune reaction, recognizing and attacking external invaders. During the evolutionary process, several immunoglobulin isotypes emerged, each having a characteristic structural organization and dedicated function. In this work, we investigate the evolution of the immunoglobulin isotypes, in order to highlight the relevant features that were preserved over time and the parts that, instead, mutated. The residues that are coupled in the evolution process are often involved in intra- or interdomain interactions, meaning that they are fundamental to maintaining the immunoglobulin fold and to ensuring interactions with other domains. The explosive growth of available sequences allows us to point out the evolutionary conserved residues and compare the biophysical properties among different animal classes and isotypes. Our study offers a general overview of the evolution of immunoglobulin isotypes and advances the knowledge of their characteristic biophysical properties, as a first step in guiding protein design from evolution.

Astrocytes express aberrant immunoglobulins as putative gatekeeper of astrocytes to neuronal progenitor conversion

Using multi-omics analyses including RNAseq, RT-PCR, RACE-PCR, and shotgun proteomic with enrichment strategies, we demonstrated that newborn rat astrocytes produce neural immunoglobulin constant and variable heavy chains as well as light chains. However, their edification is different from the ones found in B cells and they resemble aberrant immunoglobulins observed in several cancers. Moreover, the complete enzymatic V(D)J recombination complex has also been identified in astrocytes. In addition, the constant heavy chain is also present in adult rat astrocytes, whereas in primary astrocytes from human fetus we identified constant and variable kappa chains as well as the substitution lambda chains known to be involved in pre-B cells. To gather insights into the function of these neural IgGs, CRISPR-Cas9 of IgG2B constant heavy chain encoding gene (Igh6), IgG2B overexpression, proximal labeling of rat astrocytes IgG2B and targets identification through 2D gels were performed. In Igh6 KO astrocytes, overrepresentation of factors involved in hematopoietic cells, neural stem cells, and the regulation of neuritogenesis have been identified. Moreover, overexpression of IgG2B in astrocytes induces the CRTC1-CREB-BDNF signaling pathway known to be involved in gliogenesis, whereas Igh6 KO triggers the BMP/YAP1/TEAD3 pathway activated in astrocytes dedifferentiation into neural progenitors. Proximal labeling experiments revealed that IgG2B is N-glycosylated by the OST complex, addressed to vesicle membranes containing the ATPase complex, and behaves partially like CD98hc through its association with LAT1. These experiments also suggest that proximal IgG2B-LAT1 interaction occurs concomitantly with MACO-1 and C2CD2L, at the heart of a potentially novel cell signaling platform. Finally, we demonstrated that these chains are synthesized individually and associated to recognize specific targets. Indeed, intermediate filaments Eif4a2 and Pdia6 involved in astrocyte fate constitute targets for these neural IgGs. Taken together, we hypothese that neural aberrant IgG chains may act as gatekeepers of astrocytes' fate.

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.

Evolution of surrogate light chain in tetrapods and the relationship between lengths of CDR H3 and VpreB tails

In the mammalian immune system, the surrogate light chain (SLC) shapes the antibody repertoire during B cell development by serving as a checkpoint for production of functional heavy chains (HC). Structural studies indicate that tail regions of VpreB contact and cover the third complementarity-determining region of the HC (CDR H3). However, some species, particularly bovines, have CDR H3 regions that may not be compatible with this HC-SLC interaction model. With immense structural and genetic diversity in antibody repertoires across species, we evaluated the genetic origins and sequence features of surrogate light chain components. We examined tetrapod genomes for evidence of conserved gene synteny to determine the evolutionary origin of VpreB1, VpreB2, and IGLL1, as well as VpreB3 and pre-T cell receptor alpha (PTCRA) genes. We found the genes for the SLC components (VpreB1, VpreB2, and IGLL1) only in eutherian mammals. However, genes for PTCRA occurred in all amniote groups and genes for VpreB3 occurred in all tetrapod groups, and these genes were highly conserved. Additionally, we found evidence of a new VpreB gene in non-mammalian tetrapods that is similar to the VpreB2 gene of eutherian mammals, suggesting VpreB2 may have appeared earlier in tetrapod evolution and may be a precursor to traditional VpreB2 genes in higher vertebrates. Among eutherian mammals, sequence conservation between VpreB1 and VpreB2 was low for all groups except rabbits and rodents, where VpreB2 was nearly identical to VpreB1 and did not share conserved synteny with VpreB2 of other species. VpreB2 of rabbits and rodents likely represents a duplicated variant of VpreB1 and is distinct from the VpreB2 of other mammals. Thus, rabbits and rodents have two variants of VpreB1 (VpreB1-1 and VpreB1-2) but no VpreB2. Sequence analysis of VpreB tail regions indicated differences in sequence content, charge, and length; where repertoire data was available, we observed a significant relationship between VpreB2 tail length and maximum DH length. We posit that SLC components co-evolved with immunoglobulin HC to accommodate the repertoire - particularly CDR H3 length and structure, and perhaps highly unusual HC (like ultralong HC of cattle) may bypass this developmental checkpoint altogether.

Self-reactivity of CD8 T-cell clones determines their differentiation status rather than their responsiveness in infections

Mature T cells are selected for recognizing self-antigens with low to intermediate affinity in the thymus. Recently, the relative differences in self-reactivity among individual T-cell clones were appreciated as important factors regulating their fate and immune response, but the role of self-reactivity in T-cell biology is incompletely understood. We addressed the role of self-reactivity in T-cell diversity by generating an atlas of mouse peripheral CD8⁺ T cells, which revealed two unconventional populations of antigen-inexperienced T cells. In the next step, we examined the steady-state phenotype of monoclonal T cells with various levels of self-reactivity. Highly self-reactive clones preferentially differentiate into antigen-inexperienced memory-like cells, but do not form a population expressing type I interferon-induced genes, showing that these two subsets have unrelated origins. The functional comparison of naïve monoclonal CD8⁺ T cells specific to the identical model antigen did not show any correlation between the level of self-reactivity and the magnitude of the immune response.

3D structures inferred from cDNA clones identify the CD1D-Restricted γδ T cell receptor in dromedaries

The Camelidae species occupy an important immunological niche within the humoral as well as cell mediated immune response. Although recent studies have highlighted that the somatic hypermutation (SHM) shapes the T cell receptor gamma (TRG) and delta (TRD) repertoire in Camelus dromedarius, it is still unclear how γδ T cells use the TRG/TRD receptors and their respective variable V-GAMMA and V-DELTA domains to recognize antigen in an antibody-like fashion. Here we report about 3D structural analyses of the human and dromedary γδ T cell receptor. First, we have estimated the interaction energies at the interface within the human crystallized paired TRG/TRD chains and quantified interaction energies within the same human TRG/TRD chains in complex with the CD1D, an RPI-MH1-LIKE antigen presenting glycoprotein. Then, we used the human TRG/TRD-CD1D complex as template for the 3D structure of the dromedary TRG/TRD-CD1D complex and for guiding the 3D human/dromedary comparative analysis. The choice of mutated TRG alternatively combined with mutated TRD cDNA clones originating from the spleen of one single dromedary was crucial to quantify the strength of the interactions at the protein-protein interface between the paired C. dromedarius TRG and TRD V-domains and between the C. dromedarius TRG/TRD V-domains and CD1D G-domains. Interacting amino acids located in the V-domain Complementarity Determining Regions (CDR) and Framework Regions (FR) according to the IMGT unique numbering for V-domains were identified. The resulting 3D dromedary TRG V-GAMMA combined with TRD V-DELTA protein complexes allowed to deduce the most stable gamma/delta chains pairings and to propose a candidate CD1D-restricted γδ T cell receptor complex.

ATM-TCR: TCR-Epitope Binding Affinity Prediction Using a Multi-Head Self-Attention Model

TCR-epitope pair binding is the key component for T cell regulation. The ability to predict whether a given pair binds is fundamental to understanding the underlying biology of the binding mechanism as well as developing T-cell mediated immunotherapy approaches. The advent of large-scale public databases containing TCR-epitope binding pairs enabled the recent development of computational prediction methods for TCR-epitope binding. However, the number of epitopes reported along with binding TCRs is far too small, resulting in poor out-of-sample performance for unseen epitopes. In order to address this issue, we present our model ATM-TCR which uses a multi-head self-attention mechanism to capture biological contextual information and improve generalization performance. Additionally, we present a novel application of the attention map from our model to improve out-of-sample performance by demonstrating on recent SARS-CoV-2 data.

O-Glycan-Dependent Interaction between MUC1 Glycopeptide and MY.1E12 Antibody by NMR, Molecular Dynamics and Docking Simulations

Anti-mucin1 (MUC1) antibodies have been widely used for breast cancer diagnosis and treatment. This is based on the fact that MUC1 undergoes aberrant glycosylation upon cancer progression, and anti-MUC1 antibodies differentiate changes in glycan structure. MY.1E12 is a promising anti-MUC1 antibody with a distinct specificity toward MUC1 modified with an immature O-glycan (NeuAcα(2-3)Galβ(1-3)GalNAc) on a specific Thr. However, the structural basis for the interaction between MY.1E12 and MUC1 remains unclear. The aim of this study is to elucidate the mode of interaction between MY.1E12 and MUC1 O-glycopeptide by NMR, molecular dynamics (MD) and docking simulations. NMR titration using MUC1 O-glycopeptides suggests that the epitope is located within the O-linked glycan and near the O-glycosylation site. MD simulations of MUC1 glycopeptide showed that the O-glycosylation significantly limits the flexibility of the peptide backbone and side chain of the O-glycosylated Thr. Docking simulations using modeled MY.1E12 Fv and MUC1 O-glycopeptide, suggest that V_H mainly contributes to the recognition of the MUC1 peptide portion while V_L mainly binds to the O-glycan part. The V_H/V_L-shared recognition mode of this antibody may be used as a template for the rational design and development of anti-glycopeptide antibodies.

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 .

Utilizing cross-product prior knowledge to rapidly de-risk chemical liabilities in therapeutic antibody candidates

There is considerable pressure in the pharmaceutical industry to advance better molecules faster. One pervasive concern for protein-based therapeutics is the presence of potential chemical liabilities. We have developed a simple methodology for rapidly de-risking specific chemical concerns in antibody-based molecules using prior knowledge of each individual liability at a specific position in the molecule’s sequence. Our methodology hinges on the development of sequence-aligned chemical liability databases of molecules from different stages of commercialization and on sequence-aligned experimental data from prior molecules that have been developed at Amgen. This approach goes beyond the standard practice of simply flagging all instances of each motif that fall in a CDR. Instead, we de-risk motifs that are common at a specific site in commercial mAb-based molecules (and therefore did not previously pose an insurmountable barrier to commercialization) and motifs at specific sites for which we have prior experimental data indicating acceptably low levels of modification. We have used this approach successfully to identify candidates in a discovery phase program with exclusively very low risk potential chemical liabilities. Identifying these candidates in the discovery phase allowed us to bypass protein engineering and accelerate the program’s timeline by 6 months.

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.

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.

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

Systematic Engineering of Optimized Autonomous Heavy-Chain Variable Domains

Autonomous heavy-chain variable (V_H) domains are the smallest functional antibody fragments, and they possess unique features, including small size and convex paratopes, which provide enhanced targeting of concave epitopes that are difficult to access with larger conventional antibodies. However, human V_H domains have evolved to fold and function with a light chain partner, and alone, they typically suffer from low stability and high aggregation propensity. Development of autonomous human V_H domains, in which aggregation propensity is reduced without compromising antigen recognition, has proven challenging. Here, we used an autonomous human V_H domain as a scaffold to construct phage-displayed synthetic libraries in which aspartate was systematically incorporated at different paratope positions. In selections, the library yielded many anti-EphA1 receptor V_H domains, which were characterized in detail. Structural analyses of a parental anti-EphA1 V_H domain and an improved variant provided insights into the effects of aspartate and other substitutions on preventing aggregation while retaining function. Our naïve libraries and in vitro selection procedures offer a systematic approach to generating highly functional autonomous human V_H domains that resist aggregation and could be used for basic research and biomedical applications.

Improving the yield of recalcitrant Nanobodies® by simple modifications to the standard protocol

Nanobodies are single-domain antibody constructs derived from the variable regions of heavy chain only (V_HH) camelid IgGs. Their small size and single gene format make them amenable to various molecular biology applications that require a protein affinity-based approach. These features, in addition to their high solubility, allows their periplasmic expression, extraction and purification in E. coli systems with relative ease, using standardized protocols. However, some Nanobodies are recalcitrant to periplasmic expression, extraction and purification within E. coli systems. To improve their expression would require either a change in the expression host, vector or an increased scale of expression, all of which entail an increase in the complexity of their expression, and production cost. However, as shown here, specific changes in the existing standard E. coli culture protocol, aimed at reducing breakdown of selective antibiotic pressure, increasing the initial culture inoculum and improving transport to the periplasmic space, rescued the expression of several such refractory Nanobodies. The periplasmic extraction protocol was also changed to ensure efficient osmolysis, prevent both protein degradation and prevent downstream chelation of Ni²⁺ ions during IMAC purification. Adoption of this protocol will lead to an improvement of the expression of Nanobodies in general, and specifically, those that are recalcitrant.

Genomic features of humoral immunity support tolerance model in Egyptian rousette bats

Bats asymptomatically harbor many viruses that can cause severe human diseases. The Egyptian rousette bat (ERB) is the only known reservoir for Marburgviruses and Sosuga virus, making it an exceptional animal model to study antiviral mechanisms in an asymptomatic host. With this goal in mind, we constructed and annotated the immunoglobulin heavy chain locus, finding an expansion on immunoglobulin variable genes associated with protective human antibodies to different viruses. We also annotated two functional and distinct immunoglobulin epsilon genes and four distinctive functional immunoglobulin gamma genes. We described the Fc receptor repertoire in ERBs, including features that may affect activation potential, and discovered the lack of evolutionary conserved short pentraxins. These findings reinforce the hypothesis that a differential threshold of regulation and/or absence of key immune mediators may promote tolerance and decrease inflammation in ERBs.

The T Cell Receptor (TRB) Locus in Tursiops truncatus: From Sequence to Structure of the Alpha/Beta Heterodimer in the Human/Dolphin Comparison

The bottlenose dolphin (Tursiops truncatus) belongs to the Cetartiodactyla and, similarly to other cetaceans, represents the most successful mammalian colonization of the aquatic environment. Here we report a genomic, evolutionary, and expression study of T. truncatus T cell receptor beta (TRB) genes. Although the organization of the dolphin TRB locus is similar to that of the other artiodactyl species, with three in tandem D-J-C clusters located at its 3' end, its uniqueness is given by the reduction of the total length due essentially to the absence of duplications and to the deletions that have drastically reduced the number of the germline TRBV genes. We have analyzed the relevant mature transcripts from two subjects. The simultaneous availability of rearranged T cell receptor α (TRA) and TRB cDNA from the peripheral blood of one of the two specimens, and the human/dolphin amino acids multi-sequence alignments, allowed us to calculate the most likely interactions at the protein interface between the alpha/beta heterodimer in complex with major histocompatibility class I (MH1) protein. Interacting amino acids located in the complementarity-determining region according to IMGT numbering (CDR-IMGT) of the dolphin variable V-alpha and beta domains were identified. According to comparative modelization, the atom pair contact sites analysis between the human MH1 grove (G) domains and the T cell receptor (TR) V domains confirms conservation of the structure of the dolphin TR/pMH.

The Genomic Organisation of the TRA/TRD Locus Validates the Peculiar Characteristics of Dromedary δ-Chain Expression

The role of γδ T cells in vertebrate immunity is still an unsolved puzzle. Species such as humans and mice display a low percentage of these T lymphocytes (i.e., “γδ low species”) with a restricted diversity of γδ T cell receptors (TR). Conversely, artiodactyl species (i.e., “γδ high species”) account for a high proportion of γδ T cells with large γ and δ chain repertoires. The genomic organisation of the TR γ (TRG) and δ (TRD) loci has been determined in sheep and cattle, noting that a wide number of germline genes that encode for γ and δ chains characterise their genomes. Taking advantage of the current improved version of the genome assembly, we have investigated the genomic structure and gene content of the dromedary TRD locus, which, as in the other mammalian species, is nested within the TR α (TRA) genes. The most remarkable finding was the identification of a very limited number of variable germline genes (TRDV) compared to sheep and cattle, which supports our previous expression analyses for which the somatic hypermutation mechanism is able to enlarge and diversify the primary repertoire of dromedary δ chains. Furthermore, the comparison between genomic and expressed sequences reveals that D genes, up to four incorporated in a transcript, greatly contribute to the increased diversity of the dromedary δ chain antigen binding-site.

IMGT® Biocuration and Comparative Analysis of Bos taurus and Ovis aries TRA/TRD Loci

The adaptive immune response provides the vertebrate immune system with the ability to recognize and remember specific pathogens to generate immunity, and mount stronger attacks each time the pathogen is encountered. T cell receptors are the antigen receptors of the adaptive immune response expressed by T cells, which specifically recognize processed antigens, presented as peptides by the highly polymorphic major histocompatibility (MH) proteins. T cell receptors (TR) are divided into two groups, αβ and γδ, which express distinct TR containing either α and β, or γ and δ chains, respectively. The TRα locus (TRA) and TRδ locus (TRD) of bovine (Bos taurus) and the sheep (Ovis aries) have recently been described and annotated by IMGT^® biocurators. The aim of the present study is to present the results of the biocuration and to compare the genes of the TRA/TRD loci among these ruminant species based on the Homo sapiens repertoire. The comparative analysis shows similarities but also differences, including the fact that these two species have a TRA/TRD locus about three times larger than that of humans and therefore have many more genes which may demonstrate duplications and/or deletions during evolution.

The Usage of Human IGHJ Genes Follows a Particular Non-random Selection: The Recombination Signal Sequence May Affect the Usage of Human IGHJ Genes

The formation of the B cell receptor (BCR) heavy chain variable region is derived from the germline V(D)J gene rearrangement according to the “12/23” rule and the “beyond 12/23” rule. The usage frequency of each V(D)J gene in the peripheral BCR repertoires is related to the initial recombination, self-tolerance selection, and the clonal proliferative response. However, their specific differences and possible mechanisms are still unknown. We analyzed in-frame and out-of-frame BCR-H repertoires from human samples with normal physiological and various pathological conditions by high-throughput sequencing. Our results showed that IGHJ gene frequency follows a similar pattern which is previously known, where IGHJ4 is used at high frequency (>40%), IGHJ6/IGHJ3/IGHJ5 is used at medium frequencies (10∼20%), and IGH2/IGHJ1 is used at low frequency (<4%) under whether normal physiological or various pathological conditions. However, our analysis of the recombination signal sequences suggested that the conserved non-amer and heptamer and certain 23 bp spacer length may affect the initial IGHD-IGHJ recombination, which results in different frequencies of IGHJ genes among the initial BCR-H repertoire. Based on this “initial repertoire,” we recommend that re-evaluation and further investigation are needed when analyzing the significance and mechanism of IGHJ gene frequency in self-tolerance selection and the clonal proliferative response.

QSAR Implementation for HIC Retention Time Prediction of mAbs Using Fab Structure: A Comparison between Structural Representations

Monoclonal antibodies (mAbs) constitute a rapidly growing biopharmaceutical sector. However, their growth is impeded by high failure rates originating from failed clinical trials and developability issues in process development. There is, therefore, a growing need for better in silico tools to aid in risk assessment of mAb candidates to promote early-stage screening of potentially problematic mAb candidates. In this study, a quantitative structure–activity relationship (QSAR) modelling workflow was designed for the prediction of hydrophobic interaction chromatography (HIC) retention times of mAbs. Three novel descriptor sets derived from primary sequence, homology modelling, and atomistic molecular dynamics (MD) simulations were developed and assessed to determine the necessary level of structural resolution needed to accurately capture the relationship between mAb structures and HIC retention times. The results showed that descriptors derived from 3D structures obtained after MD simulations were the most suitable for HIC retention time prediction with a R² = 0.63 in an external test set. It was found that when using homology modelling, the resulting 3D structures became biased towards the used structural template. Performing an MD simulation therefore proved to be a necessary post-processing step for the mAb structures in order to relax the structures and allow them to attain a more natural conformation. Based on the results, the proposed workflow in this paper could therefore potentially contribute to aid in risk assessment of mAb candidates in early development.

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.

SETE: Sequence-based Ensemble learning approach for TCR Epitope binding prediction

Predicting the binding of T cell receptors (TCRs) to epitopes plays a vital role in the immunotherapy, because it guides the development of therapeutic vaccines and cancer treatments. Many prediction methods attempted to explain the relationship between TCR repertoires from different aspects such as the V(D)J gene locus and the biophysical features of amino acids molecules, but the extraction of these features is time consuming and the performance of these models are limited. Few studies have investigated how k-mers formed by adjacent amino acids in TCR sequences direct the epitope recognition, and the specific mechanism of TCR epitope binding is still unclear. Motivated by these, we presented SETE (Sequence-based Ensemble learning approach for TCR Epitope binding prediction), a novel model to predict the TCR epitope binding accurately. The model deconstructed the CDR3β sequence to short amino acid chains as features and learned the pattern of them between different TCR repertoires with gradient boosting decision tree algorithm. Experiments have demonstrated that SETE can be helpful in predicting the TCRs' corresponding epitopes and it outperforms other state-of-the-art methods in predicting the epitope specificity of TCR on VDJdb data set. The source codes have been uploaded at https://github.com/wonanut/SETE for academic usage only.

Evolution of the T-Cell Receptor (TR) Loci in the Adaptive Immune Response: The Tale of the TRG Locus in Mammals

T lymphocytes are the principal actors of vertebrates' cell-mediated immunity. Like B cells, they can recognize an unlimited number of foreign molecules through their antigen-specific heterodimer receptors (TRs), which consist of αβ or γδ chains. The diversity of the TRs is mainly due to the unique organization of the genes encoding the α, β, γ, and δ chains. For each chain, multi-gene families are arranged in a TR locus, and their expression is guaranteed by the somatic recombination process. A great plasticity of the gene organization within the TR loci exists among species. Marked structural differences affect the TR γ (TRG) locus. The recent sequencing of multiple whole genome provides an opportunity to examine the TR gene repertoire in a systematic and consistent fashion. In this review, we report the most recent findings on the genomic organization of TRG loci in mammalian species in order to show differences and similarities. The comparison revealed remarkable diversification of both the genomic organization and gene repertoire across species, but also unexpected evolutionary conservation, which highlights the important role of the T cells in the immune response.

Generation and Characterization of a Functional Nanobody Against Inflammatory Chemokine CXCL10, as a Novel Strategy for the Treatment of Multiple Sclerosis

BACKGROUND & OBJECTIVE

Chemokines and their receptors play a pivotal role in the pathogenesis of various autoimmune diseases such as multiple sclerosis, infectious diseases, and also in cancer metastasis via attraction of the pathogenic immune cells into the inflammation sites.

METHODS

Inflammatory chemokine CXCL10 as a T helper (Th)1-chemokine directs chemotaxis of many cell subsets especially Th1 into the central nervous system (CNS) via its receptor CXCR3 and it has been put forward as a potential therapeutic target in the treatment of multiple sclerosis. Nanobodies are the smallest intact antigen binding fragments derived from heavy chain-only antibodies occurring in camelids with unique biochemical and biophysical features which render them superior to conventional antibodies or antibody fragments. Here, we describe the generation, selection, and characterization of CXCL10-specific Nanobodies from camel immunized with CXCL10. The obtained Nanobodies displayed high affinity towards CXCL10 about 10-11-10-8 M.

RESULTS

Then a Nanobody with the highest affinity named 3Nb12 was selected and investigated as a migration inhibitor of CXCR3+ cells. Chemotaxis assay results showed that 3Nb12 blocked CXCL10- CXCR3 binding and potently inhibited chemotaxis of CXCR3-transfected HEK293T cells.

CONCLUSION

The nanobody 3Nb12 might be a promising specific and powerful blocking agent of CXCL10 function, which can be used for diagnostic, therapeutic and research purposes in MS.

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.

Comprehensive genomic analysis of the dromedary T cell receptor gamma (TRG) locus and identification of a functional TRGC5 cassette

The emergent availability in public databases of more complete genome assemblies allows us to improve genomic data obtained by classical molecular cloning. The main goal of this study was to refine the genomic map of the dromedary TRG locus by integrating our previous genomic data with the analysis of recent genomic assemblies. We identified an additional TRGC cassette, defined as a V-J-C recombination unit, located at the 5' of the locus and made up of five TRGV genes followed by three TRGJ genes and one TRGC gene. Hence, the complete dromedary TRG locus spans about 105 Kb and consists of three in tandem TRGC cassettes delimited by AMPH and STARD3NL genes at the 5' and 3' end, respectively. An expression assay carried out on peripheral blood showed the functional competency for the dromedary TRGC5 cassette and confirmed the presence of the somatic hypermutation mechanism able to enlarge the repertoire diversity of the dromedary γδ T cells.

Concerted variation of the 3' regulatory region of Ig heavy chain and Gm haplotypes across human continental populations

OBJECTIVES

The 3' regulatory region of the immunoglobulin heavy chain gene (IGH) includes the HS1.2 enhancer displaying length polymorphism with four known variants. The goal of the research was to provide an overview of this variability and of its evolutionary significance across human populations.

MATERIALS AND METHODS

We compiled published and original data on HS1.2 polymorphism in 3,100 subjects from 26 human populations. Moreover, we imputed the haplotypic arrangement of the HS1.2 region in the 1000 Genomes Project (1KGP). In this dataset, imputation could also be obtained for the G1m-G3m allotype by virtue of the precise correspondence between serological types and amino acid (and DNA) substitutions in IGHG1 and IGHG3.

RESULTS

HS1.2 variant frequencies displayed similar patterns of continental partitioning as those reported in the literature for the physically neighboring IGHG1-IGHG3 system. The 1KGP data revealed that linkage disequilibrium (LD) can explain the spread of joint HS1.2-IGHG1-IGHG3 associations across continents and within continental populations, with stronger LD out of Africa and the features of an evolutionarily stable genomic block with differential expression in lymphoblastoid cell lines.

DISCUSSION

Strong population structuring involves at least the entire 70 kb genomic region here considered, due to the tight LD which maintained HS1.2, IGHG1, and IGHG3 in nonrandom arrangements. This might be key to better understand the evolutionary path of the entire genomic region driven by immune response capabilities, during the formation of continental gene pools.

Characterization of the ferret TRB locus guided by V, D, J, and C gene expression analysis

The domestic ferret, Mustela putorius furo, is an important mammalian animal model to study human respiratory infection. However, insufficient genomic annotation hampers detailed studies of ferret T cell responses. In this study, we analyzed the published T cell receptor beta (TRB) locus and performed high-throughput sequencing (HTS) of peripheral blood of four healthy adult ferrets to identify expressed V, D, J, and C genes. The HTS data is used as a guide to manually curate the expressed V, D, J, and C genes. The ferret locus appears to be most similar to that of the dog. Like other mammalian TRB loci, the ferret TRB locus contains a library of variable genes located upstream of two D-J-C gene clusters, followed by a (in the ferret non-functional) V gene with an inverted transcriptional orientation. All TRB genes (expressed or not) reported here have been approved by the IMGT/WHO-IUIS nomenclature committee.

Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus)

Simple Summary

The study of European rabbit immunoglobulin genes has contributed decisively to the current knowledge on antibody structure and diversification. The European rabbit has also been increasingly used as an animal model for the study of many human diseases, such as syphilis, tuberculosis, and AIDS. As such, the study of its immune system genes is of crucial relevance, but the study of rabbit immunoglobulins has focused only on the IgG and IgA antibodies. In this study, we added to the knowledge of the rabbit immune system by investigating the genetic diversity of two antibodies, IgM and IgE, in wild and domestic rabbits as well as other rabbit close species. With the data obtained in this study, we showed a high similarity between the different rabbit close species studied and we pointed out important genetic differences in the wild and domestic rabbits. Our findings are a valuable tool for the management of rabbit wild populations and domestic breeds and may contribute to the identification of immunoglobulins genetic variants with greater efficiency against pathogens.

Abstract

The European rabbit (Oryctolagus cuniculus) has been an important model for immunological studies but the study of its immunoglobulins (Ig) has been restricted to its unique IgA and IgG. Here, we studied the genetic diversity of IgM and IgE in several species of leporids and performed population genetics studies on European rabbit wild populations and domestic breeds. The leporids sequencing showed that these Ig are well conserved (98% sequence similarity among leporids), For IgM the Cµ1 and Cµ4 were the most diverse and most conserved domains, respectively, while for IgE the Cε1 was the most diverse domain and Cε2 and Cε3 the most conserved domains. The differences in the pattern of most conserved and most diverse domain between the Ig isotypes are most likely related to each isotype function. The genetic population data showed contrasting results for IgM and IgE. For both Ig, as expected, a greater diversity was observed in the original species range, the Iberian Peninsula. However, unexpectedly the genetic diversity found for IgE in the domestic animals is higher than that for the French wild populations. These results will increase knowledge of the genetic diversity of leporids and wild and domestic rabbit populations and are important tools for the management of wild populations and rabbitries.

Standardized IMGT® Nomenclature of Salmonidae IGH Genes, the Paradigm of Atlantic Salmon and Rainbow Trout: From Genomics to Repertoires

In teleost fish as in mammals, humoral adaptive immunity is based on B lymphocytes expressing highly diverse immunoglobulins (IG). During B cell differentiation, IG loci are subjected to genomic rearrangements of V, D, and J genes, producing a unique antigen receptor expressed on the surface of each lymphocyte. During the course of an immune response to infections or immunizations, B cell clones specific of epitopes from the immunogen are expanded and activated, leading to production of specific antibodies. Among teleost fish, salmonids comprise key species for aquaculture. Rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) are especially important from a commercial point of view and have emerged as critical models for fish immunology. The growing interest to capture accurate and comprehensive antibody responses against common pathogens and vaccines has resulted in recent efforts to sequence the IG repertoire in these species. In this context, a unified and standardized nomenclature of salmonid IG heavy chain (IGH) genes is urgently required, to improve accuracy of annotation of adaptive immune receptor repertoire dataset generated by high-throughput sequencing (AIRRseq) and facilitate comparisons between studies and species. Interestingly, the assembly of salmonids IGH genomic sequences is challenging due to the presence of two large size duplicated IGH loci and high numbers of IG genes and pseudogenes. We used data available for Atlantic salmon to establish an IMGT standardized nomenclature of IGH genes in this species and then applied the IMGT rules to the rainbow trout IGH loci to set up a nomenclature, which takes into account the specificities of Salmonid loci. This unique, consistent nomenclature for Salmonid IGH genes was then used to construct IMGT sequence reference directories allowing accurate annotation of AIRRseq data. The complex issues raised by the genetic diversity of salmon and trout strains are discussed in the context of IG repertoire annotation.

The Camel Adaptive Immune Receptors Repertoire as a Singular Example of Structural and Functional Genomics

The adaptive immune receptors repertoire is highly plastic, with its ability to produce antigen-binding molecules and select those with high affinity for their antigen. Species have developed diverse genetic and structural strategies to create their respective repertoires required for their survival in the different environments. Camelids, until now, considered as a case of evolutionary innovation because of their only heavy-chain antibodies, represent a new mammalian model particularly useful for understanding the role of diversity in the immune system function. Here, we review the structural and functional characteristics and the current status of the genomic organization of camel immunoglobulins (IG) or antibodies, α/ß and γ/δ T cell receptors (TR), and major histocompatibility complex (MHC). In camelid humoral response, in addition to the conventional antibodies, there are IG with “only-heavy-chain” (no light chain, and two identical heavy gamma chains lacking CH1 and with a VH domain designated as VHH). The unique features of these VHH offer advantages in biotechnology and for clinical applications. The TRG and TRD rearranged variable domains of Camelus dromedarius (Arabian camel) display somatic hypermutation (SHM), increasing the intrinsic structural stability in the γ/δ heterodimer and influencing the affinity maturation to a given antigen similar to immunoglobulin genes. The SHM increases the dromedary γ/δ repertoire diversity. In Camelus genus, the general structural organization of the TRB locus is similar to that of the other artiodactyl species, with a pool of TRBV genes positioned at the 5’ end of three in tandem D-J-C clusters, followed by a single TRBV gene with an inverted transcriptional orientation located at the 3’ end. At the difference of TRG and TRD, the diversity of the TRB variable domains is not shaped by SHM and depends from the classical combinatorial and junctional diversity. The MHC locus is located on chromosome 20 in Camelus dromedarius. Cytogenetic and comparative whole genome analyses revealed the order of the three major regions “Centromere-ClassII-ClassIII-ClassI”. Unexpectedly low extent of polymorphisms and haplotypes was observed in all Old World camels despite different geographic origins.

STCRDab: the structural T-cell receptor database

The Structural T–cell Receptor Database (STCRDab; http://opig.stats.ox.ac.uk/webapps/stcrdab) is an online resource that automatically collects and curates TCR structural data from the Protein Data Bank. For each entry, the database provides annotations, such as the α/β or γ/δ chain pairings, major histocompatibility complex details, and where available, antigen binding affinities. In addition, the orientation between the variable domains and the canonical forms of the complementarity-determining region loops are also provided. Users can select, view, and download individual or bulk sets of structures based on these criteria. Where available, STCRDab also finds antibody structures that are similar to TCRs, helping users explore the relationship between TCRs and antibodies.

Comparative Analysis of the TRB Locus in the Camelus Genus

T cells can be separated into two major subsets based on the heterodimer that forms their T cell receptors. αβ T cells have receptors consisting of α and β chains, while γδ T cells are composed of γ and δ chains. αβ T cells play an essential role within the adaptive immune responses against pathogens. The recent genomic characterization of the Camelus dromedarius T cell receptor β (TRB) locus has allowed us to infer the structure of this locus from the draft genome sequences of its wild and domestic Bactrian congeners, Camelus ferus and Camelus bactrianus. The general structural organization of the wild and domestic Bactrian TRB locus is similar to that of the dromedary, with a pool of TRBV genes positioned at the 5′ end of D-J-C clusters, followed by a single TRBV gene located at the 3′ end with an inverted transcriptional orientation. Despite the fragmented nature of the assemblies, comparative genomics reveals the existence of a perfect co-linearity between the three Old World camel TRB genomic sequences, which enables the transfer of information from one sequence to another and the filling of gaps in the genomic sequences. A virtual camelid TRB locus is hypothesized with the presence of 33 TRBV genes distributed in 26 subgroups. Likewise, in the artiodactyl species, three in-tandem D-J-C clusters, each composed of one TRBD gene, six or seven TRBJ genes, and one TRBC gene, are placed at the 3′ end of the locus. As reported in the ruminant species, a group of four functional TRY genes at the 5′ end and only one gene at the 3′ end, complete the camelid TRB locus. Although the gene content is similar, differences are observed in the TRBV functional repertoire, and genes that are functional in one species are pseudogenes in the other species. Hence, variations in the functional repertoire between dromedary, wild and domestic Bactrian camels, rather than differences in the gene content, may represent the molecular basis explaining the disparity in the TRB repertoire between the Camelus species. Finally, our data contribute to the knowledge about the evolutionary history of Old World camelids.

IMGT® and 30 Years of Immunoinformatics Insight in Antibody V and C Domain Structure and Function

At the 10th Human Genome Mapping (HGM10) Workshop, in New Haven, for the first time, immunoglobulin (IG) or antibody and T cell receptor (TR) variable (V), diversity (D), joining (J), and constant (C) genes were officially recognized as 'genes', as were the conventional genes. Under these HGM auspices, IMGT®, the international ImMunoGeneTics information system®, was created in June 1989 at Montpellier (University of Montpellier and CNRS). The creation of IMGT® marked the birth of immunoinformatics, a new science, at the interface between immunogenetics and bioinformatics. The accuracy and the consistency between genes and alleles, sequences, and three-dimensional (3D) structures are based on the IMGT Scientific chart rules generated from the IMGT-ONTOLOGY axioms and concepts: IMGT standardized keywords (IDENTIFICATION), IMGT gene and allele nomenclature (CLASSIFICATION), IMGT standardized labels (DESCRIPTION), IMGT unique numbering and IMGT Collier de Perles (NUMEROTATION). These concepts provide IMGT® immunoinformatics insights for antibody V and C domain structure and function, used for the standardized description in IMGT® web resources, databases and tools, immune repertoires analysis, single cell and/or high-throughput sequencing (HTS, NGS), antibody humanization, and antibody engineering in relation with effector properties.

Role of IgG3 in Infectious Diseases

IgG3 comprises only a minor fraction of IgG and has remained relatively understudied until recent years. Key physiochemical characteristics of IgG3 include an elongated hinge region, greater molecular flexibility, extensive polymorphisms, and additional glycosylation sites not present on other IgG subclasses. These characteristics make IgG3 a uniquely potent immunoglobulin, with the potential for triggering effector functions including complement activation, antibody (Ab)-mediated phagocytosis, or Ab-mediated cellular cytotoxicity (ADCC). Recent studies underscore the importance of IgG3 effector functions against a range of pathogens and have provided approaches to overcome IgG3-associated limitations, such as allotype-dependent short Ab half-life, and excessive proinflammatory activation. Understanding the molecular and functional properties of IgG3 may facilitate the development of improved Ab-based immunotherapies and vaccines against infectious diseases.

Memory CD4+ T cells are generated in the human fetal intestine

The fetus is thought to be protected from exposure to foreign antigens, yet CD45RO⁺ T cells reside in the fetal intestine. Here we combined functional assays with mass cytometry, single-cell RNA-sequencing and high-throughput T cell antigen receptor (TCR) sequencing to characterize the CD4⁺ T cell compartment in the human fetal intestine. We identified 22 CD4⁺ T cell clusters, including naive-like, regulatory-like and memory-like subpopulations, which were confirmed and further characterized at the transcriptional level. Memory-like CD4⁺ T cells had high expression of Ki-67, indicative of cell division, and CD5, a surrogate marker of TCR avidity, and produced the cytokines IFN-γ and IL-2. Pathway analysis revealed a differentiation trajectory associated with cellular activation and proinflammatory effector functions, and TCR repertoire analysis indicated clonal expansions, distinct repertoire characteristics and interconnections between subpopulations of memory-like CD4⁺ T cells. Imaging-mass cytometry indicated that memory-like CD4⁺ T cells colocalized with antigen-presenting cells. Collectively, these results provide evidence for the generation of memory-like CD4⁺ T cells in the human fetal intestine that is consistent with exposure to foreign antigens.

Transcription factories in Igκ allelic choice and diversity

The vertebrate immune system is tasked with the challenge of responding to any pathogen the organism might encounter, and retaining memory of that pathogen in case of future infection. Recognition and memory of pathogens are encoded within the adaptive immune system and production of T and B lymphocytes with diverse antigen receptor repertoires. In B lymphocytes, diversity is generated by sequential recombination between Variable (V), Diversity (D) and Joining (J) gene segments in the immunoglobulin heavy chain gene (Igh) and subsequent V-J recombination in immunoglobulin light chain genes (Igκ followed by Igλ). However, the process by which particular V, D and J segments are selected during recombination, and stochasticity is maintained to ensure antibody repertoire diversity, is still unclear. In this review, we focus on Igκ and recent findings regarding the relationships between gene structure, the generation of diversity and allelic choice. Surprisingly, the nuclear environment in which each Igκ allele resides, including transcription factories assembled on the nuclear matrix, plays critical roles in both gene regulation and in shaping the diversity of Vκ genes accessible to recombination. These findings provide a new paradigm for understanding Igκ recombination and Vκ diversity in the context of B lymphopoiesis.

Overview of the Germline and Expressed Repertoires of the TRB Genes in Sus scrofa

The α/β T cell receptor (TR) is a complex heterodimer that recognizes antigenic peptides and binds to major histocompatibility complex (MH) molecules. Both α and β chains are encoded by different genes localized on two distinct chromosomal loci: TRA and TRB. The present study employed the recent release of the swine genome assembly to define the genomic organization of the TRB locus. According to the sequencing data, the pig TRB locus spans approximately 400 kb of genomic DNA and consists of 38 TRBV genes belonging to 24 subgroups located upstream of three in tandem TRBD-J-C clusters, which are followed by a TRBV gene in an inverted transcriptional orientation. Comparative analysis confirms that the general organization of the TRB locus is similar among mammalian species, but the number of germline TRBV genes varies greatly even between species belonging to the same order, determining the diversity and specificity of the immune response. However, sequence analysis of the TRB locus also suggests the presence of blocks of conserved homology in the genomic region across mammals. Furthermore, by analysing a public cDNA collection, we identified the usage pattern of the TRBV, TRBD, and TRBJ genes in the adult pig TRB repertoire, and we noted that the expressed TRBV repertoire seems to be broader and more diverse than the germline repertoire, in line with the presence of a high level of TRBV gene polymorphisms. Because the nucleotide differences seems to be principally concentrated in the CDR2 region, it is reasonable to presume that most T cell β-chain diversity can be related to polymorphisms in pig MH molecules. Domestic pigs represent a valuable animal model as they are even more anatomically, genetically and physiologically similar to humans than are mice. Therefore, present knowledge on the genomic organization of the pig TRB locus allows the collection of increased information on the basic aspects of the porcine immune system and contributes to filling the gaps left by rodent models.