Immunoglobulin somatic hypermutation has clinical impact in DLBCL and potential implications for immune checkpoint blockade and neoantigen-based immunotherapies

The clinical applications of immunosequencing

Technological advances in high-throughput sequencing have opened the door for the interrogation of adaptive immune responses at unprecedented scale. It is now possible to determine the sequences of antibodies or T-cell receptors produced by individual B and T cells in a sample. This capability, termed immunosequencing, has transformed the study of both infectious and non-infectious diseases by allowing the tracking of dynamic changes in B and T cell clonal populations over time. This has improved our understanding of the pathology of cancers, autoimmune diseases, and infectious diseases. However, to date there has been only limited clinical adoption of the technology. Advances over the last decade and on the horizon that reduce costs and improve interpretability could enable widespread clinical use. Many clinical applications have been proposed and, while most are still undergoing research and development, some methods relying on immunosequencing data have been implemented, the most widespread of which is the detection of measurable residual disease. Here, we review the diagnostic, prognostic, and therapeutic applications of immunosequencing for both infectious and non-infectious diseases.

Complete detection of FR1 to FR3 primer-based PCR patterns of immunoglobulin heavy chain rearrangement in the BIOMED-2 protocol is associated with poor prognosis in patients with diffuse large B-cell lymphoma

Somatic hypermutations (SHMs) in the variable region (VH) of the immunoglobulin heavy chain (IgH) gene are common in diffuse large B-cell lymphoma (DLBCL). Recently, IgH VH SHMs have become known as immunogenic neoantigens, but few studies have evaluated the prognostic impact of the frequency of VH SHMs in DLBCL. The BIOMED-2 protocol is the gold standard polymerase chain reaction (PCR) for clonality analysis in lymphoid malignancies, but can produce false negatives due to the presence of IgH VH SHMs. To overcome this problem, three primer sets were designed for the three framework regions (FR1, FR2, and FR3). We evaluated the predictive value of this PCR pattern in patients with DLBCL. To evaluate the prognostic impact of complete detection of the clonal amplifications (VHFR1-JH, VHFR2-JH, and VHFR3-JH) in the BIOMED-2 protocol, we retrospectively analyzed 301 DLBCL patients who were initially treated with anthracycline-based immunochemotherapy. Complete detection of the FR1 to FR3 primer-based IgH VH PCR patterns in the BIOMED-2 protocol was associated with low frequency of VH SHMs (p < 0.001). Patients who were positive for all these three PCRs (n = 79) were significantly associated with shorter 5-year overall survival (OS; 54.2% vs. 73.2%; p = 0.002) and progression-free survival (PFS; 34.3% vs. 59.3%; p < 0.001) compared to patients with other PCR patterns (n = 202). Specifically, the successful FR3-JH detection was associated with significantly worse OS (p < 0.001) and PFS (p < 0.001). PCR patterns of complete IgH rearrangement using the BIOMED-2 protocol are clinically meaningful indicators for prognostic stratification of DLBCL patients.

Tumor-Infiltrating Normal B Cells Revealed by Immunoglobulin Repertoire Clonotype Analysis Are Highly Prognostic and Crucial for Antitumor Immune Responses in DLBCL

PURPOSE

Tumor-infiltrating B lymphocytes (TIL-B) have demonstrated prognostic and predictive significance in solid cancers. In this study, we aimed to distinguish TIL-Bs from malignant B-cells in diffuse large B-cell lymphoma (DLBCL) and determine the clinical and biological significance.

EXPERIMENTAL DESIGN

A total of 269 patients with de novo DLBCL from the International DLBCL R-CHOP Consortium Program were studied. Ultra-deep sequencing of the immunoglobulin genes was performed to determine B-cell clonotypes. The frequencies and numbers of TIL-B clonotypes in individual repertoires were correlated with patient survival, gene expression profiling (GEP) data, and frequencies of DLBCL-infiltrating immune cells quantified by fluorescent multiplex IHC at single-cell resolution.

RESULTS

TIL-B abundance, evaluated by frequencies of normal B-cell clonotypes in the immunoglobulin repertoires, remarkably showed positive associations with significantly better survival of patients in our sequenced cohorts. DLBCLs with high versus low TIL-B abundance displayed distinct GEP signatures, increased pre-memory B-cell state and naïve CD4 T-cell state fractions, and higher CD4+ T-cell infiltration. TIL-B frequency, as a new biomarker in DLBCL, outperformed the germinal center (GC) B-cell-like/activated B-cell-like classification and TIL-T frequency. The identified TIL-B-high GEP signature, including genes upregulated during T-dependent B-cell activation and those highly expressed in normal GC B cells and T cells, showed significant favorable prognostic effects in several external validation cohorts.

CONCLUSIONS

TIL-B frequency is a significant prognostic factor in DLBCL and plays a crucial role in antitumor immune responses. This study provides novel insights into the prognostic determinants in DLBCL and TIL-B functions with important therapeutic implications.

Tutorial: integrative computational analysis of bulk RNA-sequencing data to characterize tumor immunity using RIMA

RNA-sequencing (RNA-seq) has become an increasingly cost-effective technique for molecular profiling and immune characterization of tumors. In the past decade, many computational tools have been developed to characterize tumor immunity from gene expression data. However, the analysis of large-scale RNA-seq data requires bioinformatics proficiency, large computational resources and cancer genomics and immunology knowledge. In this tutorial, we provide an overview of computational analysis of bulk RNA-seq data for immune characterization of tumors and introduce commonly used computational tools with relevance to cancer immunology and immunotherapy. These tools have diverse functions such as evaluation of expression signatures, estimation of immune infiltration, inference of the immune repertoire, prediction of immunotherapy response, neoantigen detection and microbiome quantification. We describe the RNA-seq IMmune Analysis (RIMA) pipeline integrating many of these tools to streamline RNA-seq analysis. We also developed a comprehensive and user-friendly guide in the form of a GitBook with text and video demos to assist users in analyzing bulk RNA-seq data for immune characterization at both individual sample and cohort levels by using RIMA.

The Immunology of DLBCL

Diffuse large B-cell lymphoma (DLBCL) is an aggressive malignancy and is the most common type of malignant lymphoid neoplasm. While some DLBCLs exhibit strong cell-autonomous survival and proliferation activity, others depend on interactions with non-malignant cells for their survival and proliferation. Recent next-generation sequencing studies have linked these interactions with the molecular classification of DLBCL. For example, germinal center B-cell-like DLBCL tends to show strong associations with follicular T cells and epigenetic regulation of immune recognition molecules, whereas activated B-cell-like DLBCL shows frequent genetic aberrations affecting the class I major histocompatibility complex. Single-cell technologies have also provided detailed information about cell-cell interactions and the cell composition of the microenvironment of DLBCL. Aging-related immunological deterioration, i.e., immunosenescence, also plays an important role in DLBCL pathogenesis, especially in Epstein-Barr virus-positive DLBCL. Moreover, DLBCL in "immune-privileged sites"-where multiple immune-modulating mechanisms exist-shows unique biological features, including frequent down-regulation of immune recognition molecules and an immune-tolerogenic tumor microenvironment. These advances in understanding the immunology of DLBCL may contribute to the development of novel therapies targeting immune systems.

Anti-PD-L1/PD-L2 therapeutic vaccination in untreated chronic lymphocytic leukemia patients with unmutated IgHV

Chronic lymphocytic leukemia (CLL) patients with unmutated immunoglobulin heavy chain (IgHV) are at risk of early disease progression compared to patients with mutated IgHV. As a preventive strategy, we treated 19 previously untreated CLL patients with unmutated IgHV in a phase 1/2 trial (clinicaltrials.gov, NCT03939234) exploring the efficacy and toxicity of a therapeutic cancer vaccine containing peptides derived from programmed death ligand 1 (PD-L1) and ligand 2 (PD-L2), hoping to restore immunological control of the disease. According to the International Workshop on Chronic lymphocytic Leukemia (iwCLL) response criteria, no patients obtained a response; however, during follow-up, one patient had complete normalization of the peripheral lymphocyte count and remained in biochemical remission after a follow-up time of 15 months. At the end of treatment, one patient had progressed, and 17 patients had stable disease. During follow-up with a median time of 23.5 months since inclusion, seven patients had progressed, and eight patients had stable disease. The median time to first treatment (TTFT) from diagnosis was 90.3 months with a median follow-up time of 50.1 months. This apparent favorable outcome in TTFT needs to be investigated in a randomized setting, as our population may have been biased. More than 80% of patients obtained vaccine-specific immune responses, confirming the immunogenicity of the vaccine. The vaccine was generally well tolerated with only grade I-II adverse events. Although there were some signs of clinical effects, the vaccine seems to be insufficient as monotherapy in CLL, possibly due to a high tumor burden. The efficacy of the vaccine should preferably be tested in combination with novel targeted therapies or as a consolidating treatment.

Extracellular Vesicles in Diffuse Large B Cell Lymphoma: Characterization and Diagnostic Potential

Diffuse large B cell lymphoma (DLBCL) is an aggressive B cell lymphoma characterized by a heterogeneous behavior and in need of more accurate biological characterization monitoring and prognostic tools. Extracellular vesicles are secreted by all cell types and are currently established to some extent as representatives of the cell of origin. The present study characterized and evaluated the diagnostic and prognostic potential of plasma extracellular vesicles (EVs) proteome in DLBCL by using state-of-the-art mass spectrometry. The EV proteome is strongly affected by DLBCL status, with multiple proteins uniquely identified in the plasma of DLBCL. A proof-of-concept classifier resulted in highly accurate classification with a sensitivity and specificity of 1 when tested on the holdout test data set. On the other hand, no proteins were identified to correlate with non-germinal center B-cell like (non-GCB) or GCB subtypes to a significant degree after correction for multiple testing. However, functional analysis suggested that antigen binding is regulated when comparing non-GCB and GCB. Survival analysis based on protein quantitative values and clinical parameters identified multiple EV proteins as significantly correlated to survival. In conclusion, the plasma extracellular vesicle proteome identifies DLBCL cancer patients from healthy donors and contains potential EV protein markers for prediction of survival.

Next-generation sequencing for MRD monitoring in B-lineage malignancies: from bench to bedside

Minimal residual disease (MRD) is considered the strongest relevant predictor of prognosis and an effective decision-making factor during the treatment of hematological malignancies. Remarkable breakthroughs brought about by new strategies, such as epigenetic therapy and chimeric antigen receptor-T (CAR-T) therapy, have led to considerably deeper responses in patients than ever, which presents difficulties with the widely applied gold-standard techniques of MRD monitoring. Urgent demands for novel approaches that are ultrasensitive and provide sufficient information have put a spotlight on high-throughput technologies. Recently, advances in methodology, represented by next-generation sequencing (NGS)-based clonality assays, have proven robust and suggestive in numerous high-quality studies and have been recommended by some international expert groups as disease-monitoring modalities. This review demonstrates the applicability of NGS-based clonality assessment for MRD monitoring of B-cell malignancies by summarizing the oncogenesis of neoplasms and the corresponding status of immunoglobulin (IG) rearrangements. Furthermore, we focused on the performance of NGS-based assays compared with conventional approaches and the interpretation of results, revealing directions for improvement and prospects in clinical practice.

Benchmarking computational methods for B-cell receptor reconstruction from single-cell RNA-seq data

Multiple methods have recently been developed to reconstruct full-length B-cell receptors (BCRs) from single-cell RNA sequencing (scRNA-seq) data. This need emerged from the expansion of scRNA-seq techniques, the increasing interest in antibody-based drug development and the importance of BCR repertoire changes in cancer and autoimmune disease progression. However, a comprehensive assessment of performance-influencing factors such as the sequencing depth, read length or number of somatic hypermutations (SHMs) as well as guidance regarding the choice of methodology is still lacking. In this work, we evaluated the ability of six available methods to reconstruct full-length BCRs using one simulated and three experimental SMART-seq datasets. In addition, we validated that the BCRs assembled in silico recognize their intended targets when expressed as monoclonal antibodies. We observed that methods such as BALDR, BASIC and BRACER showed the best overall performance across the tested datasets and conditions, whereas only BASIC demonstrated acceptable results on very short read libraries. Furthermore, the de novo assembly-based methods BRACER and BALDR were the most accurate in reconstructing BCRs harboring different degrees of SHMs in the variable domain, while TRUST4, MiXCR and BASIC were the fastest. Finally, we propose guidelines to select the best method based on the given data characteristics.

Recurrence Risk of Liver Cancer Post-hepatectomy Using Machine Learning and Study of Correlation With Immune Infiltration

Postoperative recurrence of liver cancer is the main obstacle to improving the survival rate of patients with liver cancer. We established an mRNA-based model to predict the risk of recurrence after hepatectomy for liver cancer and explored the relationship between immune infiltration and the risk of recurrence after hepatectomy for liver cancer. We performed a series of bioinformatics analyses on the gene expression profiles of patients with liver cancer, and selected 18 mRNAs as biomarkers for predicting the risk of recurrence of liver cancer using a machine learning method. At the same time, we evaluated the immune infiltration of the samples and conducted a joint analysis of the recurrence risk of liver cancer and found that B cell, B cell naive, T cell CD4⁺ memory resting, and T cell CD4⁺ were significantly correlated with the risk of postoperative recurrence of liver cancer. These results are helpful for early detection, intervention, and the individualized treatment of patients with liver cancer after surgical resection, and help to reveal the potential mechanism of liver cancer recurrence.

Insertion of atypical glycans into the tumor antigen-binding site identifies DLBCLs with distinct origin and behavior

Glycosylation of the surface immunoglobulin (Ig) variable region is a remarkable follicular lymphoma-associated feature rarely seen in normal B cells. Here, we define a subset of diffuse large B-cell lymphomas (DLBCLs) that acquire N-glycosylation sites selectively in the Ig complementarity-determining regions (CDRs) of the antigen-binding sites. Mass spectrometry and X-ray crystallography demonstrate how the inserted glycans are stalled at oligomannose-type structures because they are buried in the CDR loops. Acquisition of sites occurs in ∼50% of germinal-center B-cell-like DLBCL (GCB-DLBCL), mainly of the genetic EZB subtype, irrespective of IGHV-D-J use. This markedly contrasts with the activated B-cell-like DLBCL Ig, which rarely has sites in the CDR and does not seem to acquire oligomannose-type structures. Acquisition of CDR-located acceptor sites associates with mutations of epigenetic regulators and BCL2 translocations, indicating an origin shared with follicular lymphoma. Within the EZB subtype, these sites are associated with more rapid disease progression and with significant gene set enrichment of the B-cell receptor, PI3K/AKT/MTORC1 pathway, glucose metabolism, and MYC signaling pathways, particularly in the fraction devoid of MYC translocations. The oligomannose-type glycans on the lymphoma cells interact with the candidate lectin dendritic cell-specific intercellular adhesion molecule 3 grabbing non-integrin (DC-SIGN), mediating low-level signals, and lectin-expressing cells form clusters with lymphoma cells. Both clustering and signaling are inhibited by antibodies specifically targeting the DC-SIGN carbohydrate recognition domain. Oligomannosylation of the tumor Ig is a posttranslational modification that readily identifies a distinct GCB-DLBCL category with more aggressive clinical behavior, and it could be a potential precise therapeutic target via antibody-mediated inhibition of the tumor Ig interaction with DC-SIGN-expressing M2-polarized macrophages.

Genomic complexity is associated with epigenetic regulator mutations and poor prognosis in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common type of lymphoma with high mutation burdens but a low response rate to immune checkpoint inhibitors. In this study, we performed targeted next-generation sequencing and fluorescent multiplex immunohistochemistry, and investigated the clinical significance and immunological effect of mutation numbers in 424 DLBCL patients treated with standard immunochemotherapy. We found that KMT2D and TP53 nonsynonymous mutations (MUT) were significantly associated with increased nonsynonymous mutation numbers, and that high mutation numbers (MUT^high) were associated with significantly poorer clinical outcome in germinal center B-cell-like DLBCL with wild-type TP53. To understand the underlying mechanisms, we identified a gene-expression profiling signature and the association of MUT^high with decreased T cells in DLBCL patients with wild-type TP53. On the other hand, in overall cohort, MUT^high was associated with lower PD-1 expression in T cells and PD-L1 expression in macrophages, suggesting a positive role of MUT^high in immune responses. Analysis in a whole-exome sequencing dataset of 304 patients deposited by Chapuy et al. validated the correlation of MUT-KMT2D with genomic complexity and the significantly poorer survival associated with higher numbers of genomic single nucleotide variants in activated B-cell-like DLBCL with wild-type TP53. Together, these results suggest that KMT2D inactivation or epigenetic dysregulation has a role in driving DLBCL genomic instability, and that genomic complexity has adverse impact on clinical outcome in DLBCL patients with wild-type TP53 treated with standard immunochemotherapy. The oncoimmune data in this study have important implications for biomarker and therapeutic studies in DLBCL.

The tumor microenvironment of lymphomas: Insights into the potential role and modes of actions of checkpoint inhibitors

The tumor microenvironment (TME) - a term comprising non-neoplastic cells and extracellular matrix as well as various cytokines, chemokines, growth factors, and other substances in the vicinity of tumor cells - is an integrative part of most tumors including lymphomas. Interactions between lymphoma cells and the TME are vital for survival and proliferation of the former. In addition, lymphoma cells often reprogram the TME to protect them from defense mechanisms of the host's immune system. In this review, we will introduce the role of the tumor microenvironment (TME) for lymphoma cells looking at direct cell-cell interactions as well as cytokine-related communications. The immunomodulative/immunosuppressive role of the TME is more and more coming into the focus of potential new targeted therapies, and thus a special attention will be given to the interactions of immune checkpoints such as programed cell death protein 1 and L1 (PD-1/PD-L1), T-cell immunoglobulin and mucin-domain containing protein-3 (TIM-3), lymphocyte-activation gene 3 (LAG-3), and cytotoxic T-lymphocyte-associated protein-4 (CTLA4) with the TME, as well as their expression by both lymphoma cells and cells of the TME. Aspects of the TME will be discussed for indolent and aggressive B-cell lymphomas, Hodgkin lymphomas, and T-cell lymphomas. In addition, the potential influence of other immunomodulators such as lenalidomide will be briefly touched. The complex role of the TME is in the focus of new therapeutic options. In order to exploit its full therapeutic potential, however, a thorough understanding of TME biology and interaction between lymphoma cells and the TME, as well as the host's immune system and the TME is necessary.

Therapeutic vaccines for aggressive B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive B-cell lymphoma and highly heterogeneous disease. With the standard immunochemotherapy, anti-CD20 antibody rituximab (R-) plus CHOP (cyclophosphamide, doxorubicin, vincristine, and prednisone) chemotherapy, 30-40% of DLBCLs are refractory to initial immunochemotherapy or experience relapse post-therapy with poor clinical outcomes despite salvage therapies. Mechanisms underlying chemoresistance and relapse are heterogeneous across DLBCL and within individual patients, representing hurdles for targeted therapies targeting a specific oncogenic signaling pathway. In recent years, paradigm-shifting immunotherapies have shown impressive efficacy in various cancer types regardless of underlying oncogenic mechanisms. Vaccines are being developed for DLBCL to build protective immunity against relapse after first complete remission and to promote antitumor immune responses synergizing with immune checkpoint inhibitors to treat refractory/relapsed patients. This article provides a brief review of current progress in vaccine development in DLBCL and discussion on immunologic mechanisms underlying the therapeutic effectiveness and resistance.

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.

CD4+ T Cells in the Blood of MS Patients Respond to Predicted Epitopes From B cell Receptors Found in Spinal Fluid

B cells are important pathogenic players in multiple sclerosis (MS), but their exact role is not known. We have previously demonstrated that B cells from cerebrospinal fluid (CSF) of MS patients can activate T cells that specifically recognize antigenic determinants (idiotopes) from their B cell receptors (BCRs). The aim of this study was to evaluate whether in silico prediction models could identify antigenic idiotopes of immunoglobulin heavy-chain variable (IGHV) transcriptomes in MS patients. We utilized a previously assembled dataset of CSF IGHV repertoires from MS patients. To guide selection of potential antigenic idiotopes, we used in silico predicted HLA-DR affinity, endosomal processing, as well as transcript frequency from nine MS patients. Idiotopes with predicted low affinity and low likelihood of cathepsins cleavage were inert controls. Peripheral blood mononuclear cells from these patients were stimulated with the selected idiotope peptides in presence of anti-CD40 for 12 h. T cells were then labeled for activation status with anti-CD154 antibodies and CD3+CD4+ T cells phenotyped as memory (CD45RO+) or naïve (CD45RO-), with potential for brain migration (CXCR3 and/or CCR6 expression). Anti-CD14 and -CD8 were utilized to exclude monocytes and CD8+ T cells. Unstimulated cells or insulin peptides were negative controls, and EBNA-1 peptides or CD3/CD28 beads were positive controls. The mean proportion of responding memory CD4+ T cells from all nine MS patients was significantly higher for idiotope peptides with predicted high HLA-DR affinity and high likelihood of cathepsin cleavage, than toward predicted inert peptides. Responses were mainly observed toward peptides affiliated with the CDR3 region. Activated memory CD4+ T cells expressed the chemokine receptor CCR6, affiliated with a Th17 phenotype and allowing passage into the central nervous system (CNS). This in vitro study suggests that that antigenic properties of BCR idiotopes can be identified in silico using HLA affinity and endosomal processing predictions. It further indicates that MS patients have a memory T cell repertoire capable of recognizing frequent BCR idiotopes found in endogenous CSF, and that these T cells express chemokine receptors allowing them to reach the CSF B cells expressing these idiotopes.