Antigen Receptor Galaxy: A User-Friendly, Web-Based Tool for Analysis and Visualization of T and B Cell Receptor Repertoire Data

IL-7 receptor signaling drives human B-cell progenitor differentiation and expansion

Although absence of interleukin-7 (IL-7) signaling completely abrogates T and B lymphopoiesis in mice, patients with severe combined immunodeficiency caused by mutations in the IL-7 receptor α chain (IL-7Rα) still generate peripheral blood B cells. Consequently, human B lymphopoiesis has been thought to be independent of IL-7 signaling. Using flow cytometric analysis and single-cell RNA sequencing of bone marrow samples from healthy controls and patients who are IL-7Rα deficient, in combination with in vitro modeling of human B-cell differentiation, we demonstrate that IL-7R signaling plays a crucial role in human B lymphopoiesis. IL-7 drives proliferation and expansion of early B-cell progenitors but not of pre-BII large cells and has a limited role in the prevention of cell death. Furthermore, IL-7 guides cell fate decisions by enhancing the expression of BACH2, EBF1, and PAX5, which jointly orchestrate the specification and commitment of early B-cell progenitors. In line with this observation, early B-cell progenitors of patients with IL-7Rα deficiency still expressed myeloid-specific genes. Collectively, our results unveil a previously unknown role for IL-7 signaling in promoting the B-lymphoid fate and expanding early human B-cell progenitors while defining important differences between mice and humans. Our results have implications for hematopoietic stem cell transplantation strategies in patients with T- B+ severe combined immunodeficiency and provide insights into the role of IL-7R signaling in leukemogenesis.

Human IgMhiCD300a+ B Cells Are Circulating Marginal Zone Memory B Cells That Respond to Pneumococcal Polysaccharides and Their Frequency Is Decreased in People Living with HIV

CD300a is differentially expressed among B cell subsets, although its expression in immunoglobulin (Ig)M+ B cells is not well known. We identified a B cell subset expressing CD300a and high levels of IgM (IgMhiCD300a+). The results showed that IgMhiCD300a+ B cells were CD10-CD27+CD25+IgDloCD21hiCD23-CD38loCD1chi, suggesting that they are circulating marginal zone (MZ) IgM memory B cells. Regarding the immunoglobulin repertoire, IgMhiCD300a+ B cells exhibited a higher mutation rate and usage of the IgH-VDJ genes than the IgM+CD300a- counterpart. Moreover, the shorter complementarity-determining region 3 (CDR3) amino acid (AA) length from IgMhiCD300a+ B cells together with the predicted antigen experience repertoire indicates that this B cell subset has a memory phenotype. IgM memory B cells are important in T cell-independent responses. Accordingly, we demonstrate that this particular subset secretes higher amounts of IgM after stimulation with pneumococcal polysaccharides or a toll-like receptor 9 (TLR9) agonist than IgM+CD300a- cells. Finally, the frequency of IgMhiCD300a+ B cells was lower in people living with HIV-1 (PLWH) and it was inversely correlated with the years with HIV infection. Altogether, these data help to identify a memory B cell subset that contributes to T cell-independent responses to pneumococcal infections and may explain the increase in severe pneumococcal infections and the impaired responses to pneumococcal vaccination in PLWH.

IgD shapes the pre-immune naïve B cell compartment in humans

B cell maturation and immunoglobulin (Ig) repertoire selection are governed by expression of a functional B cell receptor (BCR). Naïve B cells co-express their BCR as IgM and IgD isotype. However, the role of the additionally expressed IgD on naïve B cells is not known. Here we assessed the impact of IgD on naïve B cell maturation and Ig repertoire selection in 8 individuals from 3 different families with heterozygous loss-of-function or loss-of expression mutations in IGHD. Although naïve B cells from these individuals expressed IgM on their surface, the IGHD variant in heterozygous state entailed a chimeric situation by allelic exclusion with almost half of the naïve B cell population lacking surface IgD expression. Flow cytometric analyses revealed a distinct phenotype of IgD-negative naïve B cells with decreased expression of CD19, CD20 and CD21 as well as lower BAFF-R and integrin-β7 expression. IgD-negative B cells were less responsive in vitro after engaging the IgM-BCR, TLR7/9 or CD40 pathway. Additionally, a selective disadvantage of IgD-negative B cells within the T2 transitional and mature naïve B cell compartment as well as reduced frequencies of IgM^lo/- B cells within the mature naïve B cell compartment lacking IgD were evident. RNA-Ig-seq of bulk sorted B cell populations showed an altered selection of distinct V_H segments in the IgD-negative mature naïve B cell population. We conclude that IgD expression on human naïve B cells is redundant for generation of naïve B cells in general, but further shapes the naive B cell compartment starting from T2 transitional B cells. Our observations suggest an unexpected role of IgD expression to be critical for selection of distinct Ig V_H segments into the pre-immune Ig repertoire and for the survival of IgM^lo/- naïve B cells known to be enriched in poly-/autoreactive B cell clones.

Age-dependent nasal immune responses in non-hospitalized bronchiolitis children

Bronchiolitis in children is associated with significant rates of morbidity and mortality. Many studies have been performed using samples from hospitalized bronchiolitis patients, but little is known about the immunological responses from infants suffering from mild/moderate bronchiolitis that do not require hospitalization. We have studied a collection of nasal lavage fluid (NLF) samples from outpatient bronchiolitis children as a novel strategy to unravel local humoral and cellular responses, which are not fully characterized. The children were age-stratified in three groups, two of them (GI under 2-months, GII between 2-4 months) presenting a first episode of bronchiolitis, and GIII (between 4 months and 2 years) with recurrent respiratory infections. Here we show that elevated levels of pro-inflammatory cytokines (IL1β, IL6, TNFα, IL18, IL23), regulatory cytokines (IL10, IL17A) and IFNγ were found in the three bronchiolitis cohorts. However, little or no change was observed for IL33 and MCP1, at difference to previous results from bronchiolitis hospitalized patients. Furthermore, our results show a tendency to IL1β, IL6, IL18 and TNFα increased levels in children with mild pattern of symptom severity and in those in which non RSV respiratory virus were detected compared to RSV+ samples. By contrast, no such differences were found based on gender distribution. Bronchiolitis NLFs contained more IgM, IgG1, IgG3 IgG4 and IgA than NLF from their age-matched healthy controls. NLF from bronchiolitis children predominantly contained neutrophils, and also low frequency of monocytes and few CD4⁺ and CD8⁺ T cells. NLF from infants older than 4-months contained more intermediate monocytes and B cell subsets, including naïve and memory cells. BCR repertoire analysis of NLF samples showed a biased VH1 usage in IgM repertoires, with low levels of somatic hypermutation. Strikingly, algorithmic studies of the mutation profiles, denoted antigenic selection on IgA-NLF repertoires. Our results support the use of NLF samples to analyze immune responses and may have therapeutic implications.

B-Cell Receptor Repertoire Sequencing: Deeper Digging into the Mechanisms and Clinical Aspects of Immune-mediated Diseases

B cells play an essential role in adaptive immunity and are intimately correlated with pleiotropic immune-mediated diseases. Each B cell occupies a unique B cell receptor (BCR), and all BCRs throughout our body form “BCR repertoire.” With the development of sequencing technology and coupled bioinformatics, accumulating evidence indicates that BCR repertoire largely varies under physiological and pathological conditions. Therefore, comprehensive grasp of BCR repertoire will provide new insights into the pathogenesis of immune-mediated diseases and help exploit efficient diagnostic and treatment strategies. In this review, we start with an overview of BCR repertoire and related sequencing technologies and summarize their current applications in immune-mediated diseases. We also underscore the challenges of this emerging field and propose promising future directions in advancing BCR repertoire exploration.

NGS-Based B-Cell Receptor Repertoire AnalysisRepertoire analyses in the Context of Inborn Errors of Immunity

Inborn errors of immunity (IEI) are genetic defects that can affect both the innate and the adaptive immune system. Patients with IEI usually present with recurrent infections, but many also suffer from immune dysregulation, autoimmunity, and malignancies.Inborn errors of the immune system can cause defects in the development and selection of the B-cell receptor (BCR ) repertoire. Patients with IEI can have a defect in one of the key processes of immune repertoire formation like V(D)J recombination, somatic hypermutation (SHM), class switch recombination (CSR), or (pre-)BCR signalling and proliferation. However, also other genetic defects can lead to quantitative and qualitative differences in the immune repertoire.In this chapter, we will give an overview of protocols that can be used to study the immune repertoire in patients with IEI, provide considerations to take into account before setting up experiments, and discuss analysis of the immune repertoire data using Antigen Receptor Galaxy (ARGalaxy).

A Novel AICDA Splice-Site Mutation in Two Siblings with HIGM2 Permits Somatic Hypermutation but Abrogates Mutational Targeting

Hyper-IgM syndrome type 2 (HIGM2) is a B cell intrinsic primary immunodeficiency caused by mutations in AICDA encoding activation-induced cytidine deaminase (AID) which impair immunoglobulin class switch recombination (CSR) and somatic hypermutation (SHM). Whereas autosomal-recessive AID-deficiency (AR-AID) affects both CSR and SHM, the autosomal-dominant form (AD-AID) due to C-terminal heterozygous variants completely abolishes CSR but only partially affects SHM. AR-AID patients display enhanced germinal center (GC) reactions and autoimmune manifestations, which are not present in AD-AID, suggesting that SHM but not CSR regulates GC reactions and peripheral B cell tolerance. Herein, we describe two siblings with HIGM2 due to a novel homozygous AICDA mutation (c.428-1G > T) which disrupts the splice acceptor site of exon 4 and results in the sole expression of a truncated AID variant that lacks 10 highly conserved amino acids encoded by exon 4 (AID-ΔE4a). AID-ΔE4a patients suffered from defective CSR and enhanced GC reactions and were therefore indistinguishable from other AR-AID patients. However, the AID-ΔE4a variant only partially affected SHM as observed in AD-AID patients. In addition, AID-ΔE4a but not AD-AID patients revealed impaired targeting of mutational hotspot motives and distorted mutational patterns. Hence, qualitative defects in AID function and altered SHM rather than global decreased SHM activity may account for the disease phenotype in these patients.

Long-Term Follow-Up of Newborns with 22q11 Deletion Syndrome and Low TRECs

Background

Population-based neonatal screening using T-cell receptor excision circles (TRECs) identifies infants with profound T lymphopenia, as seen in cases of severe combined immunodeficiency, and in a subgroup of infants with 22q11 deletion syndrome (22q11DS).

Purpose

To investigate the long-term prognostic value of low levels of TRECs in newborns with 22q11DS.

Methods

Subjects with 22q11DS and low TRECs at birth (22q11Low, N=10), matched subjects with 22q11DS and normal TRECs (22q11Normal, N=10), and matched healthy controls (HC, N=10) were identified. At follow-up (median age 16 years), clinical and immunological characterizations, covering lymphocyte subsets, immunoglobulins, TRECs, T-cell receptor repertoires, and relative telomere length (RTL) measurements were performed.

Results

At follow-up, the 22q11Low group had lower numbers of naïve T-helper cells, naïve T-regulatory cells, naïve cytotoxic T cells, and persistently lower TRECs compared to healthy controls. Receptor repertoires showed skewed V-gene usage for naïve T-helper cells, whereas for naïve cytotoxic T cells, shorter RTL and a trend towards higher clonality were found. Multivariate discriminant analysis revealed a clear distinction between the three groups and a skewing towards Th17 differentiation of T-helper cells, particularly in the 22q11Low individuals. Perturbations of B-cell subsets were found in both the 22q11Low and 22q11Normal group compared to the HC group, with larger proportions of naïve B cells and lower levels of memory B cells, including switched memory B cells.

Conclusions

This long-term follow-up study shows that 22q11Low individuals have persistent immunologic aberrations and increased risk for immune dysregulation, indicating the necessity of lifelong monitoring.

Clinical Implications

This study elucidates the natural history of childhood immune function in newborns with 22q11DS and low TRECs, which may facilitate the development of programs for long-term monitoring and therapeutic choices.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10875-021-01201-5.

Heavy chain sequence-based classifier for the specificity of human antibodies

Antibodies specifically bind to antigens and are an essential part of the immune system. Hence, antibodies are powerful tools in research and diagnostics. High-throughput sequencing technologies have promoted comprehensive profiling of the immune repertoire, which has resulted in large amounts of antibody sequences that remain to be further analyzed. In this study, antibodies were downloaded from IMGT/LIGM-DB and Sequence Read Archive databases. Contributing features from antibody heavy chains were formulated as numerical inputs and fed into an ensemble machine learning classifier to classify the antigen specificity of six classes of antibodies, namely anti-HIV-1, anti-influenza virus, anti-pneumococcal polysaccharide, anti-citrullinated protein, anti-tetanus toxoid and anti-hepatitis B virus. The classifier was validated using cross-validation and a testing dataset. The ensemble classifier achieved a macro-average area under the receiver operating characteristic curve (AUC) of 0.9246 from the 10-fold cross-validation, and 0.9264 for the testing dataset. Among the contributing features, the contribution of the complementarity-determining regions was 53.1% and that of framework regions was 46.9%, and the amino acid mutation rates occupied the first and second ranks among the top five contributing features. The classifier and insights provided in this study could promote the mechanistic study, isolation and utilization of potential therapeutic antibodies.

CloneRetriever: An Automated Algorithm to Identify Clonal B and T Cell Gene Rearrangements by Next-Generation Sequencing for the Diagnosis of Lymphoid Malignancies

BACKGROUND

Clonal immunoglobulin and T-cell receptor rearrangements serve as tumor-specific markers that have become mainstays of the diagnosis and monitoring of lymphoid malignancy. Next-generation sequencing (NGS) techniques targeting these loci have been successfully applied to lymphoblastic leukemia and multiple myeloma for minimal residual disease detection. However, adoption of NGS for primary diagnosis remains limited.

METHODS

We addressed the bioinformatics challenges associated with immune cell sequencing and clone detection by designing a novel web tool, CloneRetriever (CR), which uses machine-learning principles to generate clone classification schemes that are customizable, and can be applied to large datasets. CR has 2 applications-a "validation" mode to derive a clonality classifier, and a "live" mode to screen for clones by applying a validated and/or customized classifier. In this study, CR-generated multiple classifiers using 2 datasets comprising 106 annotated patient samples. A custom classifier was then applied to 36 unannotated samples.

RESULTS

The optimal classifier for clonality required clonal dominance ≥4.5× above background, read representation ≥8% of all reads, and technical replicate agreement. Depending on the dataset and analysis step, the optimal algorithm yielded sensitivities of 81%-90%, specificities of 97%-100%, areas under the curve of 91%-94%, positive predictive values of 92-100%, and negative predictive values of 88%-98%. Customization of the algorithms yielded 95%-100% concordance with gold-standard clonality determination, including rescue of indeterminate samples. Application to a set of unknowns showed concordance rates of 83%-96%.

CONCLUSIONS

CR is an out-of-the-box ready and user-friendly software designed to identify clonal rearrangements in large NGS datasets for the diagnosis of lymphoid malignancies.

Targeting B cells in the pre-phase of systemic autoimmunity globally interferes with autoimmune pathology

Systemic lupus erythematosus (SLE) is characterized by a loss of self-tolerance, systemic inflammation, and multi-organ damage. While a variety of therapeutic interventions are available, it has become clear that an early diagnosis and treatment may be key to achieve long lasting therapeutic responses and to limit irreversible organ damage. Loss of humoral tolerance including the appearance of self-reactive antibodies can be detected years before the actual onset of the clinical autoimmune disease, representing a potential early point of intervention. Not much is known, however, about how and to what extent this pre-phase of disease impacts the onset and development of subsequent autoimmunity. By targeting the B cell compartment in the pre-disease phase of a spontaneous mouse model of SLE we now show, that resetting the humoral immune system during the clinically unapparent phase of the disease globally alters immune homeostasis delaying the downstream development of systemic autoimmunity.

•

The clinically unapparent pre-phase of SLE impacts clinical disease

•

Autoreactive IgM antibodies represent a biomarker for early therapeutic intervention

•

Pre-phase B cells orchestrate clinical disease

•

Depleting pre-phase B cells diminishes disease pathology

RAPID: A Rep-Seq Dataset Analysis Platform With an Integrated Antibody Database

The antibody repertoire is a critical component of the adaptive immune system and is believed to reflect an individual’s immune history and current immune status. Delineating the antibody repertoire has advanced our understanding of humoral immunity, facilitated antibody discovery, and showed great potential for improving the diagnosis and treatment of disease. However, no tool to date has effectively integrated big Rep-seq data and prior knowledge of functional antibodies to elucidate the remarkably diverse antibody repertoire. We developed a Rep-seq dataset Analysis Platform with an Integrated antibody Database (RAPID; https://rapid.zzhlab.org/), a free and web-based tool that allows researchers to process and analyse Rep-seq datasets. RAPID consolidates 521 WHO-recognized therapeutic antibodies, 88,059 antigen- or disease-specific antibodies, and 306 million clones extracted from 2,449 human IGH Rep-seq datasets generated from individuals with 29 different health conditions. RAPID also integrates a standardized Rep-seq dataset analysis pipeline to enable users to upload and analyse their datasets. In the process, users can also select set of existing repertoires for comparison. RAPID automatically annotates clones based on integrated therapeutic and known antibodies, and users can easily query antibodies or repertoires based on sequence or optional keywords. With its powerful analysis functions and rich set of antibody and antibody repertoire information, RAPID will benefit researchers in adaptive immune studies.

Light chain skewing in autoantibodies and B-cell receptors of the citrullinated antigen-binding B-cell response in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic autoimmune disease affecting 1% of the world population. RA is associated with the presence of autoantibodies, of which anti-citrullinated protein antibodies (ACPA) are most prominent. ACPA are produced by citrullinated antigen-binding B cells that have presumably survived tolerance checkpoints. So far, it is unclear how and when such autoreactive B cells emerge. Light chain (LC) rearrangement and mutation rates can be informative with regard to selection steps during B-cell development. Therefore, we studied LC characteristics of ACPA-expressing B cells and secreted ACPA with the aim to better understand the development of this disease-specific, autoreactive B-cell response. Paired ACPA-IgG and ACPA-depleted IgG were isolated from serum (n = 87) and synovial fluid (SF, n = 21) of patients with established RA. We determined the LC composition for each fraction by ELISA using kappa(Igκ)- and lambda(Igλ) LC-specific antibodies. Cellular LC expression was determined using flow cytometry. In addition, we used a B-cell receptor (BCR)-specific PCR to obtain LC variable region sequences of citrullinated antigen- and tetanus toxoid (TT)-binding B cells. In serum, we observed an increased frequency of lambda LC in ACPA-IgG (1.64:1) compared to control IgG (2.03:1) and to the κ/λ ratio reported for healthy individuals (2:1). A similar trend towards higher frequencies of lambda LCs was observed for ACPA-IgG in SF (1.84:1). Additionally, the percentage of Igλ-expressing B cells was higher for citrullinated antigen-binding B cells (51%) compared to TT-specific (43%) and total CD19+CD20+ B cells (36%). Moreover, an increased Igλ percentage was observed in BCR-sequences derived from ACPA-expressing (49%) compared to TT-specific B cells (34%). Taken together, we report an enhanced frequency of lambda LCs in the secreted ACPA-IgG repertoire and, on the cellular level, in BCR sequences of ACPA-expressing B cells compared to control. This skewing in the autoreactive B-cell repertoire could reflect a process of active selection.

The Phenotypic Spectrum of PNKP-Associated Disease and the Absence of Immunodeficiency and Cancer Predisposition in a Dutch Cohort

BACKGROUND

We aimed to expand the number of currently known pathogenic PNKP mutations, to study the phenotypic spectrum, including radiological characteristics and genotype-phenotype correlations, and to assess whether immunodeficiency and increased cancer risk are part of the DNA repair disorder caused by mutations in the PNKP gene.

METHODS

We evaluated nine patients with PNKP mutations. A neurological history and examination was obtained. All patients had undergone neuroimaging and genetic testing as part of the prior diagnostic process. Laboratory measurements included potential biomarkers, and, in the context of a DNA repair disorder, we performed a detailed immunologic evaluation, including B cell repertoire analysis.

RESULTS

We identified three new mutations in the PNKP gene and confirm the phenotypic spectrum of PNKP-associated disease, ranging from microcephaly, seizures, and developmental delay to ataxia with oculomotor apraxia type 4. Irrespective of the phenotype, alpha-fetoprotein is a biochemical marker and increases with age and progression of the disease. On neuroimaging, (progressive) cerebellar atrophy was a universal feature. No clinical signs of immunodeficiency were present, and immunologic assessment was unremarkable. One patient developed cancer, but this was attributed to a concurrent von Hippel-Lindau mutation.

CONCLUSIONS

Immunodeficiency and cancer predisposition do not appear to be part of PNKP-associated disease, contrasting many other DNA repair disorders. Furthermore, our study illustrates that the previously described syndromes microcephaly, seizures, and developmental delay, and ataxia with oculomotor apraxia type 4, represent the extremes of an overlapping spectrum of disease. Cerebellar atrophy and elevated serum alpha-fetoprotein levels are early diagnostic findings across the entire phenotypical spectrum.

Teleost IgD+IgM- B Cells Mount Clonally Expanded and Mildly Mutated Intestinal IgD Responses in the Absence of Lymphoid Follicles

Immunoglobulin D (IgD) is an ancient antibody with dual membrane-bound and fluid-phase antigen receptor functions. The biology of secreted IgD remains elusive. Here, we demonstrate that teleost IgD⁺IgM⁻ plasmablasts constitute a major lymphocyte population in some mucosal surfaces, including the gut mucosa. Remarkably, secreted IgD binds to gut commensal bacteria, which in turn stimulate IgD gene transcription in gut B cells. Accordingly, secreted IgD from gut as well as gill mucosae, but not the spleen, show a V(D)J gene configuration consistent with microbiota-driven clonal expansion and diversification, including mild somatic hypermutation. By showing that secreted IgD establishes a mutualistic relationship with commensals, our findings suggest that secreted IgD may play an evolutionary conserved role in mucosal homeostasis.

•

IgD⁺IgM⁻ B cells constitute the main non-IgT B cell subset in rainbow trout guts

•

Gut IgD responses establish a two-way interaction with the local microbiota

•

Mucosal but not splenic IgD undergoes clonal expansion and diversification

•

Despite the lack of germinal centers, mucosal IgD is mildly mutated in rainbow trout

Methods for sequence and structural analysis of B and T cell receptor repertoires

B cell receptors (BCRs) and T cell receptors (TCRs) make up an essential network of defense molecules that, collectively, can distinguish self from non-self and facilitate destruction of antigen-bearing cells such as pathogens or tumors. The analysis of BCR and TCR repertoires plays an important role in both basic immunology as well as in biotechnology. Because the repertoires are highly diverse, specialized software methods are needed to extract meaningful information from BCR and TCR sequence data. Here, we review recent developments in bioinformatics tools for analysis of BCR and TCR repertoires, with an emphasis on those that incorporate structural features. After describing the recent sequencing technologies for immune receptor repertoires, we survey structural modeling methods for BCR and TCRs, along with methods for clustering such models. We review downstream analyses, including BCR and TCR epitope prediction, antibody-antigen docking and TCR-peptide-MHC Modeling. We also briefly discuss molecular dynamics in this context.

Polymerase δ deficiency causes syndromic immunodeficiency with replicative stress

Polymerase δ is essential for eukaryotic genome duplication and synthesizes DNA at both the leading and lagging strands. The polymerase δ complex is a heterotetramer comprising the catalytic subunit POLD1 and the accessory subunits POLD2, POLD3, and POLD4. Beyond DNA replication, the polymerase δ complex has emerged as a central element in genome maintenance. The essentiality of polymerase δ has constrained the generation of polymerase δ-knockout cell lines or model organisms and, therefore, the understanding of the complexity of its activity and the function of its accessory subunits. To our knowledge, no germline biallelic mutations affecting this complex have been reported in humans. In patients from 2 independent pedigrees, we have identified what we believe to be a novel syndrome with reduced functionality of the polymerase δ complex caused by germline biallelic mutations in POLD1 or POLD2 as the underlying etiology of a previously unknown autosomal-recessive syndrome that combines replicative stress, neurodevelopmental abnormalities, and immunodeficiency. Patients' cells showed impaired cell-cycle progression and replication-associated DNA lesions that were reversible upon overexpression of polymerase δ. The mutations affected the stability and interactions within the polymerase δ complex or its intrinsic polymerase activity. We believe our discovery of human polymerase δ deficiency identifies the central role of this complex in the prevention of replication-related DNA lesions, with particular relevance to adaptive immunity.

Evaluation of Somatic Hypermutation Status in Chronic Lymphocytic Leukemia (CLL) in the Era of Next Generation Sequencing

Somatic hypermutation (SHM) status provides an important prognostic indicator for chronic lymphocytic leukemia (CLL), a very common type of mature B-cell leukemia. Owing to the adverse prognosis associated with an unmutated immunoglobulin heavy chain variable (IGHV) status, SHM testing is performed as a standard of care in CLL. Conventionally, SHM testing has been performed using labor intensive and primarily analog Sanger sequencing method following PCR amplification of the clonal immunoglobulin heavy chain gene rearrangements in CLL cells. In comparison, recent availability of next generation sequencing (NGS) allows more versatile detection and direct identification of clonal immunoglobulin gene rearrangements in neoplastic B-cell populations. The ability to identify specific clonal IGHV signature(s) in both baseline (diagnostic) and post-treatment settings enables unique clinical applications of NGS such as determination of SHM status, minimal residual disease (MRD) monitoring, clonal heterogeneity and B cell receptor IG stereotypy. We provide a review of current practices and recommendations for SHM determination using NGS including examples of difficult cases.

Delineating Human B Cell Precursor Development With Genetically Identified PID Cases as a Model

B-cell precursors (BCP) arise from hematopoietic stem cells in bone marrow (BM). Identification and characterization of the different BCP subsets has contributed to the understanding of normal B-cell development. BCP first rearrange their immunoglobulin (Ig) heavy chain (IGH) genes to form the pre-B-cell receptor (pre-BCR) complex together with surrogate light chains. Appropriate signaling via this pre-BCR complex is followed by rearrangement of the Ig light chain genes, resulting in the formation, and selection of functional BCR molecules. Consecutive production, expression, and functional selection of the pre-BCR and BCR complexes guide the BCP differentiation process that coincides with corresponding immunophenotypic changes. We studied BCP differentiation in human BM samples from healthy controls and patients with a known genetic defect in V(D)J recombination or pre-BCR signaling to unravel normal immunophenotypic changes and to determine the effect of differentiation blocks caused by the specific genetic defects. Accordingly, we designed a 10-color antibody panel to study human BCP development in BM by flow cytometry, which allows identification of classical preB-I, preB-II, and mature B-cells as defined via BCR-related markers with further characterization by additional markers. We observed heterogeneous phenotypes associated with more than one B-cell maturation pathway, particularly for the preB-I and preB-II stages in which V(D)J recombination takes place, with asynchronous marker expression patterns. Next Generation Sequencing of complete IGH gene rearrangements in sorted BCP subsets unraveled their rearrangement status, indicating that BCP differentiation does not follow a single linear pathway. In conclusion, B-cell development in human BM is not a linear process, but a rather complex network of parallel pathways dictated by V(D)J-recombination-driven checkpoints and pre-BCR/BCR mediated-signaling occurring during B-cell production and selection. It can also be described as asynchronous, because precursor B-cells do not differentiate as full population between the different stages, but rather transit as a continuum, which seems influenced (in part) by V-D-J recombination-driven checkpoints.

sumrep: A Summary Statistic Framework for Immune Receptor Repertoire Comparison and Model Validation

The adaptive immune system generates an incredible diversity of antigen receptors for B and T cells to keep dangerous pathogens at bay. The DNA sequences coding for these receptors arise by a complex recombination process followed by a series of productivity-based filters, as well as affinity maturation for B cells, giving considerable diversity to the circulating pool of receptor sequences. Although these datasets hold considerable promise for medical and public health applications, the complex structure of the resulting adaptive immune receptor repertoire sequencing (AIRR-seq) datasets makes analysis difficult. In this paper we introduce sumrep, an R package that efficiently performs a wide variety of repertoire summaries and comparisons, and show how sumrep can be used to perform model validation. We find that summaries vary in their ability to differentiate between datasets, although many are able to distinguish between covariates such as donor, timepoint, and cell type for BCR and TCR repertoires. We show that deletion and insertion lengths resulting from V(D)J recombination tend to be more discriminative characterizations of a repertoire than summaries that describe the amino acid composition of the CDR3 region. We also find that state-of-the-art generative models excel at recapitulating gene usage and recombination statistics in a given experimental repertoire, but struggle to capture many physiochemical properties of real repertoires.

N-Glycosylation Site Analysis of Citrullinated Antigen-Specific B-Cell Receptors Indicates Alternative Selection Pathways During Autoreactive B-Cell Development

Many autoimmune diseases are hallmarked by autoreactive B and plasma cell responses that are directly or indirectly involved in disease pathogenesis. These B-cell responses show large variability between diseases, both in terms of the secreted autoantibody repertoire and the dynamics and characteristics of the underlying B-cell responses. Hence, different mechanisms have been proposed to explain the emergence of autoreactive B cells in an otherwise self-tolerant immune system. Notably, most mechanistic insights have been obtained from murine studies using models harboring genetic modifications of B and T cells. Given recent technological advances that have rendered autoreactive human B cells accessible for analysis, we here discuss the phenomenon of extensive N-glycosylation of the B-cell receptor (BCR) variable domain of a prototypic human autoreactive B-cell response and its potential role in the generation of autoimmunity. Anti-citrullinated protein antibodies (ACPA) hallmark the most disease-specific autoimmune response in Rheumatoid Arthritis (RA). Remarkably, ACPA-IgG are heavily N-glycosylated in the variable domain due to somatic mutations that generate abundant N-glycosylation consensus sequences. These sites, obtained from full-length BCR sequences of ACPA-expressing B cells from 12 ACPA-positive RA patients, were here analyzed in detail. Sites that required a single nucleotide mutation to be generated were defined as single somatic hypermutation (s-SHM) sites, whereas sites requiring multiple mutations were defined as m-SHM sites. IgG sequences of 12 healthy donors were used as control. Computational modeling of the germinal center reaction (CLONE algorithm) was used with the germline counterparts of ACPA-IgG heavy chain (HC) sequences to simulate the germinal center response. Our analyses revealed an abundance of N-glycosylation sites in ACPA-IgG HC that frequently required multiple mutations and predominated in specific positions. Based on these data, and taking into account recent insights into the dynamics of the ACPA-response during disease development, we here discuss the hypothesis that N-glycosylation sites in ACPA-IgG variable domains could lead to alternative, possibly antibody affinity-independent selection forces. Presumably, this occurs during germinal center responses allowing these B cells to escape from putative tolerance checkpoints, thereby driving autoreactive B cell development in the pathogenesis of RA.

The presence of CLL-associated stereotypic B cell receptors in the normal BCR repertoire from healthy individuals increases with age

Background

Aging is known to induce immunosenescence, resulting in alterations in both the innate and adaptive immune system. Here we evaluated the effects of aging on B cell subsets in peripheral blood of 155 immunologically healthy individuals in four age categories (range 20-95y) via multi-parameter flow cytometry. Furthermore, we studied the naive and antigen-experienced B cell receptor (BCR) repertoire of different age groups and compared it to the clonal BCR repertoire of chronic lymphocytic leukemia (CLL), a disease typically presenting in elderly individuals.

Results

Total numbers and relative frequencies of B cells were found to decline upon aging, with reductions in transitional B cells, memory cell types, and plasma blasts in the 70 + y group. The BCR repertoire of naive mature B cells and antigen-experienced B cells did not clearly alter until age 70y. Clear changes in IGHV gene usage were observed in naive mature B cells of 70 + y individuals, with a transitional pattern in the 50-70y group. IGHV gene usage of naive mature B cells of the 50-70y, but not the 70 + y, age group resembled that of both younger (50-70y) and older (70 + y) CLL patients. Additionally, CLL-associated stereotypic BCR were found as part of the healthy control BCR repertoire, with an age-associated increase in frequency of several stereotypic BCR (particularly subsets #2 and #5).

Conclusion

Composition of the peripheral B cell compartment changes with ageing, with clear reductions in non-switched and CD27 + IgG+ switched memory B cells and plasma blasts in especially the 70 + y group. The BCR repertoire is relatively stable until 70y, whereafter differences in IGHV gene usage are seen. Upon ageing, an increasing trend in the occurrence of particular CLL-associated stereotypic BCR is observed.

Electronic supplementary material

The online version of this article (10.1186/s12979-019-0163-x) contains supplementary material, which is available to authorized users.

Repertoire Sequencing of B Cells Elucidates the Role of UNG and Mismatch Repair Proteins in Somatic Hypermutation in Humans

The generation of high-affinity antibodies depends on somatic hypermutation (SHM). SHM is initiated by the activation-induced cytidine deaminase (AID), which generates uracil (U) lesions in the B-cell receptor (BCR) encoding genes. Error-prone processing of U lesions creates a typical spectrum of point mutations during SHM. The aim of this study was to determine the molecular mechanism of SHM in humans; currently available knowledge is limited by the number of mutations analyzed per patient. We collected a unique cohort of 10 well-defined patients with bi-allelic mutations in genes involved in base excision repair (BER) (UNG) or mismatch repair (MMR) (MSH2, MSH6, or PMS2) and are the first to present next-generation sequencing (NGS) data of the BCR, allowing us to study SHM extensively in humans. Analysis using ARGalaxy revealed selective skewing of SHM mutation patterns specific for each genetic defect, which are in line with the five-pathway model of SHM that was recently proposed based on mice data. However, trans-species comparison revealed differences in the role of PMS2 and MSH2 in strand targeting between mice and man. In conclusion, our results indicate a role for UNG, MSH2, MSH6, and PMS2 in the generation of SHM in humans comparable to their function in mice. However, we observed differences in strand targeting between humans and mice, emphasizing the importance of studying molecular mechanisms in a human setting. The here developed method combining NGS and ARGalaxy analysis of BCR mutation data forms the basis for efficient SHM analyses of other immune deficiencies.

BRepertoire: a user-friendly web server for analysing antibody repertoire data

Antibody repertoire analysis by high throughput sequencing is now widely used, but a persisting challenge is enabling immunologists to explore their data to discover discriminating repertoire features for their own particular investigations. Computational methods are necessary for large-scale evaluation of antibody properties. We have developed BRepertoire, a suite of user-friendly web-based software tools for large-scale statistical analyses of repertoire data. The software is able to use data preprocessed by IMGT, and performs statistical and comparative analyses with versatile plotting options. BRepertoire has been designed to operate in various modes, for example analysing sequence-specific V(D)J gene usage, discerning physico-chemical properties of the CDR regions and clustering of clonotypes. Those analyses are performed on the fly by a number of R packages and are deployed by a shiny web platform. The user can download the analysed data in different table formats and save the generated plots as image files ready for publication. We believe BRepertoire to be a versatile analytical tool that complements experimental studies of immune repertoires. To illustrate the server’s functionality, we show use cases including differential gene usage in a vaccination dataset and analysis of CDR3H properties in old and young individuals. The server is accessible under http://mabra.biomed.kcl.ac.uk/BRepertoire.

B Cell Reconstitution and Influencing Factors After Hematopoietic Stem Cell Transplantation in Children

B cell reconstitution after hematopoietic stem cell transplantation (HSCT) is variable and influenced by different patient, donor, and treatment related factors. In this review we describe B cell reconstitution after pediatric allogeneic HST, including the kinetics of reconstitution of the different B cell subsets and the development of the B cell repertoire, and discuss the influencing factors. Observational studies show important roles for stem cell source, conditioning regimen, and graft vs. host disease in B cell reconstitution. In addition, B cell recovery can play an important role in post-transplant infections and vaccine responses to encapsulated bacteria, such as pneumococcus. A substantial number of patients experience impaired B cell function and/or dependency on Ig substitution after allogeneic HSCT. The underlying mechanisms are largely unresolved. The integrated aspects of B cell recovery after HSCT, especially BCR repertoire reconstitution, are awaiting further investigation using modern techniques in order to gain more insight into B cell reconstitution and to develop strategies to improve humoral immunity after allogeneic HSCT.

IVIg-induced plasmablasts in patients with Guillain-Barré syndrome

Objective

The Guillain-Barré syndrome (GBS) is an acute, immune-mediated disease of peripheral nerves. Plasmablasts and plasma cells play a central role in GBS by producing neurotoxic antibodies. The standard treatment for GBS is high-dose intravenous immunoglobulins (IVIg), however the working mechanism is unknown and the response to treatment is highly variable. We aimed to determine whether IVIg changes the frequency of B-cell subsets in patients with GBS.

Methods

Peripheral blood mononuclear cells were isolated from 67 patients with GBS before and/or 1, 2, 4, and 12 weeks after treatment with high-dose IVIg. B-cell subset frequencies were determined by flow cytometry and related to serum immunoglobulin levels. Immunoglobulin transcripts before and after IVIg treatment were examined by next-generation sequencing. Antiglycolipid antibodies were determined by ELISA.

Results

Patients treated with IVIg demonstrated a strong increase in plasmablasts, which peaked 1 week after treatment. Flow cytometry identified a relative increase in IgG2 plasmablasts posttreatment. Within IGG sequences, dominant clones were identified which were also IGG2 and had different immunoglobulin sequences compared to pretreatment samples. High plasmablast frequencies after treatment correlated with an increase in serum IgG and IgM, suggesting endogenous production. Patients with a high number of plasmablasts started to improve earlier (P = 0.015) and were treated with a higher dose of IVIg.

Interpretation

High-dose IVIg treatment alters the distribution of B-cell subsets in the peripheral blood of GBS patients, suggesting de novo (oligo-)clonal B-cell activation. Very high numbers of plasmablasts after IVIg therapy may be a potential biomarker for fast clinical recovery.

Identification of CVID Patients With Defects in Immune Repertoire Formation or Specification

Common variable immune deficiency disorder (CVID) is the most clinically relevant cause of antibody failure. It is a highly heterogeneous disease with different underlying etiologies. CVID has been associated with a quantitative B cell defect, however, little is known about the quality of B cells present. Here, we studied the naïve and antigen selected B-cell receptor (BCR) repertoire in 33 CVID patients using next generation sequencing, to investigate B cells quality. Analysis for each individual patient revealed whether they have a defect in immune repertoire formation [V(D)J recombination] or specification (somatic hypermutation, subclass distribution, or selection). The naïve BCR repertoire was normal in most of the patients, although alterations in repertoire diversity and the junctions were found in a limited number of patients indicating possible defects in early B-cell development or V(D)J recombination in these patients. In contrast, major differences were found in the antigen selected BCR repertoire. Here, most patients (15/17) showed a reduced frequency of somatic hypermutation (SHM), changes in subclass distribution and/or minor alterations in antigen selection. Together these data show that in our CVID cohort only a small number of patients have a defect in formation of the naïve BCR repertoire, whereas the clear majority of patients have disturbances in their antigen selected repertoire, suggesting a defect in repertoire specification in the germinal centers of these patients. This highlights that CVID patients not only have a quantitative B cell defect, but that also the quality of, especially post germinal center B cells, is impaired.

ASAP - A Webserver for Immunoglobulin-Sequencing Analysis Pipeline

Reproducible and robust data on antibody repertoires are invaluable for basic and applied immunology. Next-generation sequencing (NGS) of antibody variable regions has emerged as a powerful tool in systems immunology, providing quantitative molecular information on antibody polyclonal composition. However, major computational challenges exist when analyzing antibody sequences, from error handling to hypermutation profiles and clonal expansion analyses. In this work, we developed the ASAP (A webserver for Immunoglobulin-Seq Analysis Pipeline) webserver (https://asap.tau.ac.il). The input to ASAP is a paired-end sequence dataset from one or more replicates, with or without unique molecular identifiers. These datasets can be derived from NGS of human or murine antibody variable regions. ASAP first filters and annotates the sequence reads using public or user-provided germline sequence information. The ASAP webserver next performs various calculations, including somatic hypermutation level, CDR3 lengths, V(D)J family assignments, and V(D)J combination distribution. These analyses are repeated for each replicate. ASAP provides additional information by analyzing the commonalities and differences between the repeats (“joint” analysis). For example, ASAP examines the shared variable regions and their frequency in each replicate to determine which sequences are less likely to be a result of a sample preparation derived and/or sequencing errors. Moreover, ASAP clusters the data to clones and reports the identity and prevalence of top ranking clones (clonal expansion analysis). ASAP further provides the distribution of synonymous and non-synonymous mutations within the V genes somatic hypermutations. Finally, ASAP provides means to process the data for proteomic analysis of serum/secreted antibodies by generating a variable region database for liquid chromatography high resolution tandem mass spectrometry (LC-MS/MS) interpretation. ASAP is user-friendly, free, and open to all users, with no login requirement. ASAP is applicable for researchers interested in basic questions related to B cell development and differentiation, as well as applied researchers who are interested in vaccine development and monoclonal antibody engineering. By virtue of its user-friendliness, ASAP opens the antibody analysis field to non-expert users who seek to boost their research with immune repertoire analysis.

No Overt Clinical Immunodeficiency Despite Immune Biological Abnormalities in Patients With Constitutional Mismatch Repair Deficiency

Immunoglobulin class-switch recombination (CSR) and somatic hypermutations (SHMs) are prerequisites for antibody and immunoglobulin receptor maturation and adaptive immune diversity. The mismatch repair (MMR) machinery, consisting of homologs of MutSα, MutLα, and MutSβ (MSH2/MSH6, MLH1/PMS2, and MSH2/MSH3, respectively) and other proteins, is involved in CSR, primarily acting as a backup for nonhomologous end-joining repair of activation-induced cytidine deaminase-induced DNA mismatches and, furthermore, in addition to error-prone polymerases, in the repair of SHM-induced DNA breaks. A varying degree of antibody formation defect, from IgA or selective IgG subclass deficiency to common variable immunodeficiency and hyper-IgM syndrome, has been detected in a small number of patients with constitutional mismatch repair deficiency (CMMRD) due to biallelic loss-of-function mutations in one of the MMR genes (PMS2, MSH6, MLH1, or MSH2). To elucidate the clinical relevance of a presumed primary immunodeficiency (PID) in CMMRD, we systematically collected clinical history and laboratory data of a cohort of 15 consecutive, unrelated patients (10 not previously reported) with homozygous/compound heterozygous mutations in PMS2 (n = 8), MSH6 (n = 5), and MLH1 (n = 2), most of whom manifested with typical malignancies during childhood. Detailed descriptions of their genotypes, phenotypes, and family histories are provided. Importantly, none of the patients showed any clinical warning signs of PID (infections, immune dysregulation, inflammation, failure to thrive, etc.). Furthermore, we could not detect uniform or specific patterns of laboratory abnormalities. The concentration of IgM was increased in 3 out of 12, reduced in 3 out of 12, and normal in 6 out of 12 patients, while concentrations of IgG and IgG subclasses, except IgG4, and of IgA, and specific antibody formation were normal in most. Class-switched B memory cells were reduced in 5 out of 12 patients, and in 9 out of 12 also the CD38^hiIgM⁻ plasmablasts were reduced. Furthermore, results of next generation sequencing-based analyses of antigen-selected B-cell receptor rearrangements showed a significantly reduced frequency of SHM and an increased number of rearranged immunoglobulin heavy chain (IGH) transcripts that use IGHG3, IGHG1, and IGHA1 subclasses. T cell subsets and receptor repertoires were unaffected. Together, neither clinical nor routine immunological laboratory parameters were consistently suggestive of PID in these CMMRD patients, but previously shown abnormalities in SHM and rearranged heavy chain transcripts were confirmed.

Dysregulation of B Cell Activity During Proliferative Kidney Disease in Rainbow Trout

Proliferative kidney disease (PKD) is a widespread disease caused by the endoparasite Tetracapsuloides bryosalmonae (Myxozoa: Malacosporea). Clinical disease, provoked by the proliferation of extrasporogonic parasite stages, is characterized by a chronic kidney pathology with underlying transcriptional changes indicative of altered B cell responses and dysregulated T-helper cell-like activities. Despite the relevance of PKD to European and North American salmonid aquaculture, no studies, to date, have focused on further characterizing the B cell response during the course of this disease. Thus, in this work, we have studied the behavior of diverse B cell populations in rainbow trout (Oncorhynchus mykiss) naturally infected with T. bryosalmonae at different stages of preclinical and clinical disease. Our results show a clear upregulation of all trout immunoglobulins (Igs) (IgM, IgD, and IgT) demonstrated by immunohistochemistry and Western blot analysis, suggesting the alteration of diverse B cell populations that coexist in the infected kidney. Substantial changes in IgM, IgD, and IgT repertoires were also identified throughout the course of the disease further pointing to the involvement of the three Igs in PKD through what appear to be independently regulated mechanisms. Thus, our results provide strong evidence of the involvement of IgD in the humoral response to a specific pathogen for the first time in teleosts. Nevertheless, it was IgT, a fish-specific Ig isotype thought to be specialized in mucosal immunity, which seemed to play a prevailing role in the kidney response to T. bryosalmonae. We found that IgT was the main Ig coating extrasporogonic parasite stages, IgT⁺ B cells were the main B cell subset that proliferated in the kidney with increasing kidney pathology, and IgT was the Ig for which more significant changes in repertoire were detected. Hence, although our results demonstrate a profound dysregulation of different B cell subsets during PKD, they point to a major involvement of IgT in the immune response to the parasite. These results provide further insights into the pathology of PKD that may facilitate the future development of control strategies.

VDJServer: A Cloud-Based Analysis Portal and Data Commons for Immune Repertoire Sequences and Rearrangements

Background

Recent technological advances in immune repertoire sequencing have created tremendous potential for advancing our understanding of adaptive immune response dynamics in various states of health and disease. Immune repertoire sequencing produces large, highly complex data sets, however, which require specialized methods and software tools for their effective analysis and interpretation.

Results

VDJServer is a cloud-based analysis portal for immune repertoire sequence data that provide access to a suite of tools for a complete analysis workflow, including modules for preprocessing and quality control of sequence reads, V(D)J gene segment assignment, repertoire characterization, and repertoire comparison. VDJServer also provides sophisticated visualizations for exploratory analysis. It is accessible through a standard web browser via a graphical user interface designed for use by immunologists, clinicians, and bioinformatics researchers. VDJServer provides a data commons for public sharing of repertoire sequencing data, as well as private sharing of data between users. We describe the main functionality and architecture of VDJServer and demonstrate its capabilities with use cases from cancer immunology and autoimmunity.

Conclusion

VDJServer provides a complete analysis suite for human and mouse T-cell and B-cell receptor repertoire sequencing data. The combination of its user-friendly interface and high-performance computing allows large immune repertoire sequencing projects to be analyzed with no programming or software installation required. VDJServer is a web-accessible cloud platform that provides access through a graphical user interface to a data management infrastructure, a collection of analysis tools covering all steps in an analysis, and an infrastructure for sharing data along with workflows, results, and computational provenance. VDJServer is a free, publicly available, and open-source licensed resource.

Analyzing Immunoglobulin Repertoires

Somatic assembly of T cell receptor and B cell receptor (BCR) genes produces a vast diversity of lymphocyte antigen recognition capacity. The advent of efficient high-throughput sequencing of lymphocyte antigen receptor genes has recently generated unprecedented opportunities for exploration of adaptive immune responses. With these opportunities have come significant challenges in understanding the analysis techniques that most accurately reflect underlying biological phenomena. In this regard, sample preparation and sequence analysis techniques, which have largely been borrowed and adapted from other fields, continue to evolve. Here, we review current methods and challenges of library preparation, sequencing and statistical analysis of lymphocyte receptor repertoire studies. We discuss the general steps in the process of immune repertoire generation including sample preparation, platforms available for sequencing, processing of sequencing data, measurable features of the immune repertoire, and the statistical tools that can be used for analysis and interpretation of the data. Because BCR analysis harbors additional complexities, such as immunoglobulin (Ig) (i.e., antibody) gene somatic hypermutation and class switch recombination, the emphasis of this review is on Ig/BCR sequence analysis.

Next-Generation Sequencing Analysis of the Human TCRγδ+ T-Cell Repertoire Reveals Shifts in Vγ- and Vδ-Usage in Memory Populations upon Aging

Immunological aging remodels the immune system at several levels. This has been documented in particular for the T-cell receptor (TCR)αβ+ T-cell compartment, showing reduced naive T-cell outputs and an accumulation of terminally differentiated clonally expanding effector T-cells, leading to increased proneness to autoimmunity and cancer development at older age. Even though TCRαβ+ and TCRγδ+ T-cells follow similar paths of development involving V(D)J-recombination of TCR genes in the thymus, TCRγδ+ T-cells tend to be more subjected to peripheral rather than central selection. However, the impact of aging in shaping of the peripheral TRG/TRD repertoire remains largely elusive. Next-generation sequencing analysis methods were optimized based on a spike-in method using plasmid vector DNA-samples for accurate TRG/TRD receptor diversity quantification, resulting in optimally defined primer concentrations, annealing temperatures and cycle numbers. Next, TRG/TRD repertoire diversity was evaluated during TCRγδ+ T-cell ontogeny, showing a broad, diverse repertoire in thymic and cord blood samples with Gaussian CDR3-length distributions, in contrast to the more skewed repertoire in mature circulating TCRγδ+ T-cells in adult peripheral blood. During aging the naive repertoire maintained its diversity with Gaussian CDR3-length distributions, while in the central and effector memory populations a clear shift from young (Vγ9/Vδ2 dominance) to elderly (Vγ2/Vδ1 dominance) was observed. Together with less clear Gaussian CDR3-length distributions, this would be highly suggestive of differentially heavily selected repertoires. Despite the apparent age-related shift from Vγ9/Vδ2 to Vγ2/Vδ1, no clear aging effect was observed on the Vδ2 invariant T nucleotide and canonical Vγ9-Jγ1.2 selection determinants. A more detailed look into the healthy TRG/TRD repertoire revealed known cytomegalovirus-specific TRG/TRD clonotypes in a few donors, albeit without a significant aging-effect, while Mycobacterium tuberculosis-specific clonotypes were absent. Notably, in effector subsets of elderly individuals, we could identify reported TRG and TRD receptor chains from TCRγδ+ T-cell large granular lymphocyte leukemia proliferations, which typically present in the elderly population. Collectively, our results point to relatively subtle age-related changes in the human TRG/TRD repertoire, with a clear shift in Vγ/Vδ usage in memory cells upon aging.

Immunodeficiency in Bloom's Syndrome

Bloom's syndrome (BS) is an autosomal recessive disease, caused by mutations in the BLM gene. This gene codes for BLM protein, which is a helicase involved in DNA repair. DNA repair is especially important for the development and maturation of the T and B cells. Since BLM is involved in DNA repair, we aimed to study if BLM deficiency affects T and B cell development and especially somatic hypermutation (SHM) and class switch recombination (CSR) processes. Clinical data of six BS patients was collected, and immunoglobulin serum levels were measured at different time points. In addition, we performed immune phenotyping of the B and T cells and analyzed the SHM and CSR in detail by analyzing IGHA and IGHG transcripts using next-generation sequencing. The serum immunoglobulin levels were relatively low, and patients had an increased number of infections. The absolute number of T, B, and NK cells were low but still in the normal range. Remarkably, all BS patients studied had a high percentage (20-80%) of CD4+ and CD8+ effector memory T cells. The process of SHM seems normal; however, the Ig subclass distribution was not normal, since the BS patients had more IGHG1 and IGHG3 transcripts. In conclusion, BS patients have low number of lymphocytes, but the immunodeficiency seems relatively mild since they have no severe or opportunistic infections. Most changes in the B cell development were seen in the CSR process; however, further studies are necessary to elucidate the exact role of BLM in CSR.

Human IgG2- and IgG4-expressing memory B cells display enhanced molecular and phenotypic signs of maturity and accumulate with age

The mechanisms involved in sequential immunoglobulin G (IgG) class switching are still largely unknown. Sequential IG class switching is linked to higher levels of somatic hypermutation (SHM) in vivo, but it remains unclear if these are generated temporally during an immune response or upon activation in a secondary response. We here aimed to uncouple these processes and to distinguish memory B cells from primary and secondary immune responses. SHM levels and IgG subclasses were studied with 454 pyrosequencing on blood mononuclear cells from young children and adults as models for primary and secondary immunological memory. Additional sequencing and detailed immunophenotyping with IgG subclass-specific antibodies was performed on purified IgG⁺ memory B-cell subsets. In both children and adults, SHM levels were higher in transcripts involving more downstream-located IGHG genes (esp. IGHG2 and IGHG4). In adults, SHM levels were significantly higher than in children, and downstream IGHG genes were more frequently utilized. This was associated with increased frequencies of CD27⁺IgG⁺ memory B cells, which contained higher levels of SHM, more IGHG2 usage, and higher expression levels of activation markers than CD27⁻IgG⁺ memory B cells. We conclude that secondary immunological memory accumulates with age and these memory B cells express CD27, high levels of activation markers, and carry high SHM levels and frequent usage of IGHG2. These new insights contribute to our understanding of sequential IgG subclass switching and show a potential relevance of using serum IgG2 levels or numbers of IgG2-expressing B cells as markers for efficient generation of memory responses.

CSReport: A New Computational Tool Designed for Automatic Analysis of Class Switch Recombination Junctions Sequenced by High-Throughput Sequencing

B cells ensure humoral immune responses due to the production of Ag-specific memory B cells and Ab-secreting plasma cells. In secondary lymphoid organs, Ag-driven B cell activation induces terminal maturation and Ig isotype class switch (class switch recombination [CSR]). CSR creates a virtually unique IgH locus in every B cell clone by intrachromosomal recombination between two switch (S) regions upstream of each C region gene. Amount and structural features of CSR junctions reveal valuable information about the CSR mechanism, and analysis of CSR junctions is useful in basic and clinical research studies of B cell functions. To provide an automated tool able to analyze large data sets of CSR junction sequences produced by high-throughput sequencing (HTS), we designed CSReport, a software program dedicated to support analysis of CSR recombination junctions sequenced with a HTS-based protocol (Ion Torrent technology). CSReport was assessed using simulated data sets of CSR junctions and then used for analysis of Sμ-Sα and Sμ-Sγ1 junctions from CH12F3 cells and primary murine B cells, respectively. CSReport identifies junction segment breakpoints on reference sequences and junction structure (blunt-ended junctions or junctions with insertions or microhomology). Besides the ability to analyze unprecedentedly large libraries of junction sequences, CSReport will provide a unified framework for CSR junction studies. Our results show that CSReport is an accurate tool for analysis of sequences from our HTS-based protocol for CSR junctions, thereby facilitating and accelerating their study.

High-Throughput Immunogenetics for Clinical and Research Applications in Immunohematology: Potential and Challenges

Analysis and interpretation of Ig and TCR gene rearrangements in the conventional, low-throughput way have their limitations in terms of resolution, coverage, and biases. With the advent of high-throughput, next-generation sequencing (NGS) technologies, a deeper analysis of Ig and/or TCR (IG/TR) gene rearrangements is now within reach, which impacts on all main applications of IG/TR immunogenetic analysis. To bridge the generation gap from low- to high-throughput analysis, the EuroClonality-NGS Consortium has been formed, with the main objectives to develop, standardize, and validate the entire workflow of IG/TR NGS assays for 1) clonality assessment, 2) minimal residual disease detection, and 3) repertoire analysis. This concerns the preanalytical (sample preparation, target choice), analytical (amplification, NGS), and postanalytical (immunoinformatics) phases. Here we critically discuss pitfalls and challenges of IG/TR NGS methodology and its applications in hemato-oncology and immunology.