Hanna Edelman

IgTree © : Creating Immunoglobulin variable region gene lineage trees

Lineage trees describe the microevolution of cells within an organism. They have been useful in the study of B cell affinity maturation, which is based on somatic hypermutation of immunoglobulin genes in germinal centers and selection of the resulting mutants. Our aim was to create and implement an algorithm that can generate lineage trees from immunoglobulin variable region gene sequences. The IgTree program implements the algorithm we developed, and generates lineage trees. Original sequences found in experiments are assigned to either leaves or internal nodes of the tree. Each tree node represents a single mutation separating the sequences. The mutations that separate the sequences from each other can be point mutations, deletions or insertions. The program can deal with gaps and find potential reversion mutations. The program also enumerates mutation frequencies and sequence motifs around each mutation, on a per-tree basis. The algorithm has proven useful in several studies of immunoglobulin variable region gene mutations.

Somatic hypermutation and antigen-driven selection of B cells are altered in autoimmune diseases

B cells have been found to play a critical role in the pathogenesis of several autoimmune (AI) diseases. A common feature amongst many AI diseases is the formation of ectopic germinal centers (GC) within the afflicted tissue or organ, in which activated B cells expand and undergo somatic hypermutation (SHM) and antigen-driven selection on their immunoglobulin variable region (IgV) genes. However, it is not yet clear whether these processes occurring in ectopic GCs are identical to those in normal GCs. The analysis of IgV mutations has aided in revealing many aspects concerning B cell expansion, mutation and selection in GC reactions. We have applied several mutation analysis methods, based on lineage tree construction, to a large set of data, containing IgV productive and non-productive heavy and light chain sequences from several different tissues, to examine three of the most profoundly studied AI diseases - Rheumatoid Arthritis (RA), Multiple Sclerosis (MS) and Sjögrens Syndrome (SS). We have found that RA and MS sequences exhibited normal mutation spectra and targeting motifs, but a stricter selection compared to normal controls, which was more apparent in RA. SS sequence analysis results deviated from normal controls in both mutation spectra and indications of selection, also showing differences between light and heavy chain IgV and between different tissues. The differences revealed between AI diseases and normal control mutation patterns may result from the different microenvironmental influences to which ectopic GCs are exposed, relative to those in normal secondary lymphoid tissues.

The dynamics of germinal centre selection as measured by graph-theoretical analysis of mutational lineage trees

We have developed a rigorous graph-theoretical algorithm for quantifying the shape properties of mutational lineage trees. We show that information about the dynamics of hypermutation and antigen-driven clonal selection during the humoral immune response is contained in the shape of mutational lineage trees deduced from the responding clones. Age and tissue related differences in the selection process can be studied using this method. Thus, tree shape analysis can be used as a means of elucidating humoral immune response dynamics in various situations.

Analysis of Mutational Lineage Trees from Sites of Primary and Secondary Ig Gene Diversification in Rabbits and Chickens

Lineage trees of mutated rearranged Ig V region sequences in B lymphocyte clones often serve to qualitatively illustrate claims concerning the dynamics of affinity maturation. In this study, we use a novel method for analyzing lineage tree shapes, using terms from graph theory to quantify the differences between primary and secondary diversification in rabbits and chickens. In these species, Ig gene diversification starts with rearrangement of a single (in chicken) or a few (in rabbit) VH genes. Somatic hypermutation and gene conversion contribute to primary diversification in appendix of young rabbits or in bursa of Fabricius of embryonic and young chickens and to secondary diversification during immune responses in germinal centers (GCs). We find that, at least in rabbits, primary diversification appears to occur at a constant rate in the appendix, and the type of Ag-specific selection seen in splenic GCs is absent. This supports the view that a primary repertoire is being generated within the expanding clonally related B cells in appendix of young rabbits and emphasizes the important role that gut-associated lymphoid tissues may play in early development of mammalian immune repertoires. Additionally, the data indicate a higher rate of hypermutation in rabbit and chicken GCs, such that the balance between hypermutation and selection tends more toward mutation and less toward selection in rabbit and chicken compared with murine GCs.

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.

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

Follicular lymphoma (FL), diffuse large B cell lymphoma (DLBCL) and primary central nervous system lymphoma are B cell malignancies. FL and DLBCL have a germinal center origin. We have applied mutational analyses and a novel algorithm for quantifying shape properties of mutational lineage trees to investigate the nature of the diversification, somatic hypermutation and selection processes that affect B cell clones in these malignancies and reveal whether they differ from normal responses. Lineage tree analysis demonstrated higher diversification and mutations per cell in the lymphoma clones. This was caused solely by the longer diversification times of the malignant clones, as their recent diversification processes were similar to those of normal responses, implying similar mutation frequencies. Since previous analyses of antigen-driven selection were shown to yield false positives, we performed a corrected analysis of replacement and silent mutation patterns, which revealed selection against replacement mutations in the framework regions, responsible for the structural integrity of the B cell receptor, but not for positive selection for replacements in the complementary determining regions. Most replacements, however, were neutral or conservative, suggesting that if at all selection operates in these malignancies it is for structural B cell receptor integrity but not for antigen binding.

Novel Analysis of Clonal Diversification in Blood B Cell and Bone Marrow Plasma Cell Clones in Immunoglobulin Light Chain Amyloidosis

Immunoglobulin light chain amyloidosis (AL) is characterized by a limited clonal expansion of plasma cells and amyloid formation. Here, we report restriction in the diversity of VL gene usage with a dominance of clonally related B cells in the peripheral blood (PB) isotype-specific repertoire of AL patients. A rigorous quantification of lineage trees reveals presence of intraclonal variations in the PB clones compared to the bone marrow (BM) clones, which suggests a common precursor that is still subject to somatic mutation. When compared to normal BM and PB B cells, AL clones showed significant but incomplete impairment of antigenic selection, which could not be detected by conventional R and S mutation analysis. Therefore, graphical analysis of B cell lineage trees and mathematical quantification of tree properties provide novel insights into the process of B cell clonal evolution in AL.