Restricted epitope specificity determined by variable region germline segment pairing in rodent antibody repertoires

ABSTRACT

Antibodies from B-cell clonal lineages share sequence and structural properties as well as epitope specificity. Clonally unrelated antibodies can similarly share sequence and specificity properties and are said to be convergent. Convergent antibody responses against several antigens have been described in humans and mice and include different classes of shared sequence features. In particular, some antigens and epitopes can induce convergent responses of clonally unrelated antibodies with restricted heavy (V_H) and light (V_L) chain variable region germline segment usage without similarity in the heavy chain third complementarity-determining region (CDR H3), a critical specificity determinant. Whether these V germline segment-restricted responses reflect a general epitope specificity restriction of antibodies with shared V_H/V_L pairing is not known. Here, we investigated this question by determining patterns of antigen binding competition between clonally unrelated antigen-specific rat antibodies from paired-chain deep sequencing datasets selected based solely on V_H/V_L pairing. We found that antibodies with shared V_H/V_L germline segment pairings but divergent CDR H3 sequences almost invariably have restricted epitope specificity indicated by shared binding competition patterns. This epitope restriction included 82 of 85 clonally unrelated antibodies with 13 different V_H/V_L pairings binding in 8 epitope groups in 2 antigens. The corollary that antibodies with shared V_H/V_L pairing and epitope-restricted binding can accommodate widely divergent CDR H3 sequences was confirmed by in vitro selection of variants of anti-human epidermal growth factor receptor 2 antibodies known to mediate critical antigen interactions through CDR H3. Our results show that restricted epitope specificity determined by V_H/V_L germline segment pairing is a general property of rodent antigen-specific antibodies.

KEYWORDS:

• B cell lineages
• binning
• convergence
• redundancy
• next-Gen Sequencing
• NGS

Introduction

Sequence and structural diversity in immunoglobulin variable regions enable antibodies to bind a virtually unlimited number of antigenic structures. Sequence diversity is generated somatically during the process of B cell maturation by recombination of germline-encoded heavy and light chain germline gene segments into full-length functional variable region exons.¹ Light chain variable region sequence diversity is generated by recombination of V_L and J_L gene germline segments along with relatively limited junctional diversity generated by nucleotide nibbling and incorporation between these segments. By contrast, heavy chain variable region sequence diversity is generated by recombination of three germline segments, V_H, D_H and J_H, which, along with significant nucleotide nibbling and incorporation, generates high sequence and length diversity in the third complementarity-determining region of the heavy chain (CDR H3).^2,3 The CDR H3 region is located centrally in the interface between antibody and antigen and usually provides critical contacts with the antigen. The CDR H3 region is considered to be a major determinant of antibody specificity due to its high sequence diversity and central role in antigen binding.^1,4

Immune repertoires, while highly diverse, have a significant degree of redundancy in specificities. This redundancy is mostly due to expansion of B cell clones harboring antigen-specific B cell receptors.⁵ Somatic mutation and clonal selection results in antibody clones with differing affinities for antigen and sometimes differing fine specificity for homologous antigens. However, these clonally related antibodies have essentially the same binding mode, defined as overall occluded antigen surface and binding geometry. Another level of redundancy is generated by what has been termed clonal convergence, in which antibodies from different B cell lineages share sequence elements that allow, or are presumed to allow, binding to the same antigen and epitope with similar binding modes.^6-9 Sequence features that are conserved in convergent antibodies vary and include shared V_H and/or V_L chain germline segment use, CDR H3 sequences or a combination of these. Convergent antibody responses have been observed in human and rodent immune responses to specific antigens and epitopes, including proteins, bacterial polysaccharides and viral antigens, allergens, blood group antigens, and haptens.^6,8,10-24 Antibody convergence is now a major principle behind rational human immunodeficiency virus type 1 vaccine design.²⁵

A relatively common convergence class includes groups of clonally unrelated antibodies with recurrent V_H and V_L germline segment pairing without CDR H3 sequence similarity or length conservation that bind certain antigens or epitopes.^8,10-18 These convergences have been described as V_H/V_L germline segment usage restrictions as a function of antigen or epitope specificity.^12-18 Whether these convergent responses are special cases in which certain V region germline segments are particularly suited for binding to certain epitopes independently of CDR H3 diversity or are examples of more widespread V_H/V_L germline segment convergences within antigen-specific repertoires is not clear. Here, we show that multiple epitopes in an antigen can induce antibodies with convergent V_H/V_L germline segment pairing in immunized rats. In addition, we show that shared V_H/V_L germline segment use in clonally unrelated antigen-specific antibodies is sufficient to predict shared epitope specificity with high confidence regardless of CDR H3 sequence and length, indicating that V_H/V_L germline segment pairing strongly restricts epitope specificity. Finally, we show by means of in vitro selection that the CDR H3 sequence and functional redundancy implied by widespread V_H/V_L germline segment pairing convergences can include clones with critical antigen contacts mediated by CDR H3. Our results extend our understanding of the role of V germline segment pairing and CDR H3 sequence diversity in determining epitope specificity in natural immune responses and provide insights about the organization of natural immune repertoires beyond clonal expansion.

Results

Germline segment pairing convergences are readily observed in a large anti-chicken ovalbumin antibody panel

A previously described panel of rat antibodies specific for chicken ovalbumin (OVA) generated by single-cell B-cell receptor sequencing (scBCR-seq) of antigen-sorted IgM^neg B cells pooled from 3 immunized rats was used in this study.²⁶ This panel includes 93 OVA-specific clones from different clonal lineages selected without regard for germline segment use, CDR sequence or length. As previously reported, the monovalent equilibrium dissociation constant (K_D) of the clones by surface plasmon resonance (SPR) ranged from less than 10 pM to 30 nM.²⁶ The scBCR-seq dataset includes 670 additional uncharacterized lineages, each lineage defined here as clones with the same heavy and light chain V and J germline segments, the same CDR H3 length, and at least 80% nucleotide identity within CDR H3 (IMGT® CDR H3 definition²⁷) between clones. We initially searched for correlations between sequence features and epitope specificity assessed by cross-competition binding experiments in an SPR format. In this assay, binding competition between two antibodies is determined by testing the binding of a “probe” antibody to immobilized antigen bound to a “capture” antibody (Suppl. Fig. 1). A subset of 18 clones showed cross-competition with several antibodies in different epitope groups (Suppl. Fig. 2), probably caused by recognition of alternative OVA conformational states.²⁸ These antibodies did not share V_H/V_L pairing with any other clones or amongst themselves and were excluded from further analysis. Six antibodies did not show detectable binding in this assay and were also excluded (data not shown).

Competition data was obtained for 71 of the clones, with at least 28 epitope groups identified for clones in competition assays (Suppl. Fig. 2). Eighteen of these clones could be classified in 7 groups of antibodies with recurring V_H/V_L germline segment pair use, ranging from 2 to 4 clones per group (Figure 1 and Suppl. Fig. 2, color-coded OVA clones). Within 6 of the 7 V_H/V_L pair groups, all clones cross-competed with each other in 4 epitope groups and, with some exceptions, behaved similarly in cross-competitions with other clones in the panel (Figure 1, Suppl. Fig. 2), indicating similar epitope specificity within each group. The exceptions were within the range of uncertainty for this technique, which may occasionally disagree for nearly identical clonally related antibodies, usually observed as patterns discordant between reciprocal competitions (e.g., clones OVA-80 and OVA-96, OVA-70 and OVA-13). Only one pair of clones with shared V_H/V_L germline segment pairing, clones OVA-36 and OVA-46, had cross-competition patterns with other antibodies clearly different from each other (Figure 1 and Suppl. Fig. 2, black background). We refer to clones with shared V_H/V_L pairing and epitope specificity as “parallel lineages” for simplicity.

Figure 1. Competitive binding between anti-OVA antibodies. Capture and probe antibodies refer to antibodies in solid phase used to capture OVA and antibodies in solution phase in the assay used to determine cross-competition, respectively. Red and red/black cross-hatched squares indicate no binding of the probe antibody. Green squares indicate binding of the probe antibody. Yellow squares indicate weak binding of probe antibody. OVA-36 and 46, highlighted in black, share IGHV5-25 and IGKV5S2 germline segment use but have distinct binding competition patterns. Groups of squares highlighted with white lines indicate competition of clones with the same V_H/V_L germline segment pairing. Names of antibodies within each V_H/V_L germline segment group are color-coded as in Suppl. Figure 2. Only antibodies with competition data as both capture and probe reagents as well as OVA-36 and 46 are shown, with the full dataset shown in Suppl. Figure 2

Sequences of antibodies within the parallel lineages were analyzed in detail to confirm different clonal origins. As expected, since these clones were originally selected to represent unique B-cell lineages, no significant CDR H3 sequence similarity was observed between clones within groups, with pairwise amino acid identities between clones with the same CDR H3 length below 67% (Figure 2). The independent clonal origin of the antibodies with the same CDR H3 length within groups were also confirmed by junctional DNA analysis (Suppl. Fig. 3). Our results show that convergences in the form of parallel lineages are readily apparent in a sufficiently large antibody panel.

Figure 2. V_H and V_L germline segments and CDR H3 sequences of anti-OVA clones in parallel lineages. Amino acid identities between 2 or more clones within the same parallel lineage with the same CDR H3 length are shown in black and gray backgrounds. Randomly selected clones that were and were not part of parallel lineages in the initial antibody panel are highlighted in orange and yellow, respectively. Clones selected based on V_H/V_L germline segment pair use are highlighted in blue

V_H/V_L germline segment pairing restricts epitope specificity of anti-OVA antibodies

The observation that most clones with shared V_H/V_L germline segment pairing recognize the same epitope suggests epitope specificity is restricted by V_H/V_L germline segment use. However, although unlikely, it is possible that the epitope restriction observed in the analyzed panel of 93 clones is due to chance. To test whether epitope restriction is determined by V_H/V_L germline segment pairing, additional clonally unrelated antibodies were sampled from the scBCR-seq dataset based solely on V_H and V_L germline segment use. Two classes of clones were selected. The first class includes clones with the same V_H/V_L germline segment pairings as clones in 5 of the 6 parallel lineages identified above. A second class of clones had the same V_H/V_L pairings as clones that were not part of parallel lineages in the initial panel (Figure 2, orange and yellow backgrounds, respectively, and Suppl. Fig. 2). These included clones in 2 epitope groups (2 and 14) not recognized by antibodies from parallel lineages. The only information used for clone selection from the scBCR-seq dataset was V_H and V_L germline segment use and independent clonal lineage. A total of 51 clonally unrelated OVA-specific antibodies covering 6 epitope groups and 10 V_H/V_L germline segment combinations were selected (Figure 2, blue background). Clones within the same parallel lineage with identical CDR H3 length had pairwise CDR H3 amino acid sequence identity of less than 72% (Suppl. Fig. 4). Clones with the same V_H/V_L germline segment pairing, CDR H3 length and J_H and J_L germline segments were confirmed to belong to different lineages by CDR H3 sequence identity and junctional analysis (Figure 2 and Suppl. Fig. 3). Affinities of selected clones ranged from 30 pM to 78 nM, within the same range as the clones picked randomly (Suppl. Fig. 5).

Forty-one of the clones were tested in two-way competitions with clones in other parallel lineages, as well as with clones selected to differentiate between similar epitope groups (Figure 3). The other 10 clones were tested in one- or two-way competitions with a subset of the clones in the same or other parallel lineages (Suppl. Fig. 6A and B). All clones within each of the putative parallel lineages showed restricted epitope specificity indicated by identical or nearly identical cross-competition patterns with antibodies in different epitope groups. As above, the observed discrepancies were limited to one-way patterns, typical of noise in panels of this size. Of the 56 antibodies within putative parallel lineages (51 selected plus 5 clones from the original panel not initially known to be in parallel lineages), only 5 clones showed patterns in two-way competitions that differed from other clones within the same parallel lineages. Four of these clones, OVA-100, OVA-105, OVA-104 and OVA-117, showed differences limited to 4 cross-competition data points relative to other members of their parallel lineages in two-way competitions with antibodies in other parallel lineages (Suppl. Fig. 7). However, clone OVA-120 showed patterns of cross-competitions different from other clones in its putative parallel lineage, although many of the discrepancies were only observed in one-way competition patterns (Figure 3). The atypical cross-competition pattern of clone OVA-120 was not related to polyspecificity; score of OVA-120 in a baculovirus virus particle ELISA (BV ELISA) to detect polyspecificity was 1.82, lower than the score of 5 for polyspecific antibodies.²⁹ Combined with the reference antibodies in the initial panel, 67 of the 74 (91%) anti-OVA antibodies with repeated V_H/V_L germline segment pairs showed virtually identical cross-competition patterns within each parallel lineage, 4 clones (5%) showed minor intra-group differences and only 3 clones (4%) showed competition patterns incongruent with V_H/V_L pairing.

Figure 3. Binding competitions between anti-OVA antibodies in parallel lineages. Red and red/black cross-hatched squares indicate no binding of the probe antibody. Green squares indicate binding of the probe antibody. Yellow and gray squares indicate weak and indeterminate binding of probe antibody, respectively. Squares with diagonal bars indicate discordant one-way competition results within parallel lineages. Clones are highlighted in orange, yellow and blue as in Figure 2. Clones with white background are not part of parallel lineages and were used to discriminate between epitope groups. Clone OVA-164 is a clonal variant of OVA-64. Competition patterns for reference antibodies of parallel lineages 8, 9 and 10 are shown in Suppl. Fig. 6B

V_H/V_L germline segment pairing restricts epitope specificity of antibodies against hen egg-white lysozyme

We extended our findings to a second antigen, hen egg-white lysozyme (HEL). An anti-HEL scBCR-seq dataset was obtained by similar methods used to derive the anti-OVA dataset. Significant nonspecific binding of labeled HEL was observed in sorting of B cells prior to sequencing (data not shown), resulting in only 20 of 96 selected clones showing binding to HEL by ELISA (Suppl. Fig. 8 and Suppl. File 1). Fourteen of these antibodies, along with 6 previously described anti-HEL mouse clones with known epitopes,^8,30-33 were tested in SPR-based competition assays. Reference mouse clones competed with each other as expected (Figure 4(a,b)). The HEL-specific rat clones derived from the scBCR-seq dataset bound in 11 overlapping epitope groups (Figure 4(a)). No recurring V_H/V_L germline segment pairs were observed in this relatively small panel of antibodies (Suppl. Table 1).

Figure 4. Epitope mapping of anti-HEL antibodies. (a) Competition patterns of anti-HEL clones tested. Red and red/black cross-hatched squares indicate no binding of the probe antibody. Green squares indicate binding of the probe antibody. Mouse antibodies of known epitope specificity are highlighted in shades of blue, orange, red and purple. Anti-HEL clones selected based on V_H/V_L pairing are shown with asterisks. White lines and dotted lines indicate cross-competition with other clones in parallel lineages and mouse anti-HEL antibodies. (b) Composite structure of HEL bound to reference mouse antibodies (PDB 1YQV, 1FDL, 1FBI, 1MLC, 1NDG and 1JHL). Variable regions of antibodies are color-coded as in panel a. (c) Graphical representation of epitope groups of antibodies in parallel lineages, mouse anti-HEL antibodies and additional selected rat anti-HEL antibodies. Antibodies and parallel lineages colored as in panels a and b

Two antibodies with different V_H/V_L germline segment pairs, Lyz-59 and Lyz-80, were selected to predict additional candidates binding the same epitopes based solely on V_H/V_L pairing data. Both antibodies compete with F9.13.7 and Hy-HEL-8. While Lyz-80 competes with both D11.15 and D1.3, which have closely overlapping epitopes, Lyz-59 does not compete with D11.15 (Figure 4(a–c)). Thus, subtle differences in competition between these two antibodies provide a more precise prediction of their cross-competition properties. A total of 5 and 4 lineages with HEL-specific antibodies with the same V_H/V_L germline segment pairing as clones Lyz-59 and Lyz-80, respectively, were identified, and one clone was selected from each lineage (Suppl. Table 1, Suppl. File 1). CDR H3 lengths ranged from 14 to 20 and 10 to 15 residues for clones in the Lyz-59 and Lyz-80 groups, respectively (Suppl. Table 1). Amino acid identity between 3 pairs of clones with the same CDR H3 length (Lyz-097/098, Lyz-100/101 and Lyz-80/111) ranged from 33% to 53%. These were also confirmed to belong to unique B cell lineages by junctional analysis (data not shown). Clones Lyz-80 and Lyz-111 shared an amino acid motif (TMGI) encoded by a rat D_H segment but shifted by one position, and some clones in the Lyz-59 group shared some similarities created by junctional variations of IGHD1-12 and IGHJ4, with IGHD1-12 in different registers (Suppl. Fig. 3 and Suppl. Table 1). All 9 clones within each parallel lineage showed identical patterns of cross-competition as the Lyz-59 and Lyz-80 clones (Figure 4(a,c)). This included identical competition patterns of the Lyz-59 group of clones with the reference D1.3 and D11.15 mouse antibodies with closely overlapping epitopes (Figure 4(b,c)). Our results indicate that restricted epitope specificity of antibodies sharing V_H/V_L germline segment use occurs for different antigens.

Synthetic parallel lineages of clones mediating critical interactions with antigen through CDR H3

The fact that convergences are readily detectable in large antibody panels in the form of parallel lineages, and the apparent epitope restriction of antibodies with shared V_H/V_L use, suggest that CDR H3 in parallel lineage clones may not play a critical role in antigen binding, allowing a wide range of possible CDR H3 sequences compatible with epitope binding. Alternatively, the diverse CDR H3 sequences may represent structurally distinct solutions for critical antigen interactions with a given epitope and binding mode. The latter possibility would imply that antibodies known to mediate critical antigen interactions through CDR H3 may accept distinct CDR H3 sequences for epitope binding while retaining epitope specificity.

We tested this with two humanized anti-human epidermal growth factor receptor 2 (Her2) antibodies binding different epitopes, hu4D5 and hu2C4, that are known to mediate critical antigen interactions through CDR H3.^34-37 The CDR H3 surface area buried in complex structures is 187 Ų for hu4D5 and 201 Ų for hu2C4. In hu2C4, mutation of residues Asn-108, Leu-109, Pro-111 or Tyr-114 (numbering according to the International ImMunoGeneTics information system®, IMGT®²⁷) to Ala individually results in reductions of about 100-fold in binding.³⁵ In hu4D5, the Trp-107 and Tyr-113 residue side-chains contact Her2 through hydrophobic interactions.³⁶ Both residues make large contributions to the free energy of binding, with ΔΔG values of 5.9 and 5.6 kcal/mol when individually mutated to Ala.³⁴

Variants of hu4D5 and hu2C4 were obtained from phage display libraries with randomized CDR H3 sequences, followed by testing of epitope specificity and affinity of variants with diverse CDR H3 sequences. Two sets of libraries were made for each clone. One design randomized CDR H3 sequences with NNK codons whereas the other used a mixture of trinucleotide codons (Suppl. Fig. 6). Both libraries were designed to retain the CDR H3 length of each parental clone. Libraries were panned on immobilized Her2 under low-stringency conditions to maximize sequence diversity in the output clones. Clones retaining antigen binding in ELISA and SPR after reformatting to IgG were selected for sequence analysis and functional characterization.

Despite the differences between randomization schemes and the relatively limited sequence space coverage in most libraries (Suppl. Fig. 9), results were comparable between the two libraries for each antibody and are thus presented in combination. Clones were grouped by sequence types, arbitrarily prioritizing patterns similar to the parental CDR H3 sequence that may reflect similar antigen interactions (Figure 5, Suppl. Fig. 10, Suppl. File 1). Two major consensus sequence groups were identified for hu4D5-derived clones. One group of 301 unique clones (hu4D5-g1) retained at least some of the critical contact residues of the parental antibody, substituting His or Tyr for the critical Trp in position 107 in some NNK-derived clones and all clones from the trinucleotide-based library (Figure 5(a)). Group hu4D5-g2 included 62 unique clones with strong consensus features within the group that did not retain recognizable similarity with the parental sequence (Figure 5(a)). An additional 98 unique clones in group hu4D5-g3 had limited consensus patterns (Figure 5(a)). A group of 13 clones in group hu4D5-g4 with CDR H3 sequences were one residue shorter relative to the parental clone (Figure 5(a)). Three major groups of consensus sequences were obtained for randomized hu2C4 clones. A major group of 79 unique clones (group hu2C4-g1) retained Ile/Leu in position 108, perhaps analogous to the Leu-108 residue making direct antigen contact in hu2C4,³⁷ but otherwise without recognizable similarity with the parental clone except perhaps Gly in position 109 (Figure 5(b)). Group hu2C4-g2 with 7 unique sequences had no apparent sequence similarity with the parental hu2C4 clone (Figure 5(b)). In addition, group hu2C4-g3 included 72 unique clones resembling the group hu2C4-g1 consensus, but with a CDR H3 length of 11 rather than 12 amino acid residues as in wild-type hu2C4 and group hu2C4-g1 clones (Figure 5(b)).

Figure 5. CDR H3 consensus sequences of unique anti-Her2 hu4D5 (a) and hu2C4 (b) variant clones. Parental CDR H3 sequence regions randomized in libraries are shown above the sequence logos, highlighting residues contacting Her2 through side-chains in red and through main-chain in blue. The number of unique clones for each group are shown

A total of 63 variants with higher affinity in SPR assays representing different consensus groups and unique, non-consensus patterns derived from hu4D5 and hu2C4 were expressed as human IgG1 and tested in competition assays with each of the parental antibodies. Not unexpectedly, as all clones retain non-CDR H3 residues that were optimized for binding in the parental antibodies, all variants competed only with the cognate parental antibody in SPR-based competition assays (data not shown), confirmed by competition assays with a subset of higher affinity clones from different sequence groups (Suppl. Fig. 11), indicating that epitope specificity remained unchanged in the variants relative to parental antibodies. The 63 variants were screened for affinity to Her2 and the best variants of each major consensus or unique group except the wild type-like hu4D5-g1 group were selected for detailed affinity analysis. Despite the almost complete absence of CDR H3 sequence similarity with parental antibodies, the best variants of many groups had affinities that were comparable to the parental antibodies (Table 1). Remarkably, these included hu4D5 variants that have small residues instead of Trp in position 107 and mostly non-aromatic residues in place of Tyr-113, the two CDR H3 residues that make critical direct hydrophobic contacts with Her2 through side-chains and that are critical for binding. In addition, the strongly selected hu2C4 variants with shorter CDR H3, which are a minor fraction of the library caused by inefficiencies in trinucleotide incorporation, were among the best hu2C4 variants, with affinities similar to the parental clone. The results indicate that significant CDR H3 sequence diversity is allowed in antibodies known to mediate critical interactions with antigen through this region while retaining epitope specificity and robust antigen binding.

Table 1. CDR H3 sequence and binding kinetics of best anti-Her2 variants in different consensus groups

Clone^a	CDR H3^b													ka (M^−1s^−1)	kd (s^−1)	KD (nM)
hu4D5	A	R	W	G	G	D	G	F	Y	A	M	D	Y	7.2 × 10^5 [1.06]^c	1.1 × 10^−4 [1.04]	0.13 [1.10]
hu4D5-g2.1	A	R	S	H	K	N	V	M	G	L	V	D	Y	1.6 × 10^6 [1.12]	3.5 × 10^−4 [1.09]	0.27 [1.13]
hu4D5-g2.2	A	R	S	H	A	L	V	Y	G	L	I	D	Y	1.8 × 10^6 [1.21]	4.3 × 10^−4 [1.08]	0.33 [1.27]
hu4D5-g3.1	A	R	G	H	Y	S	S	G	H	P	F	D	Y	9.6 × 10^5 [1.09]	3.2 × 10^−4 [1.03]	0.39 [1.10]
hu4D5-g3.2	A	R	G	I	G	N	G	A	E	G	M	D	Y	6.6 × 10^5 [1.05]	1.1 × 10^−3 [1.11]	1.8 [1.17]
hu2C4	A	R	N	L	G	P	–	S	F	Y	F	D	Y	5.5 × 10^4 [1.04]	2.2 × 10^−4 [1.04]	4.6 [1.14]
hu2C4-g3.1	A	R	L	M	S	Y	–	–	G	Y	F	D	Y	6.3 × 10^4 [1.17]	4.0 × 10^−4 [1.01]	6.4 [1.16]
hu2C4-g3.2	A	R	L	L	A	Y	–	–	G	Y	L	D	Y	5.8 × 10^4 [1.15]	5.7 × 10^−4 [1.09]	9.9 [1.07]
hu2C4-g1.1	A	R	A	I	A	G	–	G	V	L	L	D	Y	4.8 × 10^4 [1.16]	1.2 × 10^−3 [1.03]	25 [1.13]
hu2C4-g1.2	A	R	A	I	T	N	–	E	V	L	L	D	Y	4.9 × 10^4 [1.10]	1.4 × 10^−3 [1.01]	28 [1.07]
hu2C4-g2.1	A	R	S	S	T	P	–	P	F	F	M	D	Y	3.9 × 10^4 [1.15]	9.6 × 10^−4 [1.03]	25 [1.11]
hu2C4-g2.2	A	R	S	S	R	P	–	P	W	F	F	D	Y	6.5 × 10^4 [1.23]	1.5 × 10^−3 [1.01]	23 [1.25]

^a Clone name suffixes indicate groups according to Figure 5.

^b IMGT® positions 105 to 117 shown. Underlined boldface residues contact antigen through side-chains in parental antibodies. Highlighted hu4D5 residues are IMGT® W107 and Y113 (Kabat W95 and Y100a). Highlighted hu2C4 residues are IMGT® L108, P110 and Y114 (Kabat L96, P98 and Y100a). Residues shown in italics were not randomized. Dashes introduced to conform all clones to the IMGT® numbering system.

^c Values shown indicate geometric means and geometric standard deviation factors (in brackets) of experiments performed four times (one n = 1 and one n = 3 experiment).

Antibodies with shared V_H/V_L germline segment use and different antigen specificities have different binding modes

A possible basis for the restricted epitope specificity of antibodies with shared V_H/V_L germline segment pairs is that these antibodies bind antigen through stereotyped interactions, using a defined subset of residues in the binding interface. We addressed this question by analyzing antigen-binding interfaces in complex structures of antibodies with shared V_H/V_L germline segment use but different antigen or epitope specificities, including two previously described convergent human antibodies specific for the respiratory syncytial virus F (RSV F) glycoprotein and a human antibody with the same IGHV3-21/IGLV1-40 germline segment pairing that binds the unrelated proprotein convertase subtilisin-kexin type 9 (PCSK9).^11,38-41 All three antibodies have similar CDR H3 lengths and occur in a germline or near-germline configuration (Suppl. Fig. 12), allowing direct comparison of contacts by germline residues in the binding interface. Given the role of CDRs in antigen binding, several non-CDR H3 interface residues are shared between the anti-RSV F and anti-PCSK9 antibodies as determined by buried surface area per position (Suppl. Fig. 13). However, while the interface profiles of both anti-RSV F antibodies are nearly identical outside CDR H3 and the junctional region of light chain CDR3, significant differences are observed between the anti-RSV F and anti-PCSK9 buried surface area profiles outside CDR H3, especially in the light chain and CDR H2. Thus, the overall binding mode is not necessarily conserved between antibodies with shared V_H/V_L germline segment pairs binding different epitopes.

Discussion

Numerous reports have described convergent responses in human and mouse antibody repertoires.^6,8,10-24 One particular convergence class has been described as a restriction in V_H and V_L germline segment usage determined by antigen or epitope specificity.^12-18 Collectively, these reports show that a number of epitopes or antigens can induce convergent immune responses with recurring germline segments. Repertoire deep sequencing studies suggest that convergences in CDR H3 sequences with or without V_H germline segment recurrence, which are different from parallel lineages, may be prevalent.^21,42-46 However, with two exceptions,^21,42 there is no data to confirm the functional significance of these convergences. For V_H/V_L convergences, descriptions have been based on incidental findings of antibody characterization rather than systematic searches for these convergences.^10-23 Thus, while convergences have been described with many antigens, no report has described the extent by which different parallel lineage convergences can occur for a single antigen or, to our knowledge, predicted this convergence class to be pervasive in immune repertoires. In addition, whether V_H/V_L germline segment pairing also has a restrictive role on epitope specificity has not been previously addressed.

Here, we show, based on systematic analysis of convergences in antibody repertoires derived by high-throughput paired-chain sequencing, that V_H/V_L pairing parallel lineage convergences are common within a single antigen and that V_H/V_L pairing restricts epitope specificity. Of note, epitope restriction by germline segment use is not a logical consequence of convergent responses, structural principles or biological processes. Epitope restriction of antibodies with recurrent V_H/V_L pairing is evident because V_H/V_L convergences are prevalent while, at the same time, antibodies with the same V_H/V_L germline segment binding different epitopes in the same antigen are rare. The epitope restrictions would be unremarkable in the absence of pervasive parallel lineage responses. The widespread occurrence of parallel lineage convergences in natural repertoires led us to investigate whether highly redundant CDR H3 sequences are compatible with epitope binding in the context of a V_H/V_L germline segment pair. We show here by means of phage display that functional CDR H3 redundancy includes clones known to mediate critical interactions with antigen through CDR H3, consistent with the observed widespread occurrence of parallel lineages and possibly other types of convergences.

Shared V_H/V_L germline segment pairing, which may not seem very informative due to the relatively limited number of different possible pairings, is highly predictive of functional redundancy within antigen-specific repertoires. In total, 78 (92%) of 85 clonally unrelated antigen-specific antibodies with 13 different V_H/V_L germline segment pairings binding in 8 epitope groups in 2 antigens had restricted epitope specificity as determined by competition assays. Four additional antibodies had only minor differences in cross-competition patterns relative to other antibodies in the same parallel lineages, bringing the agreement between V_H/V_L pairing and epitope specificity to 96%. Importantly, for 64 of these antibodies the epitope restriction was identified prospectively, with clones selected based only on V_H/V_L germline segment pairing. Of the 64 antibodies, 34 had 7 different V_H/V_L pairs not known a priori to belong to epitope-restricted groups. In addition, for both antigens, epitope restriction as a function of V_H/V_L use was observed in the context of a relatively large diversity of epitope specificities in the repertoires. Thus, our dataset represents clones with shared V_H/V_L germline segment pairing that were selected randomly, except for the availability of repeated V_H/V_L pairs in the sequence sets. Based on this, we expect our results showing restricted epitope specificity in clonally unrelated antibodies sharing V_H/V_L germline segment pairing to be largely representative of unsampled antigen-specific clones in the datasets described here and in antigen-specific repertoires for other antigens.

Epitope specificity determination by cross-competition is a relatively low resolution epitope mapping technique. Antibodies that cross-compete can have significantly different overlapping epitopes or cross-compete through constant region bulk. However, the virtually indistinguishable cross-competition patterns of a relatively large number of antibodies with recurring V_H/V_L germline segment pairs to several epitopes indicates a close epitope similarity among antibodies in parallel lineages. A virtually identical binding geometry for clones within parallel lineages is the most parsimonious interpretation for the generally restricted epitope specificity of antibodies with the same V_H/V_L germline segment pair. Details in the interactions may differ among clones, especially in CDR H3, which can vary both in sequence and length between antibodies in parallel lineages and may cover slightly different epitope areas. However, overall binding geometry and many of the interactions would be expected to be similar or identical in most or all cases.

The generalized epitope restriction of antibodies with the same V_H/V_L pair that we observed would not be congruent with binding of these antibodies to partially overlapping epitopes or with significantly different rotational angles relative to the antigen, as these would be effectively as different structurally as non-overlapping epitopes. If that were the case, a high prevalence of non-overlapping epitopes would be expected within sets of antibodies with the same V_H/V_L pairing, unlike what is observed here. Consistent with our interpretation of conserved binding geometries in parallel lineages, virtually identical binding geometries have been observed in structures of antibodies with the same V_H/V_L germline segment pairs in complex with protein or peptide antigens.^8,11,47,48 It should be pointed out that conserved binding geometry may not necessarily extend to constant regions. This is because angles between variable and constant regions may differ with framework sequences encoded by J_H germline segments,^49,50 which do not necessarily recur within parallel lineages.

An unexpected implication of our results showing restricted epitope specificity determined by V_H/V_L germline segment pairing is that the number of distinct epitopes recognized by antibodies from immune repertoires generated in vivo with a particular V_H/V_L pair is apparently low, probably less than one epitope on average per small- to medium-size antigen with multiple epitope groups (HEL, 14 kDa; OVA, 43 kDa). This is unexpected, as no structural principle requires V_H/V_L pairs to have significantly restricted epitope specificities in the context of the highly variable CDR H3 region. A more general specificity preference based on light chain CDR 1 length, which correlates with binding of antibodies to different antigen classes due to different overall shape complementarity requirements, was previously described.⁵¹ However, the restricted epitope specificity based on V_H/V_L germline segment use shown here includes antibodies binding the same antigen class, globular proteins, with similar overall epitope shapes. Thus, epitope restriction based on V_H/V_L germline segment use is not expected based on that precedent.

As currently understood, V_H/V_L germline segment pairs recur in antibodies specific for different antigens, and thus have a broad range of potential specificities when combined with different CDR H3 sequences. In fact, it has been shown that CDR H3 diversity in the context of limited V_H/V_L germline segment diversity is sufficient for a relatively broad range of antigen specificities,⁴ although the number of unique antibodies reported for each antigen is low and epitope coverage in these cases is not known. It is possible that epitope restriction may be only partially valid with larger or more complex antigens such as whole viruses and cells. That is, clones with recurrent V_H/V_L germline segment pairs may have more than one distinct epitope specificity in more complex antigens. However, our data suggest that these would be restricted to a relatively limited set of epitopes. One point to note is that the epitope restriction by V_H/V_L germline segment pairing probably happens only in the context of naïve repertoires without somatic diversification prior to antigen stimulation, as it occurs in rodents and humans. Somatic mutation after antigen stimulation is constrained to retain epitope specificity except in cases of secondary responses to cross-reactive epitopes in different antigens, which is unlikely to occur within an antigen. However, in species in which the naïve repertoire undergoes somatic diversification of V germline segment regions prior to antigen stimulation, such as rabbits, the spectrum of potential epitope specificities of V_H/V_L germline segment pairs in a single antigen may be expanded. This is in fact observed with antibodies obtained from libraries with V regions randomized, or “somatically mutated”, in vitro.⁵²

A second unexpected implication of the apparent widespread occurrence of parallel lineages in repertoires is the functional redundancy of CDR H3. A certain level of sequence similarity is observed between the CDR H3 of clones in some parallel lineages, as might be expected from their shared epitope specificity. Many of these similarities involve lower diversity sequences contributed by J_H germline segments as well as the region contributed by V_H (IMGT® positions 105 and 106), which are invariant in the absence of nucleotide nibbling or somatic mutation and could be expected with relatively high frequency by chance. In addition, these similarities usually involve only subsets of clones within each parallel lineage, and the positions of amino acid identities in CDR H3 vary in different pairwise comparisons. Thus, the significance of many of these sequence similarities is evident only with the benefit of data indicating functional similarities of a group of clones. Clones in this category include anti-OVA parallel lineages 7 and 10 and subsets of other parallel lineages, as well as some clones in both anti-HEL parallel lineages (Figure 2 and Suppl. Table 1). However, the level of pairwise amino acid identity between clones with the same CDR H3 length is generally low (Suppl. Fig. 4).

The widespread occurrence of restricted epitope specificity in parallel lineages in the context of CDR H3 regions of different lengths and sequences implies that CDR H3 in parallel lineages may not be as important energetically for binding or, alternatively, may provide the necessary free energy of binding and epitope specificity with a range of distinct sequences and lengths. Rather than surveying the binding contribution of CDR H3 in parallel lineage clones, we tested whether the well-characterized hu4D5 and hu2C4 anti-Her2 antibodies known to mediate energetically important interactions with antigen through CDR H3 would accept widely distinct CDR H3 sequences while retaining epitope specificity and robust binding. We found that, contrary to a recently stated expectation that complete randomization of hu4D5 CDR H3 would be unlikely to yield a significant number of Her2 binders,⁵³ CDR H3 sequences that differ significantly from the parental antibodies can accommodate strong binding to the same epitope as determined by binding competition assays. As CDR H3 in these synthetic parallel lineage clones must provide much of the free energy of binding that is contributed by the corresponding parental CDR H3 regions, our results indicate a significant degree of sequence and functional redundancy in CDR H3. By extension, our results indicate that natural parallel lineages may include clones with CDR H3 regions with major energetic contributions to binding, explaining the relatively high frequency of these convergences in repertoires.

Parallel lineages imply a major role of V-germline segment residues in the energetics of binding. In fact, non-CDR H3 regions of natural antibodies have been shown to contribute a large fraction of binding residues and free energy of binding, although the relative contribution of germline segment and somatically mutated residues has not been determined.^54,55 Previous descriptions of structures of convergent antibodies in complex with antigen have uncovered some shared mechanisms of binding based on specific germline segment features. This is exemplified by the binding of convergent antibodies against influenza A hemagglutinin and the Staphylococcus aureus virulence factor IsdB with IGHV1-69 or IGHV4-39 germline segments, which make use of hydrophobic germline residues in CDR H2 for antigen binding.^22,56 However, other convergent clones against similar influenza A hemagglutinin epitopes using other germline segments lack that signature or any other sequence and structural signatures in common between different convergent groups,²³ suggesting different structural mechanisms for different convergent groups. Given the variety of germline segments and heavy and light chain subgroups involved in the convergences shown here, it is unlikely that these share a common structural basis, except for CDR H3 sequence and functional diversity. CDR H3 redundancy was previously examined in specific cases of convergences. In some cases, conservation of specific or chemically similar residues in certain CDR H3 positions was observed.^11,23,56 However, in other cases no conserved CDR H3 signatures were observed,^11,22 similar to our results with the hu4D5-g2 and hu2C4-g2 group clones and many of the parallel lineages described here.

One difficulty of determining the structural basis for restricted epitope specificity is that the absence of epitope specificities in the same antigen cannot be directly addressed by structural studies. Thus, the structural basis for this restriction may only be inferred. One possibility is that a given V_H/V_L pair can only interact with different epitopes through specific stereotyped interactions. Our analysis of the binding modes of different antibodies with shared V_H/V_L germline segment use and different antigen specificities indicates that this is not necessarily the case (Suppl. Fig. 13), which is not unexpected. Determining whether subtle antigen interactions are recurrently used by antibodies with the same V_H/V_L pairing with different epitopes will require systematic large-scale analyses of complex structures of analogous sets of antibodies with shared V_H/V_L pairing binding different antigens. However, it is possible that the generally restricted epitope specificity of V_H/V_L pairs within antigens may not be based on stereotypic antigen binding interactions and that this restriction may be best described statistically. Thus, epitope restriction by V_H/V_L germline segment pairing within antigens is a repertoire feature that emerges from large-scale functional analyses that, while structurally determined, may not be readily apparent from a strictly structural perspective.

Whether any V_H/V_L pairs in a set of antigen-specific clones can be potentially part of parallel lineages is not clear. About half of the anti-OVA scBCR-seq dataset, an estimated 96% of which bind OVA,²⁶ is composed of lineages with unique V_H/V_L germline segment pairs (Suppl. Fig. 14). Whether these are unique due to undersampling, relatively low frequency of compatible CDR H3 sequences, immune regulation or stochastic factors cannot be easily determined. The chemistry of the antigen interactions mediated by CDR H3 may affect the prevalence of CDR H3 sequences meeting the structural requirements for binding and hence parallel lineage formation. However, CDR H3 length does not seem to be a major factor in parallel lineage convergences within species. The CDR H3 lengths of the parallel lineage clones we characterized largely reflect the CDR H3 lengths of clones in the anti-OVA repertoire, which are slightly longer on average than the CDR H3 length of naïve rat IgG repertoires (Suppl. Fig. 15). Similarly, CDR H3 lengths of previously described human and mouse antibodies with convergent parallel lineages range from 15 to 18 and 9 to 17 residues, respectively,^{10-12,15,17,41} slightly longer than the average CDR H3 length of each species.⁵⁷ Thus, parallel lineage convergences do not seem to be generally associated with readily apparent CDR H3 length biases within species. Whether different species have different propensities for parallel lineage convergences and germline segment-restricted epitope specificity due to differences in average CDR H3 length remains to be determined.

Our results extend our understanding of the role of V_H/V_L pairing and CDR H3 diversity in determining epitope specificity. The CDR H3 functional redundancy that results in widespread parallel lineage convergences, combined with the highly predictive epitope restriction of V_H/V_L germline segment pairings, reveal a logical organization of immune repertoires beyond classical B cell lineages that is determined by structural rather than biological factors, with implications for vaccine development and antibody discovery and engineering.

Materials and methods

Selection of anti-OVA and anti-HEL clones from scBCR-seq datasets

A dataset of anti-chicken ovalbumin clones derived by scBCR-seq with paired heavy and light chain variable regions was previously described.²⁶ A subset of 93 clones selected from the scBCR-seq dataset were shown to be specific for OVA, numbered here as “OVA” followed by a two-digit number. Additional clones, indicated by “OVA” followed by a three-digit number, were selected from different B-cell lineages in scBCR-seq datasets using V_H and V_L germline segment use data and prioritizing clones within lineages with full-length J_H sequence information. Lineage definition for the anti-OVA panel was the same as previously described, using same V_H and V_L germline segments and at least 80% nucleotide identity within CDR H3, ignoring J germline segment identity, to include clones within the same lineage.²⁶ A similar definition was used for the anti-HEL panel, using 80% amino acid identity within CDR H3 instead as the threshold cutoff. The junctional nucleotide sequences of selected clones within the same lineage and 50% or greater amino acid sequence identity were further analyzed for V_H, D_H, and J_H use and boundaries using an in-house tool and a composite database of IMGT®²⁷ and in-house determined reference germline segments sequences as previously described.²⁶ Clones with the same D and J segments in the same location and similar boundaries were excluded from the dataset as likely clonal variants modified by somatic mutation. All clones were expressed as rat/human chimeric immunoglobulin G after DNA synthesis as previously described.²⁶

Sequencing of rat anti-HEL repertoires

A dataset of B-cell receptor sequences enriched for HEL-specific clones with paired heavy and light chain information was obtained as previously described.²⁶ Briefly, three Sprague Dawley (SD) rats (Charles River, Hollister, CA) were immunized with 100 µg HEL (Sigma-Aldrich, St. Louis, MO) in Complete Freund’s adjuvant (BD, Franklin Lakes, NJ) followed by bi-weekly boosts of 50 µg OVA in incomplete Freund’s adjuvant divided in three sites (intraperitoneally, subcutaneously at base of tail, at nape of neck and in both hocks). Cells from multiple lymph nodes from the 3 rats were pooled, depleted of IgM^pos B cells, stained with anti-rat IgM (Clone G53-238, BD Biosciences, San Jose, CA) conjugated to PE Cy7 and HEL conjugated to Alexa Fluor 647 (Thermo Fisher Scientific, Waltham, MA) and sorted for IgM^neg/HEL^pos cells in a FACSAriaIII sorter (BD, Franklin Lakes, NJ). Due to visible nonspecific binding of labeled antigen in control T cells, gating for cell selection was set at a high threshold to enrich for antigen-specific B cells. A total of 49,000 IgM^neg/HEL^pos events were detected from the three rats and processed for scBCR-seq. Sequences were parsed for germline segments and processed for paired V_H/V_L sequences and clones grouped by putative lineages as described below.

Selection of anti-Her2 CDR H3 variants by phage display

The variable regions of anti-Her2 hu4D5 and hu2C4 were cloned in a monovalent antigen-binding fragment (Fab) phage display vector.⁵⁸ The CDR H3 region of both clones was replaced by site-directed mutagenesis with stop codons and two sets of libraries with randomized CDR H3 sequences were produced for each clone by oligonucleotide-directed mutagenesis as previously described.⁵² CDR H3 randomization was achieved by use of NNK (IUPAC code) or trinucleotide codon mixes. The trinucleotide-based randomization was designed to partly mimic human CDR H3 diversity with position and length-specific codon mixtures (Suppl. Fig. 9). Clones were selected on Her2 directly immobilized on ELISA plates (hu4D5) or on biotinylated Her2 immobilized on streptavidin-coated ELISA plates (hu2C4). Libraries were panned for up to 5 rounds as previously described, selecting clones from rounds 3 and 5 that were confirmed to bind Her2 in phage-ELISA for Sanger sequencing.³⁵ Sequence logos were generated using the WebLogo server.⁵⁹

Sequence coverage of the anti-Her2 variant libraries was calculated based on library randomization scheme, library colony-forming units (cfu) achieved and fraction of successfully mutated clones in the library determined by Sanger sequencing of a library clone sample. Briefly, because the expected frequencies of different sequences in the library are unequal, the sequence coverage was calculated by generating j (10⁴ in this case) sequences in a Monte Carlo simulation using the probabilities of each amino acid being yielded by a given codon and calculating the expected frequency P of each sequence i. This generates a set in which each simulated sequence i represents L/j sequences of identical or very similar expected frequencies, where L is the number of successfully mutated (non-template) clones in the library. The average number of unique sequences C in a library size L is approximated as:

(1) $C\cong \sum _{i=1}^j\frac{1-e^{-L\times P_i}}{j\times P_i}$(1)

The contribution to C by clones with stop or Cys codons, which are rarely selected by phage display and effectively defective, is defined as zero here. The estimated percentage F of the total theoretical diversity V of clones without stop or Cys codons achieved in a library (functional sequence coverage) is therefore approximately

(2) $F\cong \frac{C}{V}\times 100$(2)

Note that if P_i is constant, as in an idealized library in which each sequence variant has the same expected frequency P_i = 1/V, equation (1)(1) $C\cong \sum _{i=1}^j\frac{1-e^{-L\times P_i}}{j\times P_i}$(1) simplifies to the previously described equation⁶⁰

(3) $C=V\left(1-e^{-\frac{L}{V}}\right)$(3)

Expression and purification of antibodies and Fabs

DNA fragments encoding variable regions of anti-OVA and anti-HEL scBCR-seq sequences were obtained by DNA synthesis and cloned into human immunoglobulin class 1 (IgG1)/kappa mammalian expression vectors as previously described.²⁶ The heavy chain variable region of selected anti-Her2 clones were subcloned into human IgG1 and flag-tagged Fab mammalian expression vectors. Heavy and light chain expression plasmids were co-transfected into Expi293 cells at 1 ml scale and purified in batch mode by protein A chromatography as previously described.^61,62 Fabs of hu4D5 variants were transiently expressed in 30 ml Expi293 cells for 7 days and purified by anti-flag-tag and size-exclusion chromatography. Fabs were confirmed to be monomeric by analytical size-exclusion chromatography.

Affinity and binding competition assays

Binding of clones to OVA or HEL was confirmed by ELISA and affinities determined by SPR as previously described.²⁶ Purified hu2C4 variant IgG1 was tested for binding to Her2 by SPR at 25°C using a BIAcore T200 apparatus (GE Life Sciences, Piscataway, NJ), immobilizing IgG in protein A chips and using monomeric Her2 as analyte. The hu4D5 variants were tested by SPR in BIAcore T200 apparatus with Her2 immobilized on CM5 chips, using Fabs as analytes at 25°C. Binding competition assays were performed in Carterra IBIS and LSA instruments as previously described.⁶³ Binding competition results were analyzed in a blinded manner without reference to antibody sequence groups using the Carterra Epitope Binning software following standard vendor guidelines, previously described in detail.^64,65

Polyspecificity assay

Polyspecificity of antibody OVA-120 was performed in a BV ELISA assay as previously reported,²⁹ using IgG purified by protein A and size-exclusion chromatography in the assay.

Determination of antigen binding interfaces

Buried surface area of antibody regions and residues in the interface with antigen (excluding crystal contacts) were determined using the PDBePISA server at http://www.ebi.ac.uk/pdbe/pisa.⁶⁶ PDB files with accession numbers 1N8Z (hu4D5), 1S78 (hu2C4), 6APD (anti-RSV F ADI-19425), 5U68 (anti-RSV F MPE8) and 3H42, 4K8R and 5OCA (anti-PCSK9) were used in the analyses.

Abbreviations

CDR H3	=	Heavy chain third complementarity determining region
cfu	=	Colony-forming units
HEL	=	Hen egg-white lysozyme
IgG1	=	Human immunoglobulin class 1
IMGT®	=	The International ImMunoGeneTics information system®
OVA	=	Chicken ovalbumin
PCSK9	=	Proprotein convertase subtilisin-kexin type 9
RSV F	=	Respiratory syncytial virus F glycoprotein
scBCR-seq	=	Single-cell B-cell receptor sequencing
SD	=	Sprague Dawley
SPR	=	Surface plasmon resonance
VH	=	Heavy chain variable region germline segment
VL	=	Light chain variable region germline segment

Disclosure of potential conflicts of interest

All authors except A.A. are current or former employees of Genentech and may be Roche shareholders. A.A. is an employee of SciGenom Labs. S.S. is affiliated with SciGenom Research Foundation.

Acknowledgments

We thank Bijay S. Jaiswal and Jeremy Stinson for project support and Jack Bevers III for assistance with Carterra-based epitope competition assays.

Supplementary material

Supplemental data for this article can be accessed on the publisher’s website.