Ovarian cancer is detectable from peripheral blood using machine learning over T-cell receptor repertoires.

The extraordinary diversity of T cells and B cells is critical for body maintenance. This diversity has an important role in protecting against tumor formation. In humans, the T-cell receptor (TCR) repertoire is generated through a striking stochastic process called V(D)J recombination, in which different gene segments are assembled and modified, leading to extensive variety. In ovarian cancer (OC), an unfortunate 80% of cases are detected late, leading to poor survival outcomes. However, when detected early, approximately 94% of patients live longer than 5 years after diagnosis. Thus, early detection is critical for patient survival. To determine whether the TCR repertoire obtained from peripheral blood is associated with tumor status, we collected blood samples from 85 women with or without OC and obtained TCR information. We then used machine learning to learn the characteristics of samples and to finally predict, over a set of unseen samples, whether the person is with or without OC. We successfully stratified the two groups, thereby associating the peripheral blood TCR repertoire with the formation of OC tumors. A careful study of the origin of the set of T cells most informative for the signature indicated the involvement of a specific invariant natural killer T (iNKT) clone and a specific mucosal-associated invariant T (MAIT) clone. Our findings here support the proposition that tumor-relevant signal is maintained by the immune system and is coded in the T-cell repertoire available in peripheral blood. It is also possible that the immune system detects tumors early enough for repertoire technologies to inform us near the beginning of tumor formation. Although such detection is made by the immune system, we might be able to identify it, using repertoire data from peripheral blood, to offer a pragmatic way to search for early signs of cancer with minimal patient burden, possibly with enhanced sensitivity.

The endogenous Mtv8 locus and the immunoglobulin repertoire.

The vast diversity of mammalian adaptive antigen receptors allows for robust and efficient immune responses against a wide number of pathogens. The antigen receptor repertoire is built during the recombination of B and T cell receptor (BCR, TCR) loci and hypermutation of BCR loci. V(D)J recombination rearranges these antigen receptor loci, which are organized as an array of separate V, (D), and J gene segments. Transcription activation at the recombining locus leads to changes in the local three-dimensional architecture, which subsequently contributes to which gene segments are utilized for recombination. The endogenous retrovirus (ERV) mouse mammary tumor provirus 8 (Mtv8) resides on mouse chromosome 6 interposed within the large array of light chain kappa V gene segments. As ERVs contribute to changes in genomic architecture by driving high levels of transcription of neighboring genes, it was suggested that Mtv8 could influence the BCR repertoire. We generated Mtv8-deficient mice to determine if the ERV influences V(D)J recombination to test this possibility. We find that Mtv8 does not influence the BCR repertoire.

MicroRNA, miR-501 regulate the V(D)J recombination in B cells.

The stringent regulation of RAGs (Recombination activating genes), the site-specific endonuclease responsible for V(D)J recombination, is important to prevent genomic rearrangements and chromosomal translocations in lymphoid cells. In the present study, we identify a microRNA, miR-501, which can regulate the expression of RAG1 in lymphoid cells. Overexpression of the pre-miRNA construct led to the generation of mature miRNAs and a concomitant reduction in RAG1 expression, whereas inhibition using anti-miRs resulted in its enhanced expression. The direct interaction of the 3'UTR of miR-501 with RAG1 was confirmed by the reporter assay. Importantly, overexpression of miRNAs led to inhibition of V(D)J recombination in B cells, revealing their impact on the physiological function of RAGs. Of interest is the inverse correlation observed for miR-501 with RAG1 in various leukemia patients and lymphoid cell lines, suggesting its possible use in cancer therapy. Thus, our results reveal the regulation of RAG1 by miR-501-3p in B cells and thus V(D)J recombination and its possible implications on immunoglobulin leukemogenesis.

Novel bimodal TRBD1-TRBD2 rearrangements with dual or absent D-region contribute to TRB V-(D)-J combinatorial diversity.

T-cell receptor (TR) diversity of the variable domains is generated by recombination of both the alpha (TRA) and beta (TRB) chains. The textbook process of TRB chain production starts with TRBD and TRBJ gene rearrangement, followed by the rearrangement of a TRBV gene to the partially rearranged D-J gene. Unsuccessful V-D-J TRB rearrangements lead to apoptosis of the cell. Here, we performed deep sequencing of the poorly explored pool of partial TRBD1-TRBD2 rearrangements in T-cell genomic DNA. We reconstructed full repertoires of human partial TRBD1-TRBD2 rearrangements using novel sequencing and validated them by detecting V-D-J recombination-specific byproducts: excision circles containing the recombination signal (RS) joint 5'D2-RS - 3'D1-RS. Identified rearrangements were in compliance with the classical 12/23 rule, common for humans, rats, and mice and contained typical V-D-J recombination footprints. Interestingly, we detected a bimodal distribution of D-D junctions indicating two active recombination sites producing long and short D-D rearrangements. Long TRB D-D rearrangements with two D-regions are coding joints D1-D2 remaining classically on the chromosome. The short TRB D-D rearrangements with no D-region are signal joints, the coding joint D1-D2 being excised from the chromosome. They both contribute to the TRB V-(D)-J combinatorial diversity. Indeed, short D-D rearrangements may be followed by direct V-J2 recombination. Long D-D rearrangements may recombine further with J2 and V genes forming partial D1-D2-J2 and then complete V-D1-D2-J2 rearrangement. Productive TRB V-D1-D2-J2 chains are present and expressed in thousands of clones of human antigen-experienced memory T cells proving their capacity for antigen recognition and actual participation in the immune response.

The immunogenicity of human-origin therapeutic antibodies are associated with V gene usage.

Therapeutic antibodies can elicit unwanted immune responses in a subset of patients, which leads to the production of anti-drug antibodies (ADA). Some of these ADAs have been reported to effect the pharmacokinetics, efficacy and/or safety of the therapeutic antibodies. The sequence diversity of antibodies are generated by VDJ recombination and mutagenesis. While the antibody generation process can create a large candidate pool for identifying high-affinity antibodies, it also could produce sequences that are foreign to the human immune system. However, it is not clear how VDJ recombination and mutagenesis impact the clinical ADA rate of therapeutic antibodies. In this study, we identified a positive correlation between the clinical ADA rate and the number of introduced mutations in the antibody sequences. We also found that the use of rare V alleles in human-origin antibody therapeutics is associated with higher risk of immunogenicity. The results suggest that antibody engineering projects should start with frameworks that contain commonly used V alleles and prioritize antibody candidates with low number of mutations to reduce the risk of immunogenicity.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes.

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Lymphocyte pathway analysis using naturally lymphocyte-deficient fish.

Comparative phylogenetic analyses are of potential value to establish the essential components of genetic networks underlying physiological traits. For species that naturally lack particular lymphocyte lineages, we show here that this strategy readily distinguishes trait-specific actors from pleiotropic components of the genetic network governing lymphocyte differentiation. Previously, three of the four members of the DNA polymerase X family have been implicated in the junctional diversification process during the somatic assembly of antigen receptors. Our phylogenetic analysis indicates that the presence of terminal deoxynucleotidyl transferase is strictly associated with the facility of V(D)J recombination, whereas PolL and PolM genes are retained even in species lacking Rag-mediated somatic diversification of antigen receptor genes.

Three-way contact analysis characterizes the higher order organization of the Tcra locus.

The generation of highly diverse antigen receptors in T and B lymphocytes relies on V(D)J recombination. The enhancer Eα has been implicated in regulating the accessibility of Vα and Jα genes through long-range interactions during rearrangements of the T-cell antigen receptor gene Tcra. However, direct evidence for Eα physically mediating the interaction of Vα and Jα genes is still lacking. In this study, we utilized the 3C-HTGTS assay, a chromatin interaction technique based on 3C, to analyze the higher order chromatin structure of the Tcra locus. Our analysis revealed the presence of sufficient information in the 3C-HTGTS data to detect multiway contacts. Three-way contact analysis of the Tcra locus demonstrated the co-occurrence of the proximal Jα genes, Vα genes and Eα in CD4+CD8+ double-positive thymocytes. Notably, the INT2-TEAp loop emerged as a prominent structure likely to be responsible for bringing the proximal Jα genes and the Vα genes into proximity. Moreover, the enhancer Eα utilizes this loop to establish physical proximity with the proximal Vα gene region. This study provides insights into the higher order chromatin structure of the Tcra locus, shedding light on the spatial organization of chromatin and its impact on V(D)J recombination.

A novel PRKDC mutation caused B lymphocytes V(D)J rearrangement disorder in the SLE-DAH like symptoms patient.

BACKGROUND: Analyzed the clinical features and treatment process of the patient suffering from immunodeficiency with systemic lupus erythematosus(SLE)-like syndrome in a novel mutation of PRKDC. CASE PRESENTATION: The patient had multiple positive auto-antibodies, chest CT and bronchoscopy showed Diffuse alveolar hemorrhage(DAH), and psychiatric symptoms showed brain atrophy by magnetic resonance imaging (MRI). Whole exon sequencing showed that novel complex heterozygous mutations of PRKDC gene (C. 1777 - 710_1777-709INSA (IVS16/IC16), C.1337T > A(p.Phe446Tyr). The mature B cell (CD19 + CD27 + CD38 dimIgD IgM-) were absent. The treatment of high-dose methylprednisolone (MP) and cyclophosphamide(CTX) can quickly relieve the symptoms of the patient. CONCLUSION: We described the case of an infant immunodeficiency with SLE like-syndrome, which may cause by PRKDC mutation, treated successfully with high-dose MP and CTX.

E protein binding at the Tcra enhancer promotes Tcra repertoire diversity.

V(D)J recombination of antigen receptor loci is a highly developmentally regulated process. During T lymphocyte development, recombination of the Tcra gene occurs in CD4+CD8+ double positive (DP) thymocytes and requires the Tcra enhancer (Eα). E proteins are known regulators of DP thymocyte development and have three identified binding sites in Eα. To understand the contribution of E proteins to Eα function, mutants lacking one or two of the respective binding sites were generated. The double-binding site mutant displayed a partial block at the positive selection stage of αß T cell development. Further investigation revealed loss of germline transcription within the Tcra locus at the Jα array, along with dysregulated primary and impaired secondary Vα-Jα rearrangement. Eα E protein binding increases Tcra locus accessibility and regulates TCRα recombination, thus directly promoting Tcra repertoire diversity.

Origin and evolutionary malleability of T cell receptor α diversity.

Lymphocytes of vertebrate adaptive immune systems acquired the capability to assemble, from split genes in the germline, billions of functional antigen receptors1-3. These receptors show specificity; unlike the broadly tuned receptors of the innate system, antibodies (Ig) expressed by B cells, for instance, can accurately distinguish between the two enantiomers of organic acids4, whereas T cell receptors (TCRs) reliably recognize single amino acid replacements in their peptide antigens5. In developing lymphocytes, antigen receptor genes are assembled from a comparatively small set of germline-encoded genetic elements in a process referred to as V(D)J recombination6,7. Potential self-reactivity of some antigen receptors arising from the quasi-random somatic diversification is suppressed by several robust control mechanisms8-12. For decades, scientists have puzzled over the evolutionary origin of somatically diversifying antigen receptors13-16. It has remained unclear how, at the inception of this mechanism, immunologically beneficial expanded receptor diversity was traded against the emerging risk of destructive self-recognition. Here we explore the hypothesis that in early vertebrates, sequence microhomologies marking the ends of recombining elements became the crucial targets of selection determining the outcome of non-homologous end joining-based repair of DNA double-strand breaks generated during RAG-mediated recombination. We find that, across the main clades of jawed vertebrates, TCRα repertoire diversity is best explained by species-specific extents of such sequence microhomologies. Thus, selection of germline sequence composition of rearranging elements emerges as a major factor determining the degree of diversity of somatically generated antigen receptors.

Contribution of the IGCR1 regulatory element and the 3'Igh CTCF-binding elements to regulation of Igh V(D)J recombination.

Immunoglobulin heavy chain variable region exons are assembled in progenitor-B cells, from VH, D, and JH gene segments located in separate clusters across the Igh locus. RAG endonuclease initiates V(D)J recombination from a JH-based recombination center (RC). Cohesin-mediated extrusion of upstream chromatin past RC-bound RAG presents Ds for joining to JHs to form a DJH-RC. Igh has a provocative number and organization of CTCF-binding elements (CBEs) that can impede loop extrusion. Thus, Igh has two divergently oriented CBEs (CBE1 and CBE2) in the IGCR1 element between the VH and D/JH domains, over 100 CBEs across the VH domain convergent to CBE1, and 10 clustered 3'Igh-CBEs convergent to CBE2 and VH CBEs. IGCR1 CBEs segregate D/JH and VH domains by impeding loop extrusion-mediated RAG-scanning. Downregulation of WAPL, a cohesin unloader, in progenitor-B cells neutralizes CBEs, allowing DJH-RC-bound RAG to scan the VH domain and perform VH-to-DJH rearrangements. To elucidate potential roles of IGCR1-based CBEs and 3'Igh-CBEs in regulating RAG-scanning and elucidate the mechanism of the ordered transition from D-to-JH to VH-to-DJH recombination, we tested effects of inverting and/or deleting IGCR1 or 3'Igh-CBEs in mice and/or progenitor-B cell lines. These studies revealed that normal IGCR1 CBE orientation augments RAG-scanning impediment activity and suggest that 3'Igh-CBEs reinforce ability of the RC to function as a dynamic loop extrusion impediment to promote optimal RAG scanning activity. Finally, our findings indicate that ordered V(D)J recombination can be explained by a gradual WAPL downregulation mechanism in progenitor-B cells as opposed to a strict developmental switch.

Statistical inference reveals the role of length, GC content, and local sequence in V(D)J nucleotide trimming.

To appropriately defend against a wide array of pathogens, humans somatically generate highly diverse repertoires of B cell and T cell receptors (BCRs and TCRs) through a random process called V(D)J recombination. Receptor diversity is achieved during this process through both the combinatorial assembly of V(D)J-genes and the junctional deletion and insertion of nucleotides. While the Artemis protein is often regarded as the main nuclease involved in V(D)J recombination, the exact mechanism of nucleotide trimming is not understood. Using a previously published TCRß repertoire sequencing data set, we have designed a flexible probabilistic model of nucleotide trimming that allows us to explore various mechanistically interpretable sequence-level features. We show that local sequence context, length, and GC nucleotide content in both directions of the wider sequence, together, can most accurately predict the trimming probabilities of a given V-gene sequence. Because GC nucleotide content is predictive of sequence-breathing, this model provides quantitative statistical evidence regarding the extent to which double-stranded DNA may need to be able to breathe for trimming to occur. We also see evidence of a sequence motif that appears to get preferentially trimmed, independent of GC-content-related effects. Further, we find that the inferred coefficients from this model provide accurate prediction for V- and J-gene sequences from other adaptive immune receptor loci. These results refine our understanding of how the Artemis nuclease may function to trim nucleotides during V(D)J recombination and provide another step toward understanding how V(D)J recombination generates diverse receptors and supports a powerful, unique immune response in healthy humans.

Igh and Igk loci use different folding principles for V gene recombination due to distinct chromosomal architectures of pro-B and pre-B cells.

Extended loop extrusion across the immunoglobulin heavy-chain (Igh) locus facilitates VH-DJH recombination following downregulation of the cohesin-release factor Wapl by Pax5, resulting in global changes in the chromosomal architecture of pro-B cells. Here, we demonstrate that chromatin looping and VK-JK recombination at the Igk locus were insensitive to Wapl upregulation in pre-B cells. Notably, the Wapl protein was expressed at a 2.2-fold higher level in pre-B cells compared with pro-B cells, which resulted in a distinct chromosomal architecture with normal loop sizes in pre-B cells. High-resolution chromosomal contact analysis of the Igk locus identified multiple internal loops, which likely juxtapose VK and JK elements to facilitate VK-JK recombination. The higher Wapl expression in Igµ-transgenic pre-B cells prevented extended loop extrusion at the Igh locus, leading to recombination of only the 6 most 3' proximal VH genes and likely to allelic exclusion of all other VH genes in pre-B cells. These results suggest that pro-B and pre-B cells with their distinct chromosomal architectures use different chromatin folding principles for V gene recombination, thereby enabling allelic exclusion at the Igh locus, when the Igk locus is recombined.

The DNA binding domain and the C-terminal region of DNA Ligase IV specify its role in V(D)J recombination.

DNA Ligase IV is responsible for the repair of DNA double-strand breaks (DSB), including DSBs that are generated during V(D)J recombination. Like other DNA ligases, Ligase IV contains a catalytic core with three subdomains-the DNA binding (DBD), the nucleotidyltransferase (NTD), and the oligonucleotide/oligosaccharide-fold subdomain (OBD). Ligase IV also has a unique C-terminal region that includes two BRCT domains, a nuclear localization signal sequence and a stretch of amino acid that participate in its interaction with XRCC4. Out of the three mammalian ligases, Ligase IV is the only ligase that participates in and is required for V(D)J recombination. Identification of the minimal domains within DNA Ligase IV that contribute to V(D)J recombination has remained unresolved. The interaction of the Ligase IV DNA binding domain with Artemis, and the interaction of its C-terminal region with XRCC4, suggest that both of these regions that also interact with the Ku70/80 heterodimer are important and might be sufficient for mediating participation of DNA Ligase IV in V(D)J recombination. This hypothesis was investigated by generating chimeric ligase proteins by swapping domains, and testing their ability to rescue V(D)J recombination in Ligase IV-deficient cells. We demonstrate that a fusion protein containing Ligase I NTD and OBDs flanked by DNA Ligase IV DBD and C-terminal region is sufficient to support V(D)J recombination. This chimeric protein, which we named Ligase 37, complemented formation of coding and signal joints. Coding joints generated with Ligase 37 were shorter than those observed with wild type DNA Ligase IV. The shorter length was due to increased nucleotide deletions and decreased nucleotide insertions. Additionally, overexpression of Ligase 37 in a mouse pro-B cell line supported a shift towards shorter coding joints. Our findings demonstrate that the ability of DNA Ligase IV to participate in V(D)J recombination is in large part mediated by its DBD and C-terminal region.

Locus architecture and RAG scanning determine antibody diversity.

Diverse mammalian antibody repertoires are produced via distant genomic contacts involving immunoglobulin Igh variable (V), diversity (D), and joining (J) gene segments and result in V(D)J recombination. How such interactions determine V gene usage remains unclear. The recombination-activating gene (RAG) chromatin scanning model posits that RAG recombinase bound to the recombination center (RC) linearly tracks along chromatin by means of cohesin-mediated loop extrusion; a proposition supported by cohesin depletion studies. A mechanistic role for chromatin loop extrusion has also been implicated for Igh locus contraction. In this opinion, we provide perspective on how loop extrusion interfaces with the 3D conformation of the Igh locus and newly identified enhancers that regionally regulate VH gene usage during V(D)J recombination, shaping the preselected repertoire.

AIRR-C IG Reference Sets: curated sets of immunoglobulin heavy and light chain germline genes.

Introduction: Analysis of an individual's immunoglobulin (IG) gene repertoire requires the use of high-quality germline gene reference sets. When sets only contain alleles supported by strong evidence, AIRR sequencing (AIRR-seq) data analysis is more accurate and studies of the evolution of IG genes, their allelic variants and the expressed immune repertoire is therefore facilitated. Methods: The Adaptive Immune Receptor Repertoire Community (AIRR-C) IG Reference Sets have been developed by including only human IG heavy and light chain alleles that have been confirmed by evidence from multiple high-quality sources. To further improve AIRR-seq analysis, some alleles have been extended to deal with short 3' or 5' truncations that can lead them to be overlooked by alignment utilities. To avoid other challenges for analysis programs, exact paralogs (e.g. IGHV1-69*01 and IGHV1-69D*01) are only represented once in each set, though alternative sequence names are noted in accompanying metadata. Results and discussion: The Reference Sets include less than half the previously recognised IG alleles (e.g. just 198 IGHV sequences), and also include a number of novel alleles: 8 IGHV alleles, 2 IGKV alleles and 5 IGLV alleles. Despite their smaller sizes, erroneous calls were eliminated, and excellent coverage was achieved when a set of repertoires comprising over 4 million V(D)J rearrangements from 99 individuals were analyzed using the Sets. The version-tracked AIRR-C IG Reference Sets are freely available at the OGRDB website (https://ogrdb.airr-community.org/germline_sets/Human) and will be regularly updated to include newly observed and previously reported sequences that can be confirmed by new high-quality data.

Single-Cell RNA Sequencing for the Detection of Clonotypic V(D)J Rearrangements in Multiple Myeloma.

Multiple myeloma (MM) has a highly heterogeneous genetic background, which complicates its molecular tracking over time. Nevertheless, each MM patient's malignant plasma cells (PCs) share unique V(D)J rearranged sequences at immunoglobulin loci, which represent ideal disease biomarkers. Because the tumor-specific V(D)J sequence is highly expressed in bulk RNA in MM patients, we wondered whether it can be identified by single-cell RNA sequencing (scRNA-seq). To this end we analyzed CD138+ cells purified from bone marrow aspirates of 19 samples with PC dyscrasias by both a standard method based on bulk DNA and by an implementation of the standard 10x Genomics protocol to detect expressed V(D)J sequences. A dominant clonotype was easily identified in each sample, accounting on average for 83.65% of V(D)J-rearranged cells. Compared with standard methods, scRNA-seq analysis proved highly concordant and even more effective in identifying clonal productive rearrangements, by-passing limitations related to the misannealing of consensus primers in hypermutated regions. We next validated its accuracy to track 5 clonal cells with absolute sensitivity in a virtual sample containing 3180 polyclonal cells. This shows that single-cell V(D)J analysis may be used to find rare clonal cells, laying the foundations for functional single-cell dissection of minimal residual disease.

An updated definition of V(D)J recombination signal sequences revealed by high-throughput recombination assays.

In the adaptive immune system, V(D)J recombination initiates the production of a diverse antigen receptor repertoire in developing B and T cells. Recombination activating proteins, RAG1 and RAG2 (RAG1/2), catalyze V(D)J recombination by cleaving adjacent to recombination signal sequences (RSSs) that flank antigen receptor gene segments. Previous studies defined the consensus RSS as containing conserved heptamer and nonamer sequences separated by a less conserved 12 or 23 base-pair spacer sequence. However, many RSSs deviate from the consensus sequence. Here, we developed a cell-based, massively parallel assay to evaluate V(D)J recombination activity on thousands of RSSs where the 12-RSS heptamer and adjoining spacer region contained randomized sequences. While the consensus heptamer sequence (CACAGTG) was marginally preferred, V(D)J recombination was highly active on a wide range of non-consensus sequences. Select purine/pyrimidine motifs that may accommodate heptamer unwinding in the RAG1/2 active site were generally preferred. In addition, while different coding flanks and nonamer sequences affected recombination efficiency, the relative dependency on the purine/pyrimidine motifs in the RSS heptamer remained unchanged. Our results suggest RAG1/2 specificity for RSS heptamers is primarily dictated by DNA structural features dependent on purine/pyrimidine pattern, and to a lesser extent, RAG:RSS base-specific interactions.

Functional antibodies exhibit light chain coherence.

The vertebrate adaptive immune system modifies the genome of individual B cells to encode antibodies that bind particular antigens1. In most mammals, antibodies are composed of heavy and light chains that are generated sequentially by recombination of V, D (for heavy chains), J and C gene segments. Each chain contains three complementarity-determining regions (CDR1-CDR3), which contribute to antigen specificity. Certain heavy and light chains are preferred for particular antigens2-22. Here we consider pairs of B cells that share the same heavy chain V gene and CDRH3 amino acid sequence and were isolated from different donors, also known as public clonotypes23,24. We show that for naive antibodies (those not yet adapted to antigens), the probability that they use the same light chain V gene is around 10%, whereas for memory (functional) antibodies, it is around 80%, even if only one cell per clonotype is used. This property of functional antibodies is a phenomenon that we call light chain coherence. We also observe this phenomenon when similar heavy chains recur within a donor. Thus, although naive antibodies seem to recur by chance, the recurrence of functional antibodies reveals surprising constraint and determinism in the processes of V(D)J recombination and immune selection. For most functional antibodies, the heavy chain determines the light chain.