High-Quality Library Preparation for NGS-Based Immunoglobulin Germline Gene Inference and Repertoire Expression Analysis

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.

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

Spatial transcriptomics of B cell and T cell receptors reveals lymphocyte clonal dynamics

The spatial distribution of lymphocyte clones within tissues is critical to their development, selection, and expansion. We have developed spatial transcriptomics of variable, diversity, and joining (VDJ) sequences (Spatial VDJ), a method that maps B cell and T cell receptor sequences in human tissue sections. Spatial VDJ captures lymphocyte clones that match canonical B and T cell distributions and amplifies clonal sequences confirmed by orthogonal methods. We found spatial congruency between paired receptor chains, developed a computational framework to predict receptor pairs, and linked the expansion of distinct B cell clones to different tumor-associated gene expression programs. Spatial VDJ delineates B cell clonal diversity and lineage trajectories within their anatomical niche. Thus, Spatial VDJ captures lymphocyte spatial clonal architecture across tissues, providing a platform to harness clonal sequences for therapy.

Multimodal, broadly neutralizing antibodies against SARS-CoV-2 identified by high-throughput native pairing of BCRs from bulk B cells

TruAB Discovery is an approach that integrates cellular immunology, high-throughput immunosequencing, bioinformatics, and computational biology in order to discover naturally occurring human antibodies for prophylactic or therapeutic use. We adapted our previously described pairSEQ technology to pair B cell receptor heavy and light chains of SARS-CoV-2 spike protein-binding antibodies derived from enriched antigen-specific memory B cells and bulk antibody-secreting cells. We identified approximately 60,000 productive, in-frame, paired antibody sequences, from which 2,093 antibodies were selected for functional evaluation based on abundance, isotype and patterns of somatic hypermutation. The exceptionally diverse antibodies included RBD-binders with broad neutralizing activity against SARS-CoV-2 variants, and S2-binders with broad specificity against betacoronaviruses and the ability to block membrane fusion. A subset of these RBD- and S2-binding antibodies demonstrated robust protection against challenge in hamster and mouse models. This high-throughput approach can accelerate discovery of diverse, multifunctional antibodies against any target of interest.

Probabilities of developing HIV-1 bNAb sequence features in uninfected and chronically infected individuals

HIV-1 broadly neutralizing antibodies (bNAbs) are able to suppress viremia and prevent infection. Their induction by vaccination is therefore a major goal. However, in contrast to antibodies that neutralize other pathogens, HIV-1-specific bNAbs frequently carry uncommon molecular characteristics that might prevent their induction. Here, we perform unbiased sequence analyses of B cell receptor repertoires from 57 uninfected and 46 chronically HIV-1- or HCV-infected individuals and learn probabilistic models to predict the likelihood of bNAb development. We formally show that lower probabilities for bNAbs are predictive of higher HIV-1 neutralization activity. Moreover, ranking bNAbs by their probabilities allows to identify highly potent antibodies with superior generation probabilities as preferential targets for vaccination approaches. Importantly, we find equal bNAb probabilities across infected and uninfected individuals. This implies that chronic infection is not a prerequisite for the generation of bNAbs, fostering the hope that HIV-1 vaccines can induce bNAb development in uninfected people.

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.

IGHV allele similarity clustering improves genotype inference from adaptive immune receptor repertoire sequencing data

In adaptive immune receptor repertoire analysis, determining the germline variable (V) allele associated with each T- and B-cell receptor sequence is a crucial step. This process is highly impacted by allele annotations. Aligning sequences, assigning them to specific germline alleles, and inferring individual genotypes are challenging when the repertoire is highly mutated, or sequence reads do not cover the whole V region. Here, we propose an alternative naming scheme for the V alleles, as well as a novel method to infer individual genotypes. We demonstrate the strengths of the two by comparing their outcomes to other genotype inference methods. We validate the genotype approach with independent genomic long-read data. The naming scheme is compatible with current annotation tools and pipelines. Analysis results can be converted from the proposed naming scheme to the nomenclature determined by the International Union of Immunological Societies (IUIS). Both the naming scheme and the genotype procedure are implemented in a freely available R package (PIgLET https://bitbucket.org/yaarilab/piglet). To allow researchers to further explore the approach on real data and to adapt it for their uses, we also created an interactive website (https://yaarilab.github.io/IGHV_reference_book).

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.

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.

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.

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.

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 spike-specific 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.

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.

The evolutionary and functional significance of germline immunoglobulin gene variation

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

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.

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.

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.

Antibody repertoire sequencing analysis

High-throughput sequencing for B cell receptor (BCR) repertoire provides useful insights for the adaptive immune system. With the continuous development of the BCR-seq technology, many efforts have been made to develop methods for analyzing the ever-increasing BCR repertoire data. In this review, we comprehensively outline different BCR repertoire library preparation protocols and summarize three major steps of BCR-seq data analysis, i. e., V(D)J sequence annotation, clonal phylogenetic inference, and BCR repertoire profiling and mining. Different from other reviews in this field, we emphasize background intuition and the statistical principle of each method to help biologists better understand it. Finally, we discuss data mining problems for BCR-seq data and with a highlight on recently emerging multiple-sample analysis.

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.

Individualized VDJ recombination predisposes the available Ig sequence space

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

Novel Allele Detection Tool Benchmark and Application With Antibody Repertoire Sequencing Dataset

Detailed knowledge of the diverse immunoglobulin germline genes is critical for the study of humoral immunity. Hundreds of alleles have been discovered by analyzing antibody repertoire sequencing (Rep-seq or Ig-seq) data via multiple novel allele detection tools (NADTs). However, the performance of these NADTs through antibody sequences with intrinsic somatic hypermutations (SHMs) is unclear. Here, we developed a tool to simulate repertoires by integrating the full spectrum features of an antibody repertoire such as germline gene usage, junctional modification, position-specific SHM and clonal expansion based on 2152 high-quality datasets. We then systematically evaluated these NADTs using both simulated and genuine Ig-seq datasets. Finally, we applied these NADTs to 687 Ig-seq datasets and identified 43 novel allele candidates (NACs) using defined criteria. Twenty-five alleles were validated through findings of other sources. In addition to the NACs detected, our simulation tool, the results of our comparison, and the streamline of this process may benefit further humoral immunity studies via Ig-seq.

Germline polymorphisms and alternative splicing of human immunoglobulin light chain genes

Inference of germline polymorphisms in immunoglobulin genes from B cell receptor repertoires is complicated by somatic hypermutations, sequencing/PCR errors, and by varying length of reference alleles. The light chain inference is particularly challenging owing to large gene duplications and absence of D genes. We analyzed the light chain cDNA sequences from naïve B cell receptor repertoires from 100 individuals. We optimized light chain allele inference by tweaking parameters of the TIgGER functions, extending the germline reference sequences, and establishing mismatch frequency patterns at polymorphic positions to filter out false-positive candidates. We identified 48 previously unreported variants of light chain variable genes. We selected 14 variants for validation and successfully validated 11 by Sanger sequencing. Clustering of light chain 5'UTR, L-PART1, and L-PART2 revealed partial intron retention in 11 kappa and 9 lambda V alleles. Our results provide insight into germline variation in human light chain immunoglobulin loci.

Human immune repertoire in hepatitis B virus infection

Hepatitis B virus (HBV) infection is a public health threat that affects 257 million people worldwide and can progress to liver cirrhosis, liver failure, and hepatocellular carcinoma. The HBV antigen- induced adaptive immune response plays an important role in HBV clearance. Immune repertoire sequencing (IRS) has been used to investigate the molecular mechanisms behind the immune system, find novel ways to treat HBV infection, and evaluate the genetic responses and immune characteristics of individuals infected by HBV or immunized by HBV vaccine. This review summarizes the human immune repertoire analysis methodology, and the application of the IRS in the prediction of HBV infection progression, treatment, and vaccination.

Amplification of a minimally biased antibody repertoire for in vitro display using a universal primer-based amplification method

Immune hosts are valuable sources for antibody discovery. To construct in vitro display antibody libraries from immune repertoires, singleplex or multiplex PCR amplification were employed using primers targeting multiple immunoglobulin genes. However, during this process, the B cell receptor repertoire is distorted due to interactions between multiple target genes and primers. To minimize this alternation, we devised a new method for harvesting immunoglobulin genes and tested its performance in rabbit variable heavy chain (V_H) and variable kappa light chain (V_K) genes. Double-stranded cDNA was synthesized using primers containing V/J gene-specific regions and universal sequence parts for in vitro display. V_H and V_K gene libraries were obtained through subsequent PCR amplification using primers with universal sequences. Next-generation sequencing analysis confirmed that universal PCR libraries had more diverse V_H and V_K clonotypes, and a less biased clonal distribution, than conventional singleplex or multiplex gene-specific PCR libraries.

Rapid development of neutralizing and diagnostic SARS-COV-2 mouse monoclonal antibodies

The need for high-affinity, SARS-CoV-2-specific monoclonal antibodies (mAbs) is critical in the face of the global COVID-19 pandemic, as such reagents can have important diagnostic, research, and therapeutic applications. Of greatest interest is the ~ 300 amino acid receptor binding domain (RBD) within the S1 subunit of the spike protein because of its key interaction with the human angiotensin converting enzyme 2 (hACE2) receptor present on many cell types, especially lung epithelial cells. We report here the development and functional characterization of 29 nM-affinity mouse SARS-CoV-2 mAbs created by an accelerated immunization and hybridoma screening process. Differing functions, including binding of diverse protein epitopes, viral neutralization, impact on RBD-hACE2 binding, and immunohistochemical staining of infected lung tissue, were correlated with variable gene usage and sequence.

Polyfunctional Tier 2-Neutralizing Antibodies Cloned following HIV-1 Env Macaque Immunization Mirror Native Antibodies in a Human Donor

Vaccine efforts to combat HIV are challenged by the global diversity of viral strains and shielding of neutralization epitopes on the viral envelope glycoprotein trimer. Even so, the isolation of broadly neutralizing Abs from infected individuals suggests the potential for eliciting protective Abs through vaccination. This study reports a panel of 58 mAbs cloned from a rhesus macaque (Macaca mulatta) immunized with envelope glycoprotein immunogens curated from an HIV-1 clade C-infected volunteer. Twenty mAbs showed neutralizing activity, and the strongest neutralizer displayed 92% breadth with a median IC50 of 1.35 μg/ml against a 13-virus panel. Neutralizing mAbs predominantly targeted linear epitopes in the V3 region in the cradle orientation (V3C) with others targeting the V3 ladle orientation (V3L), the CD4 binding site (CD4bs), C1, C4, or gp41. Nonneutralizing mAbs bound C1, C5, or undetermined conformational epitopes. Neutralization potency strongly correlated with the magnitude of binding to infected primary macaque splenocytes and to the level of Ab-dependent cellular cytotoxicity, but did not predict the degree of Ab-dependent cellular phagocytosis. Using an individualized germline gene database, mAbs were traced to 23 of 72 functional IgHV alleles. Neutralizing V3C Abs displayed minimal nucleotide somatic hypermutation in the H chain V region (3.77%), indicating that relatively little affinity maturation was needed to achieve in-clade neutralization breadth. Overall, this study underscores the polyfunctional nature of vaccine-elicited tier 2-neutralizing V3 Abs and demonstrates partial reproduction of the human donor's humoral immune response through nonhuman primate vaccination.

Poorly Expressed Alleles of Several Human Immunoglobulin Heavy Chain Variable Genes are Common in the Human Population

Extensive diversity has been identified in the human heavy chain immunoglobulin locus, including allelic variation, gene duplication, and insertion/deletion events. Several genes have been suggested to be deleted in many haplotypes. Such findings have commonly been based on inference of the germline repertoire from data sets covering antibody heavy chain encoding transcripts. The inference process operates under conditions that may limit identification of genes transcribed at low levels. The presence of rare transcripts that would indicate the existence of poorly expressed alleles in haplotypes that otherwise appear to have deleted these genes has been assessed in the present study. Alleles IGHV1-2*05, IGHV1-3*02, IGHV4-4*01, and IGHV7-4-1*01 were all identified as being expressed from multiple haplotypes, but only at low levels, haplotypes that by inference often appeared not to express these genes at all. These genes are thus not as commonly deleted as previously thought. An assessment of the 5' untranslated region (up to and including the TATA-box), the signal peptide-encoding part of the gene, and the 3'-heptamer suggests that the alleles have no or minimal sequence difference in these regions in comparison to highly expressed alleles. This suggest that they may be able to participate in immunoglobulin gene rearrangement, transcription and translation. However, all four poorly expressed alleles harbor unusual sequence variants within their coding region that may compromise the functionality of the encoded products, thereby limiting their incorporation into the immunoglobulin repertoire. Transcripts based on IGHV7-4-1*01 that had undergone somatic hypermutation and class switch had mutated the codon that encoded the unusual residue in framework region 3 (cysteine 92; located far from the antigen binding site). This finding further supports the poor compatibility of this unusual residue in a fully functional protein product. Indications of a linkage disequilibrium were identified as IGHV1-2*05 and IGHV4-4*01 co-localized to the same haplotypes. Furthermore, transcripts of two of the poorly expressed alleles (IGHV1-3*02 and IGHV4-4*01) mostly do not encode in-frame, functional products, suggesting that these alleles might be essentially non-functional. It is proposed that the functionality status of immunoglobulin genes should also include assessment of their ability to encode functional protein products.

Rhesus and cynomolgus macaque immunoglobulin heavy-chain genotyping yields comprehensive databases of germline VDJ alleles

Definition of the specific germline immunoglobulin (Ig) alleles present in an individual is a critical first step to delineate the ontogeny and evolution of antigen-specific antibody responses. Rhesus and cynomolgus macaques are important animal models for pre-clinical studies, with four main sub-groups being used: Indian- and Chinese-origin rhesus macaques and Mauritian and Indonesian cynomolgus macaques. We applied the (Ig) gene inference tool IgDiscover and performed extensive Sanger sequencing-based genomic validation to define germline VDJ alleles in these 4 sub-groups, comprising 45 macaques in total. There was allelic overlap between Chinese- and Indian-origin rhesus macaques and also between the two macaque species, which is consistent with substantial admixture. The island-restricted Mauritian cynomolgus population displayed the lowest number of alleles of the sub-groups, yet maintained high individual allelic diversity. These comprehensive databases of germline IGH alleles for rhesus and cynomolgus macaques provide a resource toward the study of B cell responses in these important pre-clinical models.

Immunoglobulin germline gene variation and its impact on human disease

Immunoglobulins (Ig) play an important role in the immune system both when expressed as antigen receptors on the cell surface of B cells and as antibodies secreted into extracellular fluids. The advent of high-throughput sequencing methods has enabled the investigation of human Ig repertoires at unprecedented depth. This has led to the discovery of many previously unreported germline Ig alleles. Moreover, it is becoming clear that convergent and stereotypic antibody responses are common where different individuals recognise defined antigenic epitopes with the use of the same Ig V genes. Thus, germline V gene variation is increasingly being linked to the differential capacity of generating an effective immune response, which might lead to varying disease susceptibility. Here, we review recent evidence of how germline variation in Ig genes impacts the Ig repertoire and its subsequent effects on the adaptive immune response in vaccination, infection, and autoimmunity.

Emerging Strategies for Therapeutic Antibody Discovery from Human B Cells

Monoclonal antibodies from human sources are being increasingly recognized as valuable options in many therapeutic areas. These antibodies can show exquisite specificity and high potency while maintaining a desirable safety profile, having been matured and tolerized within human patients. However, the discovery of these antibodies presents important challenges, since the B cells encoding therapeutic antibodies can be rare in a typical blood draw and are short-lived ex vivo. Furthermore, the unique pairing of V_H and V_L domains in each B cell contributes to specificity and function; therefore, maintaining antibody chain pairing presents a throughput limitation. This work will review the various approaches aimed at addressing these challenges with an eye to next-generation methods for high-throughput discovery from the human B-cell repertoire.

Polymorphisms in human immunoglobulin heavy chain variable genes and their upstream regions

Germline variations in immunoglobulin genes influence the repertoire of B cell receptors and antibodies, and such polymorphisms may impact disease susceptibility. However, the knowledge of the genomic variation of the immunoglobulin loci is scarce. Here, we report 25 potential novel germline IGHV alleles as inferred from rearranged naïve B cell cDNA repertoires of 98 individuals. Thirteen novel alleles were selected for validation, out of which ten were successfully confirmed by targeted amplification and Sanger sequencing of non-B cell DNA. Moreover, we detected a high degree of variability upstream of the V-REGION in the 5′UTR, L-PART1 and L-PART2 sequences, and found that identical V-REGION alleles can differ in upstream sequences. Thus, we have identified a large genetic variation not only in the V-REGION but also in the upstream sequences of IGHV genes. Our findings provide a new perspective for annotating immunoglobulin repertoire sequencing data.

Systematic Profiling of Full-Length Ig and TCR Repertoire Diversity in Rhesus Macaque through Long Read Transcriptome Sequencing

The diversity of Ig and TCR repertoires is a focal point of immunological studies. Rhesus macaques (Macaca mulatta) are key for modeling human immune responses, placing critical importance on the accurate annotation and quantification of their Ig and TCR repertoires. However, because of incomplete reference resources, the coverage and accuracy of the traditional targeted amplification strategies for profiling rhesus Ig and TCR repertoires are largely unknown. In this study, using long read sequencing, we sequenced four Indian-origin rhesus macaque tissues and obtained high-quality, full-length sequences for over 6000 unique Ig and TCR transcripts, without the need for sequence assembly. We constructed, to our knowledge, the first complete reference set for the constant regions of all known isotypes and chain types of rhesus Ig and TCR repertoires. We show that sequence diversity exists across the entire variable regions of rhesus Ig and TCR transcripts. Consequently, existing strategies using targeted amplification of rearranged variable regions comprised of V(D)J gene segments miss a significant fraction (27-53% and 42-49%) of rhesus Ig/TCR diversity. To overcome these limitations, we designed new rhesus-specific assays that remove the need for primers conventionally targeting variable regions and allow single cell level Ig and TCR repertoire analysis. Our improved approach will enable future studies to fully capture rhesus Ig and TCR repertoire diversity and is applicable for improving annotations in any model organism.

Analysis of a SARS-CoV-2-Infected Individual Reveals Development of Potent Neutralizing Antibodies with Limited Somatic Mutation

Antibody responses develop following SARS-CoV-2 infection, but little is known about their epitope specificities, clonality, binding affinities, epitopes, and neutralizing activity. We isolated B cells specific for the SARS-CoV-2 envelope glycoprotein spike (S) from a COVID-19-infected subject 21 days after the onset of clinical disease. 45 S-specific monoclonal antibodies were generated. They had undergone minimal somatic mutation with limited clonal expansion, and three bound the receptor-binding domain (RBD). Two antibodies neutralized SARS-CoV-2. The most potent antibody bound the RBD and prevented binding to the ACE2 receptor, while the other bound outside the RBD. Thus, most anti-S antibodies that were generated in this patient during the first weeks of COVID-19 infection were non-neutralizing and target epitopes outside the RBD. Antibodies that disrupt the SARS-CoV-2 S-ACE2 interaction can potently neutralize the virus without undergoing extensive maturation. Such antibodies have potential preventive and/or therapeutic potential and can serve as templates for vaccine design.

•

Early B cell responses to SARS-CoV-2 spike protein are analyzed from a COVID-19 patient

•

Most antibodies target non-neutralizing epitopes outside the RBD

•

A potent neutralizing mAb blocks the interaction of the S protein with ACE2

•

Neutralizing antibodies are minimally mutated

Characterization of neutralizing antibodies from a SARS-CoV-2 infected individual

B cells specific for the SARS-CoV-2 S envelope glycoprotein spike were isolated from a COVID-19-infected subject using a stabilized spike-derived ectodomain (S2P) twenty-one days post-infection. Forty-four S2P-specific monoclonal antibodies were generated, three of which bound to the receptor binding domain (RBD). The antibodies were minimally mutated from germline and were derived from different B cell lineages. Only two antibodies displayed neutralizing activity against SARS-CoV-2 pseudo-virus. The most potent antibody bound the RBD in a manner that prevented binding to the ACE2 receptor, while the other bound outside the RBD. Our study indicates that the majority of antibodies against the viral envelope spike that were generated during the first weeks of COVID-19 infection are non-neutralizing and target epitopes outside the RBD. Antibodies that disrupt the SARS-CoV-2 spike-ACE2 interaction can potently neutralize the virus without undergoing extensive maturation. Such antibodies have potential preventive/therapeutic potential and can serve as templates for vaccine-design.

VRC34-Antibody Lineage Development Reveals How a Required Rare Mutation Shapes the Maturation of a Broad HIV-Neutralizing Lineage

Rare mutations have been proposed to restrict the development of broadly neutralizing antibodies against HIV-1, but this has not been explicitly demonstrated. We hypothesized that such rare mutations might be identified by comparing broadly neutralizing and non-broadly neutralizing branches of an antibody-developmental tree. Because sequences of antibodies isolated from the fusion peptide (FP)-targeting VRC34-antibody lineage suggested it might be suitable for such rare mutation analysis, we carried out next-generation sequencing (NGS) on B cell transcripts from donor N123, the source of the VRC34 lineage, and functionally and structurally characterized inferred intermediates along broadly neutralizing and poorly neutralizing developmental branches. The broadly neutralizing VRC34.01 branch required the rare heavy-chain mutation Y33P to bind FP, whereas the early bifurcated VRC34.05 branch did not require this rare mutation and evolved less breadth. Our results demonstrate how a required rare mutation can restrict development and shape the maturation of a broad HIV-1-neutralizing antibody lineage.

Understanding the human antibody repertoire

The origins of the various elements in the human antibody repertoire have been and still are subject to considerable uncertainty. Uncertainty in respect of whether the various elements have always served a specific defense function or whether they were co-opted from other organismal roles to form a crude naïve repertoire that then became more complex as combinatorial mechanisms were added. Estimates of the current size of the human antibody naïve repertoire are also widely debated with numbers anywhere from 10 million members, based on experimentally derived numbers, to in excess of one thousand trillion members or more, based on the different sequences derived from theoretical combinatorial calculations. There are questions that are relevant at both ends of this number spectrum. At the lower bound it could be questioned whether this is an insufficient repertoire size to counter all the potential antigen-bearing pathogens. At the upper bound the question is rather simpler: How can any individual interrogate such an astronomical number of antibody-bearing B cells in a timeframe that is meaningful? This review evaluates the evolutionary aspects of the adaptive immune system, the calculations that lead to the large repertoire estimates, some of the experimental evidence pointing to a more restricted repertoire whose variation appears to derive from convergent 'structure and specificity features', and includes a theoretical model that seems to support it. Finally, a solution that may reconcile the size difference anomaly, which is still a hot subject of debate, is suggested.

A linear-amplification VDJ-seq technique for quantification of immunoglobulin and T cell receptor diversity

The V(D)J recombination is essential for generating a highly diverse repertoire of antigen receptors expressed on T and B lymphocytes. Here, we developed a linear-amplification VDJ-seq technique for quantifying V(D)J recombination of antigen receptor genes. This technique takes advantage of linear amplification using in vitro transcription and reverse transcription to avoid bias generated by the PCR amplification of low copy number of target DNA. The unrearranged alleles are removed by in vitro cleavage with the CRISPR-Cas9 system. The linear-amplification VDJ-seq assay was applied in quantification of the Vκ-Jκ recombination of the mouse Igκ gene with Jκ capture primers. The Jκ genes were detected in 95.86% of clean reads with more than half containing the Vκ gene, indicating high specificity of capturing and amplification. We also applied this approach to quantify the usage of Jα within the Trav12 gene family of the Tcra gene.