Codon insertion and deletion functions as a somatic diversification mechanism in human antibody repertoires

Fast clonal family inference from large-scale B cell repertoire sequencing data

Advances in high-throughput sequencing technologies have facilitated the large-scale characterization of B cell receptor (BCR) repertoires. However, the vast amount and high diversity of the BCR sequences pose challenges for efficient and biologically meaningful analysis. Here, we introduce fastBCR, an efficient computational approach for inferring B cell clonal families from massive BCR heavy chain sequences. We demonstrate that fastBCR substantially reduces the running time while ensuring high accuracy on simulated datasets with diverse numbers of B cell lineages and varying mutation rates. We apply fastBCR to real BCR sequencing data from peripheral blood samples of COVID-19 patients, showing that the inferred clonal families display disease-associated features, as well as corresponding antigen-binding specificity and affinity. Overall, our results demonstrate the advantages of fastBCR for analyzing BCR repertoire data, which will facilitate the identification of disease-associated antibodies and improve our understanding of the B cell immune response.

DNA repair mechanisms that promote insertion-deletion events during immunoglobulin gene diversification

Insertions and deletions (indels) are low-frequency deleterious genomic DNA alterations. Despite their rarity, indels are common, and insertions leading to long complementarity-determining region 3 (CDR3) are vital for antigen-binding functions in broadly neutralizing and polyreactive antibodies targeting viruses. Because of challenges in detecting indels, the mechanism that generates indels during immunoglobulin diversification processes remains poorly understood. We carried out ultra-deep profiling of indels and systematically dissected the underlying mechanisms using passenger-immunoglobulin mouse models. We found that activation-induced cytidine deaminase-dependent ±1-base pair (bp) indels are the most prevalent indel events, biasing deleterious outcomes, whereas longer in-frame indels, especially insertions that can extend the CDR3 length, are rare outcomes. The ±1-bp indels are channeled by base excision repair, but longer indels require additional DNA-processing factors. Ectopic expression of a DNA exonuclease or perturbation of the balance of DNA polymerases can increase the frequency of longer indels, thus paving the way for models that can generate antibodies with long CDR3. Our study reveals the mechanisms that generate beneficial and deleterious indels during the process of antibody somatic hypermutation and has implications in understanding the detrimental genomic alterations in various conditions, including tumorigenesis.

Learning the statistics and landscape of somatic mutation-induced insertions and deletions in antibodies

Affinity maturation is crucial for improving the binding affinity of antibodies to antigens. This process is mainly driven by point substitutions caused by somatic hypermutations of the immunoglobulin gene. It also includes deletions and insertions of genomic material known as indels. While the landscape of point substitutions has been extensively studied, a detailed statistical description of indels is still lacking. Here we present a probabilistic inference tool to learn the statistics of indels from repertoire sequencing data, which overcomes the pitfalls and biases of standard annotation methods. The model includes antibody-specific maturation ages to account for variable mutational loads in the repertoire. After validation on synthetic data, we applied our tool to a large dataset of human immunoglobulin heavy chains. The inferred model allows us to identify universal statistical features of indels in heavy chains. We report distinct insertion and deletion hotspots, and show that the distribution of lengths of indels follows a geometric distribution, which puts constraints on future mechanistic models of the hypermutation process.

Affinity maturation of B cell receptors is an important mechanism by which our body designs neutralizing antibodies to defend us against pathogens, including broadly neutralizing antibodies, which target a wide range of variants of the same pathogen. Such antibodies often contain key insertions and deletions to the germline gene, or “indels”, which are caused by somatic hypermutations. However, the mechanism, frequency and role of these indels are still elusive. We designed a computational method based on a probabilistic framework to infer the characteristics of this mutational process from high-throughput antibody sequencing experiments. Applied to human data, our approach provides a comprehensive quantitative description of insertions and deletions, opening avenues for better understanding the process of affinity maturation and the design of vaccines for eliciting a broad antibody response.

In vitro evolution of antibody affinity via insertional scanning mutagenesis of an entire antibody variable region

We report a systematic combinatorial exploration of affinity enhancement of antibodies by insertions and deletions (InDels). Transposon-based introduction of InDels via the method TRIAD (transposition-based random insertion and deletion mutagenesis) was used to generate large libraries with random in-frame InDels across the entire single-chain variable fragment gene that were further recombined and screened by ribosome display. Knowledge of potential insertion points from TRIAD libraries formed the basis of exploration of length and sequence diversity of novel insertions by insertional-scanning mutagenesis (InScaM). An overall 256-fold affinity improvement of an anti-IL-13 antibody BAK1 as a result of InDel mutagenesis and combination with known point mutations validates this approach, and suggests that the results of this InDel mutagenesis and conventional exploration of point mutations can synergize to generate antibodies with higher affinity.

An ultralong CDRH2 in HCV neutralizing antibody demonstrates structural plasticity of antibodies against E2 glycoprotein

A vaccine protective against diverse HCV variants is needed to control the HCV epidemic. Structures of E2 complexes with front layer-specific broadly neutralizing antibodies (bNAbs) isolated from HCV-infected individuals, revealed a disulfide bond-containing CDRH3 that adopts straight (individuals who clear infection) or bent (individuals with chronic infection) conformation. To investigate whether a straight versus bent disulfide bond-containing CDRH3 is specific to particular HCV-infected individuals, we solved a crystal structure of the HCV E2 ectodomain in complex with AR3X, a bNAb with an unusually long CDRH2 that was isolated from the chronically-infected individual from whom the bent CDRH3 bNAbs were derived. The structure revealed that AR3X utilizes both its ultralong CDRH2 and a disulfide motif-containing straight CDRH3 to recognize the E2 front layer. These results demonstrate that both the straight and bent CDRH3 classes of HCV bNAb can be elicited in a single individual, revealing a structural plasticity of VH1-69-derived bNAbs.

A proposed reverse transcription mechanism for (CAG)n and similar expandable repeats that cause neurological and other diseases

The mechanism of (CAG)n repeat generation, and related expandable repeat diseases in non-dividing cells, is currently understood in terms of a DNA template-based DNA repair synthesis process involving hairpin stabilized slippage, local error-prone repair via MutSβ (MSH2-MSH3) hairpin protective stabilization, then nascent strand extension by DNA polymerases-β and -δ. We advance a very similar slipped hairpin-stabilized model involving MSH2-MSH3 with two key differences: the copying template may also be the nascent pre-mRNA with the repair pathway being mediated by the Y-family error-prone enzymes DNA polymerase-η and DNA polymerase-κ acting as reverse transcriptases. We argue that both DNA-based and RNA-based mechanisms could well be activated in affected non-dividing brain cells in vivo. Here, we compare the advantages of the RNA/RT-based model proposed by us as an adjunct to previously proposed models. In brief, our model depends upon dysregulated innate and adaptive immunity cascades involving AID/APOBEC and ADAR deaminases that are known to be involved in normal locus-specific immunoglobulin somatic hypermutation, cancer progression and somatic mutations at many off-target non-immunoglobulin sites across the genome: we explain how these processes could also play an active role in repeat expansion diseases at RNA polymerase II-transcribed genes.

Cis- and trans-factors affecting AID targeting and mutagenic outcomes in antibody diversification

Antigen receptor diversification is a hallmark of adaptive immunity which allows specificity of the receptor to particular antigen. B cell receptor (BCR) or its secreted form, antibody, is diversified through antigen-independent and antigen-dependent mechanisms. During B cell development in bone marrow, BCR is diversified via V(D)J recombination mediated by RAG endonuclease. Upon stimulation by antigen, B cell undergo somatic hypermutation (SHM) to allow affinity maturation and class switch recombination (CSR) to change the effector function of the antibody. Both SHM and CSR are initiated by activation-induced cytidine deaminase (AID). Repair of AID-initiated lesions through different DNA repair pathways results in diverse mutagenic outcomes. Here, we focus on discussing cis- and trans-factors that target AID to its substrates and factors that affect different outcomes of AID-initiated lesions. The knowledge of mechanisms that govern AID targeting and outcomes could be harnessed to elicit rare functional antibodies and develop ex vivo antibody diversification approaches with diversifying base editors.

Deciphering evolution of immune recognition in antibodies

Background

Antibody, the primary effector molecule of the immune system, evolves after initial encounter with the antigen from a precursor form to a mature one to effectively deal with the antigen. Antibodies of a lineage diverge through antigen-directed isolated pathways of maturation to exhibit distinct recognition potential. In the context of evolution in immune recognition, diversity of antigen cannot be ignored. While there are reports on antibody lineage, structural perspective with respect to diverse recognition potential in a lineage has never been studied. Hence, it is crucial to evaluate how maturation leads to topological tailoring within a lineage enabling them to interact with significantly distinct antigens.

Results

A data-driven approach was undertaken for the study. Global experimental mouse and human antibody-antigen complex structures from PDB were compiled into a coherent database of germline-linked antibodies bound with distinct antigens. Structural analysis of all lineages showed variations in CDRs of both H and L chains. Observations of conformational adaptation made from analysis of static structures were further evaluated by characterizing dynamics of interaction in two lineages, mouse V_H1–84 and human V_H5–51. Sequence and structure analysis of the lineages explained that somatic mutations altered the geometries of individual antibodies with common structural constraints in some CDRs. Additionally, conformational landscape obtained from molecular dynamics simulations revealed that incoming pathogen led to further conformational divergence in the paratope (as observed across datasets) even while maintaining similar overall backbone topology. MM-GB/SA analysis showed binding energies to be in physiological range. Results of the study are coherent with experimental observations.

Conclusions

The findings of this study highlight basic structural principles shaping the molecular evolution of a lineage for significantly diverse antigens. Antibodies of a lineage follow different developmental pathways while preserving the imprint of the germline. From the study, it can be generalized that structural diversification of the paratope is an outcome of natural selection of a conformation from an available ensemble, which is further optimized for antigen interaction. The study establishes that starting from a common lineage, antibodies can mature to recognize a wide range of antigens. This hypothesis can be further tested and validated experimentally.

Electronic supplementary material

The online version of this article (10.1186/s12900-018-0096-1) contains supplementary material, which is available to authorized users.

Diversity in the Cow Ultralong CDR H3 Antibody Repertoire

Typical antibodies found in humans and mice usually have short CDR H3s and generally flat binding surfaces. However, cows possess a subset of antibodies with ultralong CDR H3s that can range up to 70 amino acids and form a unique “stalk and knob” structure, with the knob protruding far out of the antibody surface, where it has the potential to bind antigens with concave epitopes. Activation-induced cytidine deaminase (AID) has a proven role in diversifying antibody repertoires in humoral immunity, and it has been found to induce somatic hypermutation in bovine immunoglobulin genes both before and after contact with antigen. Due to limited use of variable and diversity genes in the V(D)J recombination events that produce ultralong CDR H3 antibodies in cows, the diversity in the bovine ultralong antibody repertoire has been proposed to rely on AID-induced mutations targeted to the IGHD8-2 gene that encodes the entire knob region. In this review, we discuss the genetics, structures, and diversity of bovine ultralong antibodies, as well as the role of AID in creating a diverse antibody repertoire.

Immunogenetic factors driving formation of ultralong VH CDR3 in Bos taurus antibodies

The antibody repertoire of Bos taurus is characterized by a subset of variable heavy (VH) chain regions with ultralong third complementarity determining regions (CDR3) which, compared to other species, can provide a potent response to challenging antigens like HIV env. These unusual CDR3 can range to over seventy highly diverse amino acids in length and form unique β-ribbon 'stalk' and disulfide bonded 'knob' structures, far from the typical antigen binding site. The genetic components and processes for forming these unusual cattle antibody VH CDR3 are not well understood. Here we analyze sequences of Bos taurus antibody VH domains and find that the subset with ultralong CDR3 exclusively uses a single variable gene, IGHV1-7 (VHBUL) rearranged to the longest diversity gene, IGHD8-2. An eight nucleotide duplication at the 3' end of IGHV1-7 encodes a longer V-region producing an extended F β-strand that contributes to the stalk in a rearranged CDR3. A low amino acid variability was observed in CDR1 and CDR2, suggesting that antigen binding for this subset most likely only depends on the CDR3. Importantly a novel, potentially AID mediated, deletional diversification mechanism of the B. taurus VH ultralong CDR3 knob was discovered, in which interior codons of the IGHD8-2 region are removed while maintaining integral structural components of the knob and descending strand of the stalk in place. These deletions serve to further diversify cysteine positions, and thus disulfide bonded loops. Hence, both germline and somatic genetic factors and processes appear to be involved in diversification of this structurally unusual cattle VH ultralong CDR3 repertoire.

Nucleotide insertions and deletions complement point mutations to massively expand the diversity created by somatic hypermutation of antibodies

During somatic hypermutation (SHM), deamination of cytidine by activation-induced cytidine deaminase and subsequent DNA repair generates mutations within immunoglobulin V-regions. Nucleotide insertions and deletions (indels) have recently been shown to be critical for the evolution of antibody binding. Affinity maturation of 53 antibodies using in vitro SHM in a non-B cell context was compared with mutation patterns observed for SHM in vivo. The origin and frequency of indels seen during in vitro maturation were similar to that in vivo. Indels are localized to CDRs, and secondary mutations within insertions further optimize antigen binding. Structural determination of an antibody matured in vitro and comparison with human-derived antibodies containing insertions reveal conserved patterns of antibody maturation. These findings indicate that activation-induced cytidine deaminase acting on V-region sequences is sufficient to initiate authentic formation of indels in vitro and in vivo and that point mutations, indel formation, and clonal selection form a robust tripartite system for antibody evolution.

Immunologic Basis for Long HCDR3s in Broadly Neutralizing Antibodies Against HIV-1

A large number of potent broadly neutralizing antibodies (bnAbs) against HIV-1 have been reported in recent years, raising hope for the possibility of an effective vaccine based on epitopes recognized by these protective antibodies. However, many of these bnAbs contain the long heavy chain complementarity-determining region 3 (HCDR3), which is viewed as an obstacle to the development of an HIV-1 vaccine targeting the bnAb responses. This mini-review summarizes the current literature and discusses the different potential immunologic mechanisms for generating long HCDR3, including D–D fusion, VH replacement, long N region addition, and skewed D–J gene usage, among which potential VH replacement products appear to be significant contributors. VH replacement occurs through recombinase activated gene-mediated secondary recombination and contributes to the diversified naïve B cell repertoire. During VH replacement, a short stretch of nucleotides from previously rearranged VH genes remains within the newly formed HCDR3, thus elongating its length. Accumulating evidence suggests that long HCDR3s are present in significant numbers in the human mature naïve B cell repertoire and are primarily generated by recombination during B cell development. These new observations indicate that long HCDR3s, though low in frequency, are a normal feature of the human antibody naïve repertoire and they appear to be selected to target conserved epitopes located in deep, partially obscured regions of the HIV-1 envelope trimer. Therefore, the presence of long HCDR3 sequences should not necessarily be viewed as an obstacle to the development of an HIV-1 vaccine based upon bnAb responses.

JEC. Secondary mechanisms of diversification in the human antibody repertoire

V(D)J recombination and somatic hypermutation (SHM) are the primary mechanisms for diversification of the human antibody repertoire. These mechanisms allow for rapid humoral immune responses to a wide range of pathogenic challenges. V(D)J recombination efficiently generate a virtually limitless diversity through random recombination of variable (V), diversity (D), and joining (J) genes with diverse non-templated junctions between the selected gene segments. Following antigen stimulation, affinity maturation by SHM produces antibodies with refined specificity mediated by mutations typically focused in complementarity determining regions (CDRs), which form the bulk of the antigen recognition site. While V(D)J recombination and SHM are responsible for much of the diversity of the antibody repertoire, there are several secondary mechanisms that, while less frequent, make substantial contributions to antibody diversity including V(DD)J recombination (or D-D fusion), SHM-associated insertions and deletions, and affinity maturation and antigen contact by non-CDR regions of the antibody. In addition to enhanced diversity, these mechanisms allow the production of antibodies that are critical to response to a variety of viral and bacterial pathogens but that would be difficult to generate using only the primary mechanisms of diversification.

Insertion/Deletion polymorphisms do play any role in G6PD deficiency individuals in the Kingdom of the Saudi Arabia

Glucose-6-phosphate dehydrogenase (G6PD) is an enzyme in the pentose phosphate pathway (PPP) that plays an important role in protecting cells from oxidative damage by producing NADPH and reduced glutathione. G6PD deficiency is considered one of the most common genetic disorders present in the X chromosome and is the most common of enzymopathic red blood cell disorder. Angiotensin converting enzyme (ACE) plays an essential role in two physiological systems, one leading to the production of angiotensin II and the other to the degradation of bradykinin. Most studies focused on an insertion/deletion (I/D) polymorphism in intron 16 of the ACE gene as a marker for a functional polymorphism. The α(2B)-adrenergic receptor gene (α(2B)AR) is a three-amino acid deletion (12Glu9) polymorphism is located on chromosome 2. (Glu(9)/Glu(9)) of this polymorphism has been first time studies in G6PD individuals. We have selected 39 G6PD deficiency male individuals and PCR was carried out with the I/D polymorphisms. ACE I/D polymorphism study was carried out in G6PD individuals and showed strong association with DD genotypes and D alleles OR=39.38, p<0.0001 (95% CI=8.80-176.1) and OR=38.58, p<0.0001 (95% CI=13.21-112.6). Another gene of α(2B)AR I/D polymorphism study cannot show any association in DD genotype OR-0.6882,p=0.9388 (95% CI=0.2035-2.327) and with D allele OR-0.9614,p=0.9388 (95% CI=0.3482-2.653). Our study shows that G6PD deficiency is showing strong association in DD genotype and D allele of ACE gene and α(2B)AR gene have not shown any important role and one of the reason could be the low sample size.

Location and length distribution of somatic hypermutation-associated DNA insertions and deletions reveals regions of antibody structural plasticity

Following the initial diversity generated by V(D)J recombination, somatic hypermutation is the principal mechanism for producing further antibody repertoire diversity in antigen-experienced B cells. While somatic hypermutation typically results in single nucleotide substitutions, the infrequent incorporation of genetic insertions and deletions has also been associated with the somatic hypermutation process. We used high throughput antibody sequencing to determine the sequence of thousands of antibody genes containing somatic hypermutation-associated insertions and deletions (SHA indels), which revealed significant differences between the location of SHA indels and somatic mutations. Further, we identified a cluster of insertions and deletions in the antibody framework 3 region which corresponds to the hypervariable region 4 (HV4) in T cell receptors. We propose that this HV4-like region, identified by SHA indel analysis, represents a region of under-appreciated affinity maturation potential. Finally, through analysis of both location and length distribution of SHA indels, we have determined regions of structural plasticity within the antibody protein.

Human peripheral blood antibodies with long HCDR3s are established primarily at original recombination using a limited subset of germline genes

A number of antibodies that efficiently neutralize microbial targets contain long heavy chain complementarity determining region 3 (HCDR3) loops. For HIV, several of the most broad and potently neutralizing antibodies have exceptionally long HCDR3s. Two broad potently neutralizing HIV-specific antibodies, PG9 and PG16, exhibit secondary structure. Two other long HCDR3 antibodies, 2F5 and 4E10, protect against mucosal challenge with SHIV. Induction of such long HCDR3 antibodies may be critical to the design of an effective vaccine strategy for HIV and other pathogens, however it is unclear at present how to induce such antibodies. Here, we present genetic evidence that human peripheral blood antibodies containing long HCDR3s are not primarily generated by insertions introduced during the somatic hypermutation process. Instead, they are typically formed by processes occurring as part of the original recombination event. Thus, the response of B cells encoding antibodies with long HCDR3s results from selection of unusual clones from the naïve repertoire rather than through accumulation of insertions. These antibodies typically use a small subset of D and J gene segments that are particularly suited to encoding long HCDR3s, resulting in the incorporation of highly conserved genetic elements in the majority of antibody sequences encoding long HCDR3s.

Coupling mammalian cell surface display with somatic hypermutation for the discovery and maturation of human antibodies

A novel approach has been developed for the isolation and maturation of human antibodies that replicates key features of the adaptive immune system by coupling in vitro somatic hypermutation (SHM) with mammalian cell display. SHM is dependent on the action of the B cell specific enzyme, activation-induced cytidine deaminase (AID), and can be replicated in non-B cells through expression of recombinant AID. A library of human antibodies, based on germline V-gene segments with recombined human regions was used to isolate low-affinity antibodies to human β nerve growth factor (hβNGF). These antibodies, initially naïve to SHM, were subjected to AID-directed SHM in vitro and selected using the same mammalian cell display system, as illustrated by the maturation of one of the antibodies to low pM K(D). This approach overcomes many of the previous limitations of mammalian cell display, enabling direct selection and maturation of antibodies as full-length, glycosylated IgGs.

Mechanism of lethal toxin neutralization by a human monoclonal antibody specific for the PA(20) region of Bacillus anthracis protective antigen

The primary immunogenic component of the currently approved anthrax vaccine is the protective antigen (PA) unit of the binary toxin system. PA-specific antibodies neutralize anthrax toxins and protect against infection. Recent research has determined that in humans, only antibodies specific for particular determinants are capable of effecting toxin neutralization, and that the neutralizing epitopes recognized by these antibodies are distributed throughout the PA monomer. The mechanisms by which the majority of these epitopes effect neutralization remain unknown. In this report we investigate the process by which a human monoclonal antibody specific for the amino-terminal domain of PA neutralizes lethal toxin in an in vitro assay of cytotoxicity, and find that it neutralizes LT by blocking the requisite cleavage of the amino-terminal 20 kD portion of the molecule (PA₂₀) from the remainder of the PA monomer. We also demonstrate that the epitope recognized by this human monoclonal does not encompass the ₁₆₆RKKR₁₆₉ furin recognition sequence in domain 1 of PA.

Domain specificity of the human antibody response to Bacillus anthracis protective antigen

Protective antigen (PA) is the cell surface recognition moiety of the Bacillus anthracis A-B toxin system, and the active immunogenic component in the currently licensed human anthrax vaccine (BioThrax, or AVA). The serum antibody response to the PA protein is polyclonal and complex both in terms of the antibody combining sites utilized to bind PA and the PA-associated epitopes recognized. We have cloned, sequenced, and expressed a large panel of PA-specific human monoclonal antibodies from seven AVA-immunized donors. Dot blots, Western blots, and radiolabeled antigen capture assays employing both proteolytic fragments of PA and engineered PA sub-domain fusion proteins were used to determine the region (domain) of the PA monomer to which each of the cloned human antibodies bound. The domain specificity of the isolated monoclonals was highly biased towards the amino-terminal 20kDa fragment of PA (PA(20)), with the majority (62%) of independently arising antibody clones reacting with determinants located on this PA fragment. A similar bias in domain specificity was also demonstrated in the serum response of AVA-vaccinated donors. Since PA(20) is cleaved from the remainder of the monomer rapidly following cell surface binding and has no known role in the intoxication process, the immunodominance of PA(20)-associated epitopes may directly affect the efficacy of PA-based anthrax vaccines.

Paratope diversity in the human antibody response to Bacillus anthracis protective antigen

The active component of the licensed human anthrax vaccine (BioThrax, or AVA) is a Bacillus anthracis toxin known as protective antigen (PA). Second generation anthrax vaccines currently under development are also based on a recombinant form of PA. Since the current and future anthrax vaccines are based on this toxin, it is important that the immunobiology of this protein in vaccinated humans be understood in detail. We have isolated and analyzed the PA-specific antibody repertoire from an AVA-vaccinated individual. When examined at the clonal level, we find an antibody response that is complex in terms of the combinatorial elements and immunoglobulin variable genes employed. All PA-specific antibodies had undergone somatic hypermutation and class switch recombination, both signs of affinity maturation. Although the antigenic epitopes recognized by the response were distributed throughout the PA monomer, the majority of antibodies arising in this individual following vaccination recognize determinants located on the amino-terminal (PA20) sub-domain of the molecule. This latter finding may have implications for the rational design of future PA-based anthrax vaccines.