Ig gene diversification and selection in follicular lymphoma, diffuse large B cell lymphoma and primary central nervous system lymphoma revealed by lineage tree and mutation analyses

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.

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.

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.

Grade I, II and III Follicular Lymphomas Express Ig VH Genes with Different Patterns of Somatic Mutation

Follicular lymphoma (FL) is an indolent, B-cell, non-Hodgkin’s lymphoma with varying cytological appearance and clinical behavior. The genetic hallmark of FL is the t(14;18) translocation, and as a germinal center derived entity it is also characterized by somatic hypermutation of the immunoglobulin heavy chain (IgH) gene. In an attempt to correlate this molecular signature with the cytological grading of FL, we have analyzed the IgH variable (IgV_H), regions in all cytological grades of FL. Four FL cases showing t(14;18) translocation were classified into grade I-III categories according to the current WHO guidelines. The IgV_H gene segments were PCR-amplified, sequenced, and compared to their respective germline IgV_H sequences. The neoplastic cells of grade I and II FLs revealed clonally related, but highly divergent IgV_H gene sequences indicating the ongoing nature of somatic hypermutation. Grade III FL also showed extensive presence of somatic hypermutation, but these mutations were not associated with intraclonal divergence. Thus, these results suggest that grade I-II and grade III FL may represent different biological entities. The presence of ongoing somatic hypermutation of IgV_H sequences in grade I and II FLs is compatible with direct follicular origin of these tumor cells, contrasting the homogenous, stable clones of grade III FL resembling a post-follicular stage of B-cell development. Our findings demonstrate that contrary to the three tiered cytological grading, molecular features of IgH genes classify FL into two distinct subcategories. These studies also suggest that with progression FL gains post-follicular–like molecular features and becomes independent of the germinal center microenvironment.

Chemotherapy Curability in Leukemia, Lymphoma, Germ Cell Tumors and Gestational Malignancies: A Reflection of the Unique Physiology of Their Cells of Origin

Cytotoxic DNA damaging chemotherapy brings clinical benefits in the treatment of many metastatic malignancies. However routine curative treatment remains restricted to a small number of malignancies including acute leukemia, high grade lymphoma, germ cell tumors, gestational malignancies and some of the rare childhood cancers. The detailed explanation for this dramatic divergence in outcomes remains to be elucidated. However, we have previously argued that there is a strong correlation between presence of the unique genetic events of immunoglobulin gene variable/diversity/joining (VDJ) recombination, somatic hypermutation (SHM), meiosis, nuclear fusion and gastrulation occurring in cells of origin of these malignancies and their high sensitivity to DNA damaging chemotherapy. In this study we have reviewed some of the basic physiological information relating to the specialized activity and sensitivity to DNA damage mediated apoptosis of normal cells undergoing these processes. In each of unique genetic events there are dramatic changes in apoptotic sensitivity. In VDJ recombination and somatic hypermutation over 95% of the cells involved undergo apoptosis, whilst in meiosis and nuclear fusion there are dramatic short term increases in the apoptotic sensitivity to DNA damage. It is apparent that each of the malignancies arising during these processes retains some of the unique phenotype associated with it. The impact of the physiological differences is most clearly seen in the two non-mutational malignancies. Gestational choriocarcinoma which arises shortly after nuclear fusion is routinely curable with chemotherapy whilst CIMP-positive ependymomas which is not linked to any of the unique genetic events is highly resistant. A similar pattern is found in a pair of malignancies driven by a single driver mutation. Infantile acute lymphoblastic leukemia (ALL) arises in a cell undergoing the early stages of VDJ recombination and has a 40% cure rate in contrast pediatric rhabdoid malignancy which is not linked to a unique genetic event responds very poorly to chemotherapy treatment. The physiological changes occurring in cancer cells at the time of the malignant transformation appear to have a major impact on the subsequent sensitivity to chemotherapy and curability. New therapies that impact on these pathways may be of therapeutic value.

Immunoglobulin gene rearrangement IGHV3-48 is a predictive marker of histological transformation into aggressive lymphoma in follicular lymphomas

Follicular lymphoma (FL) is a heterogeneous disease whose pathogenesis remains partially unknown. Around 20% of FL patients experience early progression or treatment-refractory disease and 2–3% of patients per year experience histological transformation (HT) into a more aggressive lymphoma (tFL). Here, we evaluate the immunoglobulin heavy chain variable (IGHV) gene usage and mutational status in 187 FL cases to assess its impact on clinical outcome and histological transformation. The IGHV gene repertoire was remarkably biased in FL. The IGHV4-34 (14%), IGHV3-23 (14%), IGHV3-48 (10%), IGHV3-30 (9%) and IGHV3-21 (7%) genes accounted for more than half of the whole cohort. IGHV3-48 was overrepresented in cases of tFL (19%) compared with non-transformed FL at 5 years (5%, P = 0.05). Patients with the IGHV3-48 gene were significantly more likely to have had HT after 10 years than those who used other genes (71% vs. 25%, P < 0.05), irrespective of the therapy they received. Moreover, IGHV3-30 was also overrepresented in cases of FL (9%) and tFL (13%) compared with diffuse large B-cell lymphoma in which it was nearly absent. In conclusion, our results indicate a role for antigen selection in the development of FL, while the use of IGHV3-48 could help predict histological transformation.

Computational Evaluation of B-Cell Clone Sizes in Bulk Populations

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

From genetics to the clinic: a translational perspective on follicular lymphoma

Follicular lymphoma (FL) is the most frequent indolent B cell lymphoma and is still considered to be incurable. In recent years, whole-exome sequencing studies of large cohorts of patients have greatly improved our knowledge of the FL mutational landscape. Moreover, the prolonged evolution of this disease has enabled some insights regarding the early pre-lymphoma lesions as well as the clonal evolution after treatment, allowing an evolutionary perspective on lymphomagenesis. Deciphering the earliest initiating lesions and identifying the molecular alterations leading to disease progression currently represent important goals; accomplishing these could help identify the most relevant targets for precision therapy.

Altered N-Linked Glycosylation in Follicular Lymphoma and Chronic Lymphocytic Leukemia: Involvement in Pathogenesis and Potential Therapeutic Targeting

B-cell antigen receptor (BCR) expression is indispensable for survival of most B-cell malignancies. In follicular lymphoma (FL), N-linked glycosylation sites are introduced in the immunoglobulin (Ig) variable region genes. Oligosaccharides added to the acquired sites are unusually of the high-mannose type. These glycans interact with mannose-specific lectins, especially with dendritic cell-specific intercellular adhesion molecule-3-grabbing non-integrin (DC-SIGN). Lectin binding to FL triggers persistent activating signals, suggesting that lectins within the tumor microenvironment promote cell survival and proliferation. Insertion of N-glycosylation sites in Ig variable region genes has been detected in other germinal center-associated lymphomas, specifically in subsets of diffuse large B-cell lymphomas and Burkitt’s lymphomas, suggesting involvement of altered glycans in pathogenesis of these malignancies as well. Furthermore, the BCR in chronic lymphocytic leukemia (CLL) carries high-mannose oligosaccharides, albeit in the heavy chain constant rather than variable region. The high expression level of the unique glycoform, particularly in the more aggressive unmutated CLL subset, suggests a functional significance for this glycan in CLL. As lectin interaction with the BCR is critical for FL and probably for some other lymphomas, targeting this interaction is considered to be an interesting therapeutic strategy. Reagents for blockade of lectin–BCR interaction may include antibodies against high-mannose glycans and mannose-based oligosaccharide mimics or non-carbohydrate glycomimetics. Moreover, as this interaction triggers signaling pathways similar to those demonstrated for BCR engagement by antigen, BCR signal transduction inhibitors may emerge as effective therapeutics for lectin-driven malignancies.

Next generation sequencing of the clonal IGH rearrangement detects ongoing mutations and interfollicular trafficking in in situ follicular neoplasia

Follicular lymphoma (FL) is characterized genetically by a significant intraclonal diversity of rearranged immunoglobulin heavy chain (IGH) genes and a substantial cell migration activity (follicular trafficking). Recently, in situ follicular neoplasia (ISFN), characterized by accumulations of immunohistochemically strongly BCL2-positive, t(14;18)+ clonal B cells confined to germinal centers in reactive lymph nodes, has been identified as a precursor lesion of FL with low risk of progression to manifest FL. The extent of ongoing somatic hypermutation of rearranged IGH genes and interfollicular trafficking in ISFN is not known. In this study we performed an in depth analysis of clonal evolution and cell migration patterns in a case of pure ISFN involving multiple lymph nodes. Using laser microdissection and next generation sequencing (NGS) we documented significant intraclonal diversity of the rearranged IGH gene and extensive interfollicular migration between germinal centers of the same lymph node as well as between different lymph nodes. Furthermore, we identified N-glycosylation motifs characteristic for FL in the CDR3 region.

Oncogenic events rather than antigen selection pressure may be the main driving forces for relapse in diffuse large B-cell lymphomas

Little is known on the phylogenetic relationship between diagnostic and relapse clones of diffuse large B-cell lymphoma (DLBCL). We applied high throughput sequencing (HTS) of the VDJ locus of Immunoglobulin heavy chain (IGHV) on 14 DLBCL patients with serial samples, including tumor biopsies and/or peripheral blood mononuclear cells (PBMC). Phylogenetic data were consolidated with targeted sequencing and cytogenetics. Phylogeny clearly showed that DLBCL relapse could occur according either an early or a late divergent mode. These two modes of divergence were independent from the elapsed time between diagnosis and relapse. We found no significant features for antigen selection pressure in complementary determining region both at diagnosis and relapse for 9/12 pairs and a conserved negative selection pressure for the three remaining cases. Targeted HTS and conventional cytogenetics revealed a branched vs. linear evolution for 5/5 IGHV early divergent cases, but unexpected such "oncogenetic" branched evolution could be found in at least 2/7 IGHV late divergent cases. Thus, if BCR signaling is mandatory for DLBCL emergence, oncogenetic events under chemotherapy selection pressure may be the main driving forces at relapse. Finally, circulating subclones with divergent IGHV somatic hypermutations patterns from initial biopsy could be detected in PBMC at diagnosis for 4/6 patients and, for two of them, at least one was similar to the ones found at relapse. This study highlights that oncogenetic intraclonal diversity of DLBCL should be evaluated beyond the scope a single biopsy and represents a rationale for future investigations using peripheral blood for lymphoid malignancies genotyping. Am. J. Hematol. 92:68-76, 2017. © 2016 Wiley Periodicals, Inc.

Targeting B-cell receptor signaling in leukemia and lymphoma: how and why?

B-lymphocytes are dependent on B-cell receptor (BCR) signaling for the constant maintenance of their physiological function, and in many B-cell malignancies this signaling pathway is prone to aberrant activation. This understanding has led to an ever-increasing interest in the signaling networks activated following ligation of the BCR in both normal and malignant cells, and has been critical in establishing an array of small molecule inhibitors targeting BCR-induced signaling. By dissecting how different malignancies signal through BCR, researchers are contributing to the design of more customized therapeutics which have greater efficacy and lower toxicity than previous therapies. This allows clinicians access to an array of approaches to best treat patients whose malignancies have BCR signaling as a driver of pathogenesis.

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.

DC-SIGN-expressing macrophages trigger activation of mannosylated IgM B-cell receptor in follicular lymphoma

Follicular lymphoma (FL) results from the accumulation of malignant germinal center (GC) B cells leading to the development of an indolent and largely incurable disease. FL cells remain highly dependent on B-cell receptor (BCR) signaling and on a specific cell microenvironment, including T cells, macrophages, and stromal cells. Importantly, FL BCR is characterized by a selective pressure to retain surface immunoglobulin M (IgM) BCR despite an active class-switch recombination process, and by the introduction, in BCR variable regions, of N-glycosylation acceptor sites harboring unusual high-mannose oligosaccharides. However, the relevance of these 2 FL BCR features for lymphomagenesis remains unclear. In this study, we demonstrated that IgM(+) FL B cells activated a stronger BCR signaling network than IgG(+) FL B cells and normal GC B cells. BCR expression level and phosphatase activity could both contribute to such heterogeneity. Moreover, we underlined that a subset of IgM(+) FL samples, displaying highly mannosylated BCR, efficiently bound dendritic cell-specific intercellular adhesion molecule-3-grabbing nonintegrin (DC-SIGN), which could in turn trigger delayed but long-lasting BCR aggregation and activation. Interestingly, DC-SIGN was found within the FL cell niche in situ. Finally, M2 macrophages induced a DC-SIGN-dependent adhesion of highly mannosylated IgM(+) FL B cells and triggered BCR-associated kinase activation. Interestingly, pharmacologic BCR inhibitors abolished such crosstalk between macrophages and FL B cells. Altogether, our data support an important role for DC-SIGN-expressing infiltrating cells in the biology of FL and suggest that they could represent interesting therapeutic targets.

Lectins from opportunistic bacteria interact with acquired variable-region glycans of surface immunoglobulin in follicular lymphoma

B-cell antigen receptor (BCR) expression is a key feature of most B-cell lymphomas, but the mechanisms of BCR signal induction and the involvement of autoantigen recognition remain unclear. In follicular lymphoma (FL) B cells, BCR expression is retained despite a chromosomal translocation that links the antiapoptotic gene BCL2 to the regulatory elements of immunoglobulin genes, thereby disrupting 1 heavy-chain allele. A remarkable feature of FL-BCRs is the acquisition of potential N-glycosylation sites during somatic hypermutation. The introduced glycans carry mannose termini, which create potential novel binding sites for mannose-specific lectins. Here, we investigated the effect of N-linked variable-region glycosylation for BCR interaction with cognate antigen and with lectins of different origins. N-glycans were found to severely impair BCR specificity and affinity to the initial cognate antigen. In addition, we found that lectins from Pseudomonas aeruginosa and Burkholderia cenocepacia bind and stimulate FL cells. Human exposure to these bacteria can occur by contact with soil and water. In addition, they represent opportunistic pathogens in susceptible hosts. Understanding the role of bacterial lectins might elucidate the pathogenesis of FL and establish novel therapeutic approaches.

Automated analysis of high-throughput B-cell sequencing data reveals a high frequency of novel immunoglobulin V gene segment alleles

Individual variation in germline and expressed B-cell immunoglobulin (Ig) repertoires has been associated with aging, disease susceptibility, and differential response to infection and vaccination. Repertoire properties can now be studied at large-scale through next-generation sequencing of rearranged Ig genes. Accurate analysis of these repertoire-sequencing (Rep-Seq) data requires identifying the germline variable (V), diversity (D), and joining (J) gene segments used by each Ig sequence. Current V(D)J assignment methods work by aligning sequences to a database of known germline V(D)J segment alleles. However, existing databases are likely to be incomplete and novel polymorphisms are hard to differentiate from the frequent occurrence of somatic hypermutations in Ig sequences. Here we develop a Tool for Ig Genotype Elucidation via Rep-Seq (TIgGER). TIgGER analyzes mutation patterns in Rep-Seq data to identify novel V segment alleles, and also constructs a personalized germline database containing the specific set of alleles carried by a subject. This information is then used to improve the initial V segment assignments from existing tools, like IMGT/HighV-QUEST. The application of TIgGER to Rep-Seq data from seven subjects identified 11 novel V segment alleles, including at least one in every subject examined. These novel alleles constituted 13% of the total number of unique alleles in these subjects, and impacted 3% of V(D)J segment assignments. These results reinforce the highly polymorphic nature of human Ig V genes, and suggest that many novel alleles remain to be discovered. The integration of TIgGER into Rep-Seq processing pipelines will increase the accuracy of V segment assignments, thus improving B-cell repertoire analyses.

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

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

Immunoglobulin gene repertoire diversification and selection in the stomach - from gastritis to gastric lymphomas

Chronic gastritis is characterized by gastric mucosal inflammation due to autoimmune responses or infection, frequently with Helicobacter pylori. Gastritis with H. pylori background can cause gastric mucosa-associated lymphoid tissue lymphoma (MALT-L), which sometimes further transforms into diffuse large B-cell lymphoma (DLBCL). However, gastric DLBCL can also be initiated de novo. The mechanisms underlying transformation into DLBCL are not completely understood. We analyzed immunoglobulin repertoires and clonal trees to investigate whether and how immunoglobulin gene repertoires, clonal diversification, and selection in gastritis, gastric MALT-L, and DLBCL differ from each other and from normal responses. The two gastritis types (positive or negative for H. pylori) had similarly diverse repertoires. MALT-L dominant clones (defined as the largest clones in each sample) presented higher diversification and longer mutational histories compared with all other conditions. DLBCL dominant clones displayed lower clonal diversification, suggesting the transforming events are triggered by similar responses in different patients. These results are surprising, as we expected to find similarities between the dominant clones of gastritis and MALT-L and between those of MALT-L and DLBCL.

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.

Integrating B cell lineage information into statistical tests for detecting selection in Ig sequences

Detecting selection in B cell Ig sequences is critical to understanding affinity maturation and can provide insights into Ag-driven selection in normal and pathologic immune responses. The most common sequence-based methods for detecting selection analyze the ratio of replacement and silent mutations using a binomial statistical analysis. However, these approaches have been criticized for low sensitivity. An alternative method is based on the analysis of lineage trees constructed from sets of clonally related Ig sequences. Several tree shape measures have been proposed as indicators of selection that can be statistically compared across cohorts. However, we show that tree shape analysis is confounded by underlying experimental factors that are difficult to control for in practice, including the sequencing depth and number of generations in each clone. Thus, although lineage tree shapes may reflect selection, their analysis alone is an unreliable measure of in vivo selection. To usefully capture the information provided by lineage trees, we propose a new method that applies the binomial statistical framework to mutations identified based on lineage tree structure. This hybrid method is able to detect selection with increased sensitivity in both simulated and experimental data sets. We anticipate that this approach will be especially useful in the analysis of large-scale Ig sequencing data sets generated by high-throughput sequencing technologies.

B-cell receptor signaling as a driver of lymphoma development and evolution

The B-cell receptor (BCR) is essential for normal B-cell development and maturation. In an increasing number of B-cell malignancies, BCR signaling is implicated as a pivotal pathway in tumorigenesis. Mechanisms of BCR activation are quite diverse and range from chronic antigenic drive by microbial or viral antigens to autostimulation of B-cells by self-antigens to activating mutations in intracellular components of the BCR pathway. Hepatitis C virus infection can lead to the development of splenic marginal zone lymphoma, while Helicobacter pylori infection is associated with the development of mucosa-associated lymphoid tissue lymphomas. In some of these cases, successful treatment of the infection removes the inciting antigen and results in resolution of the lymphoma. Chronic lymphocytic leukemia has been recognized for decades as a malignancy of auto-reactive B-cells and its clinical course is in part determined by the differential response of the malignant cells to BCR activation. In a number of B-cell malignancies, activating mutations in signal transduction components of the BCR pathway have been identified; prominent examples are activated B-cell-like (ABC) diffuse large B-cell lymphomas (DLBCL) that carry mutations in CD79B and CARD11 and display chronic active BCR signaling resulting in constitutive activation of the NF-κB pathway. Despite considerable heterogeneity in biology and clinical course, many mature B-cell malignancies are highly sensitive to kinase inhibitors that disrupt BCR signaling. Thus, targeted therapy through inhibition of BCR signaling is emerging as a new treatment paradigm for many B-cell malignancies. Here, we review the role of the BCR in the pathogenesis of B-cell malignancies and summarize clinical results of the emerging class of kinase inhibitors that target this pathway.

Automated analysis of immunoglobulin genes from high-throughput sequencing: life without a template

Background

Immunoglobulin (that is, antibody) and T cell receptor genes are created through somatic gene rearrangement from gene segment libraries. Immunoglobulin genes are further diversified by somatic hypermutation and selection during the immune response. Studying the repertoires of these genes yields valuable insights into immune system function in infections, aging, autoimmune diseases and cancers. The introduction of high throughput sequencing has generated unprecedented amounts of repertoire and mutation data from immunoglobulin genes. However, common analysis programs are not appropriate for pre-processing and analyzing these data due to the lack of a template or reference for the whole gene.

Results

We present here the automated analysis pipeline we created for this purpose, which integrates various software packages of our own development and others’, and demonstrate its performance.

Conclusions

Our analysis pipeline presented here is highly modular, and makes it possible to analyze the data resulting from high-throughput sequencing of immunoglobulin genes, in spite of the lack of a template gene. An executable version of the Automation program (and its source code) is freely available for downloading from our website: http://immsilico2.lnx.biu.ac.il/Software.html.

Somatic hypermutation analysis in follicular lymphoma provides evidence suggesting bidirectional cell migration between lymph node and bone marrow during disease progression and relapse

In follicular lymphoma, somatic hypermutation of the immunoglobulin heavy chain genes facilitates the identification of different lymphoma cell clones, and the construction of genealogical trees. To investigate the dissemination of lymphoma cells, and the role of bone marrow in disease progression, we simultaneously analyzed the somatic hypermutation patterns of lymph node and bone marrow specimens taken from three patients at onset and relapse of their disease. Immunoglobulin heavy chain genes were amplified by polymerase chain reaction, cloned and sequenced. Mutational pedigrees were constructed in a hierarchical order. When direct transition of one mutation pattern into that of a successor clones was not feasible, hypothetical predecessor clones were created, and a probability measurement calculation was introduced. Eighty-five sequenced clones were generated. The average mutation rates were 13.45% for the lymph node specimens, and 9.78% for the bone marrow ones. Forty-two hypothetical predecessor clones were introduced into inter-compartment pedigrees. The genealogical trees showed that early lymphoma clones with a low mutational load quickly migrate from lymph nodes into the bone marrow. Bi-directional lymphoma cell migration was detectable between the two compartments. In one case of follicular lymphoma, a clone identical to the initial lymph node clone was detected 2 years later in the bone marrow. The newly introduced algorithm allows the evaluation of both time and direction of follicular lymphoma cell migration. We found evidence that follicular lymphoma originates in the lymph node, and infiltrates the bone marrow early in the course of the disease. Moreover, inter-compartment migration between lymph nodes and bone marrow occurs in both directions.

Targeting pathological B cell receptor signalling in lymphoid malignancies

Signalling through the B cell receptor (BCR) is central to the development and maintenance of B cells. In light of the numerous proliferative and survival pathways activated downstream of the BCR, it comes as no surprise that malignant B cells would co-opt this receptor to promote their own growth and survival. However, direct evidence for BCR signalling in human lymphoma has only come to light recently. Roles for antigen-dependent and antigen-independent, or tonic, BCR signalling have now been described for several different lymphoma subtypes. Furthermore, correlative data implicate antigen-dependent BCR signalling in many other forms of lymphoma. A host of therapeutic agents targeting effectors of the BCR signalling pathway are now in clinical trials and have shown initial success against multiple forms of lymphoma.

Automated cleaning and pre-processing of immunoglobulin gene sequences from high-throughput sequencing

High-throughput sequencing (HTS) yields tens of thousands to millions of sequences that require a large amount of pre-processing work to clean various artifacts. Such cleaning cannot be performed manually. Existing programs are not suitable for immunoglobulin (Ig) genes, which are variable and often highly mutated. This paper describes Ig High-Throughput Sequencing Cleaner (Ig-HTS-Cleaner), a program containing a simple cleaning procedure that successfully deals with pre-processing of Ig sequences derived from HTS, and Ig Insertion—Deletion Identifier (Ig-Indel-Identifier), a program for identifying legitimate and artifact insertions and/or deletions (indels). Our programs were designed for analyzing Ig gene sequences obtained by 454 sequencing, but they are applicable to all types of sequences and sequencing platforms. Ig-HTS-Cleaner and Ig-Indel-Identifier have been implemented in Java and saved as executable JAR files, supported on Linux and MS Windows. No special requirements are needed in order to run the programs, except for correctly constructing the input files as explained in the text. The programs' performance has been tested and validated on real and simulated data sets.

Pathogenesis of follicular lymphoma

The hallmark t(14;18)(q32;q21) in follicular lymphoma (FL) results in constitutive overexpression of the BCL2 protein, allowing B cells to abrogate the default germinal center apoptotic program. Most tumors are characterized by recurrent secondary genetic alterations including genomic gains, losses, and mutations, some providing a growth advantage, including alterations in MLL2, EPHA7, TNFRSF14, and EZH2. The sequence in which these events occur and how they contribute to progression and ultimately to transformation is unclear. Lastly, crosstalk between neoplastic B cells and non-neoplastic immune and stromal cells in the microenvironment plays an important role in sustaining tumor cell growth, cultivating immune privilege, and promoting transformation.

Self-antigen recognition by follicular lymphoma B-cell receptors

Follicular lymphoma is a monoclonal B-cell malignancy with each patient's tumor expressing a unique cell surface immunoglobulin (Ig), or B-cell receptor (BCR), that can potentially recognize antigens and/or transduce signals into the tumor cell. Here we evaluated the reactivity of tumor derived Igs for human tissue antigens. Self-reactivity was observed in 26% of tumor Igs (25 of 98). For one follicular lymphoma patient, the recognized self-antigen was identified as myoferlin. This patient's tumor cells bound recombinant myoferlin in proportion to their level of BCR expression, and the binding to myoferlin was preserved despite ongoing somatic hypermutation of Ig variable regions. Furthermore, BCR-mediated signaling was induced after culture of tumor cells with myoferlin. These results suggest that antigen stimulation may provide survival signals to tumor cells and that there is a selective pressure to preserve antigen recognition as the tumor evolves.

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.

Early steps of follicular lymphoma pathogenesis

Follicular lymphoma (FL) pathogenesis is a complex and fascinating multi-hit process, escalating along successive derailments of the distinctive molecular and cellular mechanisms paving B-cell differentiation and activation. This progressive subversion of B-cell receptor diversification mechanisms and B-cell homeostasis likely occurs during a protracted preclinical phase of asymptomatic growth, in which premalignant clones already disseminate and establish "niches" in secondary lymphoid organs. Following FL diagnosis, a parallel indolent behavior is observed in most patients, slowly progressing over a period of many years, to eventually generate a highly refractory (and in some case transform into an aggressive subtype of) lymphoma. Novel insights in human germinal center B-cell biology recently allowed a more comprehensive understanding of the various illegitimate events sequentially involved in the premalignant progression phases. In this review, we will discuss how these new data have modified our perception of early FL pathogenesis, the new questions and challenges it opened up, and how this knowledge could impact on innovative programs of early detection, follow-up, and patient management.

Pathogenesis of human B cell lymphomas

The mechanisms that drive normal B cell differentiation and activation are frequently subverted by B cell lymphomas for their unlimited growth and survival. B cells are particularly prone to malignant transformation because the machinery used for antibody diversification can cause chromosomal translocations and oncogenic mutations. The advent of functional and structural genomics has greatly accelerated our understanding of oncogenic mechanisms in lymphomagenesis. The signaling pathways that normal B cells utilize to sense antigens are frequently derailed in B cell malignancies, leading to constitutive activation of prosurvival pathways. These malignancies co-opt transcriptional regulatory systems that characterize their normal B cell counterparts and frequently alter epigenetic regulators of chromatin structure and gene expression. These mechanistic insights are ushering in an era of targeted therapies for these cancers based on the principles of pathogenesis.