Individual variation in the germline Ig gene repertoire inferred from variable region gene rearrangements

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.

The multilevel extensive diversity across the cynomolgus macaque captured by ultra-deep adaptive immune receptor repertoire sequencing

The extent to which AIRRs differ among and within individuals remains elusive. Via ultra-deep repertoire sequencing of 22 and 25 tissues in three cynomolgus macaques, respectively, we identified 84 and 114 novel IGHV and TRBV alleles, confirming 72 (85.71%) and 100 (87.72%) of them. The heterogeneous V gene usage patterns were influenced, in turn, by genetics, isotype (for BCRs only), tissue group, and tissue. A higher proportion of intragroup shared clones in the intestinal tissues than those in other tissues suggests a close intra-intestinal adaptive immunity network. Significantly higher mutation burdens in the public clones and the inter-tissue shared IgM and IgD clones indicate that they might target the shared antigens. This study reveals the extensive heterogeneity of the AIRRs at various levels and has broad fundamental and clinical implications. The data generated here will serve as an invaluable resource for future studies on adaptive immunity in health and diseases.

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.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study reveals the importance of light-chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrate that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to the HA stem, are incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light-chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

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.

Single-Step IGHV Next-Generation Sequencing Detects Clonality and Somatic Hypermutation in Lymphoid Malignancies: A Phase III Diagnostic Accuracy Study

BACKGROUND

Multiplex PCR based on consensus primers followed by capillary electrophoresis and Sanger sequencing are considered as the gold standard method for the evaluation of clonality and somatic hypermutation in lymphoid malignancies. As an alternative, the next-generation sequencing (NGS) of immune receptor genes has recently been proposed as a solution, due to being highly effective and sensitive. Here, we designed a phase III diagnostic accuracy study intended to compare the current gold standard methods versus the first commercially available NGS approaches for testing immunoglobulin heavy chain gene rearrangements.

METHODS

We assessed IGH rearrangements in 68 samples by means of both the NGS approach (LymphoTrack® IGH assay, and LymphoTrack® IGH somatic hypermutation assay, run on Illumina MiSeq) and capillary electrophoresis/Sanger sequencing to assess clonality and somatic hypermutations (SHM).

RESULTS

In comparison to the routine capillary-based analysis, the NGS clonality assay had an overall diagnostic accuracy of 96% (63/66 cases). Other studied criteria included sensitivity (95%), specificity (100%), positive predictive value (100%) and negative predictive value (75%). In discrepant cases, the NGS results were confirmed by a different set of primers that provided coverage of the IGH leader sequence. Furthermore, there was excellent agreement of the SHM determination with both the LymphoTrack® FR1 and leader assays when compared to the Sanger sequencing analysis (84%), with NGS able to assess the SHM rate even in cases where the conventional approach failed.

CONCLUSION

Overall, conventional Sanger sequencing and next-generation-sequencing-based clonality and somatic hypermutation analyses gave comparable results. For future use in a routine diagnostic workflow, NGS-based approaches should be evaluated prospectively and an analysis of cost-effectiveness should be performed.

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

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

Variation in the antibody response has been linked to differential outcomes in disease, and suboptimal vaccine and therapeutic responsiveness, the determinants of which have not been fully elucidated. Countering models that presume antibodies are generated largely by stochastic processes, we demonstrate that polymorphisms within the immunoglobulin heavy chain locus (IGH) impact the naive and antigen-experienced antibody repertoire, indicating that genetics predisposes individuals to mount qualitatively and quantitatively different antibody responses. We pair recently developed long-read genomic sequencing methods with antibody repertoire profiling to comprehensively resolve IGH genetic variation, including novel structural variants, single nucleotide variants, and genes and alleles. We show that IGH germline variants determine the presence and frequency of antibody genes in the expressed repertoire, including those enriched in functional elements linked to V(D)J recombination, and overlapping disease-associated variants. These results illuminate the power of leveraging IGH genetics to better understand the regulation, function, and dynamics of the antibody response in disease.

Stringent and complex sequence constraints of an IGHV1-69 broadly neutralizing antibody to influenza HA stem

IGHV1-69 is frequently utilized by broadly neutralizing influenza antibodies to the hemagglutinin (HA) stem. These IGHV1-69 HA stem antibodies have diverse complementarity-determining region (CDR) H3 sequences. Besides, their light chains have minimal to no contact with the epitope. Consequently, sequence determinants that confer IGHV1-69 antibodies with HA stem specificity remain largely elusive. Using high-throughput experiments, this study revealed the importance of light chain sequence for the IGHV1-69 HA stem antibody CR9114, which is the broadest influenza antibody known to date. Moreover, we demonstrated that the CDR H3 sequences from many other IGHV1-69 antibodies, including those to HA stem, were incompatible with CR9114. Along with mutagenesis and structural analysis, our results indicate that light chain and CDR H3 sequences coordinately determine the HA stem specificity of IGHV1-69 antibodies. Overall, this work provides molecular insights into broadly neutralizing antibody responses to influenza virus, which have important implications for universal influenza vaccine development.

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.

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.

Comparing T cell receptor repertoires using optimal transport

The complexity of entire T cell receptor (TCR) repertoires makes their comparison a difficult but important task. Current methods of TCR repertoire comparison can incur a high loss of distributional information by considering overly simplistic sequence- or repertoire-level characteristics. Optimal transport methods form a suitable approach for such comparison given some distance or metric between values in the sample space, with appealing theoretical and computational properties. In this paper we introduce a nonparametric approach to comparing empirical TCR repertoires that applies the Sinkhorn distance, a fast, contemporary optimal transport method, and a recently-created distance between TCRs called TCRdist. We show that our methods identify meaningful differences between samples from distinct TCR distributions for several case studies, and compete with more complicated methods despite minimal modeling assumptions and a simpler pipeline.

Mutational patterns along different evolution paths of follicular lymphoma

Follicular lymphoma (FL) is an indolent disease, characterized by a median life expectancy of 18-20 years and by intermittent periods of relapse and remission. FL frequently transforms into the more aggressive diffuse large B cell lymphoma (t-FL). In previous studies, the analysis of immunoglobulin heavy chain variable region (IgHV) genes in sequential biopsies from the same patient revealed two different patterns of tumor clonal evolution: direct evolution, through acquisition of additional IgHV mutations over time, or divergent evolution, in which lymphoma clones from serial biopsies independently develop from a less-mutated common progenitor cell (CPC). Our goal in this study was to characterize the somatic hypermutation (SHM) patterns of IgHV genes in sequential FL samples from the same patients, and address the question of whether the mutation mechanisms (SHM targeting, DNA repair or both), or selection forces acting on the tumor clones, were different in FL samples compared to healthy control samples, or in late relapsed/transformed FL samples compared to earlier ones. Our analysis revealed differences in the distribution of mutations from each of the nucleotides when tumor and non-tumor clones were compared, while FL and transformed FL (t-FL) tumor clones displayed similar mutation distributions. Lineage tree measurements suggested that either initial clone affinity or selection thresholds were lower in FL samples compared to controls, but similar between FL and t-FL samples. Finally, we observed that both FL and t-FL tumor clones tend to accumulate larger numbers of potential N-glycosylation sites due to the introduction of new SHM. Taken together, these results suggest that transformation into t-FL, in contrast to initial FL development, is not associated with any major changes in DNA targeting or repair, or the selection threshold of the tumor clone.

Potent monoclonal antibodies neutralize Omicron sublineages and other SARS-CoV-2 variants

The emergence and global spread of the SARS-CoV-2 Omicron variants, which carry an unprecedented number of mutations, raise serious concerns due to the reduced efficacy of current vaccines and resistance to therapeutic antibodies. Here, we report the generation and characterization of two potent human monoclonal antibodies, NA8 and NE12, against the receptor-binding domain of the SARS-CoV-2 spike protein. NA8 interacts with a highly conserved region and has a breadth of neutralization with picomolar potency against the Beta variant and the Omicron BA.1 and BA.2 sublineages and nanomolar potency against BA.2.12.1 and BA.4. Combination of NA8 and NE12 retains potent neutralizing activity against the major SARS-CoV-2 variants of concern. Cryo-EM analysis provides the structural basis for the broad and complementary neutralizing activity of these two antibodies. We confirm the in vivo protective and therapeutic efficacies of NA8 and NE12 in the hamster model. These results show that broad and potent human antibodies can overcome the continuous immune escape of evolving SARS-CoV-2 variants.

Characterization of prefusion-F-specific antibodies elicited by natural infection with human metapneumovirus

Human metapneumovirus (hMPV) is a major cause of acute respiratory infections in infants and older adults, for which no vaccines or therapeutics are available. The viral fusion (F) glycoprotein is required for entry and is the primary target of neutralizing antibodies; however, little is known about the humoral immune response generated from natural infection. Here, using prefusion-stabilized F proteins to interrogate memory B cells from two older adults, we obtain over 700 paired non-IgM antibody sequences representing 563 clonotypes, indicative of a highly polyclonal response. Characterization of 136 monoclonal antibodies reveals broad recognition of the protein surface, with potently neutralizing antibodies targeting each antigenic site. Cryo-EM studies further reveal two non-canonical sites and the molecular basis for recognition of the apex of hMPV F by two prefusion-specific neutralizing antibodies. Collectively, these results provide insight into the humoral response to hMPV infection in older adults and will help guide vaccine development.

Immunoglobulin genes, reproductive isolation and vertebrate speciation

Reproductive isolation drives the formation of new species, and many genes contribute to this through Dobzhansky–Muller incompatibilities (DMIs). These incompatibilities occur when gene divergence affects loci encoding interacting products such as receptors and their ligands. We suggest here that the nature of vertebrate immunoglobulin (IG) genes must make them prone to DMIs. The genes of these complex loci form functional genes through the process of recombination, giving rise to a repertoire of heterodimeric receptors of incredible diversity. This repertoire, within individuals and within species, must defend against pathogens but must also avoid pathogenic self‐reactivity. We suggest that this avoidance of autoimmunity is only achieved through a coordination of evolution between heavy‐ and light‐chain genes, and between these genes and the rest of the genome. Without coordinated evolution, the hybrid offspring of two diverging populations will carry a heavy burden of DMIs, resulting in a loss of fitness. Critical incompatibilities could manifest as incompatibilities between a mother and her divergent offspring. During fetal development, biochemical differences between the parents of hybrid offspring could make their offspring a target of the maternal immune system. This hypothesis was conceived in the light of recent insights into the population genetics of IG genes. This has suggested that antibody genes are probably as susceptible to evolutionary forces as other parts of the genome. Further repertoire studies in human and nonhuman species should now help determine whether antibody genes have been part of the evolutionary forces that drive the development of species.

B cells in the balance: Offsetting self-reactivity avoidance with protection against foreign

Antibodies are theoretically limitless in their diversity and specificity to foreign antigens; however they are constrained by the need to avoid binding to self. Germinal centers (GC) allow diversification and maturation of the antibody response towards the foreign antigen. While self-tolerance mechanisms controlling self-reactivity during B cell maturation are well recognized, the mechanisms by which GCs balance self-tolerance and foreign binding especially in the face of cross-reactivity between self and foreign, remain much less well defined. In this review we explore the extent to which GC self-tolerance restricts affinity maturation. We present studies suggesting that the outcome is situationally dependent, affected by affinity and avidity to self-antigen, and the extent to which self-binding and foreign-binding are interdependent. While auto-reactive GC B cells can mutate away from self while maturing towards the foreign antigen, if no mutational trajectories allow for self-reactive redemption, self-tolerance prevails and GC responses to the foreign pathogen are restricted, except when self-tolerance checkpoints are relaxed. Finally, we consider whether polyreactivity is subject to the same level of restriction in GC responses, especially if polyreactivity is linked to an increase in foreign protection, as occurs in certain broadly neutralizing antibodies. Overall, the outcomes for GC B cells that bind self-antigen can range from redemption, transient relaxation in self-tolerance or restriction of the antibody response to the foreign pathogen.

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.

SARS-CoV-2 Beta variant infection elicits potent lineage-specific and cross-reactive antibodies

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Beta variant of concern (VOC) resists neutralization by major classes of antibodies from COVID-19 patients and vaccinated individuals. In this study, serum of Beta-infected patients revealed reduced cross-neutralization of wild-type virus. From these patients, we isolated Beta-specific and cross-reactive receptor-binding domain (RBD) antibodies. The Beta-specificity results from recruitment of VOC-specific clonotypes and accommodation of mutations present in Beta and Omicron into a major antibody class that is normally sensitive to these mutations. The Beta-elicited cross-reactive antibodies share genetic and structural features with wild type-elicited antibodies, including a public VH1-58 clonotype that targets the RBD ridge. These findings advance our understanding of the antibody response to SARS-CoV-2 shaped by antigenic drift, with implications for design of next-generation vaccines and therapeutics.

Extremely potent monoclonal antibodies neutralize Omicron and other SARS-CoV-2 variants

The ongoing coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has triggered a devastating global health, social and economic crisis. The RNA nature and broad circulation of this virus facilitate the accumulation of mutations, leading to the continuous emergence of variants of concern with increased transmissibility or pathogenicity ¹ . This poses a major challenge to the effectiveness of current vaccines and therapeutic antibodies ^{1, 2} . Thus, there is an urgent need for effective therapeutic and preventive measures with a broad spectrum of action, especially against variants with an unparalleled number of mutations such as the recently emerged Omicron variant, which is rapidly spreading across the globe ³ . Here, we used combinatorial antibody phage-display libraries from convalescent COVID-19 patients to generate monoclonal antibodies against the receptor-binding domain of the SARS-CoV-2 spike protein with ultrapotent neutralizing activity. One such antibody, NE12, neutralizes an early isolate, the WA-1 strain, as well as the Alpha and Delta variants with half-maximal inhibitory concentrations at picomolar level. A second antibody, NA8, has an unusual breadth of neutralization, with picomolar activity against both the Beta and Omicron variants. The prophylactic and therapeutic efficacy of NE12 and NA8 was confirmed in preclinical studies in the golden Syrian hamster model. Analysis by cryo-EM illustrated the structural basis for the neutralization properties of NE12 and NA8. Potent and broadly neutralizing antibodies against conserved regions of the SARS-CoV-2 spike protein may play a key role against future variants of concern that evade immune control.

T cell receptor beta germline variability is revealed by inference from repertoire data

BACKGROUND

T and B cell receptor (TCR, BCR) repertoires constitute the foundation of adaptive immunity. Adaptive immune receptor repertoire sequencing (AIRR-seq) is a common approach to study immune system dynamics. Understanding the genetic factors influencing the composition and dynamics of these repertoires is of major scientific and clinical importance. The chromosomal loci encoding for the variable regions of TCRs and BCRs are challenging to decipher due to repetitive elements and undocumented structural variants.

METHODS

To confront this challenge, AIRR-seq-based methods have recently been developed for B cells, enabling genotype and haplotype inference and discovery of undocumented alleles. However, this approach relies on complete coverage of the receptors' variable regions, whereas most T cell studies sequence a small fraction of that region. Here, we adapted a B cell pipeline for undocumented alleles, genotype, and haplotype inference for full and partial AIRR-seq TCR data sets. The pipeline also deals with gene assignment ambiguities, which is especially important in the analysis of data sets of partial sequences.

RESULTS

From the full and partial AIRR-seq TCR data sets, we identified 39 undocumented polymorphisms in T cell receptor Beta V (TRBV) and 31 undocumented 5 ' UTR sequences. A subset of these inferences was also observed using independent genomic approaches. We found that a single nucleotide polymorphism differentiating between the two documented T cell receptor Beta D2 (TRBD2) alleles is strongly associated with dramatic changes in the expressed repertoire.

CONCLUSIONS

We reveal a rich picture of germline variability and demonstrate how a single nucleotide polymorphism dramatically affects the composition of the whole repertoire. Our findings provide a basis for annotation of TCR repertoires for future basic and clinical studies.

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.

Computational Inference, Validation, and Analysis of 5'UTR-Leader Sequences of Alleles of Immunoglobulin Heavy Chain Variable Genes

Upstream and downstream sequences of immunoglobulin genes may affect the expression of such genes. However, these sequences are rarely studied or characterized in most studies of immunoglobulin repertoires. Inference from large, rearranged immunoglobulin transcriptome data sets offers an opportunity to define the upstream regions (5'-untranslated regions and leader sequences). We have now established a new data pre-processing procedure to eliminate artifacts caused by a 5'-RACE library generation process, reanalyzed a previously studied data set defining human immunoglobulin heavy chain genes, and identified novel upstream regions, as well as previously identified upstream regions that may have been identified in error. Upstream sequences were also identified for a set of previously uncharacterized germline gene alleles. Several novel upstream region variants were validated, for instance by their segregation to a single haplotype in heterozygotic subjects. SNPs representing several sequence variants were identified from population data. Finally, based on the outcomes of the analysis, we define a set of testable hypotheses with respect to the placement of particular alleles in complex IGHV locus haplotypes, and discuss the evolutionary relatedness of particular heavy chain variable genes based on sequences of their upstream regions.

Developmental and comparative immunology single-cell transcriptome analysis of the B-cell repertoire reveals the usage of immunoglobulins in the gray short-tailed opossum (Monodelphis domestica)

B-cells are key to humoral immunity, are found in multiple lymphoid organs, and have the unique ability to mediate the production of antigen-specific antibodies in the presence of pathogens. The marsupial immunoglobulin (Ig) heavy (H) chain locus encodes four constant region isotypes, IgA, IgG, IgM and IgE, but no IgD, and there are two light (L) chain isotypes, lambda (Igλ) and kappa (Igκ). To gain an understanding of the marsupial humoral immune system, B-cell transcriptomes generated by single-cell RNA sequencing from gray short-tailed opossum (Monodelphis domestica) splenocytes, and peripheral blood mononuclear cells were analysed. The cells used were from a single unimmunized animal and the majority of B-cells were transcribing IgM heavy chains. The ratio of Ig light chain use was roughly 2:1, Igλ:Igκ in this individual. This was not predicted due to Igκ being the more complex of the two L chain loci. The variable (V) gene segment pairs used in individual B-cells confirm greater diversity provided by the L chain V. This study is the first to report on using single cell analysis to investigate Ig repertoires in a marsupial and confirms a number of prior hypothesis, as well as revealing some surprises.

Vaccine genetics of IGHV1-2 VRC01-class broadly neutralizing antibody precursor naïve human B cells

A successful HIV vaccine eliciting broadly neutralizing antibodies (bnAbs) must overcome the hurdle of being able to activate naive precursor B cells encoding features within their germline B cell receptors (BCR) that allow recognition of broadly neutralizing epitopes. Knowledge of whether bnAb precursor B cells are circulating at sufficient frequencies within individuals in communities heavily impacted by HIV may be important. Using a germline-targeting eOD-GT8 immunogen and high-throughput droplet-based single-cell BCR sequencing, we demonstrate that large numbers of paired BCR sequences from multiple donors can be efficiently screened to elucidate precursor frequencies of rare, naive VRC01-class B cells. Further, we analyzed IGHV1-2 allelic usage among three different cohorts; we find that IGHV1-2 alleles traditionally thought to be incompatible with VRC01-class responses are relatively common in various human populations and that germline variation within IGHV1-2 associates with gene usage frequencies in the naive BCR repertoire.

Development of a Prognostic Nomogram for Acute Myeloid Leukemia on IGHD Gene Family

Purpose

Acute myeloid leukaemia (AML) is a common haematological disease in adults. The overall survival (OS) remains unsatisfactory. It is critical to identify potential prognostic biomarkers and develop a nomogram that predicts overall survival in patients with AML.

Patients and Methods

We used gene expression dataset and clinical data from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) to identify differential expression analysis, survival analysis, and prognostic value of IGHD gene family (IGHDs) in AML patients. A risk score model was built through Lasso analysis and multivariate Cox regression. We also developed a nomogram and evaluated its accuracy with Harrell’s Harmony Index (C-index) and calibration curve. Last, the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database was used for external validation.

Results

IGHD1-20 mRNA expression level was an independent prognostic factor for patients with AML by multivariate analysis. After Lasso analysis and multivariate Cox regression, we constructed a 3-gene model (IGHD1-1, IGHD1-20, IGHD3-16) associated with OS in AML. Risk score and age were validated as independent risk factors for prognosis and were used to build a nomogram. The C index and calibration curve results show that its ability to predict 1-year, 3-year and 5-year overall survival is accurate.

Conclusion

The mRNA level of IGHDs was increased in AML patients. IGHD1-20 was an independent risk factor for OS in AML patients. The IGHDs risk model (IGHD1-1, IGHD1-20, IGHD3-16) relates to the OS of AML patients. The nomogram, including risk score and age, can conveniently and effectively predict the overall survival rate of patients.

Structural and functional ramifications of antigenic drift in recent SARS-CoV-2 variants

Neutralizing antibodies (nAbs) elicited against the receptor binding site (RBS) of the spike protein of wild-type severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are generally less effective against recent variants of concern. RBS residues Glu484, Lys417, and Asn501 are mutated in variants first described in South Africa (B.1.351) and Brazil (P.1). We analyzed their effects on angiotensin-converting enzyme 2 binding, as well as the effects of two of these mutations (K417N and E484K) on nAbs isolated from COVID-19 patients. Binding and neutralization of the two most frequently elicited antibody families (IGHV3-53/3-66 and IGHV1-2), which can both bind the RBS in alternative binding modes, are abrogated by K417N, E484K, or both. These effects can be structurally explained by their extensive interactions with RBS nAbs. However, nAbs to the more conserved, cross-neutralizing CR3022 and S309 sites were largely unaffected. The results have implications for next-generation vaccines and antibody therapies.

Population matched (pm) germline allelic variants of immunoglobulin (IG) loci: Relevance in infectious diseases and vaccination studies in human populations

Immunoglobulin (IG) loci harbor inter-individual allelic variants in many different germline IG variable, diversity and joining genes of the IG heavy (IGH), kappa (IGK) and lambda (IGL) loci, which together form the genetic basis of the highly diverse antigen-specific B-cell receptors. These allelic variants can be shared between or be specific to human populations. The current immunogenetics resources gather the germline alleles, however, lack the population specificity of the alleles which poses limitations for disease-association studies related to immune responses in different human populations. Therefore, we systematically identified germline alleles from 26 different human populations around the world, profiled by "1000 Genomes" data. We identified 409 IGHV, 179 IGKV, and 199 IGLV germline alleles supported by at least seven haplotypes. The diversity of germline alleles is the highest in Africans. Remarkably, the variants in the identified novel alleles show strikingly conserved patterns, the same as found in other IG databases, suggesting over-time evolutionary selection processes. We could relate the genetic variants to population-specific immune responses, e.g. IGHV1-69 for flu in Africans. The population matched IG (pmIG) resource will enhance our understanding of the SHM-related B-cell receptor selection processes in (infectious) diseases and vaccination within and between different human populations.

Human transitional and IgMlow mature naïve B cells preserve permissive B-cell receptors

The level of immunoglobulin M (IgM) displayed on the surface of peripheral blood B cells exhibits a broad dynamic range and has been associated with both development and selection. To determine whether IgM surface expression associates with distinct immunoglobulin heavy‐chain (IGH) repertoire properties, we performed deep IgM sequencing of peripheral blood transitional and mature naïve B cells in the upper and lower quartiles of surface IgM expression for 12 healthy donors. Mature naïve B cells within the lowest quartile for surface IgM expression displayed more diverse IGH features including increased complementarity‐determining region 3 length, IGHJ6 segment usage and aromatic amino acids compared with mature naïve B cells with high surface IgM. There were no differences between IGH repertoires for transitional B cells with high or low surface IgM. These findings suggest that a selection checkpoint during progression of transitional to mature naïve B cells reduces the breadth of the IGH repertoire among high surface IgM B cells but that diversity is preserved in B cells expressing low levels of surface IgM.

Prevalent, protective, and convergent IgG recognition of SARS-CoV-2 non-RBD spike epitopes

The molecular composition and binding epitopes of the immunoglobulin G (IgG) antibodies that circulate in blood plasma after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are unknown. Proteomic deconvolution of the IgG repertoire to the spike glycoprotein in convalescent subjects revealed that the response is directed predominantly (>80%) against epitopes residing outside the receptor binding domain (RBD). In one subject, just four IgG lineages accounted for 93.5% of the response, including an amino (N)-terminal domain (NTD)-directed antibody that was protective against lethal viral challenge. Genetic, structural, and functional characterization of a multidonor class of "public" antibodies revealed an NTD epitope that is recurrently mutated among emerging SARS-CoV-2 variants of concern. These data show that "public" NTD-directed and other non-RBD plasma antibodies are prevalent and have implications for SARS-CoV-2 protection and antibody escape.

Trace Reconstruction Problems in Computational Biology

The problem of reconstructing a string from its error-prone copies, the trace reconstruction problem, was introduced by Vladimir Levenshtein two decades ago. While there has been considerable theoretical work on trace reconstruction, practical solutions have only recently started to emerge in the context of two rapidly developing research areas: immunogenomics and DNA data storage. In immunogenomics, traces correspond to mutated copies of genes, with mutations generated naturally by the adaptive immune system. In DNA data storage, traces correspond to noisy copies of DNA molecules that encode digital data, with errors being artifacts of the data retrieval process. In this paper, we introduce several new trace generation models and open questions relevant to trace reconstruction for immunogenomics and DNA data storage, survey theoretical results on trace reconstruction, and highlight their connections to computational biology. Throughout, we discuss the applicability and shortcomings of known solutions and suggest future research directions.

Diversity in immunogenomics: the value and the challenge

Immunogenomics studies have been largely limited to individuals of European ancestry, restricting the ability to identify variation in human adaptive immune responses across populations. Inclusion of a greater diversity of individuals in immunogenomics studies will substantially enhance our understanding of human immunology.

Large-scale analysis of 2,152 Ig-seq datasets reveals key features of B cell biology and the antibody repertoire

Antibody repertoire sequencing enables researchers to acquire millions of B cell receptors and investigate these molecules at the single-nucleotide level. This power and resolution in studying humoral responses have led to its wide applications. However, most of these studies were conducted with a limited number of samples. Given the extraordinary diversity, assessment of these key features with a large sample set is demanded. Thus, we collect and systematically analyze 2,152 high-quality heavy-chain antibody repertoires. Our study reveals that 52 core variable genes universally contribute to more than 99% of each individual's repertoire; a distal interspersed preferences characterize V gene recombination; the number of public clones between two repertoires follows a linear model, and the positive selection dominates at RGYW motif in somatic hypermutations. Thus, this population-level analysis resolves some critical features of the antibody repertoire and may have significant value to the large cadre of scientists.

A single donor is sufficient to produce a highly functional in vitro antibody library

Antibody complementarity determining region diversity has been considered to be the most important metric for the production of a functional antibody library. Generally, the greater the antibody library diversity, the greater the probability of selecting a diverse array of high affinity leads. According to this paradigm, the primary means of elevating library diversity has been by increasing the number of donors. In the present study we explored the possibility of creating an in vitro antibody library from a single healthy individual, showing that the number of lymphocytes, rather than the number of donors, is the key criterion in the production of a diverse and functional antibody library. We describe the construction of a high-quality phage display library comprising 5 × 10⁹ human antibodies by applying an efficient B cell extraction protocol from a single donor and a targeted V-gene amplification strategy favoring specific antibody families for their improved developability profiles. Each step of the library generation process was followed and validated by next generation sequencing to monitor the library quality and diversity. The functionality of the library was tested using several therapeutically relevant targets for which a vast number of different antibodies with desired biophysical properties were obtained.

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.

Structural and functional ramifications of antigenic drift in recent SARS-CoV-2 variants

The protective efficacy of neutralizing antibodies (nAbs) elicited during natural infection with SARS-CoV-2 and by vaccination based on its spike protein has been compromised with emergence of the recent SARS-CoV-2 variants. Residues E484 and K417 in the receptor-binding site (RBS) are both mutated in lineages first described in South Africa (B.1.351) and Brazil (B.1.1.28.1). The nAbs isolated from SARS-CoV-2 patients are preferentially encoded by certain heavy-chain germline genes and the two most frequently elicited antibody families (IGHV3-53/3-66 and IGHV1-2) can each bind the RBS in two different binding modes. However, their binding and neutralization are abrogated by either the E484K or K417N mutation, whereas nAbs to the cross-reactive CR3022 and S309 sites are largely unaffected. This structural and functional analysis illustrates why mutations at E484 and K417 adversely affect major classes of nAbs to SARS-CoV-2 with consequences for next-generation COVID-19 vaccines.

Recognition of the SARS-CoV-2 receptor binding domain by neutralizing antibodies

Immediately from the outset of the COVID-19 pandemic, researchers from diverse biomedical and biological disciplines have united to study the novel pandemic virus, SARS-CoV-2. The antibody response to SARS-CoV-2 has been a major focus of COVID-19 research due to its clinical relevance and importance in vaccine and therapeutic development. Isolation and characterization of antibodies to SARS-CoV-2 have been accumulating at an unprecedented pace. Most of the SARS-CoV-2 neutralizing antibodies to date target the spike (S) protein receptor binding domain (RBD), which engages the host receptor ACE2 for viral entry. Here we review the binding sites and molecular features of monoclonal antibodies that target the SARS-CoV-2 RBD, including a few that also cross-neutralize SARS-CoV.

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.

Current strategies for detecting functional convergence across B-cell receptor repertoires

Convergence across B-cell receptor (BCR) and antibody repertoires has become instrumental in prioritizing candidates in recent rapid therapeutic antibody discovery campaigns. It has also increased our understanding of the immune system, providing evidence for the preferential selection of BCRs to particular (immunodominant) epitopes post vaccination/infection. These important implications for both drug discovery and immunology mean that it is essential to consider the optimal way to combine experimental and computational technology when probing BCR repertoires for convergence signatures. Here, we discuss the theoretical basis for observing BCR repertoire functional convergence and explore factors of study design that can impact functional signal. We also review the computational arsenal available to detect antibodies with similar functional properties, highlighting opportunities enabled by recent clustering algorithms that exploit structural similarities between BCRs. Finally, we suggest future areas of development that should increase the power of BCR repertoire functional clustering.

Prevalent, protective, and convergent IgG recognition of SARS-CoV-2 non-RBD spike epitopes in COVID-19 convalescent plasma

Although humoral immunity is essential for control of SARS-CoV-2, the molecular composition, binding epitopes and effector functions of the immunoglobulin G (IgG) antibodies that circulate in blood plasma following infection are unknown. Proteomic deconvolution of the circulating IgG repertoire (Ig-Seq ¹ ) to the spike ectodomain (S-ECD ² ) in four convalescent study subjects revealed that the plasma response is oligoclonal and directed predominantly (>80%) to S-ECD epitopes that lie outside the receptor binding domain (RBD). When comparing antibodies directed to either the RBD, the N-terminal domain (NTD) or the S2 subunit (S2) in one subject, just four IgG lineages (1 anti-S2, 2 anti-NTD and 1 anti-RBD) accounted for 93.5% of the repertoire. Although the anti-RBD and one of the anti-NTD antibodies were equally potently neutralizing in vitro , we nonetheless found that the anti-NTD antibody was sufficient for protection to lethal viral challenge, either alone or in combination as a cocktail where it dominated the effect of the other plasma antibodies. We identified in vivo protective plasma anti-NTD antibodies in 3/4 subjects analyzed and discovered a shared class of antibodies targeting the NTD that utilize unmutated or near-germline IGHV1-24, the most electronegative IGHV gene in the human genome. Structural analysis revealed that binding to NTD is dominated by interactions with the heavy chain, accounting for 89% of the entire interfacial area, with germline residues uniquely encoded by IGHV1-24 contributing 20% (149 Å ² ). Together with recent reports of germline IGHV1-24 antibodies isolated by B-cell cloning ^3,4 our data reveal a class of shared IgG antibodies that are readily observed in convalescent plasma and underscore the role of NTD-directed antibodies in protection against SARS-CoV-2 infection.

A Novel Framework for Characterizing Genomic Haplotype Diversity in the Human Immunoglobulin Heavy Chain Locus

An incomplete ascertainment of genetic variation within the highly polymorphic immunoglobulin heavy chain locus (IGH) has hindered our ability to define genetic factors that influence antibody-mediated processes. Due to locus complexity, standard high-throughput approaches have failed to accurately and comprehensively capture IGH polymorphism. As a result, the locus has only been fully characterized two times, severely limiting our knowledge of human IGH diversity. Here, we combine targeted long-read sequencing with a novel bioinformatics tool, IGenotyper, to fully characterize IGH variation in a haplotype-specific manner. We apply this approach to eight human samples, including a haploid cell line and two mother-father-child trios, and demonstrate the ability to generate high-quality assemblies (>98% complete and >99% accurate), genotypes, and gene annotations, identifying 2 novel structural variants and 15 novel IGH alleles. We show multiplexing allows for scaling of the approach without impacting data quality, and that our genotype call sets are more accurate than short-read (>35% increase in true positives and >97% decrease in false-positives) and array/imputation-based datasets. This framework establishes a desperately needed foundation for leveraging IG genomic data to study population-level variation in antibody-mediated immunity, critical for bettering our understanding of disease risk, and responses to vaccines and therapeutics.

Structural basis of a shared antibody response to SARS-CoV-2

Molecular understanding of neutralizing antibody responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) could accelerate vaccine design and drug discovery. We analyzed 294 anti-SARS-CoV-2 antibodies and found that immunoglobulin G heavy-chain variable region 3-53 (IGHV3-53) is the most frequently used IGHV gene for targeting the receptor-binding domain (RBD) of the spike protein. Co-crystal structures of two IGHV3-53-neutralizing antibodies with RBD, with or without Fab CR3022, at 2.33- to 3.20-angstrom resolution revealed that the germline-encoded residues dominate recognition of the angiotensin I converting enzyme 2 (ACE2)-binding site. This binding mode limits the IGHV3-53 antibodies to short complementarity-determining region H3 loops but accommodates light-chain diversity. These IGHV3-53 antibodies show minimal affinity maturation and high potency, which is promising for vaccine design. Knowledge of these structural motifs and binding mode should facilitate the design of antigens that elicit this type of neutralizing response.

Structural basis of a public antibody response to SARS-CoV-2

Molecular-level understanding of human neutralizing antibody responses to SARS-CoV-2 could accelerate vaccine design and facilitate drug discovery. We analyzed 294 SARS-CoV-2 antibodies and found that IGHV3-53 is the most frequently used IGHV gene for targeting the receptor binding domain (RBD) of the spike (S) protein. We determined crystal structures of two IGHV3-53 neutralizing antibodies +/- Fab CR3022 ranging from 2.33 to 3.11 Å resolution. The germline-encoded residues of IGHV3-53 dominate binding to the ACE2 binding site epitope with no overlap with the CR3022 epitope. Moreover, IGHV3-53 is used in combination with a very short CDR H3 and different light chains. Overall, IGHV3-53 represents a versatile public VH in neutralizing SARS-CoV-2 antibodies, where their specific germline features and minimal affinity maturation provide important insights for vaccine design and assessing outcomes.

Inter- and intraspecies comparison of phylogenetic fingerprints and sequence diversity of immunoglobulin variable genes

Protection and neutralization of a vast array of pathogens is accomplished by the tremendous diversity of the B cell receptor (BCR) repertoire. For jawed vertebrates, this diversity is initiated via the somatic recombination of immunoglobulin (Ig) germline elements. While it is clear that the number of these germline segments differs from species to species, the extent of cross-species sequence diversity remains largely uncharacterized. Here we use extensive computational and statistical methods to investigate the sequence diversity and evolutionary relationship between Ig variable (V), diversity (D), and joining (J) germline segments across nine commonly studied species ranging from zebrafish to human. Metrics such as guanine-cytosine (GC) content showed low redundancy across Ig germline genes within a given species. Other comparisons, including amino acid motifs, evolutionary selection, and sequence diversity, revealed species-specific properties. Additionally, we showed that the germline-encoded diversity differs across antibody (recombined V-D-J) repertoires of various B cell subsets. To facilitate future comparative immunogenomics analysis, we created VDJgermlines, an R package that contains the germline sequences from multiple species. Our study informs strategies for the humanization and engineering of therapeutic antibodies.

Dynamics of heavy chain junctional length biases in antibody repertoires

Antibody variable domain sequence diversity is generated by recombination of germline segments. The third complementarity-determining region of the heavy chain (CDR H3) is the region of highest sequence diversity and is formed by the joining of heavy chain V_H, D_H and J_H germline segments combined with random nucleotide trimming and additions between these segments. We show that CDR H3 and junctional segment length distributions are biased in human antibody repertoires as a function of V_H, V_L and J_H germline segment utilization. Most length biases are apparent in the naive and antigen experienced B cell compartments but not in nonproductive recombination products, indicating B cell selection as a major driver of these biases. Our findings reveal biases in the antibody CDR H3 diversity landscape shaped by V_H, V_L, and J_H germline segment use during naive and antigen-experienced repertoire selection.

Sankar et al. investigate the junctional length biases (determining antibody binding potential) as a function of germline gene usage in antibody repertoires. They show that CDR H3 and junction length are biased by V_H, V_L, and J_H germline segment usage and these biases are apparent in both naive and antigen-experienced repertoires but not in non-productive repertoires.

Automated analysis of immunosequencing datasets reveals novel immunoglobulin D genes across diverse species

Immunoglobulin genes are formed through V(D)J recombination, which joins the variable (V), diversity (D), and joining (J) germline genes. Since variations in germline genes have been linked to various diseases, personalized immunogenomics focuses on finding alleles of germline genes across various patients. Although reconstruction of V and J genes is a well-studied problem, the more challenging task of reconstructing D genes remained open until the IgScout algorithm was developed in 2019. In this work, we address limitations of IgScout by developing a probabilistic MINING-D algorithm for D gene reconstruction, apply it to hundreds of immunosequencing datasets from multiple species, and validate the newly inferred D genes by analyzing diverse whole genome sequencing datasets and haplotyping heterozygous V genes.

Antibodies provide specific binding to an enormous range of antigens and represent a key component of the adaptive immune system. Immunosequencing has emerged as a method of choice for generating millions of reads that sample antibody repertoires and provides insights into monitoring immune response to disease and vaccination. Most of the previous immunogenomics studies rely on the reference germline genes in the immunoglobulin locus rather than the germline genes in a specific patient. This approach is deficient since the set of known germline genes is incomplete (particularly for non-European humans and non-human species) and contains alleles that resulted from sequencing and annotation errors. The problem of de novo inference of diversity (D) genes from immunosequencing data remained open until the IgScout algorithm was developed in 2019. We address limitations of IgScout by developing a probabilistic MINING-D algorithm for D gene reconstruction and infer multiple D genes across multiple species that are not present in standard databases.

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

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