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.

Maintenance of caecal homeostasis by diverse adaptive immune cells in the rhesus macaque.

Objectives: The caecum bridges the small and large intestine and plays a front-line role in discriminating gastrointestinal antigens. Although dysregulated in acute and chronic conditions, the tissue is often overlooked immunologically. Methods: To address this issue, we applied single-cell transcriptomic-V(D)J sequencing to FACS-isolated CD45+ caecal patch/lamina propria leukocytes from a healthy (5-year-old) female rhesus macaque ex vivo and coupled these data to VDJ deep sequencing reads from haematopoietic tissues. Results: We found caecal NK cells and ILC3s to co-exist with a spectrum of effector T cells partially derived from SOX4 + recent thymic emigrants. Tolerogenic Vγ8Vδ1-T cells, plastic CD4+ T helper cells and GZMK + EOMES + and TMIGD2 + tissue-resident memory CD8+ T cells were present and differed metabolically. An IL13 + GATA3 + Th2 subset expressing eicosanoid pathway enzymes was accompanied by IL1RL1 + GATA3 + regulatory T cells and a minor proportion of IgE+ plasma cells (PCs), illustrating tightly regulated type 2 immunity devoid of ILC2s. In terms of B lymphocyte lineages, caecal patch antigen-presenting memory B cells sat alongside germinal centre cells undergoing somatic hypermutation and differentiation into IGF1 + PCs. Prototypic gene expression signatures decreased across PC clusters, and notably, expanded IgA clonotypes could be traced in VDJ deep sequencing reads from additional compartments, including the bone marrow, supporting that these cells contribute a steady stream of systemic antibodies. Conclusions: The data advance our understanding of caecal immunological function, revealing processes involved in barrier maintenance and molecular networks relevant to disease.

Structural basis of broad SARS-CoV-2 cross-neutralization by affinity-matured public antibodies.

Descendants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant now account for almost all SARS-CoV-2 infections. The Omicron variant and its sublineages have spike glycoproteins that are highly diverged from the pandemic founder and first-generation vaccine strain, resulting in significant evasion from monoclonal antibody therapeutics and vaccines. Understanding how commonly elicited antibodies can broaden to cross-neutralize escape variants is crucial. We isolate IGHV3-53, using "public" monoclonal antibodies (mAbs) from an individual 7 months post infection with the ancestral virus and identify antibodies that exhibit potent and broad cross-neutralization, extending to the BA.1, BA.2, and BA.4/BA.5 sublineages of Omicron. Deep mutational scanning reveals these mAbs' high resistance to viral escape. Structural analysis via cryoelectron microscopy of a representative broadly neutralizing antibody, CAB-A17, in complex with the Omicron BA.1 spike highlights the structural underpinnings of this broad neutralization. By reintroducing somatic hypermutations into a germline-reverted CAB-A17, we delineate the role of affinity maturation in the development of cross-neutralization by a public class of antibodies.

Human immunoglobulin gene allelic variation impacts germline-targeting vaccine priming.

Vaccine priming immunogens that activate germline precursors for broadly neutralizing antibodies (bnAbs) have promise for development of precision vaccines against major human pathogens. In a clinical trial of the eOD-GT8 60mer germline-targeting immunogen, higher frequencies of vaccine-induced VRC01-class bnAb-precursor B cells were observed in the high dose compared to the low dose group. Through immunoglobulin heavy chain variable (IGHV) genotyping, statistical modeling, quantification of IGHV1-2 allele usage and B cell frequencies in the naive repertoire for each trial participant, and antibody affinity analyses, we found that the difference between dose groups in VRC01-class response frequency was best explained by IGHV1-2 genotype rather than dose and was most likely due to differences in IGHV1-2 B cell frequencies for different genotypes. The results demonstrate the need to define population-level immunoglobulin allelic variations when designing germline-targeting immunogens and evaluating them in clinical trials.

Multivalent antigen display on nanoparticle immunogens increases B cell clonotype diversity and neutralization breadth to pneumoviruses.

Nanoparticles for multivalent display and delivery of vaccine antigens have emerged as a promising avenue for enhancing B cell responses to protein subunit vaccines. Here, we evaluated B cell responses in rhesus macaques immunized with prefusion-stabilized respiratory syncytial virus (RSV) F glycoprotein trimer compared with nanoparticles displaying 10 or 20 copies of the same antigen. We show that multivalent display skews antibody specificities and drives epitope-focusing of responding B cells. Antibody cloning and repertoire sequencing revealed that focusing was driven by the expansion of clonally distinct B cells through recruitment of diverse precursors. We identified two antibody lineages that developed either ultrapotent neutralization or pneumovirus cross-neutralization from precursor B cells with low initial affinity for the RSV-F immunogen. This suggests that increased avidity by multivalent display facilitates the activation and recruitment of these cells. Diversification of the B cell response by multivalent nanoparticle immunogens has broad implications for vaccine design.

AIRR community curation and standardised representation for immunoglobulin and T cell receptor germline sets.

Analysis of an individual's immunoglobulin or T cell receptor gene repertoire can provide important insights into immune function. High-quality analysis of adaptive immune receptor repertoire sequencing data depends upon accurate and relatively complete germline sets, but current sets are known to be incomplete. Established processes for the review and systematic naming of receptor germline genes and alleles require specific evidence and data types, but the discovery landscape is rapidly changing. To exploit the potential of emerging data, and to provide the field with improved state-of-the-art germline sets, an intermediate approach is needed that will allow the rapid publication of consolidated sets derived from these emerging sources. These sets must use a consistent naming scheme and allow refinement and consolidation into genes as new information emerges. Name changes should be minimised, but, where changes occur, the naming history of a sequence must be traceable. Here we outline the current issues and opportunities for the curation of germline IG/TR genes and present a forward-looking data model for building out more robust germline sets that can dovetail with current established processes. We describe interoperability standards for germline sets, and an approach to transparency based on principles of findability, accessibility, interoperability, and reusability.

Unmodified rabies mRNA vaccine elicits high cross-neutralizing antibody titers and diverse B cell memory responses.

Licensed rabies virus vaccines based on whole inactivated virus are effective in humans. However, there is a lack of detailed investigations of the elicited immune response, and whether responses can be improved using novel vaccine platforms. Here we show that two doses of a lipid nanoparticle-formulated unmodified mRNA vaccine encoding the rabies virus glycoprotein (RABV-G) induces higher levels of RABV-G specific plasmablasts and T cells in blood, and plasma cells in the bone marrow compared to two doses of Rabipur in non-human primates. The mRNA vaccine also generates higher RABV-G binding and neutralizing antibody titers than Rabipur, while the degree of somatic hypermutation and clonal diversity of the response are similar for the two vaccines. The higher overall antibody titers induced by the mRNA vaccine translates into improved cross-neutralization of related lyssavirus strains, suggesting that this platform has potential for the development of a broadly protective vaccine against these viruses.

Vaccination of SARS-CoV-2-infected individuals expands a broad range of clonally diverse affinity-matured B cell lineages.

Vaccination of SARS-CoV-2 convalescent individuals generates broad and potent antibody responses. Here, we isolate 459 spike-specific monoclonal antibodies (mAbs) from two individuals who were infected with the index variant of SARS-CoV-2 and later boosted with mRNA-1273. We characterize mAb genetic features by sequence assignments to the donors' personal immunoglobulin genotypes and assess antibody neutralizing activities against index SARS-CoV-2, Beta, Delta, and Omicron variants. The mAbs used a broad range of immunoglobulin heavy chain (IGH) V genes in the response to all sub-determinants of the spike examined, with similar characteristics observed in both donors. IGH repertoire sequencing and B cell lineage tracing at longitudinal time points reveals extensive evolution of SARS-CoV-2 spike-binding antibodies from acute infection until vaccination five months later. These results demonstrate that highly polyclonal repertoires of affinity-matured memory B cells are efficiently recalled by vaccination, providing a basis for the potent antibody responses observed in convalescent persons following vaccination.

Adaptive immune receptor genotyping using the corecount program.

We present a new Rep-Seq analysis tool called corecount, for analyzing genotypic variation in immunoglobulin (IG) and T cell receptor (TCR) genes. corecount is highly efficient at identifying V alleles, including those that are infrequently used in expressed repertoires and those that contain 3' end variation that are otherwise refractory to reliable identification during germline inference from expressed libraries. Furthermore, corecount facilitates accurate D and J gene genotyping. The output is highly reproducible and facilitates the comparison of genotypes from multiple individuals, such as those from clinical cohorts. Here, we applied corecount to the genotypic analysis of IgM libraries from 16 individuals. To demonstrate the accuracy of corecount, we Sanger sequenced all the heavy chain IG alleles (65 IGHV, 27 IGHD and 7 IGHJ) from one individual from whom we also produced two independent IgM Rep-seq datasets. Genomic analysis revealed that 5 known IGHV and 2 IGHJ sequences are truncated in current reference databases. This dataset of genomically validated alleles and IgM libraries from the same individual provides a useful resource for benchmarking other bioinformatic programs that involve V, D and J assignments and germline inference, and may facilitate the development of AIRR-Seq analysis tools that can take benefit from the availability of more comprehensive reference databases.

Frequent use of IGHV3-30-3 in SARS-CoV-2 neutralizing antibody responses.

The antibody response to SARS-CoV-2 shows biased immunoglobulin heavy chain variable (IGHV) gene usage, allowing definition of genetic signatures for some classes of neutralizing antibodies. We investigated IGHV gene usage frequencies by sorting spike-specific single memory B cells from individuals infected with SARS-CoV-2 early in the pandemic. From two study participants and 703 spikespecific B cells, the most used genes were IGHV1-69, IGHV3-30-3, and IGHV3-30. Here, we focused on the IGHV3-30 group of genes and an IGHV3-30-3-using ultrapotent neutralizing monoclonal antibody, CAB-F52, which displayed broad neutralizing activity also in its germline-reverted form. IGHV3-30-3 is encoded by a region of the IGH locus that is highly variable at both the allelic and structural levels. Using personalized IG genotyping, we found that 4 of 14 study participants lacked the IGHV3-30-3 gene on both chromosomes, raising the question if other, highly similar IGHV genes could substitute for IGHV3-30-3 in persons lacking this gene. In the context of CAB-F52, we found that none of the tested IGHV3-33 alleles, but several IGHV3-30 alleles could substitute for IGHV3-30-3, suggesting functional redundancy between the highly homologous IGHV3-30 and IGHV3-30-3 genes for this antibody.

Human immunoglobulin gene allelic variation impacts germline-targeting vaccine priming.

Vaccine priming immunogens that activate germline precursors for broadly neutralizing antibodies (bnAbs) have promise for development of precision vaccines against major human pathogens. In a clinical trial of the eOD-GT8 60mer germline-targeting immunogen, higher frequencies of vaccine-induced VRC01-class bnAb-precursor B cells were observed in the high dose compared to the low dose group. Through immunoglobulin heavy chain variable (IGHV) genotyping, statistical modeling, quantification of IGHV1-2 allele usage and B cell frequencies in the naive repertoire for each trial participant, and antibody affinity analyses, we found that the difference between dose groups in VRC01-class response frequency was best explained by IGHV1-2 genotype rather than dose and was most likely due to differences in IGHV1-2 B cell frequencies for different genotypes. The results demonstrate the need to define population-level immunoglobulin allelic variations when designing germline-targeting immunogens and evaluating them in clinical trials. One-Sentence Summary: Human genetic variation can modulate the strength of vaccine-induced broadly neutralizing antibody precursor B cell responses.

Analysis of IGH allele content in a sample group of rheumatoid arthritis patients demonstrates unrevealed population heterogeneity.

Immunoglobulin heavy chain (IGH) germline gene variations influence the B cell receptor repertoire, with resulting biological consequences such as shaping our response to infections and altering disease susceptibilities. However, the lack of information on polymorphism frequencies in the IGH loci at the population level makes association studies challenging. Here, we genotyped a pilot group of 30 individuals with rheumatoid arthritis (RA) to examine IGH allele content and frequencies in this group. Eight novel IGHV alleles and one novel IGHJ allele were identified in the study. 15 cases were haplotypable using heterozygous IGHJ6 or IGHD anchors. One variant, IGHV4-34*01_S0742, was found in three out of 30 cases and included a single nucleotide change resulting in a non-canonical recombination signal sequence (RSS) heptamer. This variant allele, shown by haplotype analysis to be non-expressed, was also found in three out of 30 healthy controls and matched a single nucleotide polymorphism (SNP) described in the 1000 Genomes Project (1KGP) collection with frequencies that varied between population groups. Our finding of previously unreported alleles in a relatively small group of individuals with RA illustrates the need for baseline information about IG allelic frequencies in targeted study groups in preparation for future analysis of these genes in disease association studies.

Archaic humans have contributed to large-scale variation in modern human T cell receptor genes.

Human T cell receptors (TCRs) are critical for mediating immune responses to pathogens and tumors and regulating self-antigen recognition. Yet, variations in the genes encoding TCRs remain insufficiently defined. Detailed analysis of expressed TCR alpha, beta, gamma, and delta genes in 45 donors from four human populations-African, East Asian, South Asian, and European-revealed 175 additional TCR variable and junctional alleles. Most of these contained coding changes and were present at widely differing frequencies in the populations, a finding confirmed using DNA samples from the 1000 Genomes Project. Importantly, we identified three Neanderthal-derived, introgressed TCR regions including a highly divergent TRGV4 variant, which mediated altered butyrophilin-like molecule 3 (BTNL3) ligand reactivity and was frequent in all modern Eurasian population groups. Our results demonstrate remarkable variation in TCR genes in both individuals and populations, providing a strong incentive for including allelic variation in studies of TCR function in human biology.

Immunoglobulin germline gene polymorphisms influence the function of SARS-CoV-2 neutralizing antibodies.

The human immunoglobulin heavy-chain (IGH) locus is exceptionally polymorphic, with high levels of allelic and structural variation. Thus, germline IGH genotypes are personal, which may influence responses to infection and vaccination. For an improved understanding of inter-individual differences in antibody responses, we isolated SARS-CoV-2 spike-specific monoclonal antibodies from convalescent health care workers, focusing on the IGHV1-69 gene, which has the highest level of allelic variation of all IGHV genes. The IGHV1-69∗20-using CAB-I47 antibody and two similar antibodies isolated from an independent donor were critically dependent on allele usage. Neutralization was retained when reverting the V region to the germline IGHV1-69∗20 allele but lost when reverting to other IGHV1-69 alleles. Structural data confirmed that two germline-encoded polymorphisms, R50 and F55, in the IGHV1-69 gene were required for high-affinity receptor-binding domain interaction. These results demonstrate that polymorphisms in IGH genes can influence the function of SARS-CoV-2 neutralizing antibodies.

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.

Vaccination induces HIV broadly neutralizing antibody precursors in humans.

Broadly neutralizing antibodies (bnAbs) can protect against HIV infection but have not been induced by human vaccination. A key barrier to bnAb induction is vaccine priming of rare bnAb-precursor B cells. In a randomized, double-blind, placebo-controlled phase 1 clinical trial, the HIV vaccine-priming candidate eOD-GT8 60mer adjuvanted with AS01B had a favorable safety profile and induced VRC01-class bnAb precursors in 97% of vaccine recipients with median frequencies reaching 0.1% among immunoglobulin G B cells in blood. bnAb precursors shared properties with bnAbs and gained somatic hypermutation and affinity with the boost. The results establish clinical proof of concept for germline-targeting vaccine priming, support development of boosting regimens to induce bnAbs, and encourage application of the germline-targeting strategy to other targets in HIV and other pathogens.

Primary antibody response after influenza virus infection is first dominated by low-mutated HA-stem antibodies followed by higher-mutated HA-head antibodies.

Several studies have shown that the first encounter with influenza virus shapes the immune response to future infections or vaccinations. However, a detailed analysis of the primary antibody response is lacking as this is difficult to study in humans. It is therefore not known what the frequency and dynamics of the strain-specific hemagglutinin (HA) head- and stem-directed antibody responses are directly after primary influenza virus infection. Here, sera of twelve H1N1pdm2009 influenza virus-infected cynomolgus macaques were evaluated for HA-head and HA-stem domain antibody responses. We observed an early induction of HA-stem antibody responses, which was already decreased by day 56. In contrast, responses against the HA-head domain were low early after infection and increased at later timepoint. The HA-specific B cell repertoires in each animal showed diverse VH-gene usage with preferred VH-gene and JH-gene family usage for HA-head or HA-stem B cells but a highly diverse allelic variation within the VH-usage. HA-head B cells had shorter CDRH3s and higher VH-gene somatic hyper mutation levels relative to HA-stem B cells. In conclusion, our data suggest that HA-stem antibodies are the first to react to the infection while HA-head antibodies show a delayed response, but a greater propensity to enter the germinal center and undergo affinity maturation.

A BALB/c IGHV Reference Set, Defined by Haplotype Analysis of Long-Read VDJ-C Sequences From F1 (BALB/c x C57BL/6) Mice.

The immunoglobulin genes of inbred mouse strains that are commonly used in models of antibody-mediated human diseases are poorly characterized. This compromises data analysis. To infer the immunoglobulin genes of BALB/c mice, we used long-read SMRT sequencing to amplify VDJ-C sequences from F1 (BALB/c x C57BL/6) hybrid animals. Strain variations were identified in the Ighm and Ighg2b genes, and analysis of VDJ rearrangements led to the inference of 278 germline IGHV alleles. 169 alleles are not present in the C57BL/6 genome reference sequence. To establish a set of expressed BALB/c IGHV germline gene sequences, we computationally retrieved IGHV haplotypes from the IgM dataset. Haplotyping led to the confirmation of 162 BALB/c IGHV gene sequences. A musIGHV398 pseudogene variant also appears to be present in the BALB/cByJ substrain, while a functional musIGHV398 gene is highly expressed in the BALB/cJ substrain. Only four of the BALB/c alleles were also observed in the C57BL/6 haplotype. The full set of inferred BALB/c sequences has been used to establish a BALB/c IGHV reference set, hosted at https://ogrdb.airr-community.org. We assessed whether assemblies from the Mouse Genome Project (MGP) are suitable for the determination of the genes of the IGH loci. Only 37 (43.5%) of the 85 confirmed IMGT-named BALB/c IGHV and 33 (42.9%) of the 77 confirmed non-IMGT IGHV were found in a search of the MGP BALB/cJ genome assembly. This suggests that current MGP assemblies are unsuitable for the comprehensive documentation of germline IGHVs and more efforts will be needed to establish strain-specific reference sets.

Addressing IGHV Gene Structural Diversity Enhances Immunoglobulin Repertoire Analysis: Lessons From Rhesus Macaque.

The accurate germline gene assignment and assessment of somatic hypermutation in antibodies induced by immunization or infection are important in immunological studies. Here, we illustrate issues specific to the construction of comprehensive immunoglobulin (IG) germline gene reference databases for outbred animal species using rhesus macaques, a frequently used non-human primate model, as a model test case. We demonstrate that the genotypic variation found in macaque germline inference studies is reflected in similar levels of gene diversity in genomic assemblies. We show that the high frequency of IG heavy chain V (IGHV) region structural and gene copy number variation between subjects means that individual animals lack genes that are present in other animals. Therefore, gene databases compiled from a single or too few animals will inevitably result in inaccurate gene assignment and erroneous SHM level assessment for those genes it lacks. We demonstrate this by assigning a test macaque IgG library to the KIMDB, a database compiled of germline IGHV sequences from 27 rhesus macaques, and, alternatively, to the IMGT rhesus macaque database, based on IGHV genes inferred primarily from the genomic sequence of the rheMac10 reference assembly, supplemented with 10 genes from the Mmul_051212 assembly. We found that the use of a gene-restricted database led to overestimations of SHM by up to 5% due to misassignments. The principles described in the current study provide a model for the creation of comprehensive immunoglobulin reference databases from outbred species to ensure accurate gene assignment, lineage tracing and SHM calculations.