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

Immunomolecular and reactivity landscapes of gut IgA subclasses in homeostasis and inflammatory bowel disease

The human gut includes plasma cells (PCs) expressing immunoglobulin A1 (IgA1) or IgA2, two structurally distinct IgA subclasses with elusive regulation, function, and reactivity. We show here that intestinal IgA1+ and IgA2+ PCs co-emerged early in life, comparably accumulated somatic mutations, and were enriched within short-lived CD19+ and long-lived CD19- PC subsets, respectively. IgA2+ PCs were extensively clonally related to IgA1+ PCs and a subset of them presumably emerged from IgA1+ precursors. Of note, secretory IgA1 (SIgA1) and SIgA2 dually coated a large fraction of mucus-embedded bacteria, including Akkermansia muciniphila. Disruption of homeostasis by inflammatory bowel disease (IBD) was associated with an increase in actively proliferating IgA1+ plasmablasts, a depletion in long-lived IgA2+ PCs, and increased SIgA1+SIgA2+ gut microbiota. Such increase featured enhanced IgA1 reactivity to pathobionts, including Escherichia coli, combined with depletion of beneficial A. muciniphila. Thus, gut IgA1 and IgA2 emerge from clonally related PCs and show unique changes in both frequency and reactivity in IBD.

B cell phylogenetics in the single cell era

The widespread availability of single-cell RNA sequencing (scRNA-seq) has led to the development of new methods for understanding immune responses. Single-cell transcriptome data can now be paired with B cell receptor (BCR) sequences. However, RNA from BCRs cannot be analyzed like most other genes because BCRs are genetically diverse within individuals. In humans, BCRs are shaped through recombination followed by mutation and selection for antigen binding. As these processes co-occur with cell division, B cells can be studied using phylogenetic trees representing the mutations within a clone. B cell trees can link experimental timepoints, tissues, or cellular subtypes. Here, we review the current state and potential of how B cell phylogenetics can be combined with single-cell data to understand immune responses.

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.