Models for antigen receptor gene rearrangement: CDR3 length

Open-source antibodies as a path to enhanced research reproducibility and transparency

Antibodies are important tools with diverse uses in biomedical research. However, open access to reliable sources of well-characterized antibodies with unambiguous molecular identities remains an obstacle to research transparency and reproducibility. We propose here a community shift towards open-source antibodies, analogous to open-source computer software. The tenets of such antibodies are that 1) they are available to researchers in a ready to use form, 2) the renewable source of the antibody (e.g., hybridoma cells or plasmid) is also widely available ensuring reproducible and cost-effective access to the same antibody, and 3) the antibody sequence is publicly available. With these criteria met, the antibody can be widely used with the transparent assurance associated with a molecularly defined reagent, and the code can be edited to generate antibody variants to meet researchers' specific needs. We (the UC Davis/NIH NeuroMab Facility, the Development Studies Hybridoma Bank, and Addgene) have established a consortium to provide open-source access to a large collection of well characterized antibodies. As open-source software has benefitted both users and developers, we suggest open-source antibodies will have a similar positive impact on antibody based biomedical research. We encourage funding agencies to support initiatives to expand access to open-source antibody resources, and researchers to both utilize and to contribute to them, with a goal of enabling more reliable and cost-effective pursuit of research.

New insights into antibody structure with implications for specificity, variable region restriction and isotype choice

The mystery surrounding the mechanisms by which antibody diversity is generated was largely settled in the 1970s by the discoveries of variable gene rearrangements and somatic hypermutation. This led to the paradigm that immunoglobulins are composed of two independent domains - variable and constant - that confer specificity and effector functions, respectively. However, since these early discoveries, there have been a series of observations of communication between the variable and constant domains that affects the overall antibody structure, which suggests that immunoglobulins have a more complex, interconnected functionality than previously thought. Another unresolved issue has been the genesis of 'restricted' antibody responses, characterized by the use of only a few variable region gene segments, despite the enormous potential combinatorial diversity. In this Perspective, we place recent findings related to immunoglobulin structure and function in the context of these immunologically important, historically unsolved problems to propose a new model for how antibody specificity is achieved without autoreactivity.

Profiling the broad antibody diversity of lymphatic filariasis immune antibody repertoire by deep sequencing

Lymphatic filariasis is caused by infections of thread-like filarial worms, namely Wuchereria bancrofti, Brugia Malayi and Brugia timori. However, in-depth analysis of the antibody repertoire against Lymphatic filariasis is lacking. Using high-throughput sequencing of antibody repertoires, immunome analysis of IgG (LG) and IgM (LM) repertoires were studied. Despite significant differences between LG and LM in V(D)J gene usage, IGHV4-34, IGHV6-1, IGHD3-10 and IGHJ4 were preferred in both repertoires. The CDR3 in the LG repertoire showed a longer length than LM. Higher SHM level were observed in LG sequences and presence of oligoclonal sequences indicates the extent of clonal expansion. The prevalence of rare clonotypes in LM repertoire depicts the high clonal diversity when compared to LG repertoire. Monoclonal antibodies against closely related parasitic infections were present within the LG repertoire suggesting that immune repertoires may not be as exclusive and biased against the target infection as initially thought. The characterization of the immune repertoire can provide critical insight into the antibody response patterns in disease state, antibody generation process during infections and future antibody designs.

Characterizing adjuvants' effects at murine immunoglobulin repertoire level

Generating large-scale, high-fidelity sequencing data is challenging and, furthermore, not much has been done to characterize adjuvants' effects at the repertoire level. Thus, we introduced an IgSeq pipeline that standardized library prep protocols and data analysis functions for accurate repertoire profiling. We then studied systemically effects of CpG and Alum on the Ig heavy chain repertoire using the ovalbumin (OVA) murine model. Ig repertoires of different tissues (spleen and bone marrow) and isotypes (IgG and IgM) were examined and compared in IGHV mutation, gene usage, CDR3 length, clonal diversity, and clonal selection. We found Ig repertoires of different compartments exhibited distinguishable profiles at the non-immunized steady state, and distinctions became more pronounced upon adjuvanted immunizations. Notably, Alum and CpG effects exhibited different tissue- and isotype-preferences. The former led to increased diversity of abundant clones in bone marrow, and the latter promoted the selection of IgG clones in both tissues.

An Approach for Antigen-Agnostic Identification of Virus-Like Particle-Displayed Epitopes that Engage Specific Antibody V Gene Regions

Antibody complementarity determining regions (CDRs) participate in antigen recognition, but not all participate equally in antigen binding. Here we describe a technique for discovering strong, specific binding partners to defined motifs within the CDRs of chimeric, engineered antibodies using affinity selection and counter-selection of epitopes displayed on bacteriophage MS2-based virus-like particles (VLPs). As an example, we show how this technique can be used to identify families of VLPs that interact with antibodies displaying the CDRs encoded by the germline precursor of a broadly neutralizing monoclonal antibody against HIV-1.

Massively-multiplexed epitope mapping techniques for viral antigen discovery

Following viral infection, viral antigens bind specifically to receptors on the surface of lymphocytes thereby activating adaptive immunity in the host. An epitope, the smallest structural and functional unit of an antigen, binds specifically to an antibody or antigen receptor, to serve as key sites for the activation of adaptive immunity. The complexity and diverse range of epitopes are essential to study and map for the diagnosis of disease, the design of vaccines and for immunotherapy. Mapping the location of these specific epitopes has become a hot topic in immunology and immune therapy. Recently, epitope mapping techniques have evolved to become multiplexed, with the advent of high-throughput sequencing and techniques such as bacteriophage-display libraries and deep mutational scanning. Here, we briefly introduce the principles, advantages, and disadvantages of the latest epitope mapping techniques with examples for viral antigen discovery.

Modeling and predicting the overlap of B- and T-cell receptor repertoires in healthy and SARS-CoV-2 infected individuals

Adaptive immunity's success relies on the extraordinary diversity of protein receptors on B and T cell membranes. Despite this diversity, the existence of public receptors shared by many individuals gives hope for developing population-wide vaccines and therapeutics. Using probabilistic modeling, we show many of these public receptors are shared by chance in healthy individuals. This predictable overlap is driven not only by biases in the random generation process of receptors, as previously reported, but also by their common functional selection. However, the model underestimates sharing between repertoires of individuals infected with SARS-CoV-2, suggesting strong specific antigen-driven convergent selection. We exploit this discrepancy to identify COVID-associated receptors, which we validate against datasets of receptors with known viral specificity. We study their properties in terms of sequence features and network organization, and use them to design an accurate diagnostic tool for predicting SARS-CoV-2 status from repertoire data.

Public Immunity: Evolutionary Spandrels for Pathway-Amplifying Protective Antibodies

Humoral immunity is seeded by affinity between the B cell receptor (BCR) and cognate antigen. While the BCR is a chimeric display of diverse antigen engagement solutions, we discuss its functional activity as an ‘innate-like’ immune receptor, wherein genetically hardwired antigen complementarity can serve as reproducible templates for pathway-amplifying otherwise immunologically recessive antibody responses. We propose that the capacity for germline reactivity to new antigen emerged as a set of evolutionary spandrels or coupled traits, which can now be exploited by rational vaccine design to focus humoral immunity upon conventionally immune-subdominant antibody targets. Accordingly, we suggest that evolutionary spandrels account for the necessary but unanticipated antigen reactivity of the germline antibody repertoire.

Engineering an Antibody V Gene-Selective Vaccine

The ligand-binding surface of the B cell receptor (BCR) is formed by encoded and non-encoded antigen complementarity determining regions (CDRs). Genetically reproducible or ‘public’ antibodies can arise when the encoded CDRs play deterministic roles in antigen recognition, notably within human broadly neutralizing antibodies against HIV and influenza virus. We sought to exploit this by engineering virus-like-particle (VLP) vaccines that harbor multivalent affinity against gene-encoded moieties of the BCR antigen binding site. As proof of concept, we deployed a library of RNA bacteriophage VLPs displaying random peptides to identify a multivalent antigen that selectively triggered germline BCRs using the human V_H gene IGVH1-2*02. This VLP selectively primed IGHV1-2*02 BCRs that were present within a highly diversified germline antibody repertoire within humanized mice. Our approach thus provides methodology to generate antigens that engage specific BCR configurations of interest, in the absence of structure-based information.

An Analysis of the Effects of Spaceflight and Vaccination on Antibody Repertoire Diversity

Ab repertoire diversity plays a critical role in the host's ability to fight pathogens. CDR3 is partially responsible for Ab-Ag binding and is a significant source of diversity in the repertoire. CDR3 diversity is generated during VDJ rearrangement because of gene segment selection, gene segment trimming and splicing, and the addition of nucleotides. We analyzed the Ab repertoire diversity across multiple experiments examining the effects of spaceflight on the Ab repertoire after vaccination. Five datasets from four experiments were analyzed using rank-abundance curves and Shannon indices as measures of diversity. We discovered a trend toward lower diversity as a result of spaceflight but did not find the same decrease in our physiological model of microgravity in either the spleen or bone marrow. However, the bone marrow repertoire showed a reduction in diversity after vaccination. We also detected differences in Shannon indices between experiments and tissues. We did not detect a pattern of CDR3 usage across the experiments. Overall, we were able to find differences in the Ab repertoire diversity across experimental groups and tissues.

A Single Human VH-gene Allows for a Broad-Spectrum Antibody Response Targeting Bacterial Lipopolysaccharides in the Blood

B cell receptors (BCRs) display a combination of variable (V)-gene-encoded complementarity determining regions (CDRs) and adaptive/hypervariable CDR3 loops to engage antigens. It has long been proposed that the former tune for recognition of pathogens or groups of pathogens. To experimentally evaluate this within the human antibody repertoire, we perform immune challenges in transgenic mice that bear diverse human CDR3 and light chains but are constrained to different human V_H-genes. We find that, of six commonly deployed V_H sequences, only those CDRs encoded by IGHV1-2^∗02 enable polyclonal antibody responses against bacterial lipopolysaccharide (LPS) when introduced to the bloodstream. The LPS is from diverse strains of gram-negative bacteria, and the V_H-gene-dependent responses are directed against the non-variable and universal saccrolipid substructure of this antigen. This reveals a broad-spectrum anti-LPS response in which germline-encoded CDRs naturally hardwire the human antibody repertoire for recognition of a conserved microbial target.

Human-likeness of antibody biologics determined by back-translation and comparison with large antibody variable gene repertoires

The antibody (Ab) germline gene rearrangement of variable (V), diversity (D), and joining (J) gene segments, as well as somatic hypermutation, give rise to the human Ab variable gene sequence repertoire. It is common to characterize single nucleotide frequencies of the variable region by alignment to species-specific wildtype germline genes. The increasing application of next-generation sequencing to immune repertoire studies has led to the compilation of increasing large adaptive immunome receptor repertoire datasets. We have developed a method that maps the sequence of a target Ab onto an immunome dataset of 326 million human Ab sequences. For this purpose, we created a position- and gene-specific scoring matrix (PGSSM) and its corresponding antibody similarity score. We characterized our PGSSM score and found that it strongly correlated with the phylogenetic distance of 181,355 Ab sequences from GenBank across 20 species. The most likely human nucleotide back-translation was obtained given only PGSSMs and the amino acid sequence of an Ab achieving a nucleotide sequence recovery of 95.9% and 97.2% for human heavy and light chains, respectively. In conclusion, the scoring of our back-translation is a valuable estimate for the similarity of an Ab sequence to the natural human repertoire. As expected, Ab therapeutic molecules developed from a human source showed a higher similarity to the repertoire than engineered Abs. Thus, the PGSSM metric introduced here can be used to engineer human-like Ab therapeutics.

A Quest for a Universal Plasma-Derived Antivenom Against All Elapid Neurotoxic Snake Venoms

This review describes the research aimed at the development of universal antivenom against elapid neurotoxic snake venoms. The antivenoms produced in Thailand in the 1980s were of low potency, especially against the elapid venoms. This was thought to be due to the low immunogenicity of the α-neurotoxins, which are the most lethal toxins in these venoms. Comparisons of various α-neurotoxin conjugates and polymers, and also different immunological adjuvants, showed that the adjuvant used is the major determinant in the antibody response in horses. The potent Freund's adjuvant was not used due to its severe local side-effect in horses. Therefore, a novel immunization protocol termed 'low dose, low volume multi-site' was developed for use in horses. This immunization protocol has led to the production of highly potent monospecific antivenoms against several elapid and viperid venoms, and two potent polyspecific antivenoms, one against 4 neurotoxic and another against 3 hematotoxic venoms. The immunization protocol has also led to other improvements in antivenom production including: several fold increases in antiserum potency, a reduction in the time required to reach therapeutically useful antibody titers, a 90% reduction in the amount of venom used, and 100% of the horses responding to the immunization program. This development is partly responsible for significant decrease in the Thailand's annual snakebite death toll from a few dozens to mostly nil in recent years. Finally, a simple and novel immunization strategy, using a 'diverse toxin repertoire' composed of numerous elapid toxin fractions as immunogen, was proposed and tested. This immunization procedure has resulted in the successful production of a widely paraspecific antiserum against at least 36 neurotoxic venoms of 28 species encompassing 10 genera and from 20 countries on four continents, and possibly against all elapid venoms with α-neurotoxins as the lethal toxins. These results indicate that, with optimizations of the composition of the 'diverse toxin repertoire', the immunization scheme and antibody fractionation to increase the antivenom neutralizing potency, an effective universal antivenom against the neurotoxic elapid snakes of the world can be produced.

A single donor is sufficient to produce a highly functional in vitro antibody library

Antibody complementarity determining region diversity has been considered to be the most important metric for the production of a functional antibody library. Generally, the greater the antibody library diversity, the greater the probability of selecting a diverse array of high affinity leads. According to this paradigm, the primary means of elevating library diversity has been by increasing the number of donors. In the present study we explored the possibility of creating an in vitro antibody library from a single healthy individual, showing that the number of lymphocytes, rather than the number of donors, is the key criterion in the production of a diverse and functional antibody library. We describe the construction of a high-quality phage display library comprising 5 × 10⁹ human antibodies by applying an efficient B cell extraction protocol from a single donor and a targeted V-gene amplification strategy favoring specific antibody families for their improved developability profiles. Each step of the library generation process was followed and validated by next generation sequencing to monitor the library quality and diversity. The functionality of the library was tested using several therapeutically relevant targets for which a vast number of different antibodies with desired biophysical properties were obtained.

Antibody Focusing to Conserved Sites of Vulnerability: The Immunological Pathways for 'Universal' Influenza Vaccines

Influenza virus remains a serious public health burden due to ongoing viral evolution. Vaccination remains the best measure of prophylaxis, yet current seasonal vaccines elicit strain-specific neutralizing responses that favor the hypervariable epitopes on the virus. This necessitates yearly reformulations of seasonal vaccines, which can be limited in efficacy and also shortchange pandemic preparedness. Universal vaccine development aims to overcome these deficits by redirecting antibody responses to functionally conserved sites of viral vulnerability to enable broad coverage. However, this is challenging as such antibodies are largely immunologically silent, both following vaccination and infection. Defining and then overcoming the immunological basis for such subdominant or ‘immuno-recessive’ antibody targeting has thus become an important aspect of universal vaccine development. This, coupled with structure-guided immunogen design, has led to proof-of-concept that it is possible to rationally refocus humoral immunity upon normally ‘unseen’ broadly neutralizing antibody targets on influenza virus.

Defining best practices for tissue procurement in immuno-oncology clinical trials: consensus statement from the Society for Immunotherapy of Cancer Surgery Committee

Immunotherapy is now a cornerstone for cancer treatment, and much attention has been placed on the identification of prognostic and predictive biomarkers. The success of biomarker development is dependent on accurate and timely collection of biospecimens and high-quality processing, storage and shipping. Tumors are also increasingly used as source material for the generation of therapeutic T cells. There have been few guidelines or consensus statements on how to optimally collect and manage biospecimens and source material being used for immunotherapy and related research. The Society for Immunotherapy of Cancer Surgery Committee has brought together surgical experts from multiple subspecialty disciplines to identify best practices and to provide consensus on how best to access and manage specific tissues for immuno-oncology treatments and clinical investigation. In addition, the committee recommends early integration of surgeons and other interventional physicians with expertise in biospecimen collection, especially in clinical trials, to optimize the quality of tissue and the validity of correlative clinical studies in cancer immunotherapy.

The Adaptome as Biomarker for Assessing Cancer Immunity and Immunotherapy

In terms of diagnosing and treating diseases, our adaptive immune system is the "best doctor." It carries out these tasks with unmatched precision, with the help of both T and B cell receptors, our most diverse set of genes, distinguishing one individual from another. It does this by generating autologous extraordinary diversity in the receptors, ranging from 10¹⁵ to 10²⁵ for each chain of the rearranged receptors. By combining multiplex PCR and next-generation sequencing (NGS), we have developed high throughput methods to study adaptive immunity. The adaptome is the sum-total of expressed T and B cell receptor genes in a sample, composed of seven chains, including the alpha/beta and gamma/delta chains for T cells, and heavy/lambda or kappa chains for B cells. Immune repertoire is the sum-total of the individual clonotypes within one chain, including individual complementarity-determining regions (CDR) 3 sequences. In order to reflect the breadth and depth of the true adaptome, the following criteria assessing any method needs to be ascertained: (1) Methods need to be inclusive and quantitative; (2) Analysis should consider what questions need to be addressed and whether bulk or single cell sequencing provide the best tools for assessing the underlying biology and addressing important questions; (3) Measures of clonal diversity are key to understand the underlying structure and providence of the repertoire; and (4) Convergent evolution may allow a surprising degree of homologous or identical CDR3s associated with individual disease entities, creating hope for novel diagnostics and/or disease burden assessments. Integrating studies of the peripheral blood, lymph nodes, and tumor allows dynamic interrogation of the alterations occurring with age, treatment, and response to emergent and established therapies.

Aromatic Guanylhydrazones for the Control of Heme-Induced Antibody Polyreactivity

In a healthy immune repertoire, there exists a fraction of polyreactive antibodies that can bind to a variety of unrelated self- and foreign antigens. Apart from naturally polyreactive antibodies, in every healthy individual, there is a fraction of antibody that can gain polyreactivity upon exposure to porphyrin cofactor heme. Molecular mechanisms and biological significance of the appearance of cryptic polyreactivity are not well understood. It is believed that heme acts as an interfacial cofactor between the antibody and the newly recognized antigens. To further test this claim and gain insight into the types of interactions involved in heme binding, we herein investigated the influence of a group of aromatic guanylhydrazone molecules on the heme-induced antibody polyreactivity. From the analysis of SAR and the results of UV-vis absorbance spectroscopy, it was concluded that the most probable mechanism by which the studied molecules inhibit heme-mediated polyreactivity of the antibody is the direct binding to heme, thus preventing heme from binding to antibody and/or antigen. The inhibitory capacity of the most potent compounds was substantially higher than that of chloroquine, a well-known heme binder. Some of the guanylhydrazone molecules were able to induce polyreactivity of the studied antibody themselves, possibly by a mechanism similar to heme. Results described here point to the conclusion that heme indeed must bind to an antibody to induce its polyreactivity, and that both π-stacking interactions and iron coordination contribute to the binding affinity, while certain structures, such as guanylhydrazones, can interfere with these processes.

Comparative analysis of CDR3 regions in paired human αβ CD8 T cells

The majority of human CD8 cytotoxic T lymphocytes express αβ T-cell receptors that recognize peptide-MHC class I complexes. Considerable attention has been devoted to TCR β repertoires, but study of TCR α chains has been limited. To gain a better understanding of the features of CDR3α and CDR3β in paired samples, we comprehensively analyzed 776 unique paired αβ TCR CDR3 regions in this study. We found that (I) the CDR3 length among paired αβ TCRs had a fairly narrow distribution due to random assortment of CDR3 length in alpha and beta chains; (II) nucleotide deletions among CDR3 regions were positively correlated with insertions in both α and β TCRs; (III) the CDR3 loops of both α and β chains contained an abundance of charged/polar residues and the CDR3 base regions contained a conserved motif; and (IV) the occurrence of Gly was CDR3 length- and position-dependent in both chains, whereas the frequency of Ser at positions 106 and 107 was positively correlated with CDR3 length in TCR β. Overall, the amino acids in CDR3 loop regions were significantly different between TCR α and β, which suggests a distinct role for each chain in the recognition of antigen-MHC complexes. Here, we have provided detailed information on CDR3 in paired TCRs expressed on human CD8+ T cells and established the basis of a reference set for αβ TCR repertoires in healthy humans.

Innate-like B cell subsets during immune responses: Beyond antibody production

B lymphocytes are recognized for their crucial role in the adaptive immunity since they represent the only leukocyte lineage capable of differentiating into Ab-secreting cells. However, it has been demonstrated that these lymphocytes can exert several Ab-independent functions, including engulfing and processing Ags for presentation to T cells, secreting soluble mediators, providing co-stimulatory signals, and even participating in lymphoid tissues development. Beyond that, several reports claiming the existence of multiple B cell subsets contributing directly to innate immune responses have appeared. These "innate-like" B lymphocytes, whose phenotype, development pathways, tissue distribution, and functions are in most cases notoriously different from those of conventional B cells, are crucial to early protective responses against pathogens by exerting "crossover" defensive strategies that blur the established boundaries of innate and adaptive branches of immunity. Examples of these mechanisms include the rapid secretion of the polyspecific natural Abs, increased susceptibility to innate receptors-mediated activation, cytokine secretion, downstream priming of other innate cells, usage of specific variable immunoglobulin gene-segments, and other features. As these new insights emerge, it is becoming preponderant to redefine the functionality of B cells beyond their classical adaptive-immune tasks.

Comparison of methods for phylogenetic B-cell lineage inference using time-resolved antibody repertoire simulations (AbSim)

Motivation

The evolution of antibody repertoires represents a hallmark feature of adaptive B-cell immunity. Recent advancements in high-throughput sequencing have dramatically increased the resolution to which we can measure the molecular diversity of antibody repertoires, thereby offering for the first time the possibility to capture the antigen-driven evolution of B cells. However, there does not exist a repertoire simulation framework yet that enables the comparison of commonly utilized phylogenetic methods with regard to their accuracy in inferring antibody evolution.

Results

Here, we developed AbSim, a time-resolved antibody repertoire simulation framework, which we exploited for testing the accuracy of methods for the phylogenetic reconstruction of B-cell lineages and antibody molecular evolution. AbSim enables the (i) simulation of intermediate stages of antibody sequence evolution and (ii) the modeling of immunologically relevant parameters such as duration of repertoire evolution, and the method and frequency of mutations. First, we validated that our repertoire simulation framework recreates replicates topological similarities observed in experimental sequencing data. Second, we leveraged Absim to show that current methods fail to a certain extent to predict the true phylogenetic tree correctly. Finally, we formulated simulation-validated guidelines for antibody evolution, which in the future will enable the development of accurate phylogenetic methods.

Availability and implementation

https://cran.r-project.org/web/packages/AbSim/index.html.

Contact

sai.reddy@ethz.ch.

Supplementary information

Supplementary data are available at Bioinformatics online.

Peptides of the variable IgG domain as potential biomarker candidates in primary open-angle glaucoma (POAG)

Autoantibody profiling has gained increasing interest in the research field of glaucoma promising the detection of highly specific and sensitive marker candidates for future diagnostic purposes. Recent studies demonstrated that immune responses are characterized by the expression of congruent or similar complementarity determining regions (CDR) in different individuals and could be used as molecular targets in biomarker discovery. Main objective of this study was to characterize glaucoma-specific peptides from the variable region of sera-derived immunoglobulins using liquid chromatography--mass spectrometry (LC-MS)-based quantitative proteomics. IgG was purified from sera of 13 primary open-angle glaucoma patients (POAG) and 15 controls (CTRL) and subsequently digested into Fab and Fc by papain. Fab was further purified, tryptic digested and measured by LC-MS/MS. Discovery proteomics revealed in total 75 peptides of the variable IgG domain showing significant glaucoma-related level changes (P < 0.05; log2 fold change ≥ 0.5): 6 peptides were high abundant in POAG sera, whereas 69 peptides were low abundant in comparison to CTRL group. Via accurate inclusion mass screening strategy 28 IgG V domain peptides were further validated showing significantly decreased expression levels in POAG sera. Amongst others 5 CDR1, 2 CDR2 and 1 CDR3 sequences. In addition, we observed significant shifts in the variable heavy chain family distribution and disturbed κ/λ ratios in POAG patients in contrast to CTRL. These findings strongly indicate that glaucoma is accompanied by systemic effects on antibody production and B cell maturation possibly offering new prospects for future diagnostic or therapy purposes.

Integrated analysis of B‑cell and T‑cell receptors by high‑throughput sequencing reveals conserved repertoires in IgA nephropathy

Immunoglobulin A nephropathy (IgAN) is a type of glomerular disorder associated with immune dysregulation, and understanding B‑/T‑cell receptors (BCRs/TCRs) may be valuable for the development of specific immunotherapeutic interventions. In the present study, B and T cells were isolated from IgAN patients and healthy controls, and the composition of the BCR/TCR complementarity‑determining region (CDR)3 was analyzed by multiplex polymerase chain reaction, high‑throughput sequencing and bioinformatics. The present results revealed that the BCR/TCR CDR3 clones were expressed at very low frequencies, and the composition of clone types in patients with IgAN was skewed; the majority of clones were unique, and only 12 BCR and 228 TCR CDR3 clones were public ones, of which 16 were expressed at a significantly higher frequency in patients with IgAN (P<0.001). There were also certain conserved amino acid residues between unique clones or groups, and the residues GMDV, EQY and EQF were recurring only in the IgAN group. In addition, some VDJ gene recombinations indicated great variation between groups, including 4 high‑frequency VDJ gene recombinations in the IgAN patients (P<0.001). Immune repertoires provide novel information, and conserved BCR/TCR CDR3 clones and VDJ gene recombinations with great variation may be potential therapeutic targets for IgAN patients.

Single-cell analysis reveals sexually dimorphic repertoires of Interferon-γ and IL-17A producing T cells in salivary glands of Sjögren's syndrome mice

The development of Sjögren's syndrome (SjS) is a dynamic and temporal process with a female predilection. Following the initial influx of immune cells, T cell clusters develop, accelerating the pathology in the salivary glands. Proinflammatory cytokines, IFN-γ and IL-17A, produced by T cells contribute synergistically to the disease. In this study, we examined the sexual dimorphism in cellular infiltrates of the salivary glands by using functional single-cell microengraving analysis. Using high-throughput sequencing, we investigated the clonal diversity of the T cell receptors (TCRs) of infiltrating IFN-γ and IL-17A-producing T cells in male and female SjS-susceptible (SjSs) C57BL/6.NOD-Aec1Aec2 mice. There were elevated frequencies of IFN-γ and IL-17A-producing effector T cell populations in female SjSS mice compared to male SjSS mice. MEME analysis shows high frequency and unique, sexually dimorphic motifs in the TCR hypervariable regions in the SjSS mice. Male mice selected for TRAV8/TRAJ52 (CATDLNTGANTGKLTFG) TCR genes in Th1 cells and TRBV16/(TRBD1/2)TRBJ1-7 (CGGKRRLESIFR) in Th1 and Th17 cells. Female SjSS mice selected for TRAV8/TRAJ52 (CATDLNTGANTGKLTFG), TRAV13D-2/TRAJ23 (CVYLEHHFE), and TRBV23/(TRBD2)TRBJ2-2 (CRKLHSCATCALNFL) in Th1 cells. These findings suggest that there is an elevated prevalence of pathogenic effector T cells in the glands with a sexually dimorphic selection bias of TCR repertoires.

Early Emergence of CD19-Negative Human Antibody-Secreting Cells at the Plasmablast to Plasma Cell Transition

Long-lived human plasma cells (PCs) play central roles in immunity and autoimmunity and are enriched among the subpopulation of CD19neg human PCs. However, whether human CD19neg PCs are necessarily aged cells that have gradually lost CD19 expression is not known. Assessing peripheral blood samples at steady-state and during the acute response to influenza vaccination in healthy donors, we identify the presence of phenotypic CD19neg plasmablasts, the proliferative precursor state to mature PCs, and demonstrate by ELISPOT that these are Ab-secreting cells (ASCs). During the acute response to influenza vaccination, CD19pos, CD19low, and CD19neg ASCs secrete vaccine-specific Abs and show linked IGHV repertoires. To address precursor/product relationships, we use in vitro models that mimic T-dependent and T-independent differentiation, finding that the CD19neg state can be established at the plasmablast to PC transition, that CD19neg PCs increase as a percentage of surviving PCs in vitro, and that CD19neg and CD19pos PCs can be maintained independently. These data provide proof-of-principle for the view that newly generated ASCs can acquire a mature PC phenotype that is accompanied by loss of CD19 expression at an early stage of differentiation and that aging is not an obligate requirement for a CD19neg state to be established.

Investigating the effect of AS03 adjuvant on the plasma cell repertoire following pH1N1 influenza vaccination

Influenza pandemics require rapid deployment of effective vaccines for control. Adjuvants such as AS03 improve vaccine immunogenicity, but this mechanism is poorly understood. We used high-throughput B cell receptor sequencing of plasma cells produced following AS03-adjuvanted and non-adjuvanted 2009 pandemic H1N1 vaccination, as well as pre-pandemic seasonal influenza vaccination to elucidate the effect of the adjuvant on the humoral immune response. By analyzing mutation levels, it was possible to distinguish sequences from cells that were recently activated from naïve B cells from those that were activated by memory recall. We show that the adjuvant functions through two mechanisms. First, the adjuvant stimulates increased activation of naïve B cells, thus reducing immune interference with previous vaccine responses. Second, the adjuvant is able to increase the adaptability of the recalled cells to give improved specificity to the new vaccine antigen. We thus show how AS03 enhances pH1N1 immune responses, and reduces immune interference.

Likelihood-Based Inference of B Cell Clonal Families

The human immune system depends on a highly diverse collection of antibody-making B cells. B cell receptor sequence diversity is generated by a random recombination process called “rearrangement” forming progenitor B cells, then a Darwinian process of lineage diversification and selection called “affinity maturation.” The resulting receptors can be sequenced in high throughput for research and diagnostics. Such a collection of sequences contains a mixture of various lineages, each of which may be quite numerous, or may consist of only a single member. As a step to understanding the process and result of this diversification, one may wish to reconstruct lineage membership, i.e. to cluster sampled sequences according to which came from the same rearrangement events. We call this clustering problem “clonal family inference.” In this paper we describe and validate a likelihood-based framework for clonal family inference based on a multi-hidden Markov Model (multi-HMM) framework for B cell receptor sequences. We describe an agglomerative algorithm to find a maximum likelihood clustering, two approximate algorithms with various trade-offs of speed versus accuracy, and a third, fast algorithm for finding specific lineages. We show that under simulation these algorithms greatly improve upon existing clonal family inference methods, and that they also give significantly different clusters than previous methods when applied to two real data sets.

Antibodies must recognize a great diversity of antigens to protect us from infectious disease. The binding properties of antibodies are determined by the DNA sequences of their corresponding B cell receptors (BCRs). These BCR sequences are created in naive form by VDJ recombination, which randomly selects and trims the ends of V, D, and J genes, then joins the resulting segments together with additional random nucleotides. If they pass initial screening and bind an antigen, these sequences then undergo an evolutionary process of reproduction, mutation, and selection, revising the BCR to improve binding to its cognate antigen. It has recently become possible to determine the BCR sequences resulting from this process in high throughput. Although these sequences implicitly contain a wealth of information about both antigen exposure and the process by which we learn to resist pathogens, this information can only be extracted using computer algorithms. In this paper we describe a likelihood-based statistical method to determine, given a collection of BCR sequences, which of them are derived from the same recombination events. It is based on a hidden Markov model (HMM) of VDJ rearrangement which is able to calculate likelihoods for many sequences at once.

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.

Analysis of B Cell Repertoire Dynamics Following Hepatitis B Vaccination in Humans, and Enrichment of Vaccine-specific Antibody Sequences

Generating a diverse B cell immunoglobulin repertoire is essential for protection against infection. The repertoire in humans can now be comprehensively measured by high-throughput sequencing. Using hepatitis B vaccination as a model, we determined how the total immunoglobulin sequence repertoire changes following antigen exposure in humans, and compared this to sequences from vaccine-specific sorted cells. Clonal sequence expansions were seen 7 days after vaccination, which correlated with vaccine-specific plasma cell numbers. These expansions caused an increase in mutation, and a decrease in diversity and complementarity-determining region 3 sequence length in the repertoire. We also saw an increase in sequence convergence between participants 14 and 21 days after vaccination, coinciding with an increase of vaccine-specific memory cells. These features allowed development of a model for in silico enrichment of vaccine-specific sequences from the total repertoire. Identifying antigen-specific sequences from total repertoire data could aid our understanding B cell driven immunity, and be used for disease diagnostics and vaccine evaluation.

BCR repertoire sequencing: different patterns of B-cell activation after two Meningococcal vaccines

Next-generation sequencing was used to investigate the B-cell receptor heavy chain transcript repertoire of different B-cell subsets (naive, marginal zone (MZ), immunoglobulin M (IgM) memory and IgG memory) at baseline, and of plasma cells (PCs) 7 days following administration of serogroup ACWY meningococcal polysaccharide and protein-polysaccharide conjugate vaccines. Baseline B-cell subsets could be distinguished from each other using a small number of repertoire properties (clonality, mutation from germline and complementarity-determining region 3 (CDR3) length) that were conserved between individuals. However, analyzing the CDR3 amino-acid sequence (which is particularly important for antigen binding) of the baseline subsets showed few sequences shared between individuals. In contrast, day 7 PCs demonstrated nearly 10-fold greater sequence sharing between individuals than the baseline subsets, consistent with the PCs being induced by the vaccine antigen and sharing specificity for a more limited range of epitopes. By annotating PC sequences based on IgG subclass usage and mutation, and also comparing them with the sequences of the baseline cell subsets, we were able to identify different signatures after the polysaccharide and conjugate vaccines. PCs produced after conjugate vaccination were predominantly IgG1, and most related to IgG memory cells. In contrast, after polysaccharide vaccination, the PCs were predominantly IgG2, less mutated and were equally likely to be related to MZ, IgM memory or IgG memory cells. High-throughput B-cell repertoire sequencing thus provides a unique insight into patterns of B-cell activation not possible from more conventional measures of immunogenicity.

Antibody Repertoires in Fish

As in mammals, cartilaginous and teleost fishes possess adaptive immune systems based on antigen recognition by immunoglobulins (Ig), T cell receptors (TCR), and major histocompatibility complex molecules (MHC) I and MHC II molecules. Also it is well established that fish B cells and mammalian B cells share many similarities, including Ig gene rearrangements, and production of membrane Ig and secreted Ig forms. This chapter provides an overview of the IgH and IgL chains in cartilaginous and bony fish, including their gene organizations, expression, diversity of their isotypes, and development of the primary repertoire. Furthermore, when possible, we have included summaries of key studies on immune mechanisms such as allelic exclusion, somatic hypermutation, affinity maturation, class switching, and mucosal immune responses.

Old questions, new tools: does next-generation sequencing hold the key to unraveling intestinal B-cell responses?

Analysis of the intestinal B-cell system and properties of immunoglobulin A, the main antibody isotype produced in the gut, has dominated the rise of mucosal immunology as a discipline. Seminal work established concepts describing the induction, transport, and function of mucosal antibodies. Still, open questions remain and we lack a comprehensive view of how the various sites and pathways of immunoglobulin A induction are integrated to respond to gut antigens. Next-generation sequencing (NGS) offers a novel approach to study B-cell responses, which might substantially enhance our tool box to answer key questions in the field and to take the next steps toward therapeutic exploitation of the mucosal B-cell system. In this review we discuss the potential, challenges, and emerging solutions for gut B-cell repertoire analysis by NGS.

Peptides from the variable region of specific antibodies are shared among lung cancer patients

Late diagnosis of lung cancer is still the main reason for high mortality rates in lung cancer. Lung cancer is a heterogeneous disease which induces an immune response to different tumor antigens. Several methods for searching autoantibodies have been described that are based on known purified antigen panels. The aim of our study is to find evidence that parts of the antigen-binding-domain of antibodies are shared among lung cancer patients. This was investigated by a novel approach based on sequencing antigen-binding-fragments (Fab) of immunoglobulins using proteomic techniques without the need of previously known antigen panels. From serum of 93 participants of the NELSON trial IgG was isolated and subsequently digested into Fab and Fc. Fab was purified from the digested mixture by SDS-PAGE. The Fab containing gel-bands were excised, tryptic digested and measured on a nano-LC-Orbitrap-Mass-spectrometry system. Multivariate analysis of the mass spectrometry data by linear canonical discriminant analysis combined with stepwise logistic regression resulted in a 12-antibody-peptide model which was able to distinguish lung cancer patients from controls in a high risk population with a sensitivity of 84% and specificity of 90%. With our Fab-purification combined Orbitrap-mass-spectrometry approach, we found peptides from the variable-parts of antibodies which are shared among lung cancer patients.

RAG2's acidic hinge restricts repair-pathway choice and promotes genomic stability

V(D)J recombination-associated DNA double-strand breaks (DSBs) are normally repaired by the high-fidelity classical nonhomologous end-joining (cNHEJ) machinery. Previous studies implicated the recombination-activating gene (RAG)/DNA postcleavage complex (PCC) in regulating pathway choice by preventing access to inappropriate repair mechanisms such as homologous recombination (HR) and alternative NHEJ (aNHEJ). Here, we report that RAG2's "acidic hinge," previously of unknown function, is critical for several key steps. Mutations that reduce the hinge's negative charge destabilize the PCC, disrupt pathway choice, permit repair of RAG-mediated DSBs by the translocation-prone aNHEJ machinery, and reduce genomic stability in developing lymphocytes. Structural predictions and experimental results support our hypothesis that reduced flexibility of the hinge underlies these outcomes. Furthermore, sequence variants present in the human population reduce the hinge's negative charge, permit aNHEJ, and diminish genomic integrity.

The shape of the lymphocyte receptor repertoire: lessons from the B cell receptor

Both the B cell receptor (BCR) and the T cell receptor (TCR) repertoires are generated through essentially identical processes of V(D)J recombination, exonuclease trimming of germline genes, and the random addition of non-template encoded nucleotides. The naïve TCR repertoire is constrained by thymic selection, and TCR repertoire studies have therefore focused strongly on the diversity of MHC-binding complementarity determining region (CDR) CDR3. The process of somatic point mutations has given B cell studies a major focus on variable (IGHV, IGLV, and IGKV) genes. This in turn has influenced how both the naïve and memory BCR repertoires have been studied. Diversity (D) genes are also more easily identified in BCR VDJ rearrangements than in TCR VDJ rearrangements, and this has allowed the processes and elements that contribute to the incredible diversity of the immunoglobulin heavy chain CDR3 to be analyzed in detail. This diversity can be contrasted with that of the light chain where a small number of polypeptide sequences dominate the repertoire. Biases in the use of different germline genes, in gene processing, and in the addition of non-template encoded nucleotides appear to be intrinsic to the recombination process, imparting "shape" to the repertoire of rearranged genes as a result of differences spanning many orders of magnitude in the probabilities that different BCRs will be generated. This may function to increase the precursor frequency of naïve B cells with important specificities, and the likely emergence of such B cell lineages upon antigen exposure is discussed with reference to public and private T cell clonotypes.

Mass spectrometry analyses of κ and λ fractions result in increased number of complementarity-determining region identifications

Sera from lung cancer patients contain antibodies against tumor-associated antigens. Specific amino acid sequences of the complementarity-determining regions (CDRs) in the antigen-binding fragment (Fab) of these antibodies have potential as lung cancer biomarkers. Detection and identification of CDRs by mass spectrometry can significantly be improved by reduction of the complexity of the immunoglobulin molecule. Our aim was to molecular dissect IgG into κ and λ fragments to reduce the complexity and thereby identify substantially more CDRs than by just total Fab isolation. We purified Fab, Fab-κ, Fab-λ, κ and λ light chains from serum from 10 stage I lung adenocarcinoma patients and 10 matched controls from the current and former smokers. After purification, the immunoglobulin fragments were enzymatically digested and measured by high-resolution mass spectrometry. Finally, we compared the number of CDRs identified in these immunoglobulin fragments with that in the Fab fragments. Twice as many CDRs were identified when Fab-κ, Fab-λ, κ and λ (3330) were combined than in the Fab fraction (1663) alone. The number of CDRs and κ:λ ratio was statistically similar in both cases and controls. Molecular dissection of IgG identifies significantly more CDRs, which increases the likelihood of finding lung cancer-related CDR sequences.

Divergent human populations show extensive shared IGK rearrangements in peripheral blood B cells

We have analysed the transcribed immunoglobulin kappa (IGK) repertoire of peripheral blood B cells from four individuals from two genetically distinct populations, Papua New Guinean and Australian, using high-throughput DNA sequencing. The depth of sequencing data for each individual averaged 5,548 high-quality IGK reads, and permitted genotyping of the inferred IGKV and IGKJ germline gene segments for each individual. All individuals were homozygous at each IGKJ locus and had highly similar inferred IGKV genotypes. Preferential gene usage was seen at both the IGKV and IGKJ loci, but only IGKV segment usage varied significantly between individuals. Despite the differences in IGKV gene utilisation, the rearranged IGK repertoires showed extensive identity at the amino acid level. Public rearrangements (those shared by two or more individuals) made up 60.2% of the total sequenced IGK rearrangements. The total diversity of IGK rearrangements of each individual was estimated to range from just 340 to 549 unique amino acid sequences. Thus, the repertoire of unique expressed IGK rearrangements is dramatically less than previous theoretical estimates of IGK diversity, and the majority of expressed IGK rearrangements are likely to be extensively shared in individual human beings.

Factors important in evolutionary shaping of immunoglobulin gene loci

Background

The extraordinary diversity characterizing the antibody repertoire is generated by both evolution and lymphocyte development. Much of this diversity is due to the existence of immunoglobulin (Ig) variable region gene segment libraries, which were diversified during evolution and, in higher vertebrates, are used in generating the combinatorial diversity of antibody genes. The aim of the present study was to address the following questions: What evolutionary parameters affect the size and structure of gene libraries? Are the number of genes in libraries of contemporary species, and the corresponding gene locus structure, a random result of evolutionary history, or have these properties been optimized with respect to individual or population fitness? If a larger number of genes or different genome structures do not increase the fitness, then the current structure is probably optimized.

Results

We used a simulation of variable region gene library evolution. We measured the effect of different parameters on gene library size and diversity, and the corresponding fitness. We found compensating relationships between parameters, which optimized Ig library size and diversity.

Conclusions

We conclude that contemporary species' Ig libraries have been optimized by evolution in terms of Ig sequence lengths, the number and diversity of Ig genes, and antibody-antigen affinities.

Clustering-based identification of clonally-related immunoglobulin gene sequence sets

Background

Clonal expansion of B lymphocytes coupled with somatic mutation and antigen selection allow the mammalian humoral immune system to generate highly specific immunoglobulins (IG) or antibodies against invading bacteria, viruses and toxins. The availability of high-throughput DNA sequencing methods is providing new avenues for studying this clonal expansion and identifying the factors guiding the generation of antibodies. The identification of groups of rearranged immunoglobulin gene sequences descended from the same rearrangement (clonally-related sets) in very large sets of sequences is facilitated by the availability of immunoglobulin gene sequence alignment and partitioning software that can accurately predict component germline gene, but has required painstaking visual inspection and analysis of sequences.

Results

We have developed and implemented an algorithm for identifying sets of clonally-related sequences in large human immunoglobulin heavy chain gene variable region sequence sets. The program processes sequences that have been partitioned using iHMMune-align, and uses pairwise comparisons of CDR3 sequences and similarity in IGHV and IGHJ germline gene assignments to construct a distance matrix. Agglomerative hierarchical clustering is then used to identify likely groups of clonally-related sequences. The program is available for download from http://www.cse.unsw.edu.au/~ihmmune/ClonalRelate/ClonalRelate.zip.

Conclusions

The method was evaluated on several benchmark datasets and provided a more accurate and considerably faster identification of clonally-related immunoglobulin gene sequences than visual inspection by domain experts.

Immune responses are characterized by specific shared immunoglobulin peptides that can be detected by proteomic techniques

In the adaptive immune response, immunoglobulins develop that bind specifically to the antigens to which the organism was exposed. Immunoglobulins may bind to known or unknown antigens in a variety of diseases and have been used in the past to identify novel antigens for use as a biomarker. We propose that the immunoglobulins themselves could also be used as biomarkers in antibody-mediated disease. In this proteomic study, rats were immunized with one of two purified antigens, and immunoglobulins from pre- and postimmune sera were analyzed with nano-LC coupled mass spectrometry. It was found that the two treatment groups could be distinguished based on cluster analysis of the immunoglobulin peptides from the immune sera. In addition, we identified 684 specific peptides that were differentially present in one of the two treated groups. We could find an amino acid sequence for 44% of the features in the mass spectra by combining database-driven and de novo sequencing techniques. The latter were essential for sequence identification, as the more common database-driven approach suffers from a poor representation of immunoglobulins in the available databases. Our data show that the development of immunoglobulins during an immune response is not a fully random process, but that instead selection pressures exist that favor the best binding amino acid sequences, and that this selection is shared between different animals. This finding implies that immunoglobulin peptides could indeed be a powerful and easily accessible class of biomarkers.

Recent insights into the formation of RAG-induced chromosomal translocations

Chromosomal translocations are found in many types of tumors, where they may be either a cause or a result of malignant transformation. In lymphoid neoplasms, however, it is dear that pathogenesis is initiated by any of a number of recurrent DNA rearrangements. These particular translocations typically place an oncogene under the regulatory control of an Ig or TCR gene promoter, dysregulating cell growth, differentiation, or apoptosis. Given that physiological DNA rearrangements (V(D)J and class switch recombination) are integral to lymphocyte development, it is critical to understand how genomic stability is maintained during these processes. Recent advances in our understanding of DNA damage signaling and repair have provided clues to the kinds of mechanisms that lead to V(D)J-mediated translocations. In turn, investigations into the regulation of V(D)J joining have illuminated a formerly obscure pathway of DNA repair known as alternative NHEJ, which is error-prone and frequently involved in translocations. In this chapter we consider recent advances in our understanding of the functions of the RAG proteins, RAG interactions with DNA repair pathways, damage signaling and chromosome biology, all of which shed light on how mistakes at different stages of V(D)J recombination might lead to leukemias and lymphomas.