Massive Parallel Sequencing of Full-Length B-Cell Receptor Sequences Reveals HLA-Dependent Shaping of the B-Cell Immune Repertoire
Authors: Koning Marvyn T. , van der Zeeuw Sander A.J. , Navarrete Marcelo , van Bergen Cornelis A.M. , Nteleah Valeri , Griffioen Marieke , Kiełbasa Szymon M. , Veelken Hendrik 
Published: October 2019

Abstract
Peptides of the B-cell receptor (BCR) may be presented in HLA molecules and therefore be recognized as epitopes by T cells. Bioinformatic evidence indicates that follicular lymphoma cells are selected against expression of a clonal BCR with a high cumulative predicted binding of BCR-derived peptides to the respective patient's HLA complex (Strothmeyer, Blood 2010). This observation suggests T-cell-mediated immunosurveillance against outgrowth of follicular lymphoma cells according to BCR HLA binding strength.
Here, we investigate whether this phenomenon pertains to peripheral B cells in 6 healthy donors: 2 donors homozygous for HLA A01*01 / B08*01, 2 homozygous for HLA A02*01 / B7*02, and 2 donors heterozygous for these alleles. Unbiased representation of full-length V(D)J sequences was considered essential for correct data interpretation. PCR primers annealing to conserved motifs of BCR variable regions (e.g. BIOMED-2 protocol) fail to amplify a fraction of BCR, particularly those modified by somatic hypermutation. Therefore, we developed an improved anchored PCR strategy: cDNA was synthesized from poly(A)-RNA from peripheral blood with primers that anneal to specific Ig constant regions. In the same reaction, the 3' cDNA end is extended by switching to an oligonucleotide template containing an anchor sequence (SMART technology; Clontech). Anchor-tagged cDNA was amplified with a primer annealing to the anchor in combination with a nested constant region-specific reverse primer. Dumbbell adapters were added to the termini of 250 ng of purified PCR products. Circular consensus sequencing of single molecules was performed on the PacBio platform (Pacific Biosciences). Using one SMRT PacBio cell per amplicon, separate sequence libraries were created for μ, γ, κ, and λ BCR transcripts. Sequences covered by at least five reads were selected with SMRT Portal software to obtain >95% of sequences without sequence errors as demonstrated on multiple B-cell lines. Selected sequences were analysed by HighV-QUEST software (Alamyar, Immunome Research 2012). After exclusion of non-BCR sequences and duplicate BCR transcripts, a median of 5318 (range: 670-8752) individual BCR sequences was obtained per library. Binding affinity of nonamers in in-silico-translated BCR were calculated for the 4 HLA alleles by the NetMHC 3.4 algorithm. The fractions of BCR lacking any weak HLA binding peptide (NetMHC IC50 <500nM) within a library were compared between donors positive or negative for any HLA molecule.
μ VDJ transcripts without HLA binding peptides were significantly more frequent for all HLA alleles in donors that actually express that particular allele (Table). With the exception of HLA A01*01, similar results were observed for γ transcripts. While the fraction of κ VJ transcripts without an HLA binder was overall higher in HLA A01*01 and B08*01, HLA-positive individuals had higher proportions of non-HLA binding sequences. λ transcripts were less likely to contain HLA binders with respect to HLA B07*02 and B08*01 but not to the HLA A alleles.
Analogous analyses were performed for CDR3 regions as annotated by HighV-QUEST plus six amino acids on either flank. In 10 of 16 analyses, CDR3 sequences were less likely to contain an HLA binder in HLA-positive individuals; in three analyses an opposite effect was seen (Table).
These results indicate that the peripheral BCR repertoire is shaped by HLA alleles in healthy individuals, most likely by T-cell mediated recognition of BCR peptides. Ongoing studies expand this fundamental finding with respect to the IC50 threshold, the number of nonamers, and additional HLA alleles. Our results warrant investigation of the potential role of HLA-dependent shaping of the BCR repertoire for the immune defense and the development of autoimmune disease and B-cell lymphoma. Table 1V(D)J without HLA binding peptideCDR3 without HLA binding peptideHLADonorμγκλμγΚλ A01*01Positive21%41%61%37%87%90%98%70%Negative16%42%59%38%92%92%96%65%P<0.001n.s.<0.01n.s.<0.001n.s.<0.01<0.001 A02*01Positive6%4%3%32%77%77%77%70%Negative4%1%2%32%75%69%78%78%P<0.001<0.001<0.01n.s.<0.01<0.001n.s.<0.001 B07*02Positive31%13%3%13%79%73%91%96%Negative27%8%2%6%79%69%90%98%P<0.001<0.01<0.01<0.001n.s.<0.05<0.05<0.001 B08*01Positive30%35%64%64%89%87%92%96%Negative14%28%62%61%88%82%90%93%P<0.001<0.001<0.01<0.001<0.01<0.001<0.01<0.001

Disclosures
No relevant conflicts of interest to declare.