Generation of heterogeneous rabbit anti-DNP antibodies by gene conversion and hypermutation of rearranged VL and VH genes during clonal expansion of B cells in splenic germinal centers

Immunoglobulin heavy chain variable region analysis in dairy goats

Based on the goat genome database, we have annotated the genomic organization of the goat immunoglobulin heavy chain variable region. The goat IgH locus is present on seven genome scaffolds, and contains ten V_H, three D_H and six J_H segments. After the exclusion of three shorter segments, the V_H genes were divided into two gene families based on sequence similarity. By analyzing the IgH cDNA sequences, we further identified that V_H2 (54.2%), D_H1 (61.7%) and J_H1 (60.5%) segments were most frequently utilized in the expression of the immunoglobulin variable region, and that point mutations introduced by somatic hypermutation were the major mutation present in these expressed variable region. Compared with human and horses, D_H-D_H fusion occurred at a higher frequency in goat V(D)J recombination. These results provided variable insights into goat immunoglobulin heavy chain variable region genome loci and repertoire diversity.

From rabbit antibody repertoires to rabbit monoclonal antibodies

In this review, we explain why and how rabbit monoclonal antibodies have become outstanding reagents for laboratory research and increasingly for diagnostic and therapeutic applications. Starting with the unique ontogeny of rabbit B cells that affords highly distinctive antibody repertoires rich in in vivo pruned binders of high diversity, affinity and specificity, we describe the generation of rabbit monoclonal antibodies by hybridoma technology, phage display and alternative methods, along with an account of successful humanization strategies.

Accuracy and coverage assessment of Oryctolagus cuniculus (rabbit) genes encoding immunoglobulins in the whole genome sequence assembly (OryCun2.0) and localization of the IGH locus to chromosome 20

We report on the analyses of genes encoding immunoglobulin heavy and light chains in the rabbit 6.51× whole genome assembly. This OryCun2.0 assembly confirms previous mapping of the duplicated IGK1 and IGK2 loci to chromosome 2 and the IGL lambda light chain locus to chromosome 21. The most frequently rearranged and expressed IGHV1 that is closest to IG DH and IGHJ genes encodes rabbit VHa allotypes. The partially inbred Thorbecke strain rabbit used for whole-genome sequencing was homozygous at the IGK but heterozygous with the IGHV1a1 allele in one of 79 IGHV-containing unplaced scaffolds and IGHV1a2, IGHM, IGHG, and IGHE sequences in another. Some IGKV, IGLV, and IGHA genes are also in other unplaced scaffolds. By fluorescence in situ hybridization, we assigned the previously unmapped IGH locus to the q-telomeric region of rabbit chromosome 20. An approximately 3-Mb segment of human chromosome 14 including IGH genes predicted to map to this telomeric region based on synteny analysis could not be located on assembled chromosome 20. Unplaced scaffold chrUn0053 contains some of the genes that comparative mapping predicts to be missing. We identified discrepancies between previous targeted studies and the OryCun2.0 assembly and some new BAC clones with IGH sequences that can guide other studies to further sequence and improve the OryCun2.0 assembly. Complete knowledge of gene sequences encoding variable regions of rabbit heavy, kappa, and lambda chains will lead to better understanding of how and why rabbits produce antibodies of high specificity and affinity through gene conversion and somatic hypermutation.

Gene Conversion-Like Events in the Diversification of Human Rearranged IGHV3-23*01 Gene Sequences

Gene conversion (GCV), a mechanism mediated by activation-induced cytidine deaminase (AID) is well established as a mechanism of immunoglobulin diversification in a few species. However, definitive evidence of GCV-like events in human immunoglobulin genes is scarce. The lack of evidence of GCV in human rearranged immunoglobulin gene sequences is puzzling given the presence of highly similar germline donors and the presence of all the enzymatic machinery required for GCV. In this study, we undertook a computational analysis of rearranged IGHV3-23(*)01 gene sequences from common variable immunodeficiency (CVID) patients, AID-deficient patients, and healthy individuals to survey "GCV-like" activities. We analyzed rearranged IGHV3-23(*)01 gene sequences obtained from total PBMC RNA and single-cell polymerase chain reaction of individual B cell lysates. Our search identified strong evidence of GCV-like activity. We observed that GCV-like tracts are flanked by AID hotspot motifs. Structural modeling of IGHV3-23(*)01 gene sequence revealed that hypermutable bases flanking GCV-like tracts are in the single stranded DNA (ssDNA) of stable stem-loop structures (SLSs). ssDNA is inherently fragile and also an optimal target for AID. We speculate that GCV could have been initiated by the targeting of hypermutable bases in ssDNA state in stable SLSs, plausibly by AID. We have observed that the frequency of GCV-like events is significantly higher in rearranged IGHV3-23-(*)01 sequences from healthy individuals compared to that of CVID patients. We did not observe GCV-like events in rearranged IGHV3-23-(*)01 sequences from AID-deficient patients. GCV, unlike somatic hypermutation (SHM), can result in multiple base substitutions that can alter many amino acids. The extensive changes in antibody affinity by GCV-like events would be instrumental in protecting humans against pathogens that diversify their genome by antigenic shift.

Molecular bases of genetic diversity and evolution of the immunoglobulin heavy chain variable region (IGHV) gene locus in leporids

The rabbit has long been a model for studies of the immune system. Work using rabbits contributed both to the battle against infectious diseases such as rabies and syphilis, and to our knowledge, of antibodies' structure, function, and regulated expression. With the description of rabbit Ig allotypes, the discovery of different gene segments encoding immunoglobulins became possible. This challenged the "one gene-one protein" dogma. The observation that rabbit allotypic specificities of the variable regions were present on IgM and IgG molecules also led to the hypothesis of Ig class switching. Rabbit allotypes contributed to the documentation of phenomena such as allelic exclusion and imbalance in production of allelic gene products. During the last 30 years, the rabbit Ig allotypes revealed a number of unique features, setting them apart from mice, humans, and other mammals. Here, we review the most relevant findings concerning the rabbit IGHV. Among these are the preferential usage of one VH gene in VDJ rearrangements, the existence of trans-species polymorphism in the IGHV locus revealed by serology and confirmed by sequencing IGHV genes in Lepus, the unusually large genetic distances between allelic lineages and the fact that the antibody repertoire is diversified in this species only after birth. The whole genome sequence of a rabbit, plus re-sequencing of additional strains and related genera, will allow further evolutionary investigations of antibody variation. Continued research will help define the roles that genetic, allelic, and population diversity at antibody loci may play in host-parasite interactions.

Antigen-driven selection in germinal centers as reflected by the shape characteristics of immunoglobulin gene lineage trees: a large-scale simulation study

During the immune response, the generation of memory B lymphocytes in germinal centers involves affinity maturation of the cells' antigen receptors, based on somatic hypermutation of receptor genes and antigen-driven selection of the resulting mutants. Affinity maturation is vital for immune protection, and is the basis of humoral immune learning and memory. Lineage trees of somatically hypermutated immunoglobulin genes often serve to qualitatively illustrate claims concerning the dynamics of affinity maturation in germinal centers. Here, we derive the quantitative relationships between parameters characterizing affinity maturation dynamics (proliferation, differentiation and mutation rates, initial affinity of the Ig to the antigen, and selection thresholds) and the mathematical properties of lineage trees, using a computer simulation which combines mathematical models for all mature B cell populations, stochastic models of hypermutation and selection, lineage tree generation and measurement of graphical tree characteristics. We identified seven key lineage tree properties, and found correlations of these with initial clone affinity and with the selection threshold. These two parameters were found to be the main factors affecting lineage tree shapes in both primary and secondary response trees. The results also confirm that recycling from centrocytes back to centroblasts is highly likely.

Models for antigen receptor gene rearrangement: CDR3 length

Despite the various processing steps involved in V(D)J recombination, which could potentially introduce many biases in the length distribution of complementarity determining region 3 (CDR3) segments, the observed CDR3 length distributions for complete repertoires are very close to a normal-like distribution. This raises the question of whether this distribution is simply a result of the random steps included in the process of gene rearrangement, or has been optimized during evolution. We have addressed this issue by constructing a simulation of gene rearrangement, which takes into account the DNA modification steps included in the process, namely hairpin opening, nucleotide additions, and nucleotide deletions. We found that the near-Gaussian- shape of CDR3 length distribution can only be obtained under a relatively narrow set of parameter values, and thus our model suggests that specific biases govern the rearrangement process. In both B-cell receptor (BCR) heavy chain and T-cell receptor beta chain, we obtained a Gaussian distribution using identical parameters, despite the difference in the number and the lengths of the D segments. Hence our results suggest that these parameters most likely reflect the optimal conditions under which the rearrangement process occurs. We have subsequently used the insights gained in this study to estimate the probability of occurrence of two exactly identical BCRs over the course of a human lifetime. Whereas identical rearrangements of the heavy chain are highly unlikely to occur within one human lifetime, for the light chain we found that this probability is not negligible, and hence the light chain CDR3 alone cannot serve as an indicator of B-cell clonality.

Germinal center-independent affinity maturation in tumor necrosis factor receptor 1-deficient mice

Germinal centers (GCs) have been identified as site at which the somatic mutation of immunoglobulins occurs.However, somatic mutations in immunoglobulins have also been observed in animals that normally do not harbor germinal centers. This clearly indicates that somatic mutations can occur in the absence of germinal centers. We therefore attempted to determine whether or not GCs exist in TNFR1-deficient mice, and are essential for the somatic mutation of immunoglobulins, using (4-hydroxy-3-nitropheny)acetyl-ovalbumin (NP-OVA). Both wild-type and TNFR1-deficient mice were immunized with NPOVA, and then examined with regard to the existence of GCs. No typical B-cell follicles were detected in the TNFR1-deficient mice. Cell proliferation was detected throughout all splenic tissue types, and no in vivo immunecomplex retention was observed in the TNFR1-deficient mice. All of these data strongly suggest that no GCs were formed in the TNFR1-deficient mice. Although TNFR1-deficient mice are unable to form GCs, serological analyses indicated that affinity maturation had been achieved in both the wild-type and TNFR1-deficient mice. We therefore isolated and sequenced several DNA clones from wild-type and the TNFR1-deficient mice. Eight out of 12 wild-type clones, and 11 out of 14 clones of the TNFR-1-deficient mice contained mutations at the CDR1 site. Thus, the wild-type and TNFR1-deficient mice were not extremely different with regard to types and rates of somatic mutation. Also, high-affinity antibodies were detected in both types of mice. Collectively, our data appear to show that affinity maturation may occur in TNFR1-deficient mice, which completely lack GCs.

Gene conversion in human rearranged immunoglobulin genes

Over the past 20 years, many DNA sequences have been published suggesting that all or part of the V(H) segment of a rearranged immunoglobulin gene may be replaced in vivo. Two different mechanisms appear to be operating. One of these is very similar to primary V(D)J recombination, involving the RAG proteins acting upon recombination signal sequences, and this has recently been proven to occur. Other sequences, many of which show partial V(H) replacements with no addition of untemplated nucleotides at the V(H)-V(H) joint, have been proposed to occur by an unusual RAG-mediated recombination with the formation of hybrid (coding-to-signal) joints. These appear to occur in cells already undergoing somatic hypermutation in which, some authors are convinced, RAG genes are silenced. We recently proposed that the latter type of V(H) replacement might occur by homologous recombination initiated by the activity of AID (activation-induced cytidine deaminase), which is essential for somatic hypermutation and gene conversion. The latter has been observed in other species, but not in human Ig genes, so far. In this paper, we present a new analysis of sequences published as examples of the second type of rearrangement. This not only shows that AID recognition motifs occur in recombination regions but also that some sequences show replacement of central sections by a sequence from another gene, similar to gene conversion in the immunoglobulin genes of other species. These observations support the proposal that this type of rearrangement is likely to be AID-mediated rather than RAG-mediated and is consistent with gene conversion.

CD5+ B cells are preferentially expanded in rabbit appendix: the role of CD5 in B cell development and selection

Although only a small proportion of mouse and human B cells are CD5(+), most adult rabbit B cells express CD5. However, CD5 was not detectable on the majority of B cells in neonatal appendix 1 and 3days after birth. Cell trafficking studies demonstrated that CD5(+) and CD5(-) CD62L(+) B cells from bone marrow migrated into appendix. There, CD5(+) B cells were preferentially expanded and predominated by approximately 2weeks of age. In mutant ali/ali rabbits, VHa2(+) B cells develop through gene conversion-like alteration of rearranged VH genes upstream of deleted VH1a2. Correlated appearance of individual CD5(+) germinal centers and VHa2(+) B-cells in mutant appendix suggests that CD5 binding positively selects cells with a2(+) framework regions that bind CD5. Following negative and positive selection, cells with diversified rearranged heavy- and light-chain sequences exit appendix, migrate to peripheral tissues and constitute the preimmune repertoire of CD5(+) B cells that encounter foreign antigens.

Developmental progression of immunoglobulin heavy chain diversity in sheep

In order to assess the respective impacts of combinatorial rearrangement, junctional diversification, somatic hypermutation and gene conversion in the generation of immunoglobulin heavy chain variable regions diversity, the sequences of 42 variable regions from late fetal, newborn and young sheep were determined and compared to those of adult animals. At earlier stages of development, the use of germline diversity segments appears restricted, junctional variability is already established, and somatic hypermutations are scarce. The sequence diversity in adults is much higher, which we suggest results from a higher hymermutation activity and possibly from the use of a variety of diversity segments. Altogether, this pattern is very reminiscent of the situation observed in cattle, except for the length of the third complementarity determining regions (CDR3) which are shorter in sheep than in bovine. Unlike the chicken and rabbit systems, it seems that new rearrangements continue to occur in sheep for at least several months after birth.

Immune system learning and memory quantified by graphical analysis of B-lymphocyte phylogenetic trees

The immune system learns from its encounters with pathogens and memorizes its experiences. One of the mechanisms it uses for this purpose is the intra-individual evolution of antigen receptors on B lymphocytes, achieved via hypermutation and selection of antigen receptor variable region genes during an immune response. We have developed a novel method for analyzing the graphical properties of phylogenetic trees of receptor genes which have been mutated and selected during an immune response. In the study presented here, we address the artifacts introduced by experimental methods of cell collection for DNA analysis, the meaning of each parameter measured on the tree graphs, and the differences between the dynamics of the humoral immune response in different lymphoid tissues.

Analysis of mutational lineage trees from sites of primary and secondary Ig gene diversification in rabbits and chickens

Lineage trees of mutated rearranged Ig V region sequences in B lymphocyte clones often serve to qualitatively illustrate claims concerning the dynamics of affinity maturation. In this study, we use a novel method for analyzing lineage tree shapes, using terms from graph theory to quantify the differences between primary and secondary diversification in rabbits and chickens. In these species, Ig gene diversification starts with rearrangement of a single (in chicken) or a few (in rabbit) V(H) genes. Somatic hypermutation and gene conversion contribute to primary diversification in appendix of young rabbits or in bursa of Fabricius of embryonic and young chickens and to secondary diversification during immune responses in germinal centers (GCs). We find that, at least in rabbits, primary diversification appears to occur at a constant rate in the appendix, and the type of Ag-specific selection seen in splenic GCs is absent. This supports the view that a primary repertoire is being generated within the expanding clonally related B cells in appendix of young rabbits and emphasizes the important role that gut-associated lymphoid tissues may play in early development of mammalian immune repertoires. Additionally, the data indicate a higher rate of hypermutation in rabbit and chicken GCs, such that the balance between hypermutation and selection tends more toward mutation and less toward selection in rabbit and chicken compared with murine GCs.

Rabbit immune repertoires as sources for therapeutic monoclonal antibodies: the impact of kappa allotype-correlated variation in cysteine content on antibody libraries selected by phage display

The rabbit immune repertoire has long been a rich source of diagnostic polyclonal antibodies. Now it also holds great promise as a source of therapeutic monoclonal antibodies. On the basis of phage display technology, we recently reported the first humanization of a rabbit monoclonal antibody. The allotypic diversity of rabbit immunoglobulins prompted us to compare different rabbit immune repertoires for the generation and humanization of monoclonal antibodies that bind with strong affinity to antigens involved in tumor angiogenesis. In particular, we evaluated the diversity of unselected and selected chimeric rabbit/human Fab libraries that were derived from different kappa light chain allotypes. Most rabbit light chains have an extra disulfide bridge that links the variable and constant domains in addition to the two intrachain disulfide bridges shared with mouse and human kappa light chains. Here we evaluate the impact of this increased disulfide bridge complexity on the generation and selection of chimeric rabbit/human Fab libraries. We demonstrate that rabbits with mutant bas and wild-type parental b9 allotypes are excellent sources for therapeutic monoclonal antibodies. Featured among the selected clones with b9 allotype is a rabbit/human Fab that binds with a dissociation constant of 1nM to both human and mouse Tie-2, which will facilitate its evaluation in mouse models of human cancer. Examination of 228 new rabbit antibody sequences allowed for a comprehensive comparison of the LCDR3 and HCDR3 length diversity in rabbits. This study revealed that rabbits exhibit an HCDR3 length distribution more closely related to human antibodies than mouse antibodies.

A comparison of hydraulic and laser capture microdissection methods for collection of single B cells, PCR, and sequencing of antibody VDJ

During the development of B lymphocytes, a series of gene rearrangements assemble the sequences that encode immunoglobulin heavy and light chains (VDJ). Earlier studies of VDJ sequence diversification during expansion of cells in splenic or appendix germinal centers used hydraulic micromanipulation (HM) to collect single B cells for PCR amplification of rearranged antibody heavy and light chain genes. PCR products were directly sequenced without a cloning step. Hydraulic micromanipulation is a very tedious method. Once capability to collect single cells by laser capture microdissection (LCM) was developed, we modified previous tissue staining and fixation methods so that we could collect cells from a given stained tissue section by HM and LCM and directly compare our success rates using these two methods. Cells were alkaline lysed and after two rounds of nested PCR products were recovered and directly sequenced. Because each rearrangement of genomic DNA that occurs to form the immunoglobulin heavy-chain-encoding sequence in developing B cells is unique, this system allowed us to verify our success rate in recovering single lymphocytes from tissue sections and amplifying a single allele. The methods developed have now made LCM an efficient alternative to HM for the collection of single B cells.

Distinct clonal Ig diversification patterns in young appendix compared to antigen-specific splenic clones

The young rabbit appendix is a dynamic site for primary B cell repertoire development. To study diversification patterns during clonal expansion, we collected single appendix B cells from 3- to 9-wk-old rabbits and sequenced rearranged H and L chain genes. Single cells obtained by hydraulic micromanipulation or laser capture microdissection were lysed, PCR amplified, and products directly sequenced. Gene conversion-like changes occurred in rearranged H and L chain sequences by 3-4 wk of age. Somatic mutations were found in the D regions that lack known conversion donors and probably also occurred in the V genes. A few small sets of clonally related appendix B cells were found at 3-5 wk; by 5.5 wk, some larger clones were recovered. The diversification patterns in the clones from appendix were strikingly different from those found previously in splenic germinal centers where an immunizing Ag was driving the expansion and selection process toward high affinity. Clonally related appendix B cells developed different amino acid sequences in each complementarity-determining region (CDR) including CDR3, whereas dominant clones from spleen underwent few changes in CDR3. The variety of combining sites generated by diversification within individual clones suggests that at least some clonal expansion and selection, known to require normal gut flora, may be driven through indirect effects of microbial components rather than solely by their recognition as specific foreign Ags. This diversity of combining sites within B cell clones supports the proposed role of appendix in generating the preimmune repertoire.

Characterization of monoclonal antibodies recognizing immunoglobulin kappa and lambda chains in pigs by flow cytometry

The existence of two types of the immunoglobulin (Ig) light chain in pigs was documented>30 years ago and has been confirmed by the cloning of porcine light chain genes homologous to human and murine Ig kappa (Igkappa) and Ig lambda (Iglambda). However, immunochemical reagents defining these two light chain isotypes have not been characterized. Here, we show that rabbit antisera specific for human Igkappa and Iglambda and certain anti-porcine light chain monoclonal antibodies (mAb) are useful in distinguishing light chain isotypes by flow cytometry (FCM). Porcine B cell lines L23 and L35 stained positive only with anti-human Iglambda antiserum and were negative when tested using anti-human Igkappa antiserum. While mAbs K139.3E1, 1G6 and 27.7.1 also tested positive on these cell lines, mAb 27.2.1 did not. Therefore, FCM was used to examine the hypothesis that K139.3E1, 1G6 and 27.7.1 are Iglambda-specific whereas mAb 27.2.1 recognizes the Igkappa chain in pigs. Double staining of peripheral blood mononuclear cells (PBMC) with pairs of anti-light chain mAbs and using cocktails of anti-light chain mAbs and anti-human polyclonal antiserum, confirmed this hypothesis with the exception that mAb K139.3E1 appears to recognize only a subset of Iglambda(+) B cells in most pigs. In summary, we identified two pan-specific anti-pig Iglambda mAbs, one anti-lambda mAb that recognizes a lambda-light chain subset and one anti-pig Igkappa mAb.