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

Insights into genetic aberrations and signalling pathway interactions in chronic lymphocytic leukemia: from pathogenesis to treatment strategies

Chronic lymphocytic leukemia (CLL) is prevalent in adults and is characterized by the accumulation of mature B cells in the blood, bone marrow, lymph nodes, and spleens. Recent progress in therapy and the introduction of targeted treatments [inhibitors of Bruton's tyrosine kinase (BTKi) or inhibitor of anti-apoptotic B-cell lymphoma-2 (Bcl-2i) protein (venetoclax)] in place of chemoimmunotherapy have significantly improved the outcomes of patients with CLL. These advancements have shifted the importance of traditional predictive markers, leading to a greater focus on resistance genes and reducing the significance of mutations, such as TP53 and del(17p). Despite the significant progress in CLL treatment, some patients still experience disease relapse. This is due to the substantial heterogeneity of CLL as well as the interconnected genetic resistance mechanisms and pathway adaptive resistance mechanisms to targeted therapies in CLL. Although the knowledge of the pathomechanism of CLL has expanded significantly in recent years, the precise origins of CLL and the interplay between various genetic factors remain incompletely understood, necessitating further research. This review enhances the molecular understanding of CLL by describing how BCR signalling, NF-κB PI3K/AKT, and ROR1 pathways sustain CLL cell survival, proliferation, and resistance to apoptosis. It also presents genetic and pathway-adaptive resistance mechanisms in CLL. Identifying B-cell receptor (BCR) signalling as a pivotal driver of CLL progression, the findings advocate personalized treatment strategies based on molecular profiling, emphasizing the need for further research to unravel the complex interplay between BCR signalling and its associated pathways to improve patient outcomes.

Ultrasensitive allele inference from immune repertoire sequencing data with MiXCR

Allelic variability in the adaptive immune receptor loci, which harbor the gene segments that encode B cell and T cell receptors (BCR/TCR), is of critical importance for immune responses to pathogens and vaccines. Adaptive immune receptor repertoire sequencing (AIRR-seq) has become widespread in immunology research making it the most readily available source of information about allelic diversity in immunoglobulin (IG) and T cell receptor (TR) loci. Here, we present a novel algorithm for extrasensitive and specific variable (V) and joining (J) gene allele inference, allowing the reconstruction of individual high-quality gene segment libraries. The approach can be applied for inferring allelic variants from peripheral blood lymphocyte BCR and TCR repertoire sequencing data, including hypermutated isotype-switched BCR sequences, thus allowing high-throughput novel allele discovery from a wide variety of existing data sets. The developed algorithm is a part of the MiXCR software. We demonstrate the accuracy of this approach using AIRR-seq paired with long-read genomic sequencing data, comparing it to a widely used algorithm, TIgGER. We applied the algorithm to a large set of IG heavy chain (IGH) AIRR-seq data from 450 donors of ancestrally diverse population groups, and to the largest reported full-length TCR alpha and beta chain (TRA and TRB) AIRR-seq data set, representing 134 individuals. This allowed us to assess the genetic diversity within the IGH, TRA, and TRB loci in different populations and to establish a database of alleles of V and J genes inferred from AIRR-seq data and their population frequencies with free public access through VDJ.online database.

The long and short of hyperdivergent regions

The increasing prevalence of genome sequencing and assembly has uncovered evidence of hyperdivergent genomic regions - loci with excess genetic diversity - in species across the tree of life. Hyperdivergent regions are often enriched for genes that mediate environmental responses, such as immunity, parasitism, and sensory perception. Especially in self-fertilizing species where the majority of the genome is homozygous, the existence of hyperdivergent regions might imply the historical action of evolutionary forces such as introgression and/or balancing selection. We anticipate that the application of new sequencing technologies, broader taxonomic sampling, and evolutionary modeling of hyperdivergent regions will provide insights into the mechanisms that generate and maintain genetic diversity within and between species.

Reading the repertoire: Progress in adaptive immune receptor analysis using machine learning

The adaptive immune system holds invaluable information on past and present immune responses in the form of B and T cell receptor sequences, but we are limited in our ability to decode this information. Machine learning approaches are under active investigation for a range of tasks relevant to understanding and manipulating the adaptive immune receptor repertoire, including matching receptors to the antigens they bind, generating antibodies or T cell receptors for use as therapeutics, and diagnosing disease based on patient repertoires. Progress on these tasks has the potential to substantially improve the development of vaccines, therapeutics, and diagnostics, as well as advance our understanding of fundamental immunological principles. We outline key challenges for the field, highlighting the need for software benchmarking, targeted large-scale data generation, and coordinated research efforts.

Novel polymorphic and copy number diversity in the antibody IGH locus of South African individuals

The heavy chain of an antibody is crucial for mediating antigen binding. IGHV genes, which partially encode the heavy chain of antibodies, exhibit vast genetic diversity largely through polymorphism and copy number variation (CNV). These genetic variations impact population-level expression levels. In this study, we analyzed expressed antibody transcriptomes and matched germline IGHV genes from donors from KwaZulu-Natal, South Africa. Amplicon NGS targeting germline IGHV sequences was performed on genomic DNA from 70 participants, eight of whom had matched datasets of expressed antibody transcriptomes. Germline IGHV sequencing identified 161 unique IGHV alleles, of which 32 were novel. A further 21 novel IGHV alleles were detected in the expressed transcriptomes of these donors. We also examined the datasets for CNV, uncovering gene duplications of 10 IGHV genes from germline sequencing and 33 genes in the expressed transcriptomes. Many of the IGHV gene duplications have not been described in other populations. This study expands our understanding of genetic differences in distinct populations and suggests the potential impact of genetic diversity on immune responses.

Applying phylogenetic methods for species delimitation to distinguish B-cell clonal families

The adaptive immune system generates a diverse array of B-cell receptors through the processes of V(D)J recombination and somatic hypermutation. B-cell receptors that bind to an antigen will undergo clonal expansion, creating a Darwinian evolutionary dynamic within individuals. A key step in studying these dynamics is to identify sequences derived from the same ancestral V(D)J recombination event (i.e. a clonal family). There are a number of widely used methods for accomplishing this task but a major limitation of all of them is that they rely, at least in part, on the ability to map sequences to a germline reference set. This requirement is particularly problematic in non-model systems where we often know little about the germline allelic diversity in the study population. Recognizing that delimiting B-cell clonal families is analogous to delimiting species from single locus data, we propose a novel strategy of reconstructing the phylogenetic tree of all B-cell sequences in a sample and using a popular species delimitation method, multi-rate Poisson Tree Processes (mPTP), to delimit clonal families. Using extensive simulations, we show that not only does this phylogenetically explicit approach perform well for the purpose of delimiting clonal families when no reference allele set is available, it performs similarly to state-of-the-art techniques developed specifically for B-cell data even when we have a complete reference allele set. Additionally, our analysis of an empirical dataset shows that mPTP performs similarly to leading methods in the field. These findings demonstrate the utility of using off-the-shelf phylogenetic techniques for analyzing B-cell clonal dynamics in non-model systems, and suggests that phylogenetic inference techniques may be potentially combined with mapping based approaches for even more robust inferences, even in model systems.

Thrifty wide-context models of B cell receptor somatic hypermutation

Somatic hypermutation (SHM) is the diversity-generating process in antibody affinity maturation. Probabilistic models of SHM are needed for analyzing rare mutations, for understanding the selective forces guiding affinity maturation, and for understanding the underlying biochemical process. High throughput data offers the potential to develop and fit models of SHM on relevant data sets. In this paper we model SHM using modern frameworks. We are motivated by recent work suggesting the importance of a wider context for SHM, however, assigning an independent rate to each k-mer leads to an exponential proliferation of parameters. Thus, using convolutions on 3-mer embeddings, we develop "thrifty" models of SHM that have fewer free parameters than a 5-mer model and yet have a significantly wider context. These offer a slight performance improvement over a 5-mer model. We also find that a per-site effect is not necessary to explain SHM patterns given nucleotide context. Also, the two current methods for fitting an SHM model - on out-of-frame sequence data and on synonymous mutations - produce significantly different results, and augmenting out-of-frame data with synonymous mutations does not aid out-of-sample performance.

Diverse priming outcomes under conditions of very rare precursor B cells

Rare B cells can have special pathogen-recognition features giving them the potential to make outsized contributions to protective immunity. However, rare naive B cells infrequently participate in immune responses. We investigated how germline-targeting vaccine antigen delivery and adjuvant selection affect priming of exceptionally rare BG18-like HIV broadly neutralizing antibody-precursor B cells (~1 in 50 million) in non-human primates. Only escalating dose (ED) priming immunization using the saponin adjuvant SMNP elicited detectable BG18-like cells in germinal centers (GCs). All groups had strong GC responses, but only ED+SMNP and bolus+SMNP induced BG18-like memory B cells in >50% of animals. One group had vaccine-specific GC responses equivalent to ED+SMNP, but BG18-like memory B cells were rarely detected. Following homologous boosting, BG18-like memory B cells were more frequent in a bolus priming group, but had lower somatic hypermutation and affinities. This outcome was inversely associated with post-prime antibody titers, suggesting antibody feedback can significantly influence rare precursor B cell responses.

The landscape of structural variation in aye-ayes (Daubentonia madagascariensis)

Aye-ayes (Daubentonia madagascariensis) are one of the 25 most critically endangered primate species in the world. Endemic to Madagascar, their small and highly fragmented populations make them particularly vulnerable to both genetic disease and anthropogenic environmental changes. Over the past decade, conservation genomic efforts have largely focused on inferring and monitoring population structure based on single nucleotide variants to identify and protect critical areas of genetic diversity. However, the recent release of a highly contiguous genome assembly allows, for the first time, for the study of structural genomic variation (deletions, duplications, insertions, and inversions) which are likely to impact a substantial proportion of the species' genome. Based on whole-genome, short-read sequencing data from 14 individuals, >1,000 high-confidence autosomal structural variants were detected, affecting ~240 kb of the aye-aye genome. The majority of these variants (>85%) were deletions shorter than 200 bp, consistent with the notion that longer structural mutations are often associated with strongly deleterious fitness effects. For example, two deletions longer than 850 bp located within disease-linked genes were predicted to impose substantial fitness deficits owing to a resulting frameshift and gene fusion, respectively; whereas several other major effect variants outside of coding regions are likely to impact gene regulatory landscapes. Taken together, this first glimpse into the landscape of structural variation in aye-ayes will enable future opportunities to advance our understanding of the traits impacting the fitness of this endangered species, as well as allow for enhanced evolutionary comparisons across the full primate clade.

Molecular characterization of the IgH locus and V(D)J recombination in large yellow croaker (Larimichthys crocea)

V(D)J recombination is crucial for generating a diverse repertoire of immunoglobulins. Although the V(D)J recombination process has been well characterized in mammals, this process remains largely unexplored in teleosts. In this study, we comprehensively analyzed the IgH locus of a marine fish species large yellow croaker (Larimichthys crocea), and identified 28 V, 19 D, and 8 J gene segments, following a pattern of V-Dζ-Jζ-Cζ-Dμ-Jμ-Cμ1-Cμ2. The V, D, and J gene segments are flanked by consensus recombination signal sequences, with spacer lengths similar to those observed in mammals. The V gene segments are categorized into three distinct families, and exhibited a higher sequence identity compared to those in mammals. Additionally, we designed a set of primers for the examination of the V(D)J recombination in large yellow croaker. RNA-seq analysis showed increased expression of genes related to immunoglobulin production and lymphocyte chemotaxis in IgM + B cells upon Pseudomonas plecoglossicida infection, accompanied by altered expression of V gene segments, suggesting their involvement in the response to P. plecoglossicida infection. Taken together, we identified the IgH locus and V(D)J recombination process of large yellow croaker, which contribute to the understanding of immunoglobulin production and B cell immunity in teleosts, and may provide insights into vaccine development in large yellow croaker.

Complete sequencing of ape genomes

We present haplotype-resolved reference genomes and comparative analyses of six ape species, namely: chimpanzee, bonobo, gorilla, Bornean orangutan, Sumatran orangutan, and siamang. We achieve chromosome-level contiguity with unparalleled sequence accuracy (<1 error in 500,000 base pairs), completely sequencing 215 gapless chromosomes telomere-to-telomere. We resolve challenging regions, such as the major histocompatibility complex and immunoglobulin loci, providing more in-depth evolutionary insights. Comparative analyses, including human, allow us to investigate the evolution and diversity of regions previously uncharacterized or incompletely studied without bias from mapping to the human reference. This includes newly minted gene families within lineage-specific segmental duplications, centromeric DNA, acrocentric chromosomes, and subterminal heterochromatin. This resource should serve as a definitive baseline for all future evolutionary studies of humans and our closest living ape relatives.

Enriched G4 forming repeats in the human genome are associated with robust well-coordinated transcription and reduced cancer transcriptome variation

Non-B DNA G-quadruplex (G4) structures with guanine (G) runs of 2 to 4 repeats can trigger opposing experimental transcriptional impacts. Here, we used bioinformatic algorithms to comprehensively assess correlations of steady-state RNA transcript levels with all putative G4 sequence (pG4) locations genome-wide in three mammalian genomes and in normal and tumor human tissues. The human pG4-containing gene set displays higher expression levels than the set without pG4, supporting and extending some prior observations. pG4 enrichment at transcription start sites (TSSs) in human, but not chimpanzee and mouse genomes, suggests possible positive selection pressure for pG4 at human TSS, potentially driving genome rewiring and gene expression divergence between human and chimpanzee. Comprehensive bioinformatic analyses revealed lower pG4-containing gene set variability in humans and among different pG4 genes in tumors. As G4 stabilizers are under therapeutic consideration for cancer and pathogens, such distinctions between human normal and tumor G4s along with other species merit attention. Furthermore, in germline and cancer sequences, the most mutagenic pG4 mapped to regions promoting alternative DNA structures. Overall findings establish high pG4 at TSS as a human genome attribute statistically associated with robust well-coordinated transcription and reduced cancer transcriptome variation with implications for biology, model organisms, and medicine.

An unbiased comparison of immunoglobulin sequence aligners

Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) is critical for our understanding of the adaptive immune system's dynamics in health and disease. Reliable analysis of AIRR-seq data depends on accurate rearranged immunoglobulin (Ig) sequence alignment. Various Ig sequence aligners exist, but there is no unified benchmarking standard representing the complexities of AIRR-seq data, obscuring objective comparisons of aligners across tasks. Here, we introduce GenAIRR, a modular simulation framework for generating Ig sequences alongside their ground truths. GenAIRR realistically simulates the intricacies of V(D)J recombination, somatic hypermutation, and an array of sequence corruptions. We comprehensively assessed prominent Ig sequence aligners across various metrics, unveiling unique performance characteristics for each aligner. The GenAIRR-produced datasets, combined with the proposed rigorous evaluation criteria, establish a solid basis for unbiased benchmarking of immunogenetics computational tools. It sets up the ground for further improving the crucial task of Ig sequence alignment, ultimately enhancing our understanding of adaptive immunity.

Evaluation of the Immune Response within the Tumor Microenvironment in African American and Non-Hispanic White Patients with Non-Small Cell Lung Cancer

BACKGROUND

African Americans have higher incidence and mortality from lung cancer than non-Hispanic Whites, but investigations into differences in immune response have been minimal. Therefore, we compared components of the tumor microenvironment among African Americans and non-Hispanic Whites diagnosed with non-small cell lung cancer based on PDL1 or tertiary lymphoid structure (TLS) status to identify differences of translational relevance.

METHODS

Using a cohort of 280 patients with non-small cell lung cancer from the Inflammation, Health, Ancestry, and Lung Epidemiology study (non-Hispanic White: n = 155; African American: n = 125), we evaluated PDL1 tumor proportion score (<1% vs. ≥1%) and TLS status (presence/absence), comparing differences within the tumor microenvironment based on immune cell distribution and differential expression of genes.

RESULTS

Tumors from African Americans had a higher proportion of plasma cell signatures within the tumor microenvironment than non-Hispanic Whites. In addition, gene expression patterns in African American PDL1-positive samples suggest that these tumors contained greater numbers of γδ T cells and resting dendritic cells, along with fewer CD8+ T cells after adjusting for age, sex, pack-years, stage, and histology. Investigation of differential expression of B cell/plasma cell-related genes between the two patient populations revealed that two immunoglobulin genes (IGKV2-29 and IGLL5) were associated with decreased mortality risk in African Americans.

CONCLUSIONS

In the first known race-stratified analysis of tumor microenvironment components in lung cancer based on PDL1 expression or TLS status, differences within the immune cell composition and transcriptomic signature were identified that may have therapeutic implications.

IMPACT

Future investigation of racial variation within the tumor microenvironment may help direct the use of immunotherapy.

B cell tolerance and autoimmunity: Lessons from repertoires

Adaptive immune cell function is regulated by a highly diverse receptor recombined from variable germline-encoded segments that can recognize an almost unlimited array of epitopes. While this diversity enables the recognition of any pathogen, it also poses a risk of self-recognition, leading to autoimmunity. Many layers of regulation are present during both the generation and activation of B cells to prevent this phenomenon, although they are evidently imperfect. In recent years, our ability to analyze immune repertoires at scale has drastically increased, both through advances in sequencing and single-cell analyses. Here, we review the current knowledge on B cell repertoire analyses, focusing on their implication for autoimmunity. These studies demonstrate that a failure of tolerance occurs at multiple independent checkpoints in different autoimmune contexts, particularly during B cell maturation, plasmablast differentiation, and within germinal centers. These failures are marked by distinct repertoire features that may be used to identify disease- or patient-specific therapeutic approaches.

Addressing Technical Pitfalls in Pursuit of Molecular Factors That Mediate Immunoglobulin Gene Regulation

The expressed Ab repertoire is a critical determinant of immune-related phenotypes. Ab-encoding transcripts are distinct from other expressed genes because they are transcribed from somatically rearranged gene segments. Human Abs are composed of two identical H and L chain polypeptides derived from genes in IGH locus and one of two L chain loci. The combinatorial diversity that results from Ab gene rearrangement and the pairing of different H and L chains contributes to the immense diversity of the baseline Ab repertoire. During rearrangement, Ab gene selection is mediated by factors that influence chromatin architecture, promoter/enhancer activity, and V(D)J recombination. Interindividual variation in the composition of the Ab repertoire associates with germline variation in IGH, implicating polymorphism in Ab gene regulation. Determining how IGH variants directly mediate gene regulation will require integration of these variants with other functional genomic datasets. In this study, we argue that standard approaches using short reads have limited utility for characterizing regulatory regions in IGH at haplotype resolution. Using simulated and chromatin immunoprecipitation sequencing reads, we define features of IGH that limit use of short reads and a single reference genome, namely 1) the highly duplicated nature of the DNA sequence in IGH and 2) structural polymorphisms that are frequent in the population. We demonstrate that personalized diploid references enhance performance of short-read data for characterizing mappable portions of the locus, while also showing that long-read profiling tools will ultimately be needed to fully resolve functional impacts of IGH germline variation on expressed Ab repertoires.

RAGging on recombination signal sequence strength for diffusion-mediated recombination

In this article, we discuss new insights into the distinct mechanisms for V(D)J recombination for different immunoglobulin loci. This follows the recent revelation that recombination signal sequences (RSS) within the IGKV locus have evolved to be more efficient mediators of recombination activating gene (RAG) recombination compared to the same elements in the IGH locus. This difference in RSS strength is proposed to be driven by different molecular mechanisms for RAG-mediated recombination between the two loci.

Assessing Assembly Errors in Immunoglobulin Loci: A Comprehensive Evaluation of Long-read Genome Assemblies Across Vertebrates

Long-read sequencing technologies have revolutionized genome assembly producing near-complete chromosome assemblies for numerous organisms, which are invaluable to research in many fields. However, regions with complex repetitive structure continue to represent a challenge for genome assembly algorithms, particularly in areas with high heterozygosity. Robust and comprehensive solutions for the assessment of assembly accuracy and completeness in these regions do not exist. In this study we focus on the assembly of biomedically important antibody-encoding immunoglobulin (IG) loci, which are characterized by complex duplications and repeat structures. High-quality full-length assemblies for these loci are critical for resolving haplotype-level annotations of IG genes, without which, functional and evolutionary studies of antibody immunity across vertebrates are not tractable. To address these challenges, we developed a pipeline, "CloseRead", that generates multiple assembly verification metrics for analysis and visualization. These metrics expand upon those of existing quality assessment tools and specifically target complex and highly heterozygous regions. Using CloseRead, we systematically assessed the accuracy and completeness of IG loci in publicly available assemblies of 74 vertebrate species, identifying problematic regions. We also demonstrated that inspecting assembly graphs for problematic regions can both identify the root cause of assembly errors and illuminate solutions for improving erroneous assemblies. For a subset of species, we were able to correct assembly errors through targeted reassembly. Together, our analysis demonstrated the utility of assembly assessment in improving the completeness and accuracy of IG loci across species.

Resolving haplotype variation and complex genetic architecture in the human immunoglobulin kappa chain locus in individuals of diverse ancestry

Immunoglobulins (IGs), critical components of the human immune system, are composed of heavy and light protein chains encoded at three genomic loci. The IG Kappa (IGK) chain locus consists of two large, inverted segmental duplications. The complexity of the IG loci has hindered use of standard high-throughput methods for characterizing genetic variation within these regions. To overcome these limitations, we use long-read sequencing to create haplotype-resolved IGK assemblies in an ancestrally diverse cohort (n = 36), representing the first comprehensive description of IGK haplotype variation. We identify extensive locus polymorphism, including novel single nucleotide variants (SNVs) and novel structural variants harboring functional IGKV genes. Among 47 functional IGKV genes, we identify 145 alleles, 67 of which were not previously curated. We report inter-population differences in allele frequencies for 10 IGKV genes, including alleles unique to specific populations within this dataset. We identify haplotypes carrying signatures of gene conversion that associate with SNV enrichment in the IGK distal region, and a haplotype with an inversion spanning the proximal and distal regions. These data provide a critical resource of curated genomic reference information from diverse ancestries, laying a foundation for advancing our understanding of population-level genetic variation in the IGK locus.

Computational mining of B cell receptor repertoires reveals antigen-specific and convergent responses to Ebola vaccination

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

Single-cell immune repertoire analysis

Single-cell T cell and B cell antigen receptor-sequencing data analysis can potentially perform in-depth assessments of adaptive immune cells that inform on understanding immune cell development to tracking clonal expansion in disease and therapy. However, it has been extremely challenging to analyze and interpret T cells and B cells and their adaptive immune receptor repertoires at the single-cell level due to not only the complexity of the data but also the underlying biology. In this Review, we delve into the computational breakthroughs that have transformed the analysis of single-cell T cell and B cell antigen receptor-sequencing data.

Coeliac disease: the paradox of diagnosing a food hypersensitivity disorder with autoantibodies

Serum antibodies to the autoantigen transglutaminase 2 (TG2) are increasingly harnessed to diagnose coeliac disease. Diagnostic guidelines for children give recommendation for a no-biopsy-based diagnosis through detection of high amounts of IgA anti-TG2 antibodies in serum with confirmation of positivity in a separate blood sample by characteristic autoantibody-staining of tissue. While measurement of IgA anti-TG2 also is important in the diagnostic workup of adults, the adult guidelines still mandate examination of gut biopsies. This requirement might well change in the future, as might the necessity for confirming autoantibody positivity by tissue staining. The key role of autoantibody serology for diagnosis of coeliac disease is paradoxical. Coeliac disease was considered, and still can be considered, a food intolerance disorder where autoantibodies at face value are out of place. The immunological mechanisms underlying the formation of autoantibodies in response to gluten exposure have been dissected. This review presents the current insights demonstrating that the autoantibodies in coeliac disease are intimately integrated in the maladapted immune response to gluten.

Adaptive immune receptor germline gene variation

Recognition of antigens by T cell receptors (TCRs) and B cell receptors (BCRs) is a key step in lymphocyte activation. T and B cells mediate adaptive immune responses, which protect us against infections and provide immunological memory, and also, in some instances, drive pathogenic responses in autoimmune diseases. TCRs and BCRs are encoded within loci that are known to be genetically diverse. However, the extent and functional impact of this variation, both in humans and model animals used in immunological research, remain largely unknown. Experimental and genetic evidence has demonstrated that the complementarity determining regions 1 and 2 (HCDR1 and HCDR2), encoded by the variable (V) region of TCRs and BCRs, also often make critical contacts with the targeted antigen. Thus, knowledge about allelic variation in the genes encoding TCRs and BCRs is critically important for understanding adaptive immune responses in outbred populations and to define responder and non-responder phenotypes.

Innovations and trends in antibody repertoire analysis

Monoclonal antibodies have revolutionized the treatment of human diseases, which has made them the fastest-growing class of therapeutics, with global sales expected to reach $346.6 billion USD by 2028. Advances in antibody engineering and development have led to the creation of increasingly sophisticated antibody-based therapeutics (e.g. bispecific antibodies and chimeric antigen receptor T cells). However, approaches for antibody discovery have remained comparatively grounded in conventional yet reliable in vitro assays. Breakthrough developments in high-throughput single B-cell sequencing and immunoglobulin proteomic serology, however, have enabled the identification of high-affinity antibodies directly from endogenous B cells or circulating immunoglobulin produced in vivo. Moreover, advances in artificial intelligence offer vast potential for antibody discovery and design with large-scale repertoire datasets positioned as the optimal source of training data for such applications. We highlight advances and recent trends in how these technologies are being applied to antibody repertoire analysis.

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.

AIRR-C IG Reference Sets: curated sets of immunoglobulin heavy and light chain germline genes

Introduction

Analysis of an individual's immunoglobulin (IG) gene repertoire requires the use of high-quality germline gene reference sets. When sets only contain alleles supported by strong evidence, AIRR sequencing (AIRR-seq) data analysis is more accurate and studies of the evolution of IG genes, their allelic variants and the expressed immune repertoire is therefore facilitated.

Methods

The Adaptive Immune Receptor Repertoire Community (AIRR-C) IG Reference Sets have been developed by including only human IG heavy and light chain alleles that have been confirmed by evidence from multiple high-quality sources. To further improve AIRR-seq analysis, some alleles have been extended to deal with short 3' or 5' truncations that can lead them to be overlooked by alignment utilities. To avoid other challenges for analysis programs, exact paralogs (e.g. IGHV1-69*01 and IGHV1-69D*01) are only represented once in each set, though alternative sequence names are noted in accompanying metadata.

Results and discussion

The Reference Sets include less than half the previously recognised IG alleles (e.g. just 198 IGHV sequences), and also include a number of novel alleles: 8 IGHV alleles, 2 IGKV alleles and 5 IGLV alleles. Despite their smaller sizes, erroneous calls were eliminated, and excellent coverage was achieved when a set of repertoires comprising over 4 million V(D)J rearrangements from 99 individuals were analyzed using the Sets. The version-tracked AIRR-C IG Reference Sets are freely available at the OGRDB website (https://ogrdb.airr-community.org/germline_sets/Human) and will be regularly updated to include newly observed and previously reported sequences that can be confirmed by new high-quality data.

Ultrasensitive allele inference from immune repertoire sequencing data with MiXCR

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

Widespread impact of immunoglobulin V-gene allelic polymorphisms on antibody reactivity

The ability of the human immune system to generate antibodies to any given antigen can be strongly influenced by immunoglobulin V-gene allelic polymorphisms. However, previous studies have provided only limited examples. Therefore, the prevalence of this phenomenon has been unclear. By analyzing >1,000 publicly available antibody-antigen structures, we show that many V-gene allelic polymorphisms in antibody paratopes are determinants for antibody binding activity. Biolayer interferometry experiments further demonstrate that paratope allelic polymorphisms on both heavy and light chains often abolish antibody binding. We also illustrate the importance of minor V-gene allelic polymorphisms with low frequency in several broadly neutralizing antibodies to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza virus. Overall, this study not only highlights the pervasive impact of V-gene allelic polymorphisms on antibody binding but also provides mechanistic insights into the variability of antibody repertoires across individuals, which in turn have important implications for vaccine development and antibody discovery.

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

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