Diversity, cellular origin and autoreactivity of antibody-secreting cell population expansions in acute systemic lupus erythematosus.

Acute systemic lupus erythematosus (SLE) courses with surges of antibody-secreting cells (ASCs) whose origin, diversity and contribution to serum autoantibodies remain unknown. Here, deep sequencing, proteomic profiling of autoantibodies and single-cell analysis demonstrated highly diversified ASCs punctuated by clones expressing the variable heavy-chain region VH4-34 that produced dominant serum autoantibodies. A fraction of ASC clones contained autoantibodies without mutation, a finding consistent with differentiation outside the germinal centers. A substantial ASC segment was derived from a distinct subset of newly activated naive cells of considerable clonality that persisted in the circulation for several months. Thus, selection of SLE autoreactivities occurred during polyclonal activation, with prolonged recruitment of recently activated naive B cells. Our findings shed light on the pathogenesis of SLE, help explain the benefit of agents that target B cells and should facilitate the design of future therapies.

Diversification of memory B cells drives the continuous adaptation of secretory antibodies to gut microbiota.

Secretory immunoglobulin A (SIgA) shields the gut epithelium from luminal antigens and contributes to host-microbe symbiosis. However, how antibody responses are regulated to achieve sustained host-microbe interactions is unknown. We found that mice and humans exhibited longitudinal persistence of clonally related B cells in the IgA repertoire despite major changes in the microbiota during antibiotic treatment or infection. Memory B cells recirculated between inductive compartments and were clonally related to plasma cells in gut and mammary glands. Our findings suggest that continuous diversification of memory B cells constitutes a central process for establishing symbiotic host-microbe interactions and offer an explanation of how maternal antibodies are optimized throughout life to protect the newborn.

MHC Class I-Dependent Shaping of the NK Cell Ly49 Receptor Repertoire Takes Place Early during Maturation in the Bone Marrow.

MHC class I (MHC I) expression in the host influences NK cells in a process termed education. The result of this education is reflected in the responsiveness of NK cells at the level of individual cells as well as in the repertoire of inhibitory MHC I-specific receptors at the NK cell system level. The presence of MHC I molecules in the host environment gives rise to a skewed receptor repertoire in spleen NK cells where subsets expressing few (one or two) inhibitory receptors are expanded whereas subsets with many (three or more) receptors are contracted. It is not known whether this MHC I-dependent skewing is imposed during development or after maturation of NK cells. In this study, we tested the hypothesis that the NK cell receptor repertoire is shaped already early during NK cell development in the bone marrow. We used mice with a repertoire imposed by a single MHC I allele, as well as a C57BL/6 mutant strain with exaggerated repertoire skewing, to investigate Ly49 receptor repertoires at different stages of NK cell differentiation. Our results show that NK cell inhibitory receptor repertoire skewing can indeed be observed in the bone marrow, even during the earliest developmental steps where Ly49 receptors are expressed. This may partly be accounted for by selective proliferation of certain NK cell subsets, but other mechanisms must also be involved. We propose a model for how repertoire skewing is established during a developmental phase in the bone marrow, based on sequential receptor expression as well as selective proliferation.

Murine peripheral NK-cell populations originate from site-specific immature NK cells more than from BM-derived NK cells.

Murine NK cells can be divided by the expression of two cell surface markers, CD27 and Mac-1 (a.k.a. CD11b), into four separate subsets. These subsets suggest a linear development model: CD27(-) Mac-1(-) → CD27(+) Mac-1(-) → CD27(+) Mac-1(+) → CD27(-) Mac-1(+) . Here, we used a combination of BrdU labeling experiments and mathematical modeling to gain insights regarding NK-cell development in mouse bone marrow (BM), spleen and liver. The modeling results that best fit the experimental data show that the majority of NK cells already express CD27 upon entering the NK-cell developmental pathway. Additionally, only a small fraction of NK cells exit the BM to other sites, suggesting that peripheral NK-cell populations originate from site-specific immature NK cells more than from BM-derived mature NK cells.

Aging affects B-cell antigen receptor repertoire diversity in primary and secondary lymphoid tissues.

The elderly immune system is characterized by reduced responses to infections and vaccines, and an increase in the incidence of autoimmune diseases and cancer. Age-related deficits in the immune system may be caused by peripheral homeostatic pressures that limit bone marrow B-cell production or migration to the peripheral lymphoid tissues. Studies of peripheral blood B-cell receptor spectratypes have shown that those of the elderly are characterized by reduced diversity, which is correlated with poor health status. In the present study, we performed for the first time high-throughput sequencing of immunoglobulin genes from archived biopsy samples of primary and secondary lymphoid tissues in old (74 ± 7 years old, range 61-89) versus young (24 ± 5 years old, range 18-45) individuals, analyzed repertoire diversities and compared these to results in peripheral blood. We found reduced repertoire diversity in peripheral blood and lymph node repertoires from old people, while in the old spleen samples the diversity was larger than in the young. There were no differences in somatic hypermutation characteristics between age groups. These results support the hypothesis that age-related immune frailty stems from altered B-cell homeostasis leading to narrower memory B-cell repertoires, rather than changes in somatic hypermutation mechanisms.

The V(H) repertoire and clonal diversification of B cells in inflammatory myopathies.

The contribution of antigen-driven B-cell adaptive immune responses within the inflamed muscle of inflammatory myopathies (IMs) is largely unknown. In this study, we investigated the immunoglobulin V(H) gene repertoire, somatic hypermutation, clonal diversification, and selection of infiltrating B cells in muscle biopsies from IM patients (dermatomyositis and polymyositis), to determine whether B cells and/or plasma cells contribute to the associated pathologies of these diseases. The data reveal that Ig V(H) gene repertoires of muscle-infiltrating B cells deviate from the normal V(H) gene repertoire in individual patients, and differ between different types of IMs. Analysis of somatic mutations revealed clonal diversification of muscle-infiltrating B cells and evidence for a chronic B-cell response within the inflamed muscle. We conclude that muscle-infiltrating B cells undergo selection, somatic hypermutation and clonal diversification in situ during antigen-driven immune responses in patients with IMs, providing insight into the contribution of B cells to the pathological mechanisms of these disorders.

Classification of human natural killer cells based on migration behavior and cytotoxic response.

Despite intense scrutiny of the molecular interactions between natural killer (NK) and target cells, few studies have been devoted to dissection of the basic functional heterogeneity in individual NK cell behavior. Using a microchip-based, time-lapse imaging approach allowing the entire contact history of each NK cell to be recorded, in the present study, we were able to quantify how the cytotoxic response varied between individual NK cells. Strikingly, approximately half of the NK cells did not kill any target cells at all, whereas a minority of NK cells was responsible for a majority of the target cell deaths. These dynamic cytotoxicity data allowed categorization of NK cells into 5 distinct classes. A small but particularly active subclass of NK cells killed several target cells in a consecutive fashion. These "serial killers" delivered their lytic hits faster and induced faster target cell death than other NK cells. Fast, necrotic target cell death was correlated with the amount of perforin released by the NK cells. Our data are consistent with a model in which a small fraction of NK cells drives tumor elimination and inflammation.

Frequency and phenotype of B cell subpopulations in young and aged HIV-1 infected patients receiving ART.

BACKGROUND: Aged individuals respond poorly to vaccination and have a higher risk of contracting infections in comparison to younger individuals; whether age impacts on the composition and function of B cell subpopulations relevant for immune responses is still controversial. It is also not known whether increased age during HIV-1 infection further synergizes with the virus to alter B cell subpopulations. In view of the increased number of HIV-1 infected patients living to high age as a result of anti-retroviral treatment this is an important issue to clarify. RESULTS: In this work we evaluated the distribution of B cell subpopulations in young and aged healthy individuals and HIV-1 infected patients, treated and naïve to treatment. B cell populations were characterized for the expression of inhibitory molecules (PD-1 and FcRL4) and activation markers (CD25 and CD69); the capacity of B cells to respond to activation signals through up-regulation of IL-6 expression was also evaluated. Increased frequencies of activated and tissue-like memory B cells occurring during HIV-1 infection are corrected by prolonged ART therapy. Our findings also reveal that, in spite of prolonged treatment, resting memory B cells in both young and aged HIV-1 infected patients are reduced in number, and all memory B cell subsets show a low level of expression of the activation marker CD25. CONCLUSIONS: The results of our study show that resting memory B cells in ART-treated young and aged HIV-1 infected patients are reduced in number and memory B cell subsets exhibit low expression of the activation marker CD25. Aging per se in the HIV-1 infected population does not worsen impairments initiated by HIV-1 in the memory B cell populations of young individuals.

Lyn deficiency affects B-cell maturation as well as survival.

Lyn, an Src-family protein tyrosine kinase expressed in B lymphocytes, contributes to initiation of BCR signaling and is also responsible for feedback inhibition of BCR signaling. Lyn-deficient mice have a decreased number of follicular B cells and also spontaneously develop a lupus-like autoimmunity. We used flow cytometric analysis, BrdU labeling and our mathematical models of B-cell population dynamics, to analyze how Lyn deficiency impacts B-cell maturation and survival. We found that Lyn-deficient transitional 1 (T1) cells develop normally, but T2 cells develop primarily from the T1 subset in the spleen and fail to also develop directly from BM immature B cells. Lyn-deficient T2 cells either mature to the follicular B-cell type at a close to normal rate, or die in this compartment rather than access the T3 anergic subset. The ≈ 40% of WT follicular cells that were short-lived exited primarily by joining the T3 anergic subset, whereas the ≈ 15% Lyn(-/-) follicular cells that were not long lived had a high death rate and died in this compartment rather than entering the T3 subset. We hypothesize that exaggerated BCR signaling resulting from weak interactions with self-antigens is largely responsible for these alterations in Lyn-deficient B cells.

BCR CDR3 length distributions differ between blood and spleen and between old and young patients, and TCR distributions can be used to detect myelodysplastic syndrome.

Complementarity-determining region 3 (CDR3) is the most hyper-variable region in B cell receptor (BCR) and T cell receptor (TCR) genes, and the most critical structure in antigen recognition and thereby in determining the fates of developing and responding lymphocytes. There are millions of different TCR Vß chain or BCR heavy chain CDR3 sequences in human blood. Even now, when high-throughput sequencing becomes widely used, CDR3 length distributions (also called spectratypes) are still a much quicker and cheaper method of assessing repertoire diversity. However, distribution complexity and the large amount of information per sample (e.g. 32 distributions of the TCRα chain, and 24 of TCRß) calls for the use of machine learning tools for full exploration. We have examined the ability of supervised machine learning, which uses computational models to find hidden patterns in predefined biological groups, to analyze CDR3 length distributions from various sources, and distinguish between experimental groups. We found that (a) splenic BCR CDR3 length distributions are characterized by low standard deviations and few local maxima, compared to peripheral blood distributions; (b) healthy elderly people's BCR CDR3 length distributions can be distinguished from those of the young; and (c) a machine learning model based on TCR CDR3 distribution features can detect myelodysplastic syndrome with approximately 93% accuracy. Overall, we demonstrate that using supervised machine learning methods can contribute to our understanding of lymphocyte repertoire diversity.

PCR amplification and high throughput sequencing of immunoglobulin heavy chain genes from formalin-fixed paraffin-embedded human biopsies.

The use of high throughput sequencing (HTS) technologies in biomedicine is expanding in a variety of fields in recent years. The 454 system is an HTS platform that is ideally suited to characterize B cell receptor (BCR) repertoires by sequencing of immunoglobulin (Ig) genes, as it is able to sequence stretches of several hundred nucleotides. Most studies that used this platform for antibody repertoire analyses have started from fresh or frozen tissues or peripheral blood samples, and rely on starting with optimal quality DNA. In this paper we demonstrate that BCR repertoire analysis can be done using DNA from formalin-fixed paraffin-embedded (FFPE) human tissue samples. The heterogeneity of BCR repertoires we obtained confirms the plausibility of HTS of DNA from FFPE specimens. The establishment of experimental protocols and computational tools that enable sequence data analysis from the low quality DNA of FFPE tissues is important for enabling research, as it would enable the use of the rich source of preserved samples in clinical biobanks and biopsy archives.

Simulations of the NK cell immune synapse reveal that activation thresholds can be established by inhibitory receptors acting locally.

NK cell activation is regulated by a balance between activating and inhibitory signals. To address the question of how these signals are spatially integrated, we created a computer simulation of activating and inhibitory NK cell immunological synapse (NKIS) assembly, implementing either a "quantity-based" inhibition model or a "distance-based" inhibition model. The simulations mimicked the observed molecule distributions in inhibitory and activating NKIS and yielded several new insights. First, the total signal is highly influenced by activating complex dissociation rates but not by adhesion and inhibitory complex dissociation rates. Second, concerted motion of receptors in clusters significantly accelerates NKIS maturation. Third, when the potential of a cis interaction between Ly49 receptors and MHC class I on murine NK cells was added to the model, the integrated signal as a function of receptor and ligand numbers was only slightly increased, at least up to the level of 50% cis-bound Ly49 receptors reached in the model. Fourth, and perhaps most importantly, the integrated signal behavior obtained when using the distance-based inhibition signal model was closer to the experimentally observed behavior, with an inhibition radius of the order 3-10 molecules. Microscopy to visualize Vav activation in NK cells on micropatterned surfaces of activating and inhibitory strips revealed that Vav is only locally activated where activating receptors are ligated within a single NK cell contact. Taken together, these data are consistent with a model in which inhibitory receptors act locally; that is, that every bound inhibitory receptor acts on activating receptors within a certain radius around it.

Chronic B cell deficiency from birth prevents age-related alterations in the B lineage.

Aging is accompanied by a decline in B lymphopoiesis in the bone marrow and accumulation of long-lived B cells in the periphery. The mechanisms underlying these changes are unclear. To explore whether aging in the B lineage is subjected to homeostatic regulation, we used mutant mice bearing chronic B cell deficiency from birth. We show that chronic B cell deficiency from birth, resulting from impaired maturation (CD19(-/-) and CD74(-/-)) or reduced survival (baff-r(-/-)), prevents age-related changes in the B lineage. Thus, frequencies of early and late hematopoietic stem cells, B lymphopoiesis, and the rate of B cell production do not substantially change with age in these mice, as opposed to wild-type mice where kinetic experiments indicate that the output from the bone marrow is impaired. Further, we found that long-lived B cells did not accumulate and peripheral repertoire was not altered with age in these mice. Collectively, our results suggest that aging in the B lineage is not autonomously progressing but subjected to homeostatic regulation.

[B lymphocyte clonal evolution of human reactive lymph nodes revealed by lineage tree analysis].

INTRODUCTION: Hypermutation and selection processes, characterizing T-dependent B cell responses taking place in germinal centers of lymph nodes, lead to B cell receptor affinity maturation. Those immune responses lead to the development of memory B cells and plasma cells that secrete high amounts of antibody molecules. The dynamics of B cell clonal evolution during affinity maturation has significant importance in infectious and autoimmune diseases, malignancies and aging. Immunoglobulin (Ig) gene mutational Lineage tree construction by comparing variable regions of Ig-gene sequences to the Ig germline gene is an interesting approach for studying B cell cLonal evolution. Lineage tree shapes and Ig gene mutations can be evaluated not only qualitatively and intuitively, but also quantitatively, and thus reveal important information related to hypermutation and selection. AIM: In this paper we describe the experimental protocols that we used for PCR amplification of Igvariable region genes from human formalin fixed paraffin embedded reactive lymph node tissues and the subsequent bioinformatical analyses of sequencing data using Ig mutational lineage trees. RESULTS: B cell populations of three out of four reactive Lymph node tissues were composed of several clones. Most of the Ig gene mutational lineage trees were small and narrow. Significant differences were not detected by quantification of Lineage trees. SUMMARY: B lymphocyte clones that were detected in human reactive lymph node tissues represent major responding clones in normal polyclonal immune response. This result is in line with the polyclonal profile of B Lymphocyte populations that reside in reactive lymph node tissues.

B cell M-CLL clones retain selection against replacement mutations in their immunoglobulin gene framework regions.

Introduction: Chronic lymphocytic leukemia (CLL) is the most common adult leukemia, accounting for 30-40% of all adult leukemias. The dynamics of B-lymphocyte CLL clones with mutated immunoglobulin heavy chain variable region (IgHV) genes in their tumor (M-CLL) can be studied using mutational lineage trees. Methods: Here, we used lineage tree-based analyses of somatic hypermutation (SHM) and selection in M-CLL clones, comparing the dominant (presumably malignant) clones of 15 CLL patients to their non-dominant (presumably normal) B cell clones, and to those of healthy control repertoires. This type of analysis, which was never previously published in CLL, yielded the following novel insights. Results: CLL dominant clones undergo - or retain - more replacement mutations that alter amino acid properties such as charge or hydropathy. Although, as expected, CLL dominant clones undergo weaker selection for replacement mutations in the complementarity determining regions (CDRs) and against replacement mutations in the framework regions (FWRs) than non-dominant clones in the same patients or normal B cell clones in healthy controls, they surprisingly retain some of the latter selection in their FWRs. Finally, using machine learning, we show that even the non-dominant clones in CLL patients differ from healthy control clones in various features, most notably their expression of higher fractions of transition mutations. Discussion: Overall, CLL seems to be characterized by significant loosening - but not a complete loss - of the selection forces operating on B cell clones, and possibly also by changes in SHM mechanisms.

Gastric DLBCL clonal evolution as function of patient age.

Diffuse large B cell lymphoma (DLBCL) is the most common type of NHL, accounting for about 40% of NHL cases, and is one of the most aggressive lymphomas. DLBCL is widespread in individuals aged more than 50 years old, with a maximum incidence in the seventh decade, but it may also occur in younger patients. DLBCL may occur in any immune system tissue, including those around the gastrointestinal tract, and even in the stomach, though gastric DLBCL has yet to be sufficiently investigated. This study aimed to understand changes in gastric Diffuse Large B cell lymphoma (gastric DLBCL) development with age. Immunoglobulin (Ig) heavy chain variable region genes were amplified from sections of nine preserved biopsies, from patients whose age varied between 25 and 89 years, sequenced and analyzed. We show first that identification of the malignant clone based on the biopsies is much less certain than was previously assumed; and second that, contrary to expectations, the repertoire of gastric B cell clones is more diverse among the elderly DLBCL patients than among the young.

Mutational patterns along different evolution paths of follicular lymphoma.

Follicular lymphoma (FL) is an indolent disease, characterized by a median life expectancy of 18-20 years and by intermittent periods of relapse and remission. FL frequently transforms into the more aggressive diffuse large B cell lymphoma (t-FL). In previous studies, the analysis of immunoglobulin heavy chain variable region (IgHV) genes in sequential biopsies from the same patient revealed two different patterns of tumor clonal evolution: direct evolution, through acquisition of additional IgHV mutations over time, or divergent evolution, in which lymphoma clones from serial biopsies independently develop from a less-mutated common progenitor cell (CPC). Our goal in this study was to characterize the somatic hypermutation (SHM) patterns of IgHV genes in sequential FL samples from the same patients, and address the question of whether the mutation mechanisms (SHM targeting, DNA repair or both), or selection forces acting on the tumor clones, were different in FL samples compared to healthy control samples, or in late relapsed/transformed FL samples compared to earlier ones. Our analysis revealed differences in the distribution of mutations from each of the nucleotides when tumor and non-tumor clones were compared, while FL and transformed FL (t-FL) tumor clones displayed similar mutation distributions. Lineage tree measurements suggested that either initial clone affinity or selection thresholds were lower in FL samples compared to controls, but similar between FL and t-FL samples. Finally, we observed that both FL and t-FL tumor clones tend to accumulate larger numbers of potential N-glycosylation sites due to the introduction of new SHM. Taken together, these results suggest that transformation into t-FL, in contrast to initial FL development, is not associated with any major changes in DNA targeting or repair, or the selection threshold of the tumor clone.

IgTreeZ, A Toolkit for Immunoglobulin Gene Lineage Tree-Based Analysis, Reveals CDR3s Are Crucial for Selection Analysis.

Somatic hypermutation (SHM) is an important diversification mechanism that plays a part in the creation of immune memory. Immunoglobulin (Ig) variable region gene lineage trees were used over the last four decades to model SHM and the selection mechanisms operating on B cell clones. We hereby present IgTreeZ (Immunoglobulin Tree analyZer), a python-based tool that analyses many aspects of Ig gene lineage trees and their repertoires. Using simulations, we show that IgTreeZ can be reliably used for mutation and selection analyses. We used IgTreeZ on empirical data, found evidence for different mutation patterns in different B cell subpopulations, and gained insights into antigen-driven selection in corona virus disease 19 (COVID-19) patients. Most importantly, we show that including the CDR3 regions in selection analyses - which is only possible if these analyses are lineage tree-based - is crucial for obtaining correct results. Overall, we present a comprehensive lineage tree analysis tool that can reveal new biological insights into B cell repertoire dynamics.

Natural killer cell education in mice with single or multiple major histocompatibility complex class I molecules.

The ability of murine NK cells to reject cells lacking self MHC class I expression results from an in vivo education process. To study the impact of individual MHC class I alleles on this process, we generated mice expressing single MHC class I alleles (K(b), D(b), D(d), or L(d)) or combinations of two or more alleles. All single MHC class I mice rejected MHC class I-deficient cells in an NK cell-dependent way. Expression of K(b) or D(d) conveyed strong rejection of MHC class I-deficient cells, whereas the expression of D(b) or L(d) resulted in weaker responses. The educating impact of weak ligands (D(b) and L(d)) was further attenuated by the introduction of additional MHC class I alleles, whereas strong ligands (K(b) and D(d)) maintained their educating impact under such conditions. An analysis of activating and inhibitory receptors in single MHC class I mice suggested that the educating impact of a given MHC class I molecule was controlled both by the number of NK cells affected and by the strength of each MHC class I-Ly49 receptor interaction, indicating that NK cell education may be regulated by a combination of qualitative and quantitative events.

Ectopic GC in the thymus of myasthenia gravis patients show characteristics of normal GC.

Young patients with myasthenia gravis (MG) frequently have ectopic GC in their thymus. We investigated these ectopic GC by microdissection of GC B cells and analysis of their Ig gene characteristics, in comparison to normal GC. CDR3 length distribution, a measure of clonal variability, and Ig gene family usage were similar in MG and normal tonsil samples. Lineage tree analysis demonstrated similar diversification and mutations per cell compared with normal control trees. Mutations were observed in the framework regions, responsible for the structural integrity of the BCR; however, these mutations were mostly conservative or neutral, confirming that a functional BCR is conserved in MG. In the CDR, responsible for Ag binding, selection against replacement mutations was revealed. This may indicate that the MG clones analyzed are already highly Ag-specific, and therefore potential affinity-reducing replacement mutations in the CDR3 are not propagated, due to Ag-driven selection. Somatic hypermutation (SHM) targeting motifs and aa substitution preferences in MG were similar to those of normal controls. Overall, these results suggest that B cells in the ectopic GC in MG appear to undergo normal diversification and selection, in spite of the chronic nature and different environment of the response.