Novel Method for High-Throughput Full-Length IGHV-D-J Sequencing of the Immune Repertoire from Bulk B-Cells with Single-Cell Resolution

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.

Anchor Clustering for million-scale immune repertoire sequencing data

BACKGROUND

The clustering of immune repertoire data is challenging due to the computational cost associated with a very large number of pairwise sequence comparisons. To overcome this limitation, we developed Anchor Clustering, an unsupervised clustering method designed to identify similar sequences from millions of antigen receptor gene sequences. First, a Point Packing algorithm is used to identify a set of maximally spaced anchor sequences. Then, the genetic distance of the remaining sequences to all anchor sequences is calculated and transformed into distance vectors. Finally, distance vectors are clustered using unsupervised clustering. This process is repeated iteratively until the resulting clusters are small enough so that pairwise distance comparisons can be performed.

RESULTS

Our results demonstrate that Anchor Clustering is faster than existing pairwise comparison clustering methods while providing similar clustering quality. With its flexible, memory-saving strategy, Anchor Clustering is capable of clustering millions of antigen receptor gene sequences in just a few minutes.

CONCLUSIONS

This method enables the meta-analysis of immune-repertoire data from different studies and could contribute to a more comprehensive understanding of the immune repertoire data space.

CLL stereotyped B-cell receptor immunoglobulin sequences are recurrent in the B-cell repertoire of healthy individuals: Apparent lack of central and early peripheral tolerance censoring

Introduction

The leukemic cells of patients with chronic lymphocytic leukemia (CLL) are often unique, expressing remarkably similar IGHV-IGHD-IGHJ gene rearrangements, “stereotyped BCRs”. The B-cell receptors (BCRs) on CLL cells are also distinctive in often deriving from autoreactive B lymphocytes, leading to the assumption of a defect in immune tolerance.

Results

Using bulk and single-cell immunoglobulin heavy and light chain variable domain sequencing, we enumerated CLL stereotype-like IGHV-IGHD-IGHJ sequences (CLL-SLS) in B cells from cord blood (CB) and adult peripheral blood (PBMC) and bone marrow (BM of healthy donors. CLL-SLS were found at similar frequencies among CB, BM, and PBMC, suggesting that age does not influence CLL-SLS levels. Moreover, the frequencies of CLL-SLS did not differ among B lymphocytes in the BM at early stages of development, and only re-circulating marginal zone B cells contained significantly higher CLL-SLS frequencies than other mature B-cell subpopulations. Although we identified CLL-SLS corresponding to most of the CLL major stereotyped subsets, CLL-SLS frequencies did not correlate with those found in patients. Interestingly, in CB samples, half of the CLL-SLS identified were attributed to two IGHV-mutated subsets. We also found satellite CLL-SLS among the same normal samples, and they were also enriched in naïve B cells but unexpectedly, these were ~10-fold higher than standard CLL-SLS. In general, IGHV-mutated CLL-SLS subsets were enriched among antigen-experienced B-cell subpopulations, and IGHV-unmutated CLL-SLS were found mostly in antigen-inexperienced B cells. Nevertheless, CLL-SLS with an IGHV-mutation status matching that of CLL clones varied among the normal B-cell subpopulations, suggesting that specific CLL-SLS could originate from distinct subpopulations of normal B cells. Lastly, using single-cell DNA sequencing, we identified paired IGH and IGL rearrangements in normal B lymphocytes resembling those of stereotyped BCRs in CLL, although some differed from those in patients based on IG isotype or somatic mutation.

Discussion

CLL-SLS are present in normal B-lymphocyte populations at all stages of development. Thus, despite their autoreactive profile they are not deleted by central tolerance mechanisms, possibly because the level of autoreactivity is not registered as dangerous by deletion mechanisms or because editing of L-chain variable genes occurred which our experimental approach could not identify.

Nanopore sequencing of clonal IGH rearrangements in cell-free DNA as a biomarker for acute lymphoblastic leukemia

Background

Acute Lymphoblastic Leukemia (ALL) is the most common pediatric cancer, and patients with relapsed ALL have a poor prognosis. Detection of ALL blasts remaining at the end of treatment, or minimal residual disease (MRD), and spread of ALL into the central nervous system (CNS) have prognostic importance in ALL. Current methods to detect MRD and CNS disease in ALL rely on the presence of ALL blasts in patient samples. Cell-free DNA, or small fragments of DNA released by cancer cells into patient biofluids, has emerged as a robust and sensitive biomarker to assess cancer burden, although cfDNA analysis has not previously been applied to ALL.

Methods

We present a simple and rapid workflow based on NanoporeMinION sequencing of PCR amplified B cell-specific rearrangement of the (IGH) locus in cfDNA from B-ALL patient samples. A cohort of 5 pediatric B-ALL patient samples was chosen for the study based on the MRD and CNS disease status.

Results

Quantitation of IGH-variable sequences in cfDNA allowed us to detect clonal heterogeneity and track the response of individual B-ALL clones throughout treatment. cfDNA was detected in patient biofluids with clinical diagnoses of MRD and CNS disease, and leukemic clones could be detected even when diagnostic cell-count thresholds for MRD were not met. These data suggest that cfDNA assays may be useful in detecting the presence of ALL in the patient, even when blasts are not physically present in the biofluid sample.

Conclusions

The Nanopore IGH detection workflow to monitor cell-free DNA is a simple, rapid, and inexpensive assay that may ultimately serve as a valuable complement to traditional clinical diagnostic approaches for ALL.

Case Report: Sequential Development of Three Mature Lymphoid Neoplasms in a Single Patient: Clonal Relationship and Molecular Insights

Two main variants of Richter syndrome (RS) are recognized, namely, the diffuse large B-cell lymphoma (DLBCL) and the Hodgkin’s lymphoma (HL) variant. Clonal relationship, defined as an identity of the immunoglobulin heavy chain variable (IGHV) region sequence between chronic lymphocytic leukemia (CLL) and RS clones, characterizes patients with a poor prognosis. Due to method sensitivity, this categorization is performed without considering the possibility of small-size ancillary clones, sharing the same phenotype with the preexisting predominant CLL clone, but with different IGHV rearrangements. Here we describe and molecularly profile the peculiar case of a patient with a CLL-like monoclonal B-cell lymphocytosis (MBL), who sequentially developed a DLBCL, which occurred concomitantly to progression of MBL to CLL, and a subsequent HL. Based on standard IGHV clonality analysis, DLBCL was considered clonally unrelated to the concomitantly expanded CLL clone and treated as a de novo lymphoma, achieving a persistent response. Three years later, the patient further developed a clonally unrelated HL, refractory to bendamustine, which was successfully treated with brentuximab vedotin and radiotherapy, and later with pembrolizumab. We retrospectively performed additional molecular testing, by applying next-generation sequencing (NGS) of immunoglobulin repertoire (Ig-rep) techniques and a more sensitive allele-specific oligonucleotide-droplet digital PCR (ASO-ddPCR) strategy, in order to quantitatively investigate the presence of the rearranged IGHV genes in tumor specimens collected during the disease course. In this highly complex case, the application of modern and sensitive molecular technologies uncovered that DLBCL, initially considered as a de novo lymphoma, was instead the result of the transformation of a preexisting ancillary B-cell clone, which was already present at the time of first MBL diagnosis. A similar approach was also applied on the HL sample, showing its clonal unrelatedness to the previous MBL and DLBCL.

Characterizing Features of Human Circulating B Cells Carrying CLL-Like Stereotyped Immunoglobulin Rearrangements

Chronic Lymphocytic Leukemia (CLL) is characterized by the accumulation of monoclonal CD5⁺ B cells with low surface immunoglobulins (IG). About 40% of CLL clones utilize quasi-identical B cell receptors, defined as stereotyped BCR. CLL-like stereotyped-IG rearrangements are present in normal B cells as a part of the public IG repertoire. In this study, we collected details on the representation and features of CLL-like stereotyped-IG in the IGH repertoire of B-cell subpopulations purified from the peripheral blood of nine healthy donors. The B-cell subpopulations were also fractioned according to the expression of surface CD5 molecules and IG light chain, IGκ and IGλ. IG rearrangements, obtained by high throughput sequencing, were scanned for the presence of CLL-like stereotyped-IG. CLL-like stereotyped-IG did not accumulate preferentially in the CD5⁺ B cells, nor in specific B-cell subpopulations or the CD5⁺ cell fraction thereof, and their distribution was not restricted to a single IG light chain type. CLL-like stereotyped-IG shared with the corresponding CLL stereotype rearrangements the IGHV mutational status. Instead, for other features such as IGHV genes and frequency, CLL stereotyped-IGs presented a CLL-like subset specific behavior which could, or could not, be consistent with CLL stereotyped-IGs. Therefore, as opposed to the immuno-phenotype, the features of the CLL stereotyped-IG repertoire suggest a CLL stereotyped subset-specific ontogeny. Overall, these findings suggest that the immune-genotype can provide essential details in tracking and defining the CLL cell of origin.

Antibody upstream sequence diversity and its biological implications revealed by repertoire sequencing

The sequence upstream of the antibody variable region (antibody upstream sequence [AUS]) consists of a 5' untranslated region (5' UTR) and a preceding leader region. The sequence variations in AUS affect antibody engineering and PCR based antibody quantification and may also be implicated in mRNA transcription and translation. However, the diversity of AUSs remains elusive. Using 5' rapid amplification of cDNA ends and high-throughput antibody repertoire sequencing technique, we acquired full-length AUSs for human, rhesus macaque, cynomolgus macaque, mouse, and rat. We designed a bioinformatics pipeline and identified 3307 unique AUSs, corresponding to 3026 and 1457 unique sequences for 5' UTR and leader region, respectively. Comparative analysis indicated that 928 (63.69%) leader sequences are novel relative to those recorded in the international ImMunoGeneTics information system. Evolutionarily, leader sequences are more conserved than 5' UTR and seem to coevolve with their downstream V genes. Besides, single-nucleotide polymorphisms are position dependent for leader regions and may contribute to the functional reversal of the downstream V genes. Finally, the AUGs in AUSs were found to have little impact on gene expression. Taken together, our findings can facilitate primer design for capturing antibodies efficiently and provide a valuable resource for antibody engineering and molecule-level antibody studies.

A Bayesian model based computational analysis of the relationship between bisulfite accessible single-stranded DNA in chromatin and somatic hypermutation of immunoglobulin genes

The B cells in our body generate protective antibodies by introducing somatic hypermutations (SHM) into the variable region of immunoglobulin genes (IgVs). The mutations are generated by activation induced deaminase (AID) that converts cytosine to uracil in single stranded DNA (ssDNA) generated during transcription. Attempts have been made to correlate SHM with ssDNA using bisulfite to chemically convert cytosines that are accessible in the intact chromatin of mutating B cells. These studies have been complicated by using different definitions of "bisulfite accessible regions" (BARs). Recently, deep-sequencing has provided much larger datasets of such regions but computational methods are needed to enable this analysis. Here we leveraged the deep-sequencing approach with unique molecular identifiers and developed a novel Hidden Markov Model based Bayesian Segmentation algorithm to characterize the ssDNA regions in the IGHV4-34 gene of the human Ramos B cell line. Combining hierarchical clustering and our new Bayesian model, we identified recurrent BARs in certain subregions of both top and bottom strands of this gene. Using this new system, the average size of BARs is about 15 bp. We also identified potential G-quadruplex DNA structures in this gene and found that the BARs co-locate with G-quadruplex structures in the opposite strand. Using various correlation analyses, there is not a direct site-to-site relationship between the bisulfite accessible ssDNA and all sites of SHM but most of the highly AID mutated sites are within 15 bp of a BAR. In summary, we developed a novel platform to study single stranded DNA in chromatin at a base pair resolution that reveals potential relationships among BARs, SHM and G-quadruplexes. This platform could be applied to genome wide studies in the future.

Role of Dot1L and H3K79 methylation in regulating somatic hypermutation of immunoglobulin genes

Somatic hypermutation (SHM) and class-switch recombination (CSR) of the immunoglobulin (Ig) genes allow B cells to make antibodies that protect us against a wide variety of pathogens. SHM is mediated by activation-induced deaminase (AID), occurs at a million times higher frequency than other mutations in the mammalian genome, and is largely restricted to the variable (V) and switch (S) regions of Ig genes. Using the Ramos human Burkitt's lymphoma cell line, we find that H3K79me2/3 and its methyltransferase Dot1L are more abundant on the V region than on the constant (C) region, which does not undergo mutation. In primary naïve mouse B cells examined ex vivo, the H3K79me2/3 modification appears constitutively in the donor Sμ and is inducible in the recipient Sγ1 upon CSR stimulation. Knockout and inhibition of Dot1L in Ramos cells significantly reduces V region mutation and the abundance of H3K79me2/3 on the V region and is associated with a decrease of polymerase II (Pol II) and its S2 phosphorylated form at the IgH locus. Knockout of Dot1L also decreases the abundance of BRD4 and CDK9 (a subunit of the P-TEFb complex) on the V region, and this is accompanied by decreased nascent transcripts throughout the IgH gene. Treatment with JQ1 (inhibitor of BRD4) or DRB (inhibitor of CDK9) decreases SHM and the abundance of Pol II S2P at the IgH locus. Since all these factors play a role in transcription elongation, our studies reinforce the idea that the chromatin context and dynamics of transcription are critical for SHM.

Rituximab-resistant splenic memory B cells and newly engaged naive B cells fuel relapses in patients with immune thrombocytopenia

Rituximab (RTX), an antibody targeting CD20, is widely used as a first-line therapeutic strategy in B cell-mediated autoimmune diseases. However, a large proportion of patients either do not respond to the treatment or relapse during B cell reconstitution. Here, we characterize the cellular basis responsible for disease relapse in secondary lymphoid organs in humans, taking advantage of the opportunity offered by therapeutic splenectomy in patients with relapsing immune thrombocytopenia. By analyzing the B and plasma cell immunoglobulin gene repertoire at bulk and antigen-specific single-cell level, we demonstrate that relapses are associated with two responses coexisting in germinal centers and involving preexisting mutated memory B cells that survived RTX treatment and naive B cells generated upon reconstitution of the B cell compartment. To identify distinctive characteristics of the memory B cells that escaped RTX-mediated depletion, we analyzed RTX refractory patients who did not respond to treatment at the time of B cell depletion. We identified, by single-cell RNA sequencing (scRNA-seq) analysis, a population of quiescent splenic memory B cells that present a unique, yet reversible, RTX-shaped phenotype characterized by down-modulation of B cell-specific factors and expression of prosurvival genes. Our results clearly demonstrate that these RTX-resistant autoreactive memory B cells reactivate as RTX is cleared and give rise to plasma cells and further germinal center reactions. Their continued surface expression of CD19 makes them efficient targets for current anti-CD19 therapies. This study thus identifies a pathogenic contributor to autoimmune diseases that can be targeted by available therapeutic agents.

Post-Transformation IGHV-IGHD-IGHJ Mutations in Chronic Lymphocytic Leukemia B Cells: Implications for Mutational Mechanisms and Impact on Clinical Course

Analyses of IGHV gene mutations in chronic lymphocytic leukemia (CLL) have had a major impact on the prognostication and treatment of this disease. A hallmark of IGHV-mutation status is that it very rarely changes clonally over time. Nevertheless, targeted and deep DNA sequencing of IGHV-IGHD-IGHJ regions has revealed intraclonal heterogeneity. We used a DNA sequencing approach that achieves considerable depth and minimizes artefacts and amplification bias to identify IGHV-IGHD-IGHJ subclones in patients with prolonged temporal follow-up. Our findings extend previous studies, revealing intraclonal IGHV-IGHD-IGHJ diversification in almost all CLL clones. Also, they indicate that some subclones with additional IGHV-IGHD-IGHJ mutations can become a large fraction of the leukemic burden, reaching numerical criteria for monoclonal B-cell lymphocytosis. Notably, the occurrence and complexity of post-transformation IGHV-IGHD-IGHJ heterogeneity and the expansion of diversified subclones are similar among U-CLL and M-CLL patients. The molecular characteristics of the mutations present in the parental, clinically dominant CLL clone (CDC) differed from those developing post-transformation (post-CDC). Post-CDC mutations exhibit significantly lower fractions of mutations bearing signatures of activation induced deaminase (AID) and of error-prone repair by Polη, and most of the mutations were not ascribable to those enzymes. Additionally, post-CDC mutations displayed a lower percentage of nucleotide transitions compared with transversions that was also not like the action of AID. Finally, the post-CDC mutations led to significantly lower ratios of replacement to silent mutations in VH CDRs and higher ratios in VH FRs, distributions different from mutations found in normal B-cell subsets undergoing an AID-mediated process. Based on these findings, we propose that post-transformation mutations in CLL cells either reflect a dysfunctional standard somatic mutational process or point to the action of another mutational process not previously associated with IG V gene loci. If the former option is the case, post-CDC mutations could lead to a lesser dependence on antigen dependent BCR signaling and potentially a greater influence of off-target, non-IG genomic mutations. Alternatively, the latter activity could add a new stimulatory survival/growth advantage mediated by the BCR through structurally altered FRs, such as that occurring by superantigen binding and stimulation.

Large-scale analysis of 2,152 Ig-seq datasets reveals key features of B cell biology and the antibody repertoire

Antibody repertoire sequencing enables researchers to acquire millions of B cell receptors and investigate these molecules at the single-nucleotide level. This power and resolution in studying humoral responses have led to its wide applications. However, most of these studies were conducted with a limited number of samples. Given the extraordinary diversity, assessment of these key features with a large sample set is demanded. Thus, we collect and systematically analyze 2,152 high-quality heavy-chain antibody repertoires. Our study reveals that 52 core variable genes universally contribute to more than 99% of each individual's repertoire; a distal interspersed preferences characterize V gene recombination; the number of public clones between two repertoires follows a linear model, and the positive selection dominates at RGYW motif in somatic hypermutations. Thus, this population-level analysis resolves some critical features of the antibody repertoire and may have significant value to the large cadre of scientists.

Correlations in Somatic Hypermutation Between Sites in IGHV Genes Can Be Explained by Interactions Between AID and/or Polη Hotspots

The somatic hypermutation (SHM) of Immunoglobulin (Ig) genes is a key process during antibody affinity maturation in B cells. The mutagenic enzyme activation induced deaminase (AID) is required for SHM and has a preference for WRC hotspots in DNA. Error-prone repair mechanisms acting downstream of AID introduce further mutations, including DNA polymerase eta (Polη), part of the non-canonical mismatch repair pathway (ncMMR), which preferentially generates mutations at WA hotspots. Previously proposed mechanistic models lead to a variety of predictions concerning interactions between hotspots, for example, how mutations in one hotspot will affect another hotspot. Using a large, high-quality, Ig repertoire sequencing dataset, we evaluated pairwise correlations between mutations site-by-site using an unbiased measure similar to mutual information which we termed “mutational association” (MA). Interactions are dominated by relatively strong correlations between nearby sites (short-range MAs), which can be almost entirely explained by interactions between overlapping hotspots for AID and/or Polη. We also found relatively weak dependencies between almost all sites throughout each gene (longer-range MAs), although these arise mostly as a statistical consequence of high pairwise mutation frequencies. The dominant short-range interactions are also highest within the most highly mutating IGHV sub-regions, such as the complementarity determining regions (CDRs), where there is a high hotspot density. Our results suggest that the hotspot preferences for AID and Polη have themselves evolved to allow for greater interactions between AID and/or Polη induced mutations.

Polymorphisms in human immunoglobulin heavy chain variable genes and their upstream regions

Germline variations in immunoglobulin genes influence the repertoire of B cell receptors and antibodies, and such polymorphisms may impact disease susceptibility. However, the knowledge of the genomic variation of the immunoglobulin loci is scarce. Here, we report 25 potential novel germline IGHV alleles as inferred from rearranged naïve B cell cDNA repertoires of 98 individuals. Thirteen novel alleles were selected for validation, out of which ten were successfully confirmed by targeted amplification and Sanger sequencing of non-B cell DNA. Moreover, we detected a high degree of variability upstream of the V-REGION in the 5′UTR, L-PART1 and L-PART2 sequences, and found that identical V-REGION alleles can differ in upstream sequences. Thus, we have identified a large genetic variation not only in the V-REGION but also in the upstream sequences of IGHV genes. Our findings provide a new perspective for annotating immunoglobulin repertoire sequencing data.

Network Organization of Antibody Interactions in Sequence and Structure Space: the RADARS Model

Adaptive immunity in vertebrates is a complex self-organizing network of molecular interactions. While deep sequencing of the immune-receptor repertoire may reveal clonal relationships, functional interpretation of such data is hampered by the inherent limitations of converting sequence to structure to function. In this paper, a novel model of antibody interaction space and network, termed radial adjustment of system resolution, RAdial ADjustment of System Resolution (RADARS), is proposed. The model is based on the radial growth of interaction affinity of antibodies towards an infinity of directions in structure space, each direction corresponding to particular shapes of antigen epitopes. Levels of interaction affinity appear as free energy shells of the system, where hierarchical B-cell development and differentiation takes place. Equilibrium in this immunological thermodynamic system can be described by a power law distribution of antibody-free energies with an ideal network degree exponent of phi square, representing a scale-free fractal network of antibody interactions. Plasma cells are network hubs, memory B cells are nodes with intermediate degrees, and B1 cells function as nodes with minimal degree. Overall, the RADARS model implies that a finite number of antibody structures can interact with an infinite number of antigens by immunologically controlled adjustment of interaction energy distribution. Understanding quantitative network properties of the system should help the organization of sequence-derived predicted structural data.

AID Overlapping and Polη Hotspots Are Key Features of Evolutionary Variation Within the Human Antibody Heavy Chain (IGHV) Genes

Somatic hypermutation (SHM) of the immunoglobulin variable (IgV) loci is a key process in antibody affinity maturation. The enzyme activation-induced deaminase (AID), initiates SHM by creating C → U mismatches on single-stranded DNA (ssDNA). AID has preferential hotspot motif targets in the context of WRC/GYW (W = A/T, R = A/G, Y = C/T) and particularly at WGCW overlapping hotspots where hotspots appear opposite each other on both strands. Subsequent recruitment of the low-fidelity DNA repair enzyme, Polymerase eta (Polη), during mismatch repair, creates additional mutations at WA/TW sites. Although there are more than 50 functional immunoglobulin heavy chain variable (IGHV) segments in humans, the fundamental differences between these genes and their ability to respond to all possible foreign antigens is still poorly understood. To better understand this, we generated profiles of WGCW hotspots in each of the human IGHV genes and found the expected high frequency in complementarity determining regions (CDRs) that encode the antigen binding sites but also an unexpectedly high frequency of WGCW in certain framework (FW) sub-regions. Principal Components Analysis (PCA) of these overlapping AID hotspot profiles revealed that one major difference between IGHV families is the presence or absence of WGCW in a sub-region of FW3 sometimes referred to as “CDR4.” Further differences between members of each family (e.g., IGHV1) are primarily determined by their WGCW densities in CDR1. We previously suggested that the co-localization of AID overlapping and Polη hotspots was associated with high mutability of certain IGHV sub-regions, such as the CDRs. To evaluate the importance of this feature, we extended the WGCW profiles, combining them with local densities of Polη (WA) hotspots, thus describing the co-localization of both types of hotspots across all IGHV genes. We also verified that co-localization is associated with higher mutability. PCA of the co-localization profiles showed CDR1 and CDR2 as being the main contributors to variance among IGHV genes, consistent with the importance of these sub-regions in antigen binding. Our results suggest that AID overlapping (WGCW) hotspots alone or in conjunction with Polη (WA/TW) hotspots are key features of evolutionary variation between IGHV genes.

Additional B-cell malignancies in patients with chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL)

Family and migration studies suggest a genetic risk of developing chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL). We hypothesized that CLL patients have an increased risk of additional clonally unrelated B-cell malignancies. To test this, we studied 467 CLL patients (2743 person-years (PYs)) at a single institution over 17 years. The incidence rate (IR) of any additional B-cell lymphoid malignancy was 10.9 per 1000 PYs (n = 30, 6.4%). Eighteen (4%) patients had a clonally unrelated B-cell malignancy (IR = 6.6 per 1000 PYs). Standardized incidence ratios (SIRs) were used to compare the incidence of additional clonally unrelated B-cell malignancies in CLL patients to the age- and sex-matched expected rates in the USA generated from the Surveillance, Epidemiology, and End Results (SEER) database. For the subset of 13 patients having data for comparison in the SEER database, the SIR was 5.41 (95% CI = 2.9, 9.3) which is supportive of our hypothesis.

Tracing CLL-biased stereotyped immunoglobulin gene rearrangements in normal B cell subsets using a high-throughput immunogenetic approach

Background

B cell receptor Immunoglobulin (BcR IG) repertoire of Chronic Lymphocytic Leukemia (CLL) is characterized by the expression of quasi-identical BcR IG. These are observed in approximately 30% of patients, defined as stereotyped receptors and subdivided into subsets based on specific VH CDR3 aa motifs and phylogenetically related IGHV genes. Although relevant to CLL ontogeny, the distribution of CLL-biased stereotyped immunoglobulin rearrangements (CBS-IG) in normal B cells has not been so far specifically addressed using modern sequencing technologies. Here, we have investigated the presence of CBS-IG in splenic B cell subpopulations (s-BCS) and in CD5⁺ and CD5⁻ B cells from the spleen and peripheral blood (PB).

Methods

Fractionation of splenic B cells into 9 different B cell subsets and that of spleen and PB into CD5⁺ and CD5⁻ cells were carried out by FACS sorting. cDNA sequences of BcR IG gene rearrangements were obtained by NGS. Identification of amino acidic motifs typical of CLL stereotyped subsets was carried out on IGHV1-carrying gene sequences and statistical evaluation has been subsequently performed to assess stereotypes distribution.

Results

CBS-IG represented the 0.26% average of IGHV1 genes expressing sequences, were detected in all of the BCS investigated. CBS-IG were more abundant in splenic and circulating CD5⁺ B (0.57%) cells compared to CD5⁻ B cells (0.17%). In all instances, most CBS IG did not exhibit somatic hypermutation similar to CLL stereotyped receptors. However, compared to CLL, they exhibited a different CLL subset distribution and a broader utilization of the genes of the IGHV1 family.

Conclusions

CBS-IG receptors appear to represent a part of the “public” BcR repertoire in normal B cells. This repertoire is observed in all BCS excluding the hypothesis that CLL stereotyped BcR accumulate in a specific B cell subset, potentially capable of originating a leukemic clone. The different relative representation of CBS-IG in normal B cell subgroups suggests the requirement for additional selective processes before a full transformation into CLL is achieved.

PD-1/PD-L1 immune checkpoint and p53 loss facilitate tumor progression in activated B-cell diffuse large B-cell lymphomas

Refractory or relapsed diffuse large B-cell lymphoma (DLBCL) often associates with the activated B-cell-like (ABC) subtype and genetic alterations that drive constitutive NF-κB activation and impair B-cell terminal differentiation. Here, we show that DNA damage response by p53 is a central mechanism suppressing the pathogenic cooperation of IKK2ca-enforced canonical NF-κB and impaired differentiation resulting from Blimp1 loss in ABC-DLBCL lymphomagenesis. We provide evidences that the interplay between these genetic alterations and the tumor microenvironment select for additional molecular addictions that promote lymphoma progression, including aberrant coexpression of FOXP1 and the B-cell mutagenic enzyme activation-induced deaminase, and immune evasion through major histocompatibility complex class II downregulation, PD-L1 upregulation, and T-cell exhaustion. Consistently, PD-1 blockade cooperated with anti-CD20-mediated B-cell cytotoxicity, promoting extended T-cell reactivation and antitumor specificity that improved long-term overall survival in mice. Our data support a pathogenic cooperation among NF-κB-driven prosurvival, genetic instability, and immune evasion mechanisms in DLBCL and provide preclinical proof of concept for including PD-1/PD-L1 blockade in combinatorial immunotherapy for ABC-DLBCL.

High-Quality Library Preparation for NGS-Based Immunoglobulin Germline Gene Inference and Repertoire Expression Analysis

Next generation sequencing (NGS) of immunoglobulin (Ig) repertoires (Rep-seq) enables examination of the adaptive immune system at an unprecedented level. Applications include studies of expressed repertoires, gene usage, somatic hypermutation levels, Ig lineage tracing and identification of genetic variation within the Ig loci through inference methods. All these applications require starting libraries that allow the generation of sequence data with low error rate and optimal representation of the expressed repertoire. Here, we provide detailed protocols for the production of libraries suitable for human Ig germline gene inference and Ig repertoire studies. Various parameters used in the process were tested in order to demonstrate factors that are critical to obtain high quality libraries. We demonstrate an improved 5'RACE technique that reduces the length constraints of Illumina MiSeq based Rep-seq analysis but allows for the acquisition of sequences upstream of Ig V genes, useful for primer design. We then describe a 5' multiplex method for library preparation, which yields full length V(D)J sequences suitable for genotype identification and novel gene inference. We provide comprehensive sets of primers targeting IGHV, IGKV, and IGLV genes. Using the optimized protocol, we produced IgM, IgG, IgK, and IgL libraries and analyzed them using the germline inference tool IgDiscover to identify expressed germline V alleles. This process additionally uncovered three IGHV, one IGKV, and six IGLV novel alleles in a single individual, which are absent from the IMGT reference database, highlighting the need for further study of Ig genetic variation. The library generation protocols presented here enable a robust means of analyzing expressed Ig repertoires, identifying novel alleles and producing individualized germline gene databases from humans.

Inferred Allelic Variants of Immunoglobulin Receptor Genes: A System for Their Evaluation, Documentation, and Naming

Immunoglobulins or antibodies are the main effector molecules of the B-cell lineage and are encoded by hundreds of variable (V), diversity (D), and joining (J) germline genes, which recombine to generate enormous IG diversity. Recently, high-throughput adaptive immune receptor repertoire sequencing (AIRR-seq) of recombined V-(D)-J genes has offered unprecedented insights into the dynamics of IG repertoires in health and disease. Faithful biological interpretation of AIRR-seq studies depends upon the annotation of raw AIRR-seq data, using reference germline gene databases to identify the germline genes within each rearrangement. Existing reference databases are incomplete, as shown by recent AIRR-seq studies that have inferred the existence of many previously unreported polymorphisms. Completing the documentation of genetic variation in germline gene databases is therefore of crucial importance. Lymphocyte receptor genes and alleles are currently assigned by the Immunoglobulins, T cell Receptors and Major Histocompatibility Nomenclature Subcommittee of the International Union of Immunological Societies (IUIS) and managed in IMGT®, the international ImMunoGeneTics information system® (IMGT). In 2017, the IMGT Group reached agreement with a group of AIRR-seq researchers on the principles of a streamlined process for identifying and naming inferred allelic sequences, for their incorporation into IMGT®. These researchers represented the AIRR Community, a network of over 300 researchers whose objective is to promote all aspects of immunoglobulin and T-cell receptor repertoire studies, including the standardization of experimental and computational aspects of AIRR-seq data generation and analysis. The Inferred Allele Review Committee (IARC) was established by the AIRR Community to devise policies, criteria, and procedures to perform this function. Formalized evaluations of novel inferred sequences have now begun and submissions are invited via a new dedicated portal (https://ogrdb.airr-community.org). Here, we summarize recommendations developed by the IARC-focusing, to begin with, on human IGHV genes-with the goal of facilitating the acceptance of inferred allelic variants of germline IGHV genes. We believe that this initiative will improve the quality of AIRR-seq studies by facilitating the description of human IG germline gene variation, and that in time, it will expand to the documentation of TR and IG genes in many vertebrate species.

Single-Cell Genomics

BACKGROUND

Single-cell genomics is an approach to investigate cell heterogeneity and to identify new molecular features correlated with clinical outcomes. This approach allows identification of the complexity of cell diversity in a sample without the loss of information that occurs when multicellular or bulk tissue samples are analyzed.

CONTENT

The first single-cell RNA-sequencing study was published in 2009, and since then many more studies and single-cell sequencing methods have been published. These studies have had a major impact on several fields, including microbiology, neurobiology, cancer, and developmental biology. Recently, improvements in reliability and the development of commercial single-cell isolation platforms are opening the potential of this technology to the clinical laboratory.

SUMMARY

In this review we provide an overview of the current state of single-cell genomics. We describe opportunities in clinical research and medical applications.

Observed Antibody Space: A Resource for Data Mining Next-Generation Sequencing of Antibody Repertoires

Abs are immune system proteins that recognize noxious molecules for elimination. Their sequence diversity and binding versatility have made Abs the primary class of biopharmaceuticals. Recently, it has become possible to query their immense natural diversity using next-generation sequencing of Ig gene repertoires (Ig-seq). However, Ig-seq outputs are currently fragmented across repositories and tend to be presented as raw nucleotide reads, which means nontrivial effort is required to reuse the data for analysis. To address this issue, we have collected Ig-seq outputs from 55 studies, covering more than half a billion Ab sequences across diverse immune states, organisms (primarily human and mouse), and individuals. We have sorted, cleaned, annotated, translated, and numbered these sequences and make the data available via our Observed Antibody Space (OAS) resource at http://antibodymap.org The data within OAS will be regularly updated with newly released Ig-seq datasets. We believe OAS will facilitate data mining of immune repertoires for improved understanding of the immune system and development of better biotherapeutics.

Computational Evaluation of B-Cell Clone Sizes in Bulk Populations

B cell clones expand and contract during adaptive immune responses and can persist or grow uncontrollably in lymphoproliferative disorders. One way to monitor and track B cell clones is to perform large-scale sampling of bulk cell populations, amplifying, and sequencing antibody gene rearrangements by next-generation sequencing (NGS). Here, we describe a series of computational approaches for estimating B cell clone size in NGS immune repertoire profiling data of antibody heavy chain gene rearrangements. We define three different measures of B cell clone size-copy numbers, instances, and unique sequences-and show how these measures can be used to rank clones, analyze their diversity, and study their distribution within and between individuals. We provide a detailed, step-by-step procedure for performing these analyses using two different data sets of spleen samples from human organ donors. In the first data set, 19 independently generated biological replicates from a single individual are analyzed for B cell clone size, diversity and sampling sufficiency for clonal overlap analysis. In the second data set, B cell clones are compared in eight different organ donors. We comment upon frequently encountered pitfalls and offer practical advice with alternative approaches. Overall, we provide a series of pragmatic analytical approaches and show how different clone size measures can be used to study the clonal landscape in bulk B cell immune repertoire profiling data.

Evaluation of Antigen-Conjugated Fluorescent Beads to Identify Antigen-Specific B Cells

Selection of single antigen-specific B cells to identify their expressed antibodies is of considerable interest for evaluating human immune responses. Here, we present a method to identify single antibody-expressing cells using antigen-conjugated fluorescent beads. To establish this, we selected Folate Receptor alpha (FRα) as a model antigen and a mouse B cell line, expressing both the soluble and the membrane-bound forms of a human/mouse chimeric antibody (MOv18 IgG1) specific for FRα, as test antibody-expressing cells. Beads were conjugated to FRα using streptavidin/avidin-biotin bridges and used to select single cells expressing the membrane-bound form of anti-FRα. Bead-bound cells were single cell-sorted and processed for single cell RNA retrotranscription and PCR to isolate antibody heavy and light chain variable regions. Variable regions were then cloned and expressed as human IgG1/k antibodies. Like the original clone, engineered antibodies from single cells recognized native FRα. To evaluate whether antigen-coated beads could identify specific antibody-expressing cells in mixed immune cell populations, human peripheral blood mononuclear cells (PBMCs) were spiked with test antibody-expressing cells. Antigen-specific cells could comprise up to 75% of cells selected with antigen-conjugated beads when the frequency of the antigen-positive cells was 1:100 or higher. In PBMC pools, beads conjugated to recombinant antigens FRα and HER2 bound antigen-specific anti-FRα MOv18 and anti-HER2 Trastuzumab antibody-expressing cells, respectively. From melanoma patient-derived B cells selected with melanoma cell line-derived protein-coated fluorescent beads, we generated a monoclonal antibody that recognized melanoma antigen-coated beads. This approach may be further developed to facilitate analysis of B cells and their antibody profiles at the single cell level and to help unravel humoral immune repertoires.