Genome-wide mutational signatures revealed distinct developmental paths for human B cell lymphomas

Identification of ferroptosis-related genes associated with diffuse large B-cell lymphoma via bioinformatics and machine learning approaches

Ferroptosis has emerged as a critical mechanism in the development and progression of various tumors, particularly diffuse large B-cell lymphoma (DLBCL). However, the thorough characterization of ferroptosis-related genes in DLBCL remains inadequately explored. We retrieved datasets associated with DLBCL and ferroptosis gene sets from the Gene Expression Omnibus (GEO) database and the Ferroptosis Database (FerrDb), resulting in the identification of 27 differentially expressed ferroptosis-related genes (DE-FRGs) linked to DLBCL. Utilizing the LASSO and Support Vector Machine Recursive Feature Elimination (SVW-RFE) algorithms, we identified 10 genes-MT1G, MTOR, BRD4, ACO1, SAT1, PEBP1, LPIN1, ATM, SRXN1, and PRDX1-as key biomarker candidates with significant diagnostic potential. Functional enrichment analyses revealed that these biomarker genes are likely involved in regulating several critical biological pathways implicated in DLBCL pathogenesis, including immune response, oxidative phosphorylation, and cell cycle regulation. Moreover, we identified 246 potential therapeutic agents targeting these 10 biomarker genes. Concurrently, competitive endogenous RNA (ceRNA) network analysis uncovered a complex regulatory network centered on the identified biomarker genes. Additionally, CIBERSORT analysis highlighted notable alterations in the immune microenvironment of DLBCL patients. We propose a diagnostic strategy that provides novel insights into the molecular mechanisms underlying DLBCL. Nevertheless, further validation of the practical value of this strategy for DLBCL diagnosis is necessary before its clinical application.

Orientation Regulation of Class-switch Recombination in Human B Cells

We developed a linear amplification-mediated high-throughput genome-wide translocation sequencing method to profile Ig class-switch recombination (CSR) in human B cells in an unbiased and quantitative manner. This enables us to characterize CSR junctions resulting from either deletional recombination or inversion for each Ig class/subclass. Our data showed that more than 90% of CSR junctions detected in peripheral blood in healthy control subjects were due to deletional recombination. We further identified two major CSR junction signatures/patterns in human B cells. Signature 1 consists of recombination junctions resulting from both IgG and IgA switching, with a dominance of Sµ-Sγ junctions (72%) and deletional recombination (87%). Signature 2 is contributed mainly by Sµ-Sα junctions (96%), and these junctions were almost all due to deletional recombination (99%) and were characterized by longer microhomologies. CSR junctions identified in healthy individuals can be assigned to both signatures but with a dominance of signature 1, whereas almost all CSR junctions found in patients with defects in DNA-PKcs or Artemis, two classical nonhomologous end joining (c-NHEJ) factors, align with signature 2. Thus, signature 1 may represent c-NHEJ activity during CSR, whereas signature 2 is associated with microhomology-mediated alternative end joining in the absence of the studied c-NHEJ factors. Our findings suggest that in human B cells, the efficiency of the c-NHEJ machinery and the features of switch regions are crucial for the regulation of CSR orientation. Finally, our high-throughput method can also be applied to study the mechanism of rare types of recombination, such as switching to IgD and locus suicide switching.

Molecular Determinants of Sensitivity to Polatuzumab Vedotin in Diffuse Large B-Cell Lymphoma

Polatuzumab vedotin (Pola-V) is an antibody-drug conjugate directed to the CD79B subunit of the B-cell receptor (BCR). When combined with conventional immunochemotherapy, Pola-V improves outcomes in diffuse large B-cell lymphoma (DLBCL). To identify determinants of Pola-V sensitivity, we used CRISPR-Cas9 screening for genes that modulated Pola-V toxicity for lymphomas or the surface expression of its target, CD79B. Our results reveal the striking impact of CD79B glycosylation on Pola-V epitope availability on the lymphoma cell surface and on Pola-V toxicity. Genetic, pharmacological, and enzymatic approaches that remove sialic acid from N-linked glycans enhanced lymphoma killing by Pola-V. Pola-V toxicity was also modulated by KLHL6, an E3 ubiquitin ligase that is recurrently inactivated in germinal center derived lymphomas. We reveal how KLHL6 targets CD79B for degradation in normal and malignant germinal center B cells, thereby determining expression of the surface BCR complex. Our findings suggest precision medicine strategies to optimize Pola-V as a lymphoma therapeutic. Significance: These findings unravel the molecular basis of response heterogeneity to Pola-V and identify approaches that might be deployed therapeutically to enhance the efficacy of CD79B-specific tumor killing. In addition, they reveal a novel post-translational mechanism used by normal and malignant germinal center B cells to regulate expression of the BCR. See related commentary by Leveille, p. 1577 See related article by Meriranta et al.

Motive and opportunity: MYC rearrangements in high-grade B-cell lymphoma with MYC and BCL2 rearrangements (an LLMPP study)

ABSTRACT

Rearrangements that place the oncogenes MYC, BCL2, or BCL6 adjacent to superenhancers are common in mature B-cell lymphomas. Lymphomas with diffuse large B-cell lymphoma (DLBCL) or high-grade morphology with both MYC and BCL2 rearrangements are classified as high-grade B-cell lymphoma with MYC and BCL2 rearrangements ("double hit"; HGBCL-DH-BCL2) and are associated with aggressive disease and poor outcomes. Although it is established that MYC rearrangements involving immunoglobulin (IG) loci are associated with inferior outcomes relative to those involving other non-IG superenhancers, the frequency of and mechanisms driving IG vs non-IG MYC rearrangements have not been elucidated. Here, we used custom targeted capture and/or whole-genome sequencing to characterize oncogene rearrangements across 883 mature B-cell lymphomas including Burkitt lymphoma, follicular lymphoma, DLBCL, and HGBCL-DH-BCL2 tumors. We demonstrate that, although BCL2 rearrangement topology is consistent across entities, HGBCL-DH-BCL2 have distinct MYC rearrangement architecture relative to tumors with single MYC rearrangements or with both MYC and BCL6 rearrangements (HGBCL-DH-BCL6), including both a higher frequency of non-IG rearrangements and different architecture of MYC::IGH rearrangements. The distinct MYC rearrangement patterns in HGBCL-DH-BCL2 occur on the background of high levels of somatic hypermutation across MYC partner loci in HGBCL-DH-BCL2, creating more opportunity to form these rearrangements. Furthermore, because 1 IGH allele is already disrupted by the existing BCL2 rearrangement, the MYC rearrangement architecture in HGBCL-DH-BCL2 likely reflects selective pressure to preserve both BCL2 and B-cell receptor expression. These data provide new mechanistic explanations for the distinct patterns of MYC rearrangements observed across different lymphoma entities.

Taming AID mutator activity in somatic hypermutation

Activation-induced cytidine deaminase (AID) initiates somatic hypermutation (SHM) by introducing base substitutions into antibody genes, a process enabling antibody affinity maturation in immune response. How a mutator is tamed to precisely and safely generate programmed DNA lesions in a physiological process remains unsettled, as its dysregulation drives lymphomagenesis. Recent research has revealed several hidden features of AID-initiated mutagenesis: preferential activity on flexible DNA substrates, restrained activity within chromatin loop domains, unique DNA repair factors to differentially decode AID-caused lesions, and diverse consequences of aberrant deamination. Here, we depict the multifaceted regulation of AID activity with a focus on emerging concepts/factors and discuss their implications for the design of base editors (BEs) that install somatic mutations to correct deleterious genomic variants.

The RNA tether model for human chromosomal translocation fragile zones

One of the two chromosomal breakage events in recurring translocations in B cell neoplasms is often due to the recombination-activating gene complex (RAG complex) releasing DNA ends before end joining. The other break occurs in a fragile zone of 20-600 bp in a non-antigen receptor gene locus, with a more complex and intriguing set of mechanistic factors underlying such narrow fragile zones. These factors include activation-induced deaminase (AID), which acts only at regions of single-stranded DNA (ssDNA). Recent work leads to a model involving the tethering of AID to the nascent RNA as it emerges from the RNA polymerase. This mechanism may have relevance in class switch recombination (CSR) and somatic hypermutation (SHM), as well as broader relevance for other DNA enzymes.

AID in non-Hodgkin B-cell lymphomas: The consequences of on- and off-target activity

Activation induced cytidine deaminase (AID) is a key element of the adaptive immune system, required for immunoglobulin isotype switching and affinity maturation of B-cells as they undergo the germinal center (GC) reaction in peripheral lymphoid tissue. The inherent DNA damaging activity of this enzyme can also have off-target effects in B-cells, producing lymphomagenic chromosomal translocations that are characteristic features of various classes of non-Hodgkin B-cell lymphoma (B-NHL), and generating oncogenic mutations, so-called aberrant somatic hypermutation (aSHM). Additionally, AID has been found to affect gene expression through demethylation as well as altered interactions between gene regulatory elements. These changes have been most thoroughly studied in B-NHL arising from GC B-cells. Here, we describe the most common classes of GC-derived B-NHL and explore the consequences of on- and off-target AID activity in B and plasma cell neoplasms. The relationships between AID expression, including effects of infection and other exposures/agents, mutagenic activity and lymphoma biology are also discussed.

Follicular lymphoma B cells exhibit heterogeneous transcriptional states with associated somatic alterations and tumor microenvironments

Follicular lymphoma (FL) is an indolent non-Hodgkin lymphoma of germinal center origin, which presents with significant biologic and clinical heterogeneity. Using RNA-seq on B cells sorted from 87 FL biopsies, combined with machine-learning approaches, we identify 3 transcriptional states that divide the biological ontology of FL B cells into inflamed, proliferative, and chromatin-modifying states, with relationship to prior GC B cell phenotypes. When integrated with whole-exome sequencing and immune profiling, we find that each state was associated with a combination of mutations in chromatin modifiers, copy-number alterations to TNFAIP3, and T follicular helper cells (Tfh) cell interactions, or primarily by a microenvironment rich in activated T cells. Altogether, these data define FL B cell transcriptional states across a large cohort of patients, contribute to our understanding of FL heterogeneity at the tumor cell level, and provide a foundation for guiding therapeutic intervention.

5-Chloro-2'-deoxycytidine Induces a Distinctive High-Resolution Mutational Spectrum of Transition Mutations In Vivo

The biomarker 5-chlorocytosine (5ClC) appears in the DNA of inflamed tissues. Replication of a site-specific 5ClC in a viral DNA genome results in C → T mutations, which is consistent with 5ClC acting as a thymine mimic in vivo. Direct damage of nucleic acids by immune-cell-derived hypochlorous acid is one mechanism by which 5ClC could appear in the genome. A second, nonmutually exclusive mechanism involves damage of cytosine nucleosides or nucleotides in the DNA precursor pool, with subsequent utilization of the 5ClC deoxynucleotide triphosphate as a precursor for DNA synthesis. The present work characterized the mutagenic properties of 5ClC in the nucleotide pool by exposing cells to the nucleoside 5-chloro-2'-deoxycytidine (5CldC). In both Escherichia coli and mouse embryonic fibroblasts (MEFs), 5CldC in the growth media was potently mutagenic, indicating that 5CldC enters cells and likely is erroneously incorporated into the genome from the nucleotide pool. High-resolution sequencing of DNA from MEFs derived from the gptΔ C57BL/6J mouse allowed qualitative and quantitative characterization of 5CldC-induced mutations; CG → TA transitions in 5'-GC(Y)-3' contexts (Y = a pyrimidine) were dominant, while TA → CG transitions appeared at a much lower frequency. The high-resolution mutational spectrum of 5CldC revealed a notable similarity to the Catalogue of Somatic Mutations in Cancer mutational signatures SBS84 and SBS42, which appear in human lymphoid tumors and in occupationally induced cholangiocarcinomas, respectively. SBS84 is associated with the expression of activation-induced cytidine deaminase (AID), a cytosine deaminase associated with inflammation, as well as immunoglobulin gene diversification during antibody maturation. The similarity between the spectra of AID activation and 5CldC could be coincidental; however, the administration of 5CldC did induce some AID expression in MEFs, which have no inherent expression of its gene. In summary, this work shows that 5CldC induces a distinct pattern of mutations in cells. Moreover, that pattern resembles human mutational signatures induced by inflammatory processes, such as those triggered in certain malignancies.

Relapses in early-stage follicular lymphoma frequently develop via a divergent evolution from their clonally related precursor cells

Follicular lymphoma (FL) develops through a stepwise acquisition of cooperative genetic changes with t(14;18)(q32;q21)/IGH::BCL2 occurring early at the pre-B stage of B-cell development. Patients with FL typically show an indolent clinical course, remitting and relapsing with the eventual development of resistance to treatments. Interestingly, the majority of transformed FL do not progress directly from FL but originate from their clonally related lymphoma precursor (CLP) cells. To examine whether such divergent tumour evolution also underpins the relapses in patients with early-stage FL, we investigated by targeted next-generation sequencing 13 cases (stage I = 9, stage II = 4), who showed complete remission (mean: 5 years; range: 1-11.5 years) following local radiotherapy but subsequently relapsed (≥2 in 5). A clonal relationship between the diagnostic FL and relapses was confirmed in 11 cases. In six cases, common and distinct variants were seen between the paired diagnostic and relapsed lymphomas, indicating their divergent evolution from a CLP. In two cases, different B-cell clones were involved in the diagnostic and relapsed lymphomas, including one case involving two different BCL2 translocations. In the remaining five cases, the relapsed lymphoma developed via a linear progression (n = 4) or a mixed evolutionary path (n = 1). These findings may bear important implications in the routine diagnosis and management of relapsed FL. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Genetic and transcriptomic analyses of diffuse large B-cell lymphoma patients with poor outcomes within two years of diagnosis

Despite the improvements in clinical outcomes for DLBCL, a significant proportion of patients still face challenges with refractory/relapsed (R/R) disease after receiving first-line R-CHOP treatment. To further elucidate the underlying mechanism of R/R disease and to develop methods for identifying patients at risk of early disease progression, we integrated clinical, genetic and transcriptomic data derived from 2805 R-CHOP-treated patients from seven independent cohorts. Among these, 887 patients exhibited R/R disease within two years (poor outcome), and 1918 patients remained in remission at two years (good outcome). Our analysis identified four preferentially mutated genes (TP53, MYD88, SPEN, MYC) in the untreated (diagnostic) tumor samples from patients with poor outcomes. Furthermore, transcriptomic analysis revealed a distinct gene expression pattern linked to poor outcomes, affecting pathways involved in cell adhesion/migration, T-cell activation/regulation, PI3K, and NF-κB signaling. Moreover, we developed and validated a 24-gene expression score as an independent prognostic predictor for treatment outcomes. This score also demonstrated efficacy in further stratifying high-risk patients when integrated with existing genetic or cell-of-origin subtypes, including the unclassified cases in these models. Finally, based on these findings, we developed an online analysis tool ( https://lymphprog.serve.scilifelab.se/app/lymphprog ) that can be used for prognostic prediction for DLBCL patients.

Activation-induced cytidine deaminase causes recurrent splicing mutations in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma. A major mutagenic process in DLBCL is aberrant somatic hypermutation (aSHM) by activation-induced cytidine deaminase (AID), which occurs preferentially at RCH/TW sequence motifs proximal to transcription start sites. Splice sequences are highly conserved, rich in RCH/TW motifs, and recurrently mutated in DLBCL. Therefore, we hypothesized that aSHM may cause recurrent splicing mutations in DLBCL. In a meta-cohort of > 1,800 DLBCLs, we found that 77.5% of splicing mutations in 29 recurrently mutated genes followed aSHM patterns. In addition, in whole-genome sequencing (WGS) data from 153 DLBCLs, proximal mutations in splice sequences, especially in donors, were significantly enriched in RCH/TW motifs (p < 0.01). We validated this enrichment in two additional DLBCL cohorts (N > 2,000; p < 0.0001) and confirmed its absence in 12 cancer types without aSHM (N > 6,300). Comparing sequencing data from mouse models with and without AID activity showed that the splice donor sequences were the top genomic feature enriched in AID-induced mutations (p < 0.0001). Finally, we observed that most AID-related splice site mutations are clonal within a sample, indicating that aSHM may cause early loss-of-function events in lymphomagenesis. Overall, these findings support that AID causes an overrepresentation of clonal splicing mutations in DLBCL.

Analysis of immunoglobulin/T-cell receptor repertoires by high-throughput RNA sequencing reveals a continuous dynamic of positive clonal selection in follicular lymphoma

Follicular lymphoma (FL) course is highly variable, making its clinical management challenging. In this incurable and recurring pathology, the interval between relapses tends to decrease while aggressiveness increases, sometimes resulting in the transformation to higher-grade lymphoma. These evolutions are particularly difficult to anticipate, resulting from complex clonal evolutions where multiple subclones compete and thrive due to their capacity to proliferate and resist therapies. Here, to apprehend further these processes, we used a high-throughput RNA sequencing approach to address simultaneously the B-cell immunoglobulin repertoires and T-cell immunoglobulin repertoires repertoires of lymphoma cells and their lymphoid microenvironment in a large cohort of 131 FL1/2-3A patients. Our data confirm the existence of a high degree of intra-clonal heterogeneity in this pathology, resulting from ongoing somatic hyper-mutation and class switch recombination. Through the evaluation of the Simpson ecological-diversity index, we show that the contribution of the cancerous cells increases during the course of the disease to the detriment of the reactive compartment, a phenomenon accompanied by a concomitant decrease in the diversity of the tumoral population. Clonal evolution in FL thus contrasts with many tumors, where clonal heterogeneity steadily increases over time and participates in treatment evasion. In this pathology, the selection of lymphoma subclones with proliferative advantages progressively outweighs clonal diversification, ultimately leading in extreme cases to transformation to high-grade lymphoma resulting from the rapid emergence of homogeneous subpopulations.

SurvIAE: Survival prediction with Interpretable Autoencoders from Diffuse Large B-Cells Lymphoma gene expression data

BACKGROUND

In Diffuse Large B-Cell Lymphoma (DLBCL), several methodologies are emerging to derive novel biomarkers to be incorporated in the risk assessment. We realized a pipeline that relies on autoencoders (AE) and Explainable Artificial Intelligence (XAI) to stratify prognosis and derive a gene-based signature.

METHODS

AE was exploited to learn an unsupervised representation of the gene expression (GE) from three publicly available datasets, each with its own technology. Multi-layer perceptron (MLP) was used to classify prognosis from latent representation. GE data were preprocessed as normalized, scaled, and standardized. Four different AE architectures (Large, Medium, Small and Extra Small) were compared to find the most suitable for GE data. The joint AE-MLP classified patients on six different outcomes: overall survival at 12, 36, 60 months and progression-free survival (PFS) at 12, 36, 60 months. XAI techniques were used to derive a gene-based signature aimed at refining the Revised International Prognostic Index (R-IPI) risk, which was validated in a fourth independent publicly available dataset. We named our tool SurvIAE: Survival prediction with Interpretable AE.

RESULTS

From the latent space of AEs, we observed that scaled and standardized data reduced the batch effect. SurvIAE models outperformed R-IPI with Matthews Correlation Coefficient up to 0.42 vs. 0.18 for the validation-set (PFS36) and to 0.30 vs. 0.19 for the test-set (PFS60). We selected the SurvIAE-Small-PFS36 as the best model and, from its gene signature, we stratified patients in three risk groups: R-IPI Poor patients with High levels of GAB1, R-IPI Poor patients with Low levels of GAB1 or R-IPI Good/Very Good patients with Low levels of GPR132, and R-IPI Good/Very Good patients with High levels of GPR132.

CONCLUSIONS

SurvIAE showed the potential to derive a gene signature with translational purpose in DLBCL. The pipeline was made publicly available and can be reused for other pathologies.

Noncoding mutations cause super-enhancer retargeting resulting in protein synthesis dysregulation during B cell lymphoma progression

Whole-genome sequencing of longitudinal tumor pairs representing transformation of follicular lymphoma to high-grade B cell lymphoma with MYC and BCL2 rearrangements (double-hit lymphoma) identified coding and noncoding genomic alterations acquired during lymphoma progression. Many of these transformation-associated alterations recurrently and focally occur at topologically associating domain resident regulatory DNA elements, including H3K4me3 promoter marks located within H3K27ac super-enhancer clusters in B cell non-Hodgkin lymphoma. One region found to undergo recurrent alteration upon transformation overlaps a super-enhancer affecting the expression of the PAX5/ZCCHC7 gene pair. ZCCHC7 encodes a subunit of the Trf4/5-Air1/2-Mtr4 polyadenylation-like complex and demonstrated copy number gain, chromosomal translocation and enhancer retargeting-mediated transcriptional upregulation upon lymphoma transformation. Consequently, lymphoma cells demonstrate nucleolar dysregulation via altered noncoding 5.8S ribosomal RNA processing. We find that a noncoding mutation acquired during lymphoma progression affects noncoding rRNA processing, thereby rewiring protein synthesis leading to oncogenic changes in the lymphoma proteome.

Divergent evolution of metachronous follicular lymphoma and extranodal marginal zone lymphoma of mucosa-associated lymphoid tissue from a common precursor

The translocation t(14;18)(q32:q21)/IGH::BCL2 occurs at the pre-B stage of B-cell development in the bone marrow and is insufficient for malignant transformation, although it leads to the formation of in situ follicular B-cell neoplasia (ISFN). Despite that, the translocation is the genetic hallmark of follicular lymphoma (FL), it occurs infrequently in metachronous/synchronous lymphomas, including extranodal marginal zone lymphoma of mucosa-associated lymphoid tissue (EMZL), mantle cell lymphoma, and Hodgkin's lymphoma. In each of these scenarios, the two lymphomas often appear to be clonally related by analyses of IGH::BCL2 and/or rearranged IG genes. However, it remains largely unknown whether one lymphoma originates from the other or they develop independently. We studied five cases of metachronous EMZL and FL. In four cases, the two lymphomas were clonally related, as shown by identical IGH::BCL2 and/or rearranged IG genes or shared mutations. There were common and unique mutations between the paired EMZL and FL, indicating that they developed independently from a common premalignant cell population, harbouring IGH::BCL2 in three cases. Furthermore, case 1 presented with three metachronous FLs, and all of them originated from a common precursor cell population via divergent evolution. Our findings highlight the multi-malignant potential of IGH::BCL2-positive B-cells. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

The off-target effects of AID in carcinogenesis

Activation-induced cytidine deaminase (AID) plays a crucial role in promoting B cell diversification through somatic hypermutation (SHM) and class switch recombination (CSR). While AID is primarily associated with the physiological function of humoral immune response, it has also been linked to the initiation and progression of lymphomas. Abnormalities in AID have been shown to disrupt gene networks and signaling pathways in both B-cell and T-cell lineage lymphoblastic leukemia, although the full extent of its role in carcinogenesis remains unclear. This review proposes an alternative role for AID and explores its off-target effects in regulating tumorigenesis. In this review, we first provide an overview of the physiological function of AID and its regulation. AID plays a crucial role in promoting B cell diversification through SHM and CSR. We then discuss the off-target effects of AID, which includes inducing mutations of non-Igs, epigenetic modification, and the alternative role as a cofactor. We also explore the networks that keep AID in line. Furthermore, we summarize the off-target effects of AID in autoimmune diseases and hematological neoplasms. Finally, we assess the off-target effects of AID in solid tumors. The primary focus of this review is to understand how and when AID targets specific gene loci and how this affects carcinogenesis. Overall, this review aims to provide a comprehensive understanding of the physiological and off-target effects of AID, which will contribute to the development of novel therapeutic strategies for autoimmune diseases, hematological neoplasms, and solid tumors.

ADAR1-mediated RNA editing promotes B cell lymphomagenesis

Diffuse large B cell lymphoma (DLBCL) is one of the most common types of aggressive lymphoid malignancies. Here, we explore the contribution of RNA editing to DLBCL pathogenesis. We observed that DNA mutations and RNA editing events are often mutually exclusive, suggesting that tumors can modulate pathway outcomes by altering sequences at either the genomic or the transcriptomic level. RNA editing targets transcripts within known disease-driving pathways such as apoptosis, p53 and NF-κB signaling, as well as the RIG-I-like pathway. In this context, we show that ADAR1-mediated editing within MAVS transcript positively correlates with MAVS protein expression levels, associating with increased interferon/NF-κB signaling and T cell exhaustion. Finally, using targeted RNA base editing tools to restore editing within MAVS 3'UTR in ADAR1-deficient cells, we demonstrate that editing is likely to be causal to an increase in downstream signaling in the absence of activation by canonical nucleic acid receptor sensing.

Mesoscale DNA feature in antibody-coding sequence facilitates somatic hypermutation

Somatic hypermutation (SHM), initiated by activation-induced cytidine deaminase (AID), generates mutations in the antibody-coding sequence to allow affinity maturation. Why these mutations intrinsically focus on the three nonconsecutive complementarity-determining regions (CDRs) remains enigmatic. Here, we found that predisposition mutagenesis depends on the single-strand (ss) DNA substrate flexibility determined by the mesoscale sequence surrounding AID deaminase motifs. Mesoscale DNA sequences containing flexible pyrimidine-pyrimidine bases bind effectively to the positively charged surface patches of AID, resulting in preferential deamination activities. The CDR hypermutability is mimicable in in vitro deaminase assays and is evolutionarily conserved among species using SHM as a major diversification strategy. We demonstrated that mesoscale sequence alterations tune the in vivo mutability and promote mutations in an otherwise cold region in mice. Our results show a non-coding role of antibody-coding sequence in directing hypermutation, paving the way for the synthetic design of humanized animal models for optimal antibody discovery and explaining the AID mutagenesis pattern in lymphoma.

DNA repair mechanisms that promote insertion-deletion events during immunoglobulin gene diversification

Insertions and deletions (indels) are low-frequency deleterious genomic DNA alterations. Despite their rarity, indels are common, and insertions leading to long complementarity-determining region 3 (CDR3) are vital for antigen-binding functions in broadly neutralizing and polyreactive antibodies targeting viruses. Because of challenges in detecting indels, the mechanism that generates indels during immunoglobulin diversification processes remains poorly understood. We carried out ultra-deep profiling of indels and systematically dissected the underlying mechanisms using passenger-immunoglobulin mouse models. We found that activation-induced cytidine deaminase-dependent ±1-base pair (bp) indels are the most prevalent indel events, biasing deleterious outcomes, whereas longer in-frame indels, especially insertions that can extend the CDR3 length, are rare outcomes. The ±1-bp indels are channeled by base excision repair, but longer indels require additional DNA-processing factors. Ectopic expression of a DNA exonuclease or perturbation of the balance of DNA polymerases can increase the frequency of longer indels, thus paving the way for models that can generate antibodies with long CDR3. Our study reveals the mechanisms that generate beneficial and deleterious indels during the process of antibody somatic hypermutation and has implications in understanding the detrimental genomic alterations in various conditions, including tumorigenesis.

Minimal information for reporting a genomics experiment

Exome and genome sequencing has facilitated the identification of hundreds of genes and other regions that are recurrently mutated in hematologic neoplasms. The data sets from these studies theoretically provide opportunities. Quality differences between data sets can confound secondary analyses. We explore the consequences of these on the conclusions from some recent studies of B-cell lymphomas. We highlight the need for a minimum reporting standard to increase transparency in genomic research.

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.

Molecular subclusters of follicular lymphoma: a report from the United Kingdom's Haematological Malignancy Research Network

Follicular lymphoma (FL) is morphologically and clinically diverse, with mutations in epigenetic regulators alongside t(14;18) identified as disease-initiating events. Identification of additional mutational entities confirms this cancer's heterogeneity, but whether mutational data can be resolved into mechanistically distinct subsets remains an open question. Targeted sequencing was applied to an unselected population-based FL cohort (n = 548) with full clinical follow-up (n = 538), which included 96 diffuse large B-cell lymphoma (DLBCL) transformations. We investigated whether molecular subclusters of FL can be identified and whether mutational data provide predictive information relating to transformation. DNA extracted from FL samples was sequenced with a 293-gene panel representing genes frequently mutated in DLBCL and FL. Three clusters were resolved using mutational data alone, independent of translocation status: FL_aSHM, with high burden of aberrant somatic hypermutation (aSHM) targets; FL_STAT6, with high STAT6 & CREBBP mutation and low aSHM; and FL_Com, with the absence of features of other subtypes and enriched KMT2D mutation. Analysis of mutation signatures demonstrated differential enrichment of predicted mutation signatures between subgroups and a dominant preference in the FL_aSHM subgroup for G(C>T)T and G(C>T)C transitions consistent with previously defined aSHM-like patterns. Of transformed cases with paired samples, 17 of 26 had evidence of branching evolution. Poorer overall survival (OS) in the aSHM group (P = .04) was associated with older age; however, overall tumor genetics provided limited information to predict individual patient risk. Our approach identifies 3 molecular subclusters of FL linked to differences in underlying mechanistic pathways. These clusters, which may be further resolved by the inclusion of translocation status and wider mutation profiles, have implications for understanding pathogenesis as well as improving treatment strategies in the future.

Genome-wide detection of human variants that disrupt intronic branchpoints

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter.

Unifying Different Cancer Theories in a Unique Tumour Model: Chronic Inflammation and Deaminases as Meeting Points

The increase in cancer incidences shows that there is a need to better understand tumour heterogeneity to achieve efficient treatments. Interestingly, there are several common features among almost all types of cancers, with chronic inflammation induction and deaminase dysfunctions singled out. Deaminases are a family of enzymes with nucleotide-editing capacity, which are classified into two main groups: DNA-based and RNA-based. Remarkably, a close relationship between inflammation and the dysregulation of these molecules has been widely documented, which may explain the characteristic intratumor heterogeneity, both at DNA and transcriptional levels. Indeed, heterogeneity in cancer makes it difficult to establish a unique tumour progression model. Currently, there are three main cancer models—stochastic, hierarchic, and dynamic—although there is no consensus on which one better resembles cancer biology because they are usually overly simplified. Here, to accurately explain tumour progression, we propose interactions among chronic inflammation, deaminases dysregulation, intratumor genetic heterogeneity, cancer phenotypic plasticity, and even the previously proposed appearance of cancer stem-like cell populations in the edges of advanced solid tumour masses (instead of being the cells of origin of primary malignancies). The new tumour development model proposed in this study does not contradict previously accepted models and it may open up a window to interesting therapeutic approaches.

The genomic landscape of canine diffuse large B-cell lymphoma identifies distinct subtypes with clinical and therapeutic implications

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoid neoplasm in dogs and in humans. It is characterized by a remarkable degree of clinical heterogeneity that is not completely elucidated by molecular data. This poses a major barrier to understanding the disease and its response to therapy, or when treating dogs with DLBCL within clinical trials. We performed an integrated analysis of exome (n = 77) and RNA sequencing (n = 43) data in a cohort of canine DLBCL to define the genetic landscape of this tumor. A wide range of signaling pathways and cellular processes were found in common with human DLBCL, but the frequencies of the most recurrently mutated genes (TRAF3, SETD2, POT1, TP53, MYC, FBXW7, DDX3X and TBL1XR1) differed. We developed a prognostic model integrating exonic variants and clinical and transcriptomic features to predict the outcome in dogs with DLBCL. These results comprehensively define the genetic drivers of canine DLBCL and can be prospectively utilized to identify new therapeutic opportunities.

Genomic characterization of lymphomas in patients with inborn errors of immunity

Inborn errors of immunity-associated lymphomas are characterized by distinct clinical features and genetic signatures.

Both germline and somatic alterations contribute to lymphomagenesis in patients with inborn errors of immunity.

Patients with inborn errors of immunity (IEI) have a higher risk of developing cancer, especially lymphoma. However, the molecular basis for IEI-related lymphoma is complex and remains elusive. Here, we perform an in-depth analysis of lymphoma genomes derived from 23 IEI patients. We identified and validated disease-causing or -associated germline mutations in 14 of 23 patients involving ATM, BACH2, BLM, CD70, G6PD, NBN, PIK3CD, PTEN, and TNFRSF13B. Furthermore, we profiled somatic mutations in the lymphoma genome and identified 8 genes that were mutated at a significantly higher level in IEI-associated diffuse large B-cell lymphomas (DLBCLs) than in non-IEI DLBCLs, such as BRCA2, NCOR1, KLF2, FAS, CCND3, and BRWD3. The latter, BRWD3, is furthermore preferentially mutated in tumors of a subgroup of activated phosphoinositide 3-kinase δ syndrome patients. We also identified 5 genomic mutational signatures, including 2 DNA repair deficiency-related signatures, in IEI-associated lymphomas and a strikingly high number of inter- and intrachromosomal structural variants in the tumor genome of a Bloom syndrome patient. In summary, our comprehensive genomic characterization of lymphomas derived from patients with rare genetic disorders expands our understanding of lymphomagenesis and provides new insights for targeted therapy.

Distinct clinical and genetic features of hepatitis B virus-associated follicular lymphoma in Chinese patients

Hepatitis B virus (HBV) infection has been associated with an increased risk for B-cell lymphomas. We previously showed that 20% of diffuse large B-cell lymphoma (DLBCL) patients from China, an endemic area of HBV infection, have chronic HBV infection (surface antigen-positive, HBsAg+) and are characterized by distinct clinical and genetic features. Here, we showed that 24% of follicular lymphoma (FL) Chinese patients are HBsAg+. Compared with the HBsAg- FL patients, HBsAg+ patients are younger, have a higher histological grade at diagnosis, and have a higher incidence of disease progression within 24 months. Moreover, by sequencing the genomes of 109 FL tumors, we observed enhanced mutagenesis and distinct genetic profile in HBsAg+ FLs, with a unique set of preferentially mutated genes (TNFAIP3, FAS, HIST1H1C, KLF2, TP53, PIM1, TMSB4X, DUSP2, TAGAP, LYN, and SETD2) but lack of the hallmark of HBsAg- FLs (ie, IGH/BCL2 translocations and CREBBP mutations). Transcriptomic analyses further showed that HBsAg+ FLs displayed gene-expression signatures resembling the activated B-cell-like subtype of diffuse large B-cell lymphoma, involving IRF4-targeted genes and NF-κB/MYD88 signaling pathways. Finally, we identified an increased infiltration of CD8+ memory T cells, CD4+ Th1 cells, and M1 macrophages and higher T-cell exhaustion gene signature in HBsAg+ FL samples. Taken together, we present new genetic/epigenetic evidence that links chronic HBV infection to B-cell lymphomagenesis, and HBV-associated FL is likely to have a distinct cell-of-origin and represent as a separate subtype of FL. Targetable genetic/epigenetic alterations identified in tumors and their associated tumor microenvironment may provide potential novel therapeutic approaches for this subgroup of patients.

Structure of S1PR2-heterotrimeric G13 signaling complex

Sphingosine-1-phosphate (S1P) regulates immune cell trafficking, angiogenesis, and vascular function via its five receptors. Inherited mutations in S1P receptor 2 (S1PR2) occur in individuals with hearing loss, and acquired mutations in S1PR2 and Gα13 occur in a malignant lymphoma. Here, we present the cryo-electron microscopy structure of S1P-bound S1PR2 coupled to the heterotrimeric G13. Interaction between S1PR2 intracellular loop 2 (ICL2) and transmembrane helix 4 confines ICL2 to engage the α5 helix of Gα13. Transforming growth factor-α shedding assays and cell migration assays support the key roles of the residues in S1PR2-Gα13 complex assembly. The structure illuminates the mechanism of receptor disruption by disease-associated mutations. Unexpectedly, we showed that FTY720-P, an agonist of the other four S1PRs, can trigger G13 activation via S1PR2. S1PR2F274I variant can increase the activity of G13 considerably with FTY720-P and S1P, thus revealing a basis for S1PR drug selectivity.

AID overexpression leads to aggressive murine CLL and nonimmunoglobulin mutations that mirror human neoplasms

Most cancers become more dangerous by the outgrowth of malignant subclones with additional DNA mutations that favor proliferation or survival. Using chronic lymphocytic leukemia (CLL), a disease that exemplifies this process and is a model for neoplasms in general, we created transgenic mice overexpressing the enzyme activation-induced deaminase (AID), which has a normal function of inducing DNA mutations in B lymphocytes. AID not only allows normal B lymphocytes to develop more effective immunoglobulin-mediated immunity, but is also able to mutate nonimmunoglobulin genes, predisposing to cancer. In CLL, AID expression correlates with poor prognosis, suggesting a role for this enzyme in disease progression. Nevertheless, direct experimental evidence identifying the specific genes that are mutated by AID and indicating that those genes are associated with disease progression is not available. To address this point, we overexpressed Aicda in a murine model of CLL (Eμ-TCL1). Analyses of TCL1/AID mice demonstrate a role for AID in disease kinetics, CLL cell proliferation, and the development of cancer-related target mutations with canonical AID signatures in nonimmunoglobulin genes. Notably, our mouse models can accumulate mutations in the same genes that are mutated in human cancers. Moreover, some of these mutations occur at homologous positions, leading to identical or chemically similar amino acid substitutions as in human CLL and lymphoma. Together, these findings support a direct link between aberrant AID activity and CLL driver mutations that are then selected for their oncogenic effects, whereby AID promotes aggressiveness in CLL and other B-cell neoplasms.

Genome-wide CRISPR screens reveal synthetic lethal interaction between CREBBP and EP300 in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common type of aggressive lymphoid malignancy and a highly heterogeneous disease. In this study, we performed whole-genome and transcriptome sequencing, and a genome-wide CRISPR-Cas9-knockout screen to study an activated B-cell-like DLBCL cell line (RC-K8). We identified a distinct pattern of genetic essentialities in RC-K8, including a dependency on CREBBP and MDM2. The dependency on CREBBP is associated with a balanced translocation involving EP300, which results in a truncated form of the protein that lacks the critical histone acetyltransferase (HAT) domain. The synthetic lethal interaction between CREBBP and EP300 genes, two frequently mutated epigenetic modulators in B-cell lymphoma, was further validated in the previously published CRISPR-Cas9 screens and inhibitor assays. Our study suggests that integration of the unbiased functional screen results with genomic and transcriptomic data can identify both common and unique druggable vulnerabilities in DLBCL and histone acetyltransferases inhibition could be a therapeutic option for CREBBP or EP300 mutated cases.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.