Potent Neutralizing Antibodies against SARS-CoV-2 Identified by High-Throughput Single-Cell Sequencing of Convalescent Patients' B Cells.

The COVID-19 pandemic urgently needs therapeutic and prophylactic interventions. Here, we report the rapid identification of SARS-CoV-2-neutralizing antibodies by high-throughput single-cell RNA and VDJ sequencing of antigen-enriched B cells from 60 convalescent patients. From 8,558 antigen-binding IgG1+ clonotypes, 14 potent neutralizing antibodies were identified, with the most potent one, BD-368-2, exhibiting an IC50 of 1.2 and 15 ng/mL against pseudotyped and authentic SARS-CoV-2, respectively. BD-368-2 also displayed strong therapeutic and prophylactic efficacy in SARS-CoV-2-infected hACE2-transgenic mice. Additionally, the 3.8 Å cryo-EM structure of a neutralizing antibody in complex with the spike-ectodomain trimer revealed the antibody's epitope overlaps with the ACE2 binding site. Moreover, we demonstrated that SARS-CoV-2-neutralizing antibodies could be directly selected based on similarities of their predicted CDR3H structures to those of SARS-CoV-neutralizing antibodies. Altogether, we showed that human neutralizing antibodies could be efficiently discovered by high-throughput single B cell sequencing in response to pandemic infectious diseases.

Single-Cell-Based High-Throughput Ig and TCR Repertoire Sequencing Analysis in Rhesus Macaques.

Recent advancements in microfluidics and high-throughput sequencing technologies have enabled recovery of paired H and L chains of Igs and VDJ and VJ chains of TCRs from thousands of single cells simultaneously in humans and mice. Despite rhesus macaques being one of the most well-studied model organisms for the human adaptive immune response, high-throughput single-cell immune repertoire sequencing assays are not yet available due to the complexity of these polyclonal receptors. We used custom primers that capture all known rhesus macaque Ig and TCR isotypes and chains that are fully compatible with a commercial solution for single-cell immune repertoire profiling. Using these rhesus-specific assays, we sequenced Ig and TCR repertoires in >60,000 cells from cryopreserved rhesus PBMCs, splenocytes, and FACS-sorted B and T cells. We were able to recover every Ig isotype and TCR chain, measure clonal expansion in proliferating T cells, and pair Ig and TCR repertoires with gene expression profiles of the same single cells. Our results establish the ability to perform high-throughput immune repertoire analysis in rhesus macaques at the single-cell level.

Global characterization of B cell receptor repertoire in COVID-19 patients by single-cell V(D)J sequencing.

The world is facing a pandemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Adaptive immune responses are essential for SARS-CoV-2 virus clearance. Although a large body of studies have been conducted to investigate the immune mechanism in COVID-19 patients, we still lack a comprehensive understanding of the BCR repertoire in patients. In this study, we used the single-cell V(D)J sequencing to characterize the BCR repertoire across convalescent COVID-19 patients. We observed that the BCR diversity was significantly reduced in disease compared with healthy controls. And BCRs tend to skew toward different V gene segments in COVID-19 and healthy controls. The CDR3 sequences of heavy chain in clonal BCRs in patients were more convergent than that in healthy controls. In addition, we discovered increased IgG and IgA isotypes in the disease, including IgG1, IgG3 and IgA1. In all clonal BCRs, IgG isotypes had the most frequent class switch recombination events and the highest somatic hypermutation rate, especially IgG3. Moreover, we found that an IgG3 cluster from different clonal groups had the same IGHV, IGHJ and CDR3 sequences (IGHV4-4-CARLANTNQFYDSSSYLNAMDVW-IGHJ6). Overall, our study provides a comprehensive characterization of the BCR repertoire in COVID-19 patients, which contributes to the understanding of the mechanism for the immune response to SARS-CoV-2 infection.

The 3D structure of the immunoglobulin heavy-chain locus: implications for long-range genomic interactions.

The immunoglobulin heavy-chain (Igh) locus is organized into distinct regions that contain multiple variable (V(H)), diversity (D(H)), joining (J(H)) and constant (C(H)) coding elements. How the Igh locus is structured in 3D space is unknown. To probe the topography of the Igh locus, spatial distance distributions were determined between 12 genomic markers that span the entire Igh locus. Comparison of the distance distributions to computer simulations of alternative chromatin arrangements predicted that the Igh locus is organized into compartments containing clusters of loops separated by linkers. Trilateration and triple-point angle measurements indicated the mean relative 3D positions of the V(H), D(H), J(H), and C(H) elements, showed compartmentalization and striking conformational changes involving V(H) and D(H)-J(H) elements during early B cell development. In pro-B cells, the entire repertoire of V(H) regions (2 Mbp) appeared to have merged and juxtaposed to the D(H) elements, mechanistically permitting long-range genomic interactions to occur with relatively high frequency.

Public antibodies to malaria antigens generated by two LAIR1 insertion modalities.

In two previously described donors, the extracellular domain of LAIR1, a collagen-binding inhibitory receptor encoded on chromosome 19 (ref. 1), was inserted between the V and DJ segments of an antibody. This insertion generated, through somatic mutations, broadly reactive antibodies against RIFINs, a type of variant antigen expressed on the surface of Plasmodium falciparum-infected erythrocytes. To investigate how frequently such antibodies are produced in response to malaria infection, we screened plasma from two large cohorts of individuals living in malaria-endemic regions. Here we report that 5-10% of malaria-exposed individuals, but none of the European blood donors tested, have high levels of LAIR1-containing antibodies that dominate the response to infected erythrocytes without conferring enhanced protection against febrile malaria. By analysing the antibody-producing B cell clones at the protein, cDNA and gDNA levels, we characterized additional LAIR1 insertions between the V and DJ segments and discovered a second insertion modality whereby the LAIR1 exon encoding the extracellular domain and flanking intronic sequences are inserted into the switch region. By exon shuffling, this mechanism leads to the production of bispecific antibodies in which the LAIR1 domain is precisely positioned at the elbow between the VH and CH1 domains. Additionally, in one donor the genomic DNA encoding the VH and CH1 domains was deleted, leading to the production of a camel-like LAIR1-containing antibody. Sequencing of the switch regions of memory B cells from European blood donors revealed frequent templated inserts originating from transcribed genes that, in rare cases, comprised exons with orientations and frames compatible with expression. These results reveal different modalities of LAIR1 insertion that lead to public and dominant antibodies against infected erythrocytes and suggest that insertion of templated DNA represents an additional mechanism of antibody diversification that can be selected in the immune response against pathogens and exploited for B cell engineering.

Alignment free identification of clones in B cell receptor repertoires.

Following antigenic challenge, activated B cells rapidly expand and undergo somatic hypermutation, yielding groups of clonally related B cells with diversified immunoglobulin receptors. Inference of clonal relationships based on the receptor sequence is an essential step in many adaptive immune receptor repertoire sequencing studies. These relationships are typically identified by a multi-step process that involves: (i) grouping sequences based on shared V and J gene assignments, and junction lengths and (ii) clustering these sequences using a junction-based distance. However, this approach is sensitive to the initial gene assignments, which are error-prone, and fails to identify clonal relatives whose junction length has changed through accumulation of indels. Through defining a translation-invariant feature space in which we cluster the sequences, we develop an alignment free clonal identification method that does not require gene assignments and is not restricted to a fixed junction length. This alignment free approach has higher sensitivity compared to a typical junction-based distance method without loss of specificity and PPV. While the alignment free procedure identifies clones that are broadly consistent with the junction-based distance method, it also identifies clones with characteristics (multiple V or J gene assignments or junction lengths) that are not detectable with the junction-based distance method.

Orientation-specific joining of AID-initiated DNA breaks promotes antibody class switching.

During B-cell development, RAG endonuclease cleaves immunoglobulin heavy chain (IgH) V, D, and J gene segments and orchestrates their fusion as deletional events that assemble a V(D)J exon in the same transcriptional orientation as adjacent Cµ constant region exons. In mice, six additional sets of constant region exons (CHs) lie 100-200 kilobases downstream in the same transcriptional orientation as V(D)J and Cµ exons. Long repetitive switch (S) regions precede Cµ and downstream CHs. In mature B cells, class switch recombination (CSR) generates different antibody classes by replacing Cµ with a downstream CH (ref. 2). Activation-induced cytidine deaminase (AID) initiates CSR by promoting deamination lesions within Sµ and a downstream acceptor S region; these lesions are converted into DNA double-strand breaks (DSBs) by general DNA repair factors. Productive CSR must occur in a deletional orientation by joining the upstream end of an Sµ DSB to the downstream end of an acceptor S-region DSB. However, the relative frequency of deletional to inversional CSR junctions has not been measured. Thus, whether orientation-specific joining is a programmed mechanistic feature of CSR as it is for V(D)J recombination and, if so, how this is achieved is unknown. To address this question, we adapt high-throughput genome-wide translocation sequencing into a highly sensitive DSB end-joining assay and apply it to endogenous AID-initiated S-region DSBs in mouse B cells. We show that CSR is programmed to occur in a productive deletional orientation and does so via an unprecedented mechanism that involves in cis Igh organizational features in combination with frequent S-region DSBs initiated by AID. We further implicate ATM-dependent DSB-response factors in enforcing this mechanism and provide an explanation of why CSR is so reliant on the 53BP1 DSB-response factor.

Flexible ordering of antibody class switch and V(D)J joining during B-cell ontogeny.

V(D)J joining is mediated by RAG recombinase during early B-lymphocyte development in the bone marrow (BM). Activation-induced deaminase initiates isotype switching in mature B cells of secondary lymphoid structures. Previous studies questioned the strict ontological partitioning of these processes. We show that pro-B cells undergo robust switching to a subset of immunoglobulin H (IgH) isotypes. Chromatin studies reveal that in pro-B cells, the spatial organization of the Igh locus may restrict switching to this subset of isotypes. We demonstrate that in the BM, V(D)J joining and switching are interchangeably inducible, providing an explanation for the hyper-IgE phenotype of Omenn syndrome.

The T Cell Repertoires from Nickel Sensitized Joint Implant Failure Patients.

Nickel (Ni2+) is one of the most common allergens, affecting around 10-15% of the general population. As the demand for orthopedic implant surgery rises, the number of surgical revisions due to joint implant failure also increases. There is evidence that some patients develop joint failure due to an immune response to a component of the implant, and we have found that Ni2+ is an especially important cause. Hence, understanding the mechanisms by which Ni2+ allergy induces joint implant failure becomes a critical research question. The structural basis of Ni2+ activation of pathogenic T cells is still not clear. The purpose of this study was to characterize Ni2+-reactive T cell repertoires derived from the peripheral blood of joint failure patients due to Ni2+ sensitization using single-cell sequencing techniques. We stimulated the proliferation of Ni2+ -reactive T cells from two implant failure patients in vitro, and sorted them for single-cell VDJ sequencing (10× genomics). We identified 2650 productive V-J spanning pairs. Both TCR α chains and ß chains were enriched. TRBV18 usage is the highest in the P7 CD4+ population (18.1%), and TRBV5-1 usage is the highest in the P7 CD8+ population (12.1%). TRBV19 and TRBV20-1 segments are present in a high percentage of both P7 and P9 sequenced T cells. Remarkably, the alpha and beta chain combination of TRAV41-TRBV18 accounts for 13.5% of the CD4+ population of P7 patient. Compared to current Ni specific T cell repertoire studies of contact dermatitis, the Vα and Vß usages of these joint implant failure patients were different. This could be due to the different availability of self-peptides in these two different tissues. However, TRBV19 (Vß17) was among frequently used TCR ß chains, which are common in previous reports. This implies that some pathogenic T cells could be similar in Ni2+ hypersensitivities in skin and joints. The alignment of the TCR CDR3ß sequences showed a conserved glutamic acid (Glu) that could potentially interact with Ni2+. The study of these Ni2+ specific TCRs may shed light on the molecular mechanism of T cell activation by low molecular weight chemical haptens.

Biased IGH VDJ gene repertoire and clonal expansions in B cells of chronically hepatitis C virus-infected individuals.

Patients chronically infected with hepatitis C virus (HCV) frequently develop mixed cryoglobulinemia (MC), a monoclonal expansion of immunoglobulin M (IgM)+ autoreactive B cells, and also have an increased B-cell lymphoma risk. Whether HCV infection also impacts the B-cell compartment and the B-cell receptor repertoire in patients not affected by MC or lymphomas is poorly understood. Flow cytometric analysis of peripheral blood B cells of 30 MC-negative HCV-infected patients and 15 healthy controls revealed that frequencies of class-switched memory B cells were increased in the patients, whereas frequencies of transitional and naive B cells were decreased. For 22 HCV+ patients and 7 healthy controls, we performed high-throughput sequencing of immunoglobulin heavy chain VDJ rearrangements of naive, mature CD5+, IgM+ memory, and class-switched memory B cells. An increased usage of several IGHV genes, including IGHV1-69 and IGHV4-59, which are closely linked to MC and HCV-associated lymphomas, was specifically seen among IgM+ memory B cells of the patients. Moreover, many, and partly very large, expanded clones were seen predominantly among IgM+ memory B cells of all HCV-infected patients analyzed. Thus, chronic HCV infection massively disturbs the B-cell compartment even in patients without clinically detectable B-cell lymphoproliferation and generates many large B-cell clones, especially among non-class-switched memory B cells. Because B-cell clones in MC and lymphomas derive from this B-cell subset, this establishes IgM+ memory B cells as a general target of lymphoproliferation in HCV+ patients, affecting apparently all patients.

Concurrent immunoglobulin G-lambda type multiple myeloma and mixed cellularity classical Hodgkin lymphoma: A case report.

A 66-year-old man with a swollen right inguinal lymph node (LN) had pain on the lower side of the back. Computed tomography revealed bone disease in the back and swollen right inguinal LNs. Laboratory studies showed anemia and serum immunoglobulin G-lambda (IgG-λ) type monoclonal protein. The bone marrow contained 39.6% plasma cells. He was diagnosed with IgG-λ type multiple myeloma (MM). However, the pathological findings of the right inguinal LN were mixed cellular classical Hodgkin lymphoma (HL). The administration of melphalan, prednisone, and bortezomib (MPB) was started for MM; however, swelling in the right inguinal LN increased. After three cycles of MPB, the administration of doxorubicin, bleomycin, vinblastine, and dacarbazine (ABVD) was started for HL. However, HL was refractory to ABVD. Pancytopenia subsequently progressed and rapid swelling occurred in his LNs. He died 7 months after diagnosis. Multiple myeloma was diagnosed, based on the typical symptoms, although the pathological findings of the LN indicated a diagnosis of HL. We analyzed the molecular relationship between MM and HL cells using a direct sequencing method. The sequencing results demonstrated that the variable-diversity-joining (VDJ) region of the IgH gene was identified with 94.4% of IGLV3-32*01 in the bone marrow sample at diagnosis. Furthermore, clonotypic IgH sequence was identified in CD30-positive cells from the LN. These results suggested that the clonal HL cells were derived from the same source as the clonal MM cells and demonstrated that MM and HL in this patient may have originated from the same B cell progenitor.

Insights into immune system development and function from mouse T-cell repertoires.

The ability of the adaptive immune system to respond to arbitrary pathogens stems from the broad diversity of immune cell surface receptors. This diversity originates in a stochastic DNA editing process (VDJ recombination) that acts on the surface receptor gene each time a new immune cell is created from a stem cell. By analyzing T-cell receptor (TCR) sequence repertoires taken from the blood and thymus of mice of different ages, we quantify the changes in the VDJ recombination process that occur from embryo to young adult. We find a rapid increase with age in the number of random insertions and a dramatic increase in diversity. Because the blood accumulates thymic output over time, blood repertoires are mixtures of different statistical recombination processes, and we unravel the mixture statistics to obtain a picture of the time evolution of the early immune system. Sequence repertoire analysis also allows us to detect the statistical impact of selection on the output of the VDJ recombination process. The effects we find are nearly identical between thymus and blood, suggesting that our analysis mainly detects selection for proper folding of the TCR receptor protein. We further find that selection is weaker in laboratory mice than in humans and it does not affect the diversity of the repertoire.

Stability and uniqueness of clonal immunoglobulin CDR3 sequences for MRD tracking in multiple myeloma.

Minimal residual disease (MRD) tracking, by next generation sequencing of immunoglobulin sequences, is moving towards clinical implementation in multiple myeloma. However, there is only sparse information available to address whether clonal sequences remain stable for tracking over time, and to what extent light chain sequences are sufficiently unique for tracking. Here, we analyzed immunoglobulin repertoires from 905 plasma cell myeloma and healthy control samples, focusing on the third complementarity determining region (CDR3). Clonal heavy and/or light chain expression was identified in all patients at baseline, with one or more subclones related to the main clone in 3.2%. In 45 patients with 101 sequential samples, the dominant clonal CDR3 sequences remained identical over time, despite differential clonal evolution by whole exome sequencing in 49% of patients. The low frequency of subclonal CDR3 variants, and absence of evolution over time in active multiple myeloma, indicates that tumor cells at this stage are not under selective pressure to undergo antibody affinity maturation. Next, we establish somatic hypermutation and non-templated insertions as the most important determinants of light chain clonal uniqueness, identifying a potentially trackable sequence in the majority of patients. Taken together, we show that dominant clonal sequences identified at baseline are reliable biomarkers for long-term tracking of the malignant clone, including both IGH and the majority of light chain clones.

Interplay between Target Sequences and Repair Pathways Determines Distinct Outcomes of AID-Initiated Lesions.

Activation-induced deaminase (AID) functions by deaminating cytosines and causing U:G mismatches, a rate-limiting step of Ab gene diversification. However, precise mechanisms regulating AID deamination frequency remain incompletely understood. Moreover, it is not known whether different sequence contexts influence the preferential access of mismatch repair or uracil glycosylase (UNG) to AID-initiated U:G mismatches. In this study, we employed two knock-in models to directly compare the mutability of core Sµ and VDJ exon sequences and their ability to regulate AID deamination and subsequent repair process. We find that the switch (S) region is a much more efficient AID deamination target than the V region. Igh locus AID-initiated lesions are processed by error-free and error-prone repair. S region U:G mismatches are preferentially accessed by UNG, leading to more UNG-dependent deletions, enhanced by mismatch repair deficiency. V region mutation hotspots are largely determined by AID deamination. Recurrent and conserved S region motifs potentially function as spacers between AID deamination hotspots. We conclude that the pattern of mutation hotspots and DNA break generation is influenced by sequence-intrinsic properties, which regulate AID deamination and affect the preferential access of downstream repair. Our studies reveal an evolutionarily conserved role for substrate sequences in regulating Ab gene diversity and AID targeting specificity.

The ATM Kinase Restrains Joining of Both VDJ Signal and Coding Ends.

The evidence that ATM affects resolution of RAG-induced DNA double-strand breaks is profuse and unequivocal; moreover, it is clear that the RAG complex itself cooperates (in an undetermined way) with ATM to facilitate repair of these double-strand breaks by the classical nonhomologous end-joining pathway. The mechanistic basis for the cooperation between ATM and the RAG complex has not been defined, although proposed models invoke ATM and RAG2's C terminus in maintaining the RAG postcleavage complex. In this study, we show that ATM reduces the rate of both coding and signal joining in a robust episomal assay; we suggest that this is the result of increased stability of the postcleavage complex. ATM's ability to inhibit VDJ joining requires its enzymatic activity. The noncore C termini of both RAG1 and RAG2 are also required for ATM's capacity to limit signal (but not coding) joining. Moreover, potential phosphorylation targets within the C terminus of RAG2 are also required for ATM's capacity to limit signal joining. These data suggest a model whereby the RAG signal end complex is stabilized by phosphorylation of RAG2 by ATM.

CTCF-binding elements 1 and 2 in the Igh intergenic control region cooperatively regulate V(D)J recombination.

Ig heavy chain (IgH) variable region exons are assembled from V, D, and J gene segments during early B-lymphocyte differentiation. A several megabase region at the "distal" end of the mouse IgH locus (Igh) contains hundreds of V(H)s, separated by an intergenic region from Igh Ds, J(H)s, and constant region exons. Diverse primary Igh repertoires are generated by joining Vs, Ds, and Js in different combinations, with a given B cell productively assembling only one combination. The intergenic control region 1 (IGCR1) in the V(H)-to-D intergenic region regulates Igh V(D)J recombination in the contexts of developmental order, lineage specificity, and feedback from productive rearrangements. IGCR1 also diversifies IgH repertoires by balancing proximal and distal V(H) use. IGCR1 functions in all these regulatory contexts by suppressing predominant rearrangement of D-proximal V(H)s. Such IGCR1 functions were neutralized by simultaneous mutation of two CCCTC-binding factor (CTCF)-binding elements (CBE1 and CBE2) within it. However, it was unknown whether only one CBE mediates IGCR1 functions or whether both function in this context. To address these questions, we generated mice in which either IGCR1 CBE1 or CBE2 was replaced with scrambled sequences that do not bind CTCF. We found that inactivation of CBE1 or CBE2 individually led to only partial impairment of various IGCR1 functions relative to the far greater effects of inactivating both binding elements simultaneously, demonstrating that they function cooperatively to achieve full IGCR1 regulatory activity. Based on these and other findings, we propose an orientation-specific looping model for synergistic CBE1 and CBE2 functions.

Deep sequencing of the TCR-ß repertoire of human forkhead box protein 3 (FoxP3)+ and FoxP3- T cells suggests that they are completely distinct and non-overlapping.

Maintenance of peripheral tolerance requires a balance between autoreactive conventional T cells (Tconv ) and thymically derived forkhead box protein 3 (FoxP3)+ regulatory T cells (tTregs ). Considerable controversy exists regarding the similarities/differences in T cell receptor (TCR) repertoires expressed by Tconv and tTregs . We generated highly purified populations of human adult and cord blood Tconv and tTregs based on the differential expression of CD25 and CD127. The purity of the sorted populations was validated by intracellular staining for FoxP3 and Helios. We also purified an overlap group of CD4 T cells from adult donors to ensure that considerable numbers of shared clonotypes could be detected when present. We used deep sequencing of entire TCR-ß CDR3 sequences to analyse the TCR repertoire of Tconv and tTregs . Our studies suggest that both neonatal and adult human Tconv and tTreg cells are, in fact, entirely distinct CD4 T cell lineages.

IgH-V(D)J NGS-MRD measurement pre- and early post-allotransplant defines very low- and very high-risk ALL patients.

Positive detection of minimal residual disease (MRD) by multichannel flow cytometry (MFC) prior to hematopoietic cell transplantation (HCT) of patients with acute lymphoblastic leukemia (ALL) identifies patients at high risk for relapse, but many pre-HCT MFC-MRD negative patients also relapse, and the predictive power MFC-MRD early post-HCT is poor. To test whether the increased sensitivity of next-generation sequencing (NGS)-MRD better identifies pre- and post-HCT relapse risk, we performed immunoglobulin heavy chain (IgH) variable, diversity, and joining (V[D]J) DNA sequences J NGS-MRD on 56 patients with B-cell ALL enrolled in Children's Oncology Group trial ASCT0431. NGS-MRD predicted relapse and survival more accurately than MFC-MRD (P < .0001), especially in the MRD negative cohort (relapse, 0% vs 16%; P = .02; 2-year overall survival, 96% vs 77%; P = .003). Post-HCT NGS-MRD detection was better at predicting relapse than MFC-MRD (P < .0001), especially early after HCT (day 30 MFC-MRD positive relapse rate, 35%; NGS-MRD positive relapse rate, 67%; P = .004). Any post-HCT NGS positivity resulted in an increase in relapse risk by multivariate analysis (hazard ratio, 7.7; P = .05). Absence of detectable IgH-V(D)J NGS-MRD pre-HCT defines good-risk patients potentially eligible for less intense treatment approaches. Post-HCT NGS-MRD is highly predictive of relapse and survival, suggesting a role for this technique in defining patients early who would be eligible for post-HCT interventions. The trial was registered at www.clinicaltrials.gov as #NCT00382109.

Consistency of VDJ Rearrangement and Substitution Parameters Enables Accurate B Cell Receptor Sequence Annotation.

VDJ rearrangement and somatic hypermutation work together to produce antibody-coding B cell receptor (BCR) sequences for a remarkable diversity of antigens. It is now possible to sequence these BCRs in high throughput; analysis of these sequences is bringing new insight into how antibodies develop, in particular for broadly-neutralizing antibodies against HIV and influenza. A fundamental step in such sequence analysis is to annotate each base as coming from a specific one of the V, D, or J genes, or from an N-addition (a.k.a. non-templated insertion). Previous work has used simple parametric distributions to model transitions from state to state in a hidden Markov model (HMM) of VDJ recombination, and assumed that mutations occur via the same process across sites. However, codon frame and other effects have been observed to violate these parametric assumptions for such coding sequences, suggesting that a non-parametric approach to modeling the recombination process could be useful. In our paper, we find that indeed large modern data sets suggest a model using parameter-rich per-allele categorical distributions for HMM transition probabilities and per-allele-per-position mutation probabilities, and that using such a model for inference leads to significantly improved results. We present an accurate and efficient BCR sequence annotation software package using a novel HMM "factorization" strategy. This package, called partis (https://github.com/psathyrella/partis/), is built on a new general-purpose HMM compiler that can perform efficient inference given a simple text description of an HMM.

Reversing Tolerance in Isotype Switch-Competent Anti-Insulin B Lymphocytes.

Autoreactive B lymphocytes that escape central tolerance and mature in the periphery are a liability for developing autoimmunity. IgG insulin autoantibodies that predict type 1 diabetes and complicate insulin therapies indicate that mechanisms for tolerance to insulin are flawed. To examine peripheral tolerance in anti-insulin B cells, we generated C57BL/6 mice that harbor anti-insulin VDJH-125 site directed to the native IgH locus (VH125(SD)). Class switch-competent anti-insulin B cells fail to produce IgG Abs following T cell-dependent immunization of VH125(SD) mice with heterologous insulin, and they exhibit markedly impaired proliferation to anti-CD40 plus insulin in vitro. In contrast, costimulation with LPS plus insulin drives robust anti-insulin B cell proliferation. Furthermore, VH125(SD) mice produce both IgM and IgG2a anti-insulin Abs following immunization with insulin conjugated to type 1 T cell-independent Brucella abortus ring test Ag (BRT). Anti-insulin B cells undergo clonal expansion in vivo and emerge as IgM(+) and IgM(-) GL7(+)Fas(+) germinal center (GC) B cells following immunization with insulin-BRT, but not BRT alone. Analysis of Igκ genes in VH125(SD) mice immunized with insulin-BRT reveals that anti-insulin Vκ from the preimmune repertoire is selected into GCs. These data demonstrate that class switch-competent anti-insulin B cells remain functionally silent in T cell-dependent immune responses, yet these B cells are vulnerable to reversal of anergy following combined BCR/TLR engagement that promotes Ag-specific GC responses and Ab production. Environmental factors that lead to infection and inflammation could play a critical yet underappreciated role in driving loss of tolerance and promoting autoimmune disease.