STAT6 mutations enriched at diffuse large B-cell lymphoma relapse reshape the tumor microenvironment

Diffuse large B-cell lymphoma (DLBCL) relapses in approximately 40% of patients following frontline therapy. We reported that STAT6D419 mutations are enriched in relapsed/refractory DLBCL (rrDLBCL) samples, suggesting that JAK/STAT signaling plays a role in therapeutic resistance. We hypothesized that STAT6D419 mutations can improve DLBCL cell survival by reprogramming the microenvironment to sustain STAT6 activation. Thus, we investigated the role of STAT6D419 mutations on DLBCL cell growth and its microenvironment. We found that phospho-STAT6D419N was retained in the nucleus longer than phospho-STAT6WT following IL-4 stimulation, and STAT6D419N recognized a more restricted DNA-consensus sequence than STAT6WT. Upon IL-4 induction, STAT6D419N expression led to a higher magnitude of gene expression changes, but in a more selective list of gene targets compared with STATWT. The most significantly expressed genes induced by STAT6D419N were those implicated in survival, proliferation, migration, and chemotaxis, in particular CCL17. This chemokine, also known as TARC, attracts helper T-cells to the tumor microenvironment, especially in Hodgkin's lymphoma. To this end, in DLBCL, phospho-STAT6+ rrDLBCL cells had a greater proportion of infiltrating CD4+ T-cells than phospho-STAT6- tumors. Our findings suggest that STAT6D419 mutations in DLBCL lead to cell autonomous changes, enhanced signaling, and altered composition of the tumor microenvironment.

Biological heterogeneity in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is heterogeneous both in clinical outcomes and the underlying disease biology. Over the last 2 decades, several different approaches for dissecting biological heterogeneity have emerged. Gene expression profiling (GEP) stratifies DLBCL into 3 broad groups (ABC, GCB, and DZsig/MHG), each with parallels to different normal mature B cell developmental states and prognostic implications. More recently, several different genomic approaches have been developed to categorize DLBCL based on the co-occurrence of tumor somatic mutations, identifying more granular biologically unified subgroups that complement GEP-based approaches. We review the molecular approaches and clinical evidence supporting the stratification of DLBCL patients based on tumor biology. By offering a platform for subtype-guided therapy, these divisions remain a promising avenue for improving patient outcomes, especially in subgroups with inferior outcomes with current standard-of-care therapy.

Relapse Timing Is Associated With Distinct Evolutionary Dynamics in Diffuse Large B-Cell Lymphoma

PURPOSE

Diffuse large B-cell lymphoma (DLBCL) is cured in more than 60% of patients, but outcomes remain poor for patients experiencing disease progression or relapse (refractory or relapsed DLBCL [rrDLBCL]), particularly if these events occur early. Although previous studies examining cohorts of rrDLBCL have identified features that are enriched at relapse, few have directly compared serial biopsies to uncover biological and evolutionary dynamics driving rrDLBCL. Here, we sought to confirm the relationship between relapse timing and outcomes after second-line (immuno)chemotherapy and determine the evolutionary dynamics that underpin that relationship.

PATIENTS AND METHODS

Outcomes were examined in a population-based cohort of 221 patients with DLBCL who experienced progression/relapse after frontline treatment and were treated with second-line (immuno)chemotherapy with an intention-to-treat with autologous stem-cell transplantation (ASCT). Serial DLBCL biopsies from a partially overlapping cohort of 129 patients underwent molecular characterization, including whole-genome or whole-exome sequencing in 73 patients.

RESULTS

Outcomes to second-line therapy and ASCT are superior for late relapse (>2 years postdiagnosis) versus primary refractory (<9 months) or early relapse (9-24 months). Diagnostic and relapse biopsies were mostly concordant for cell-of-origin classification and genetics-based subgroup. Despite this concordance, the number of mutations exclusive to each biopsy increased with time since diagnosis, and late relapses shared few mutations with their diagnostic counterpart, demonstrating a branching evolution pattern. In patients with highly divergent tumors, many of the same genes acquired new mutations independently in each tumor, suggesting that the earliest mutations in a shared precursor cell constrain tumor evolution toward the same genetics-based subgroups at both diagnosis and relapse.

CONCLUSION

These results suggest that late relapses commonly represent genetically distinct and chemotherapy-naïve disease and have implications for optimal patient management.

Genetic subdivisions of follicular lymphoma defined by distinct coding and noncoding mutation patterns

Follicular lymphoma (FL) accounts for ∼20% of all new lymphoma cases. Increases in cytological grade are a feature of the clinical progression of this malignancy, and eventual histologic transformation (HT) to the aggressive diffuse large B-cell lymphoma (DLBCL) occurs in up to 15% of patients. Clinical or genetic features to predict the risk and timing of HT have not been described comprehensively. In this study, we analyzed whole-genome sequencing data from 423 patients to compare the protein coding and noncoding mutation landscapes of untransformed FL, transformed FL, and de novo DLBCL. This revealed 2 genetically distinct subgroups of FL, which we have named DLBCL-like (dFL) and constrained FL (cFL). Each subgroup has distinguishing mutational patterns, aberrant somatic hypermutation rates, and biological and clinical characteristics. We implemented a machine learning-derived classification approach to stratify patients with FL into cFL and dFL subgroups based on their genomic features. Using separate validation cohorts, we demonstrate that cFL status, whether assigned with this full classifier or a single-gene approximation, is associated with a reduced rate of HT. This implies distinct biological features of cFL that constrain its evolution, and we highlight the potential for this classification to predict HT from genetic features present at diagnosis.

Relapse timing is associated with distinct evolutionary dynamics in DLBCL

Diffuse large B-cell lymphoma (DLBCL) is cured in over 60% of patients, but outcomes are poor for patients with relapsed or refractory disease (rrDLBCL). Here, we performed whole genome/exome sequencing (WGS/WES) on tumors from 73 serially-biopsied patients with rrDLBCL. Based on the observation that outcomes to salvage therapy/autologous stem cell transplantation are related to time-to-relapse, we stratified patients into groups according to relapse timing to explore the relationship to genetic divergence and sensitivity to salvage immunochemotherapy. The degree of mutational divergence increased with time between biopsies, yet tumor pairs were mostly concordant for cell-of-origin, oncogene rearrangement status and genetics-based subgroup. In patients with highly divergent tumors, several genes acquired exclusive mutations independently in each tumor, which, along with concordance of genetics-based subgroups, suggests that the earliest mutations in a shared precursor cell constrain tumor evolution. These results suggest that late relapses commonly represent genetically distinct and chemotherapy-naïve disease.

Genetic subgroups inform on pathobiology in adult and pediatric Burkitt lymphoma

Burkitt lymphoma (BL) accounts for most pediatric non-Hodgkin lymphomas, being less common but significantly more lethal when diagnosed in adults. Much of the knowledge of the genetics of BL thus far has originated from the study of pediatric BL (pBL), leaving its relationship to adult BL (aBL) and other adult lymphomas not fully explored. We sought to more thoroughly identify the somatic changes that underlie lymphomagenesis in aBL and any molecular features that associate with clinical disparities within and between pBL and aBL. Through comprehensive whole-genome sequencing of 230 BL and 295 diffuse large B-cell lymphoma (DLBCL) tumors, we identified additional significantly mutated genes, including more genetic features that associate with tumor Epstein-Barr virus status, and unraveled new distinct subgroupings within BL and DLBCL with 3 predominantly comprising BLs: DGG-BL (DDX3X, GNA13, and GNAI2), IC-BL (ID3 and CCND3), and Q53-BL (quiet TP53). Each BL subgroup is characterized by combinations of common driver and noncoding mutations caused by aberrant somatic hypermutation. The largest subgroups of BL cases, IC-BL and DGG-BL, are further characterized by distinct biological and gene expression differences. IC-BL and DGG-BL and their prototypical genetic features (ID3 and TP53) had significant associations with patient outcomes that were different among aBL and pBL cohorts. These findings highlight shared pathogenesis between aBL and pBL, and establish genetic subtypes within BL that serve to delineate tumors with distinct molecular features, providing a new framework for epidemiologic, diagnostic, and therapeutic strategies.

Emerging roles for heterogeneous ribonuclear proteins in normal and malignant B cells

Heterogeneous nuclear ribonucleoproteins (hnRNPs) are among the most abundantly expressed RNA binding proteins in the cell and play major roles in all facets of RNA metabolism. hnRNPs are increasingly appreciated as essential for mammalian B cell development by regulating the carefully ordered expression of specific genes. Due to this tight regulation of the hnRNP-RNA network, it is no surprise that a growing number of genes encoding hnRNPs have been causally associated with the onset or progression of many cancers, including B cell neoplasms. Here we discuss our current understanding of hnRNP-driven regulation in normal, perturbed, and malignant B cells, and the most recent and emerging therapeutic innovations aimed at targeting the hnRNP-RNA network in lymphoma.

PeptideRanger: An R Package to Optimize Synthetic Peptide Selection for Mass Spectrometry Applications

Targeted and semitargeted mass spectrometry-based approaches are reliable methods to consistently detect and quantify low abundance proteins including proteins of clinical significance. Despite their potential, the development of targeted and semitargeted assays is time-consuming and often requires the purchase of costly libraries of synthetic peptides. To improve the efficiency of this rate-limiting step, we developed PeptideRanger, a tool to identify peptides from protein of interest with physiochemical properties that make them more likely to be suitable for mass spectrometry analysis. PeptideRanger is a flexible, extensively annotated, and intuitive R package that uses a random forest model trained on a diverse data set of thousands of MS experiments spanning a variety of sample types profiled with different chromatography setups and instruments. To support a variety of applications and to leverage rapidly growing public MS databases, PeptideRanger can readily be retrained with experiment-specific data sets and customized to prioritize and filter peptides based on selected properties.

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.

Genomic profiling for clinical decision making in lymphoid neoplasms

With the introduction of large-scale molecular profiling methods and high-throughput sequencing technologies, the genomic features of most lymphoid neoplasms have been characterized at an unprecedented scale. Although the principles for the classification and diagnosis of these disorders, founded on a multidimensional definition of disease entities, have been consolidated over the past 25 years, novel genomic data have markedly enhanced our understanding of lymphomagenesis and enriched the description of disease entities at the molecular level. Yet, the current diagnosis of lymphoid tumors is largely based on morphological assessment and immunophenotyping, with only few entities being defined by genomic criteria. This paper, which accompanies the International Consensus Classification of mature lymphoid neoplasms, will address how established assays and newly developed technologies for molecular testing already complement clinical diagnoses and provide a novel lens on disease classification. More specifically, their contributions to diagnosis refinement, risk stratification, and therapy prediction will be considered for the main categories of lymphoid neoplasms. The potential of whole-genome sequencing, circulating tumor DNA analyses, single-cell analyses, and epigenetic profiling will be discussed because these will likely become important future tools for implementing precision medicine approaches in clinical decision making for patients with lymphoid malignancies.

Abstract A12: Pervasive hypermutation of super-enhancer regions dysregulates oncogene expression in diffuse large B-cell lymphoma

Coding-genome sequencing efforts have identified several genes/pathways altered in Diffuse Large B-cell Lymphoma (DLBCL), including new potential therapeutic targets. However, the mutational contribution of the non-coding genome of DLBCL remains largely unexplored. To identify functionally relevant non-coding mutations targeting regulatory domains, we integrated enhancer (E)/super-enhancer (SE) identification, by ChIP-seq analysis of the hallmark H3K27 histone mark, with whole genome sequencing (WGS) and RNA-seq analysis in 29 DLBCL cell lines representative of the major DLBCL subtypes, along with germinal center (GC) B cells representing the normal counterpart of DLBCL. Following validation in 93 normal/tumor DLBCL biopsies, we found that active SEs are specifically hypermutated (≥3 somatic mutations/Kb) in 99% of DLBCLs, as compared to the same loci when not active as SE (in total, 159 hypermutated SEs, with at least 2/case). As evidence of oncogenic relevance, we have shown that the hypermutated SEs linked to the BCL6 and BCL2 proto-oncogenes prevent the binding and transcriptional downregulation by transcriptional repressors. Genetic correction of selected mutations using the CRISPR/Cas9 technology restored repressor DNA-binding, downregulated target gene expression, and led to the counter-selection of cells harboring corrected alleles, indicating oncogenic dependency from the SE mutations. A third recurrently hypermutated SE was the CXCR4-SE, affected in 19% of cases. CXCR4 encodes for a G-protein coupled chemokine receptor involved in cell migration within the GC, and is mutationally activated in 40% of Waldenström macrobulinemia (WM), a post-GC lymphoproliferative disease. We identified a mutational hotspot preventing the binding of the NR3C1 glucocorticoid receptor in 7.5% of DLBCL cases. Correction of the mutations in cell lines led to counter selection and reduced expression of CXCR4. Conversely, the introduction of a mutation disrupting NR3C1 binding in WT DLBCL cells increased CXCR4 expression. Together, these results are consistent with dependency of the mutant cells on the CXCR4-SE mutation, and confirm a direct link between SE mutations in this region and deregulated gene expression through escape from NR3C1-mediated suppression. In primary cases, mutations in the NR3C1 binding site of the CXCR4-SE were observed in the absence of NR3C1 coding mutations and correlated with increased CXCR4 transcript levels, as documented by RNA-seq. Together with the oncogenic role of CXCR4 in WM, these findings reveal this chemokine receptor as a potential oncogenic target in DLBCL. Collectively, these findings reveal a new layer of genetic alterations, which identifies novel mechanisms of dysregulation for known oncogenes, as well as new dysregulated genes and pathways, with implications for precision classification and therapeutic targeting of DLBCL. * Co-senior authors.

Citation Format: Elodie Bal, Rahul Kumar, Mohammad Hadigol, Antony B. Holmes, Laura K. Hilton, Jui Wan Loh, Kostiantyn Dreval, Jasper Wong, Sofija Vlasevska, Clarissa Corinaldesi, Rajesh Kumar Soni Soni, Katia Basso, Ryan D. Morin, Hossein Khiabanian, Laura Pasqualucci, Riccardo Dalla-Favera. Pervasive hypermutation of super-enhancer regions dysregulates oncogene expression in diffuse large B-cell lymphoma [abstract]. In: Proceedings of the Third AACR International Meeting: Advances in Malignant Lymphoma: Maximizing the Basic-Translational Interface for Clinical Application; 2022 Jun 23-26; Boston, MA. Philadelphia (PA): AACR; Blood Cancer Discov 2022;3(5_Suppl):Abstract nr A12.

Molecular determinants of outcomes in relapsed or refractory mantle cell lymphoma treated with ibrutinib or temsirolimus in the MCL3001 (RAY) trial

Mantle cell lymphoma (MCL) is a rare, incurable lymphoma subtype characterized by heterogeneous outcomes. To better understand the clinical behavior and response to treatment, predictive biomarkers are needed. Using residual archived material from patients enrolled in the MCL3001 (RAY) study, we performed detailed analyses of gene expression and targeted genetic sequencing. This phase III clinical trial randomized patients with relapsed or refractory MCL to treatment with either ibrutinib or temsirolimus. We confirmed the prognostic capability of the gene expression proliferation assay MCL35 in this cohort treated with novel agents; it outperformed the simplified MCL International Prognostic Index in discriminating patients with different outcomes. Regardless of treatment arm, our data demonstrated that this assay captures the risk conferred by known biological factors, including increased MYC expression, blastoid morphology, aberrations of TP53, and truncated CCND1 3' untranslated region. We showed the negative impact of BIRC3 mutations/deletions on outcomes in this cohort and identified that deletion of chromosome 8p23.3 also negatively impacts survival. Restricted to patients with deletions/alterations in TP53, ibrutinib appeared to abrogate the deleterious impact on outcome. These data illustrate the potential to perform a molecular analysis of predictive biomarkers on routine patient samples that can meaningfully inform clinical practice.

BCL10 Mutations Define Distinct Dependencies Guiding Precision Therapy for DLBCL

Activated B cell-like diffuse large B-cell lymphomas (ABC-DLBCL) have unfavorable outcomes and chronic activation of CARD11-BCL10-MALT1 (CBM) signal amplification complexes that form due to polymerization of BCL10 subunits, which is affected by recurrent somatic mutations in ABC-DLBCLs. Herein, we show that BCL10 mutants fall into at least two functionally distinct classes: missense mutations of the BCL10 CARD domain and truncation of its C-terminal tail. Truncating mutations abrogated a motif through which MALT1 inhibits BCL10 polymerization, trapping MALT1 in its activated filament-bound state. CARD missense mutations enhanced BCL10 filament formation, forming glutamine network structures that stabilize BCL10 filaments. Mutant forms of BCL10 were less dependent on upstream CARD11 activation and thus manifested resistance to BTK inhibitors, whereas BCL10 truncating but not CARD mutants were hypersensitive to MALT1 inhibitors. Therefore, BCL10 mutations are potential biomarkers for BTK inhibitor resistance in ABC-DLBCL, and further precision can be achieved by selecting therapy based on specific biochemical effects of distinct mutation classes.

SIGNIFICANCE

ABC-DLBCLs feature frequent mutations of signaling mediators that converge on the CBM complex. We use structure-function approaches to reveal that BCL10 mutations fall into two distinct biochemical classes. Both classes confer resistance to BTK inhibitors, whereas BCL10 truncations confer hyperresponsiveness to MALT1 inhibitors, providing a road map for precision therapies in ABC-DLBCLs. See related commentary by Phelan and Oellerich, p. 1844. This article is highlighted in the In This Issue feature, p. 1825.

Super-enhancer hypermutation alters oncogene expression in B cell lymphoma

Diffuse large B cell lymphoma (DLBCL) is the most common B cell non-Hodgkin lymphoma and remains incurable in around 40% of patients. Efforts to sequence the coding genome identified several genes and pathways that are altered in this disease, including potential therapeutic targets1-5. However, the non-coding genome of DLBCL remains largely unexplored. Here we show that active super-enhancers are highly and specifically hypermutated in 92% of samples from individuals with DLBCL, display signatures of activation-induced cytidine deaminase activity, and are linked to genes that encode B cell developmental regulators and oncogenes. As evidence of oncogenic relevance, we show that the hypermutated super-enhancers linked to the BCL6, BCL2 and CXCR4 proto-oncogenes prevent the binding and transcriptional downregulation of the corresponding target gene by transcriptional repressors, including BLIMP1 (targeting BCL6) and the steroid receptor NR3C1 (targeting BCL2 and CXCR4). Genetic correction of selected mutations restored repressor DNA binding, downregulated target gene expression and led to the counter-selection of cells containing corrected alleles, indicating an oncogenic dependency on the super-enhancer mutations. This pervasive super-enhancer mutational mechanism reveals a major set of genetic lesions deregulating gene expression, which expands the involvement of known oncogenes in DLBCL pathogenesis and identifies new deregulated gene targets of therapeutic relevance.

Combinatorial and Machine Learning Approaches for Improved Somatic Variant Calling From Formalin-Fixed Paraffin-Embedded Genome Sequence Data

Formalin fixation of paraffin-embedded tissue samples is a well-established method for preserving tissue and is routinely used in clinical settings. Although formalin-fixed, paraffin-embedded (FFPE) tissues are deemed crucial for research and clinical applications, the fixation process results in molecular damage to nucleic acids, thus confounding their use in genome sequence analysis. Methods to improve genomic data quality from FFPE tissues have emerged, but there remains significant room for improvement. Here, we use whole-genome sequencing (WGS) data from matched Fresh Frozen (FF) and FFPE tissue samples to optimize a sensitive and precise FFPE single nucleotide variant (SNV) calling approach. We present methods to reduce the prevalence of false-positive SNVs by applying combinatorial techniques to five publicly available variant callers. We also introduce FFPolish, a novel variant classification method that efficiently classifies FFPE-specific false-positive variants. Our combinatorial and statistical techniques improve precision and F1 scores compared to the results of publicly available tools when tested individually.

Mutated RAS-associating proteins and ERK activation in relapse/refractory diffuse large B cell lymphoma

Diffuse large B cell lymphoma (DLBCL) is successfully treated with combination immuno-chemotherapy, but relapse with resistant disease occurs in ~ 40% of patients. However, little is known regarding relapsed/refractory DLBCL (rrDLBCL) genetics and alternative therapies. Based on findings from other tumors, we hypothesized that RAS-MEK-ERK signaling would be upregulated in resistant tumors, potentially correlating with mutations in RAS, RAF, or associated proteins. We analyzed mutations and phospho-ERK levels in tumor samples from rrDLBCL patients. Unlike other tumor types, rrDLBCL is not mutated in any Ras or Raf family members, despite having increased expression of p-ERK. In paired biopsies comparing diagnostic and relapsed specimens, 33% of tumors gained p-ERK expression, suggesting a role in promoting survival. We did find mutations in several Ras-associating proteins, including GEFs, GAPs, and downstream effectors that could account for increased ERK activation. We further investigated mutations in one such protein, RASGRP4. In silico modeling indicated an increased interaction between H-Ras and mutant RASGRP4. In cell lines, mutant RASGRP4 increased basal p-ERK expression and lead to a growth advantage in colony forming assays when challenged with doxorubicin. Relapsed/refractory DLBCL is often associated with increased survival signals downstream of ERK, potentially corresponding with mutations in protein controlling RAS/MEK/ERK signaling.

Evolutionary conservation of systemic and reversible amyloid aggregation

In response to environmental stress, human cells have been shown to form reversible amyloid aggregates within the nucleus, termed amyloid bodies (A-bodies). These protective physiological structures share many of the biophysical characteristics associated with the pathological amyloids found in Alzheimer's and Parkinson's disease. Here, we show that A-bodies are evolutionarily conserved across the eukaryotic domain, with their detection in Drosophila melanogaster and Saccharomyces cerevisiae marking the first examples of these functional amyloids being induced outside of a cultured cell setting. The conditions triggering amyloidogenesis varied significantly among the species tested, with results indicating that A-body formation is a severe, but sublethal, stress response pathway that is tailored to the environmental norms of an organism. RNA-sequencing analyses demonstrate that the regulatory low-complexity long non-coding RNAs that drive A-body aggregation are both conserved and essential in human, mouse and chicken cells. Thus, the identification of these natural and reversible functional amyloids in a variety of evolutionarily diverse species highlights the physiological significance of this protein conformation, and will be informative in advancing our understanding of both functional and pathological amyloid aggregation events. This article has an associated First Person interview with the first author of the paper.

Proteotranscriptomic classification and characterization of pancreatic neuroendocrine neoplasms

Pancreatic neuroendocrine neoplasms (PNENs) are biologically and clinically heterogeneous. Here, we use a multi-omics approach to uncover the molecular factors underlying this heterogeneity. Transcriptomic analysis of 84 PNEN specimens, drawn from two cohorts, is substantiated with proteomic profiling and identifies four subgroups: Proliferative, PDX1-high, Alpha cell-like and Stromal/Mesenchymal. The Proliferative subgroup, consisting of both well- and poorly differentiated specimens, is associated with inferior overall survival probability. The PDX1-high and Alpha cell-like subgroups partially resemble previously described subtypes, and we further uncover distinctive metabolism-related features in the Alpha cell-like subgroup. The Stromal/Mesenchymal subgroup exhibits molecular characteristics of YAP1/WWTR1(TAZ) activation suggestive of Hippo signaling pathway involvement in PNENs. Whole-exome sequencing reveals subgroup-enriched mutational differences, supported by activity inference analysis, and identifies hypermorphic proto-oncogene variants in 14.3% of sequenced PNENs. Our study reveals differences in cellular signaling axes that provide potential directions for PNEN patient stratification and treatment strategies.

Molecular profiling in diffuse large B-cell lymphoma: why so many types of subtypes?

The term diffuse large B-cell lymphoma (DLBCL) includes a heterogeneous collection of biologically distinct tumours. This heterogeneity currently presents a barrier to the successful deployment of novel, biologically targeted therapies. Molecular profiling studies have recently proposed new molecular classification systems. These have the potential to resolve the biological heterogeneity of DLBCL into manageable subgroups of tumours that rely on shared oncogenic programmes. In many cases these biological programmes straddle the boundaries of our existing systems for classifying B-cell lymphomas. Here we review the findings from these major molecular profiling studies with a specific focus on those that propose new genetic subgroups of DLBCL. We highlight the areas of consensus and discordance between these studies and discuss the implications for current clinical practice and for clinical trials. Finally, we address the outstanding challenges and solutions to the introduction of genomic subtyping and precision medicine in DLBCL.

Prognostic significance of FCGR2B expression for the response of DLBCL patients to rituximab or obinutuzumab treatment

Fc γ receptor IIB (FcγRIIB) is an inhibitory molecule capable of reducing antibody immunotherapy efficacy. We hypothesized its expression could confer resistance in patients with diffuse large B-cell lymphoma (DLBCL) treated with anti-CD20 monoclonal antibody (mAb) chemoimmunotherapy, with outcomes varying depending on mAb (rituximab [R]/obinutuzumab [G]) because of different mechanisms of action. We evaluated correlates between FCGR2B messenger RNA and/or FcγRIIB protein expression and outcomes in 3 de novo DLBCL discovery cohorts treated with R plus cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) reported by Arthur, Schmitz, and Reddy, and R-CHOP/G-CHOP-treated patients in the GOYA trial (NCT01287741). In the discovery cohorts, higher FCGR2B expression was associated with significantly shorter progression-free survival (PFS; Arthur: hazard ratio [HR], 1.09; 95% confidence interval [CI], 1.01-1.19; P = .0360; Schmitz: HR, 1.13; 95% CI, 1.02-1.26; P = .0243). Similar results were observed in GOYA with R-CHOP (HR, 1.26; 95% CI, 1.00-1.58; P = .0455), but not G-CHOP (HR, 0.91; 95% CI, 0.69-1.20; P = .50). A nonsignificant trend that high FCGR2B expression favored G-CHOP over R-CHOP was observed (HR, 0.67; 95% CI, 0.44-1.02; P = .0622); however, low FCGR2B expression favored R-CHOP (HR, 1.58; 95% CI, 1.00-2.50; P = .0503). In Arthur and GOYA, FCGR2B expression was associated with tumor FcγRIIB expression; correlating with shorter PFS for R-CHOP (HR, 2.17; 95% CI, 1.04-4.50; P = .0378), but not G-CHOP (HR, 1.37; 95% CI, 0.66-2.87; P = .3997). This effect was independent of established prognostic biomarkers. High FcγRIIB/FCGR2B expression has prognostic value in R-treated patients with DLBCL and may confer differential responsiveness to R-CHOP/G-CHOP.

Mast Cell and Eosinophil Activation Are Associated With COVID-19 and TLR-Mediated Viral Inflammation: Implications for an Anti-Siglec-8 Antibody

Coronavirus disease 2019 (COVID-19) caused by SARS-CoV-2 infection represents a global health crisis. Immune cell activation via pattern recognition receptors has been implicated as a driver of the hyperinflammatory response seen in COVID-19. However, our understanding of the specific immune responses to SARS-CoV-2 remains limited. Mast cells (MCs) and eosinophils are innate immune cells that play pathogenic roles in many inflammatory responses. Here we report MC-derived proteases and eosinophil-associated mediators are elevated in COVID-19 patient sera and lung tissues. Stimulation of viral-sensing toll-like receptors in vitro and administration of synthetic viral RNA in vivo induced features of hyperinflammation, including cytokine elevation, immune cell airway infiltration, and MC-protease production—effects suppressed by an anti-Siglec-8 monoclonal antibody which selectively inhibits MCs and depletes eosinophils. Similarly, anti-Siglec-8 treatment reduced disease severity and airway inflammation in a respiratory viral infection model. These results suggest that MC and eosinophil activation are associated with COVID-19 inflammation and anti-Siglec-8 antibodies are a potential therapeutic approach for attenuating excessive inflammation during viral infections.

Abstract PO-26: Prognostic significance of Fc gamma receptor IIB expression in the response of previously untreated diffuse large B-cell lymphomas to anti-CD20 monoclonal antibodies: Differing impact of rituximab and obinutuzumab

Background: The anti-CD20 monoclonal antibody (mAb), obinutuzumab (G), has shown improved outcomes versus rituximab (R) in indolent lymphomas; however, no improvement was seen in diffuse large B-cell lymphoma (DLBCL), and the molecular basis is unclear. The inhibitory Fc gamma receptor IIB (Fc γRIIB), expressed on lymphoma cells, can impair the effects of direct targeting mAbs, such as R, by binding and internalizing them, reducing opsonization and limiting FcγR-mediated killing. We hypothesized that FcγRIIB expression on cellular effectors and/or the lymphoma confers treatment resistance in some patients (pts) and evaluated if outcomes differ for therapies involving a non-internalized mAb (G).

Methods: We evaluated correlates between FCGR2B mRNA and/or FcγRIIB protein expression and pt outcomes in two discovery cohorts of de novo DLBCL treated with R-CHOP (Arthur et al. 2018, n=372; Schmitz et al. 2018, n=234), and the phase III GOYA trial (NCT01287741; n=552), which compared R-CHOP with G-CHOP in pts with previously untreated DLBCL. FCGR2B mRNA expression was assessed by RNA-Seq, and protein expression was assessed in evaluable cohorts by immunohistochemistry, using tissue microarrays with macrophages identified by CD68. FCGR2B expression was also measured by a NanoString assay (Arthur cohort). Cox regression analyzed the impact of FcγRIIB/FCGR2B expression on progression-free survival (PFS), with univariate and multivariate models adjusted for International Prognostic Index (IPI), cell of origin (COO), and BCL2 protein expression.

Results: In the discovery cohorts, a higher FCGR2B expression was significantly associated with shorter PFS (Arthur: HR 1.09 [95% CI: 1.01–1.19], P=0.036; Schmitz: HR 1.13 [95% CI: 1.02–1.26], P=0.0243). Expression by NanoString strongly correlated with RNA-Seq, confirming the association with shorter PFS (HR 1.13 [95% CI: 1.04–1.23], P=0.0048). In GOYA, a significant association between PFS and FCGR2B was observed in the R arm (HR 1.26 [95% CI: 1.00–1.58], P=0.0455), with no prognostic effect observed for G (HR 0.91 [95% CI: 0.69–1.20], P=0.5). Pts with high FCGR2B expression appeared to benefit more from G than R (HR 0.67 [95% CI: 0.44–1.02], P=0.0622), in contrast to pts with low FCGR2B expression (HR 1.58 [95% CI: 1.00–2.50], P=0.0503). In both Arthur and GOYA cohorts, FCGR2B expression by RNA-Seq was associated with FcγRIIB on the tumor, which correlated with a shorter PFS for R (HR 2.17 [95% CI: 1.04–4.50], P=0.03), but not G (HR 1.37 [95% CI: 0.66–2.87], P=0.4). This prognostic effect on PFS was independent of established prognostic biomarkers, IPI, COO and BCL2.

Conclusion: High FcγRIIB/FCGR2B expression in pts with DLBCL has prognostic value in those treated with R and may confer differential responsiveness to R or G.

Citation Format: Laura K. Hilton, Malgorzata Nowicka, Margaret Ashton-Key, Chantal E. Hargreaves, Chern Lee, Russell Foxall, Matthew J. Carter, Stephen A. Beers, Kathleen N. Potter, Christopher R. Bolen, Christian Klein, Andrea Knapp, Farheen Mir, Matthew Rose-Zerilli, Cathy Burton, Wolfram Klapper, David W. Scott, Laurie H. Sehn, Umberto Vitolo, Maurizio Martelli, Marek Trneny, Graham W. Slack, Pedro Farinha, Jonathan C. Strefford, Mikkel Z. Oestergaard, Ryan D. Morin, Mark S. Cragg. Prognostic significance of Fc gamma receptor IIB expression in the response of previously untreated diffuse large B-cell lymphomas to anti-CD20 monoclonal antibodies: Differing impact of rituximab and obinutuzumab [abstract]. In: Proceedings of the AACR Virtual Meeting: Advances in Malignant Lymphoma; 2020 Aug 17-19. Philadelphia (PA): AACR; Blood Cancer Discov 2020;1(3_Suppl):Abstract nr PO-26.

Abstract PO-32: NFKBIZ 3′ UTR mutations confer selective growth advantage and affect drug response in diffuse large B-cell lymphoma

Introduction: The activated B cell-like (ABC) molecular subgroup of diffuse large B-cell lymphoma (DLBCL) is characterized by activation of NF-κB signaling. Recurrent mutations affecting genes such as MYD88, CD79A/B, and TNFAIP3 contribute to this in some cases, but there remain tumors with no known genetic basis for this pathway activation. This suggests that our understanding of ABC DLBCL drivers remains incomplete. Previously, NFKBIZ was shown to be amplified in 10% of ABC DLBCLs and to contribute to activation of NF-κB signaling. We recently described a novel pattern of mutations affecting the 3′ UTR of NFKBIZ resulting in an overall mutation rate of 30% (UTR or AMP) in ABC DLBCL. These NFKBIZ UTR mutations are mutually exclusive with MYD88 mutations, thus suggesting they may also lead to activation of NF-κB signaling. The NFKBIZ protein interacts with NF-κB transcription factors and is thought to regulate canonical NF-κB signaling. We hypothesized that NFKBIZ UTR mutations affect the normally rapid degradation of this mRNA by disrupting secondary structures recognized by RNA-binding proteins such as ribonucleases. The resulting elevated NFKBIZ mRNA levels would lead to accumulation of protein and may be a novel mechanism to promote cell growth and survival in ABC DLBCL.

Methods: NFKBIZ 3′ UTR mutations were introduced into a DLBCL cell line using the CRISPR-Cas9 system. A competitive growth assay with wild-type (WT) and CRISPR-mutant lines was performed to assess whether UTR mutations provide a growth advantage in culture (in vitro) and in mouse xenografts (in vivo). RNA-sequencing was then performed on WT and a subset of CRISPR-mutant lines and analyses were performed to identify genes upregulated by IκB-ζ in mutant lines. The IC50 of relevant drugs was determined by WST-1 assays after drug treatment on WT and mutant lines.

Results: Introduction of NFKBIZ mutations into a DLBCL cell line confirmed that UTR deletions lead to increased mRNA and protein levels. NFKBIZ UTR deletions give DLBCL cells a selective growth advantage over WT both in vitro and in vivo. RNA-sequencing of mutant and WT lines revealed possible transcriptional targets of NFKBIZ, including NF-κB targets and genes commonly overexpressed in ABC DLBCL. Novel candidate NFKBIZ targets were also discovered through this analysis, including CD274, the gene encoding PD-L1. Mutant cell lines had significantly higher IC50 compared to WT for the drugs ibrutinib, idelalisib, and masitinib, but not bortezomib, suggesting that NKFBIZ UTR mutations make cell lines more resistant to specific NF-κB pathway-targeted drugs.

Conclusions: This work directly establishes a role for NFKBIZ amplifications and 3′ UTR mutations in driving ABC DLBCL through NF-κB signaling. We demonstrate that these mutations can lead to overexpression of NFKBIZ and provide a selective growth advantage to cells both in vitro and in vivo. In addition, we found that these mutant lines were more resistant to some targeted lymphoma drugs but not others.

Citation Format: Sarah E. Arthur, Nicole Thomas, Jeffrey Tang, Christopher K. Rushton, Miguel Alcaide, Adèle Telenius, Shannon Healy, Anja Mottok, David W. Scott, Christian Steidl, Ryan D. Morin. NFKBIZ 3′ UTR mutations confer selective growth advantage and affect drug response in diffuse large B-cell lymphoma [abstract]. In: Proceedings of the AACR Virtual Meeting: Advances in Malignant Lymphoma; 2020 Aug 17-19. Philadelphia (PA): AACR; Blood Cancer Discov 2020;1(3_Suppl):Abstract nr PO-32.

Abstract PO-04: Noncoding mutations in mantle cell lymphoma disrupt regulation of HNRNPH1 by alternative splicing

Non-Hodgkin lymphomas (NHL) are a collection of cancers with each malignancy having distinct clinical management and prognosis. Mantle cell lymphoma (MCL) is particularly genetically heterogeneous and is considered incurable. Through a combination of exome, genome, and targeted sequencing of MCL tumors, we identified recurrent mutations in HNRNPH1 (heterogeneous nuclear ribonucleoprotein H1). These mutations are largely intronic or silent and are associated with a putative cis regulatory region involving a single exon. In RNA-seq data from matched cases, we identified variable representation of two distinct HNRNPH1 isoforms. Based on the reading frame of the affected exons, canonical splicing is predicted to produce a functional protein, while alternative splicing introduces a premature termination codon leading to nonsense-mediated decay. We observed a significantly higher proportion of canonical transcripts in MCL tumors bearing HNRNPH1 mutations, leading us to conclude that mutations in HNRNPH1 significantly alter its splicing. Furthermore, increased canonical splicing results in higher HNRNPH1 protein abundance as determined by immunohistochemical analysis of an MCL tissue microarray. HNRNPH1 mutation status and splicing ratio are associated with shorter survival of MCL patients. We developed an in vitro reporter minigene and introduced three specific mutations corresponding to the patient-identified mutations in HNRNPH1. These mutations lead to an increase in canonical splicing and translation of the HNRNPH1 minigene-derived peptide. Additionally, when HNRNPH1 is overexpressed, we observe a decrease of this peptide, implicating HNRNPH1 in the regulation of its own splicing. This is supported by preliminary data indicating that HNRNPH1 binds its own RNA. Beyond its own splicing, HNRNPH1 is likely also involved in the splicing of additional splicing factors. Overexpression of HNRNPH1 affects the splicing of SRSF3 and HNRNPDL, and many other targets are likely to be identified. This work elucidates a functional role for recurrent noncoding HNRNPH1 mutations and implicates HNRNPH1 expression and splicing of its downstream targets in lymphomagenesis. We continue to explore trans regulatory targets of HNRNPH1 using in vitro and cell-based models. While splicing is a growing field of interest in lymphoma biology, the unique pattern and consequences of these largely silent mutations specifically implicate alternative splicing as an oncogenic mechanism in MCL.

Citation Format: Krysta M. Coyle, Quratulain Qureshi, Prasath Pararajalingam, Nicole Thomas, Timothy E. Audas, Ryan D. Morin. Noncoding mutations in mantle cell lymphoma disrupt regulation of HNRNPH1 by alternative splicing [abstract]. In: Proceedings of the AACR Virtual Meeting: Advances in Malignant Lymphoma; 2020 Aug 17-19. Philadelphia (PA): AACR; Blood Cancer Discov 2020;1(3_Suppl):Abstract nr PO-04.

Abstract IA42: Detecting and quantifying mutations associated with treatment resistance in aggressive lymphomas using ctDNA

A significant proportion of diffuse large B-cell lymphoma (DLBCL) patients treated with immunochemotherapy containing rituximab (R-CHOP) exhibit either primary or acquired treatment resistance. The advancement of therapeutics in the relapse setting has likely been encumbered by our limited understanding of the molecular features that underlie resistance to R-CHOP. Unfortunately, our knowledge of DLBCL genetics is mostly limited to analyses conducted on diagnostic tissue biopsies, which have not been exposed to the selective pressures imposed by therapy. Identifying genetic alterations that contribute to treatment resistance may reveal additional treatment options and lead to biomarkers allowing patients to be paired with appropriate treatments. Genetic subgroups are gaining popularity as a new strategy to implement precision medicine in DLBCL (1). The relevance of these and other biomarkers in the relapse setting remains unclear due to limited genetic exploration of relapsed and refractory DLBCL (rrDLBCL). Progress has been limited, in part, by the requirement of tissue biopsies collected after relapse. It is well established that quantitative genomic techniques such as digital PCR and targeted sequencing can be used to determine the proportion of tumor DNA in plasma from lymphoma patients (2). With a sufficiently broad panel, sequencing affords additional opportunities including the ability to identify subclonal structure and population dynamics over time. This presentation will discuss our recent analysis of a large collection of ctDNA primarily comprising DLBCL patients on various clinical trials (3). Targeted sequencing of these samples and comparison to exome data from a meta-cohort of previously characterized untreated DLBCL biopsies revealed six genes significantly enriched for mutations upon relapse. We found both TP53 and KMT2D were mutated in the majority of rrDLBCLs, and these mutations persisted in the dominant clone following relapse, suggesting a role in primary treatment resistance. By inferring subclonal dynamics, we observed recurrent patterns of clonal expansion and contraction following rituximab-based therapy, with MS4A1 mutations representing the only example of consistent clonal expansion. MS4A1 missense mutations within the transmembrane domains led to loss of CD20 expression in vitro, and patient tumors harboring these mutations lacked CD20 protein expression. Our analysis nominates TP53 and KMT2D mutation status as novel prognostic factors that may facilitate the identification of high-risk patients prior to therapy. Moreover, we have demonstrated the potential to identify tumors with loss of CD20 surface expression stemming from MS4A1 mutations. Implementation of noninvasive assays to detect such features of acquired treatment resistance may allow timely transition to more effective treatment regimens. In certain scenarios whole-exome sequencing (WES) or whole-genome sequencing (WGS) can be successfully applied to ctDNA, thereby allowing the identification of mutations, structural variation, and copy number changes. Low-pass sequencing of shotgun libraries can also be used to ascertain course estimates of ctDNA levels as well as the copy number landscape (4). Given the importance of copy number and structural alterations in the inference of genetic subgroups, these methods may allow the exploration of these groups and their stability over time. Through a series of illustrative examples, this presentation will explore the benefits of each of these techniques in the study of tumor evolution and acquired treatment resistance in DLBCL.

References: 1. Morin RD, Scott DW. DLBCL subclassification: Divide and conquer? Blood 2020;135:1722–4. 2. Rossi D et al. The development of liquid biopsy for research and clinical practice in lymphomas: Report of the 15-ICML workshop on ctDNA. Hematol Oncol 2020;38:34–7. 3. Rushton CK et al. Genetic and evolutionary patterns of treatment resistance in relapsed B-cell lymphoma. Blood Adv 2020;4:2886–98. 4. Adalsteinsson VA et al. Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors. Nat Commun 2017;8:1324.

Citation Format: Nicole Thomas, Laura K. Hilton, Neil Michaud, Kevin Bushell, Ryan Rys, Michael Jain, Lois Shepherd, Marco A. Marra, John Kuruvilla, Michael Crump, Koren Mann, Sarit Assouline, Christian Steidl, Mark S. Cragg, David W. Scott, Nathalie Johnson, Ryan D. Morin, Christopher K. Rushton, Sarah E. Arthur, Miguel Alcaide, Matthew Cheung, Aixiang Jiang, Krysta M. Coyle, Kirstie L. S. Cleary. Detecting and quantifying mutations associated with treatment resistance in aggressive lymphomas using ctDNA [abstract]. In: Proceedings of the AACR Virtual Meeting: Advances in Malignant Lymphoma; 2020 Aug 17-19. Philadelphia (PA): AACR; Blood Cancer Discov 2020;1(3_Suppl):Abstract nr IA42.

PRPS-ST: A protocol-agnostic self-training method for gene expression-based classification of blood cancers

Gene expression classifiers are gaining increasing popularity for stratifying tumors into subgroups with distinct biological features. A fundamental limitation shared by current classifiers is the requirement for comparable training and testing data sets. Here, we describe a self-training implementation of our probability ratio-based classification prediction score method (PRPS-ST), which facilitates the porting of existing classification models to other gene expression data sets. In comparison to gold standards, we demonstrate favorable performance of PRPS-ST in gene expression-based classification of DLBCL and B-ALL using a diverse variety of gene expression data types and pre-processing methods, including in classifications with a high degree of class imbalance. Tumors classified by our method were significantly enriched for prototypical genetic features of their respective subgroups. Interestingly, this included cases that were unclassifiable by established methods, implying the potential enhanced sensitivity of PRPS-ST.

Coding and noncoding drivers of mantle cell lymphoma identified through exome and genome sequencing

Mantle cell lymphoma (MCL) is an uncommon B-cell non-Hodgkin lymphoma (NHL) that is incurable with standard therapies. The genetic drivers of this cancer have not been firmly established, and the features that contribute to differences in clinical course remain limited. To extend our understanding of the biological pathways involved in this malignancy, we performed a large-scale genomic analysis of MCL using data from 51 exomes and 34 genomes alongside previously published exome cohorts. To confirm our findings, we resequenced the genes identified in the exome cohort in 191 MCL tumors, each having clinical follow-up data. We confirmed the prognostic association of TP53 and NOTCH1 mutations. Our sequencing revealed novel recurrent noncoding mutations surrounding a single exon of the HNRNPH1gene. In RNA-seq data from 103 of these cases, MCL tumors with these mutations had a distinct imbalance of HNRNPH1 isoforms. This altered splicing of HNRNPH1 was associated with inferior outcomes in MCL and showed a significant increase in protein expression by immunohistochemistry. We describe a functional role for these recurrent noncoding mutations in disrupting an autoregulatory feedback mechanism, thereby deregulating HNRNPH1 protein expression. Taken together, these data strongly imply a role for aberrant regulation of messenger RNA processing in MCL pathobiology.

Genetic and evolutionary patterns of treatment resistance in relapsed B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) patients are typically treated with immunochemotherapy containing rituximab (rituximab, cyclophosphamide, hydroxydaunorubicin-vincristine (Oncovin), and prednisone [R-CHOP]); however, prognosis is extremely poor if R-CHOP fails. To identify genetic mechanisms contributing to primary or acquired R-CHOP resistance, we performed target-panel sequencing of 135 relapsed/refractory DLBCLs (rrDLBCLs), primarily comprising circulating tumor DNA from patients on clinical trials. Comparison with a metacohort of 1670 diagnostic DLBCLs identified 6 genes significantly enriched for mutations upon relapse. TP53 and KMT2D were mutated in the majority of rrDLBCLs, and these mutations remained clonally persistent throughout treatment in paired diagnostic-relapse samples, suggesting a role in primary treatment resistance. Nonsense and missense mutations affecting MS4A1, which encodes CD20, are exceedingly rare in diagnostic samples but show recurrent patterns of clonal expansion following rituximab-based therapy. MS4A1 missense mutations within the transmembrane domains lead to loss of CD20 in vitro, and patient tumors harboring these mutations lacked CD20 protein expression. In a time series from a patient treated with multiple rounds of therapy, tumor heterogeneity and minor MS4A1-harboring subclones contributed to rapid disease recurrence, with MS4A1 mutations as founding events for these subclones. TP53 and KMT2D mutation status, in combination with other prognostic factors, may be used to identify high-risk patients prior to R-CHOP for posttreatment monitoring. Using liquid biopsies, we show the potential to identify tumors with loss of CD20 surface expression stemming from MS4A1 mutations. Implementation of noninvasive assays to detect such features of acquired treatment resistance may allow timely transition to more effective treatment regimens.

TBL1XR1 Mutations Drive Extranodal Lymphoma by Inducing a Pro-tumorigenic Memory Fate

The most aggressive B cell lymphomas frequently manifest extranodal distribution and carry somatic mutations in the poorly characterized gene TBL1XR1. Here, we show that TBL1XR1 mutations skew the humoral immune response toward generating abnormal immature memory B cells (MB), while impairing plasma cell differentiation. At the molecular level, TBL1XR1 mutants co-opt SMRT/HDAC3 repressor complexes toward binding the MB cell transcription factor (TF) BACH2 at the expense of the germinal center (GC) TF BCL6, leading to pre-memory transcriptional reprogramming and cell-fate bias. Upon antigen recall, TBL1XR1 mutant MB cells fail to differentiate into plasma cells and instead preferentially reenter new GC reactions, providing evidence for a cyclic reentry lymphomagenesis mechanism. Ultimately, TBL1XR1 alterations lead to a striking extranodal immunoblastic lymphoma phenotype that mimics the human disease. Both human and murine lymphomas feature expanded MB-like cell populations, consistent with a MB-cell origin and delineating an unforeseen pathway for malignant transformation of the immune system.

A Probabilistic Classification Tool for Genetic Subtypes of Diffuse Large B Cell Lymphoma with Therapeutic Implications

The development of precision medicine approaches for diffuse large B cell lymphoma (DLBCL) is confounded by its pronounced genetic, phenotypic, and clinical heterogeneity. Recent multiplatform genomic studies revealed the existence of genetic subtypes of DLBCL using clustering methodologies. Here, we describe an algorithm that determines the probability that a patient's lymphoma belongs to one of seven genetic subtypes based on its genetic features. This classification reveals genetic similarities between these DLBCL subtypes and various indolent and extranodal lymphoma types, suggesting a shared pathogenesis. These genetic subtypes also have distinct gene expression profiles, immune microenvironments, and outcomes following immunochemotherapy. Functional analysis of genetic subtype models highlights distinct vulnerabilities to targeted therapy, supporting the use of this classification in precision medicine trials.

Evaluation of protocols for rRNA depletion-based RNA sequencing of nanogram inputs of mammalian total RNA

Next generation RNA-sequencing (RNA-seq) is a flexible approach that can be applied to a range of applications including global quantification of transcript expression, the characterization of RNA structure such as splicing patterns and profiling of expressed mutations. Many RNA-seq protocols require up to microgram levels of total RNA input amounts to generate high quality data, and thus remain impractical for the limited starting material amounts typically obtained from rare cell populations, such as those from early developmental stages or from laser micro-dissected clinical samples. Here, we present an assessment of the contemporary ribosomal RNA depletion-based protocols, and identify those that are suitable for inputs as low as 1–10 ng of intact total RNA and 100–500 ng of partially degraded RNA from formalin-fixed paraffin-embedded tissues.

SUBSTRA: Supervised Bayesian Patient Stratification

Motivation

Patient stratification methods are key to the vision of precision medicine. Here, we consider transcriptional data to segment the patient population into subsets relevant to a given phenotype. Whereas most existing patient stratification methods focus either on predictive performance or interpretable features, we developed a method striking a balance between these two important goals.

Results

We introduce a Bayesian method called SUBSTRA that uses regularized biclustering to identify patient subtypes and interpretable subtype-specific transcript clusters. The method iteratively re-weights feature importance to optimize phenotype prediction performance by producing more phenotype-relevant patient subtypes. We investigate the performance of SUBSTRA in finding relevant features using simulated data and successfully benchmark it against state-of-the-art unsupervised stratification methods and supervised alternatives. Moreover, SUBSTRA achieves predictive performance competitive with the supervised benchmark methods and provides interpretable transcriptional features in diverse biological settings, such as drug response prediction, cancer diagnosis, or kidney transplant rejection.

Availability and implementation

The R code of SUBSTRA is available at https://github.com/sahandk/SUBSTRA.

Supplementary information

Supplementary data are available at Bioinformatics online.

Integrative genomic analysis of matched primary and metastatic pediatric osteosarcoma

Despite being the most common childhood bone tumor, the genomic characterization of osteosarcoma remains incomplete. In particular, very few osteosarcoma metastases have been sequenced to date, critical to better understand mechanisms of progression and evolution in this tumor. We performed an integrated whole genome and exome sequencing analysis of paired primary and metastatic pediatric osteosarcoma specimens to identify recurrent genomic alterations. Sequencing of 13 osteosarcoma patients including 13 primary, 10 metastatic, and 3 locally recurring tumors revealed a highly heterogeneous mutational landscape, including cases of hypermutation and microsatellite instability positivity, but with virtually no recurrent alterations except for mutations involving the tumor suppressor genes RB1 and TP53. At the germline level, we detected alterations in multiple cancer related genes in the majority of the cohort, including those potentially disrupting DNA damage response pathways. Metastases retained only a minimal number of short variants from their corresponding primary tumors, while copy number alterations showed higher conservation. One recurrently amplified gene, KDR, was highly expressed in advanced cases and associated with poor prognosis. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Abstract 3480: TMEM30A loss-of-function mutations drive lymphomagenesis and confer therapeutically exploitable vulnerability in B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common lymphoma subtype worldwide, accounting for 40% of all non-Hodgkin lymphomas. DLBCL presents as an aggressive disease requiring immediate treatment. Although significant improvement in outcome has been achieved, ~40% of patients still experience treatment failure. Here, we characterized the recurrent genetic alterations and transcriptomic signatures in diagnostic biopsies from a population registry-based cohort of 347 patients with de novo DLBCL uniformly treated with R-CHOP. This analysis revealed bi-allelic loss of function mutations of TMEM30A that were associated with favorable treatment outcome. TMEM30A is a chaperone protein, involved in maintaining the asymmetric distribution of phosphatidylethanolamine and phosphatidylserine, an integral component of the plasma membrane and “eat-me” signal recognized by macrophages. Using TMEM30A knockout systems by CRISPR genome editing techniques, we have functionally characterized this loss-of-function mutation in representative human and mouse DLBCL cell line models. We have discovered that TMEM30A loss is associated with increased B-cell signaling following antigen stimulation, including a two-fold increase in the diffusion rate of B-cell receptor (BCR) clustering, using high resolution Single Particle Tracking (SPT) technology. In addition, we have measured three-fold increase in chemotherapeutic drug accumulation in both knockout cell lines and randomly selected patient biopsies with TMEM30A biallelic loss. This observation was validated in a xenograft mouse model, which presented improved survival and limited tumor growth following vincristine treatment in mice injected with TMEM30A null DLBCL cell lines compared with native cell lines. This phenotype explains the improved prognosis observed in DLBCL patients following R-CHOP treatment. Furthermore, we have observed over two fold higher numbers of tumor-associated macrophages in B-cell lymphoma syngeneic mouse models with Tmem30a loss-of-function, prior to any form of treatment, suggesting the existence of “hot” and primed tumors. Our data highlight a multi-faceted role for TMEM30A and plasma membrane physiology in B-cell lymphomagenesis, and characterize intrinsic and extrinsic vulnerabilities of cancer cells that can be therapeutically exploited. Characterization of these mechanisms will address a missing link in the cancer field as related insights in lymphoma will outline therapeutic approaches that can be extended to cancer therapy in general.

Citation Format: Shannon Healy, Daisuke Ennishi, Ali Bashashati, Saeed Saberi, Christoffer Hother, Anja Mottok, Fong Chun Chan, Lauren Chong, Robert Kridel, Merrill Boyle, Barbara Meissner, Tomohiro Aoki, Katsuyoshi Takata, Bruce W. Woolcock, Elena Vigano, Libin Abraham, Michael Gold, Adele Telenius, Pedro Farinha, Graham Slack, Susana Ben-Neriah, Daniel Lai, Allen W. Zhang, Sohrab Salehi, Hennady P. Shulha, Derek S. Chiu, Sara Mostafavi, Alina S. Gerrie, Diego Villa, Laurie H. Sehn, Kerry J. J. Savage, Andrew J. J. Mungall, Andrew P. Weng, Marcel Bally, Ryan D. Morin, Gabriela V. Cohen Freue, Joseph M. Connors, Marco A. Marra, Sohrab P. Shah, Randy D. Gascoyne1, David W. Scott, Christian Steidl, Ulrich Steidl. TMEM30A loss-of-function mutations drive lymphomagenesis and confer therapeutically exploitable vulnerability in B-cell lymphoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 3480.

Genome-wide discovery of somatic coding and noncoding mutations in pediatric endemic and sporadic Burkitt lymphoma

Although generally curable with intensive chemotherapy in resource-rich settings, Burkitt lymphoma (BL) remains a deadly disease in older patients and in sub-Saharan Africa. Epstein-Barr virus (EBV) positivity is a feature in more than 90% of cases in malaria-endemic regions, and up to 30% elsewhere. However, the molecular features of BL have not been comprehensively evaluated when taking into account tumor EBV status or geographic origin. Through an integrative analysis of whole-genome and transcriptome data, we show a striking genome-wide increase in aberrant somatic hypermutation in EBV-positive tumors, supporting a link between EBV and activation-induced cytidine deaminase (AICDA) activity. In addition to identifying novel candidate BL genes such as SIN3A, USP7, and CHD8, we demonstrate that EBV-positive tumors had significantly fewer driver mutations, especially among genes with roles in apoptosis. We also found immunoglobulin variable region genes that were disproportionally used to encode clonal B-cell receptors (BCRs) in the tumors. These include IGHV4-34, known to produce autoreactive antibodies, and IGKV3-20, a feature described in other B-cell malignancies but not yet in BL. Our results suggest that tumor EBV status defines a specific BL phenotype irrespective of geographic origin, with particular molecular properties and distinct pathogenic mechanisms. The novel mutation patterns identified here imply rational use of DNA-damaging chemotherapy in some patients with BL and targeted agents such as the CDK4/6 inhibitor palbociclib in others, whereas the importance of BCR signaling in BL strengthens the potential benefit of inhibitors for PI3K, Syk, and Src family kinases among these patients.

Disruption of the Gut Microbiota With Antibiotics Exacerbates Acute Vascular Rejection

BACKGROUND

The gut microbiota influences many immunological processes but how its disruption affects transplant rejection is poorly understood.

METHODS

Interposition grafting of aortic segments was used to examine vascular rejection. The gut microbiota was disrupted in graft recipients using an antibiotic cocktail (ampicillin, vancomycin, metronidazole, neomycin sulfate) in their drinking water.

RESULTS

Treatment of mice with antibiotics severely reduced total bacterial content in the intestine and disrupted the bacterial composition. Short-term treatment of mice for only the first 3 weeks of life resulted in the population of the intestine in mature mice with bacterial communities that were mildly different from untreated mice, containing slightly more Clostridia and less Bacteroides. Antibiotic disruption of the gut microbiota of graft recipients, either for their entire life or only during the first 3 weeks of life, resulted in increased medial injury of allograft arteries that is reflective of acute vascular rejection but did not affect intimal thickening reflective of transplant arteriosclerosis. Exacerbated vascular rejection resulting from disruption of the gut microbiota was related to increased infiltration of allograft arteries by neutrophils.

CONCLUSIONS

Disruption of the gut microbiota early in life results in exacerbation of immune responses that cause acute vascular rejection.

Molecular and Genetic Characterization of MHC Deficiency Identifies EZH2 as Therapeutic Target for Enhancing Immune Recognition

We performed a genomic, transcriptomic, and immunophenotypic study of 347 patients with diffuse large B-cell lymphoma (DLBCL) to uncover the molecular basis underlying acquired deficiency of MHC expression. Low MHC-II expression defines tumors originating from the centroblast-rich dark zone of the germinal center (GC) that was associated with inferior prognosis. MHC-II-deficient tumors were characterized by somatically acquired gene mutations reducing MHC-II expression and a lower amount of tumor-infiltrating lymphocytes. In particular, we demonstrated a strong enrichment of EZH2 mutations in both MHC-I- and MHC-II-negative primary lymphomas, and observed reduced MHC expression and T-cell infiltrates in murine lymphoma models expressing mutant Ezh2 ^Y641. Of clinical relevance, EZH2 inhibitors significantly restored MHC expression in EZH2-mutated human DLBCL cell lines. Hence, our findings suggest a tumor progression model of acquired immune escape in GC-derived lymphomas and pave the way for development of complementary therapeutic approaches combining immunotherapy with epigenetic reprogramming. SIGNIFICANCE: We demonstrate how MHC-deficient lymphoid tumors evolve in a cell-of-origin-specific context. Specifically, EZH2 mutations were identified as a genetic mechanism underlying acquired MHC deficiency. The paradigmatic restoration of MHC expression by EZH2 inhibitors provides the rationale for synergistic therapies combining immunotherapies with epigenetic reprogramming to enhance tumor recognition and elimination.See related commentary by Velcheti et al., p. 472.This article is highlighted in the In This Issue feature, p. 453.

Double-Hit Gene Expression Signature Defines a Distinct Subgroup of Germinal Center B-Cell-Like Diffuse Large B-Cell Lymphoma

PURPOSE

High-grade B-cell lymphoma with MYC and BCL2 and/or BCL6 rearrangements (HGBL-DH/TH) has a poor outcome after standard chemoimmunotherapy. We sought to understand the biologic underpinnings of HGBL-DH/TH with BCL2 rearrangements (HGBL-DH/TH- BCL2) and diffuse large B-cell lymphoma (DLBCL) morphology through examination of gene expression.

PATIENTS AND METHODS

We analyzed RNA sequencing data from 157 de novo germinal center B-cell-like (GCB)-DLBCLs, including 25 with HGBL-DH/TH- BCL2, to define a gene expression signature that distinguishes HGBL-DH/TH- BCL2 from other GCB-DLBCLs. To assess the genetic, molecular, and phenotypic features associated with this signature, we analyzed targeted resequencing, whole-exome sequencing, RNA sequencing, and immunohistochemistry data.

RESULTS

We developed a 104-gene double-hit signature (DHITsig) that assigned 27% of GCB-DLBCLs to the DHITsig-positive group, with only one half harboring MYC and BCL2 rearrangements (HGBL-DH/TH- BCL2). DHITsig-positive patients had inferior outcomes after rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisone immunochemotherapy compared with DHITsig-negative patients (5-year time to progression rate, 57% and 81%, respectively; P < .001), irrespective of HGBL-DH/TH- BCL2 status. The prognostic value of DHITsig was confirmed in an independent validation cohort. DHITsig-positive tumors are biologically characterized by a putative non-light zone germinal center cell of origin and a distinct mutational landscape that comprises genes associated with chromatin modification. A new NanoString assay (DLBCL90) recapitulated the prognostic significance and RNA sequencing assignments. Validating the association with HGBL-DH/TH- BCL2, 11 of 25 DHITsig-positive-transformed follicular lymphomas were classified as HGBL-DH/TH- BCL2 compared with zero of 50 in the DHITsig-negative group. Furthermore, the DHITsig was shared with the majority of B-cell lymphomas with high-grade morphology tested.

CONCLUSION

We have defined a clinically and biologically distinct subgroup of tumors within GCB-DLBCL characterized by a gene expression signature of HGBL-DH/TH- BCL2. This knowledge has been translated into an assay applicable to routinely available biopsy samples, which enables exploration of its utility to guide patient management.

Kronos: a workflow assembler for genome analytics and informatics

Background:

The field of next-generation sequencing informatics has matured to a point where algorithmic advances in sequence alignment and individual feature detection methods have stabilized. Practical and robust implementation of complex analytical workflows (where such tools are structured into “best practices” for automated analysis of next-generation sequencing datasets) still requires significant programming investment and expertise.

Results:

We present Kronos, a software platform for facilitating the development and execution of modular, auditable, and distributable bioinformatics workflows. Kronos obviates the need for explicit coding of workflows by compiling a text configuration file into executable Python applications. Making analysis modules would still require programming. The framework of each workflow includes a run manager to execute the encoded workflows locally (or on a cluster or cloud), parallelize tasks, and log all runtime events. The resulting workflows are highly modular and configurable by construction, facilitating flexible and extensible meta-applications that can be modified easily through configuration file editing. The workflows are fully encoded for ease of distribution and can be instantiated on external systems, a step toward reproducible research and comparative analyses. We introduce a framework for building Kronos components that function as shareable, modular nodes in Kronos workflows.

Conclusions:

The Kronos platform provides a standard framework for developers to implement custom tools, reuse existing tools, and contribute to the community at large. Kronos is shipped with both Docker and Amazon Web Services Machine Images. It is free, open source, and available through the Python Package Index and at https://github.com/jtaghiyar/kronos.

Enhancing knowledge discovery from cancer genomics data with Galaxy

The field of cancer genomics has demonstrated the power of massively parallel sequencing techniques to inform on the genes and specific alterations that drive tumor onset and progression. Although large comprehensive sequence data sets continue to be made increasingly available, data analysis remains an ongoing challenge, particularly for laboratories lacking dedicated resources and bioinformatics expertise. To address this, we have produced a collection of Galaxy tools that represent many popular algorithms for detecting somatic genetic alterations from cancer genome and exome data. We developed new methods for parallelization of these tools within Galaxy to accelerate runtime and have demonstrated their usability and summarized their runtimes on multiple cloud service providers. Some tools represent extensions or refinement of existing toolkits to yield visualizations suited to cohort-wide cancer genomic analysis. For example, we present Oncocircos and Oncoprintplus, which generate data-rich summaries of exome-derived somatic mutation. Workflows that integrate these to achieve data integration and visualizations are demonstrated on a cohort of 96 diffuse large B-cell lymphomas and enabled the discovery of multiple candidate lymphoma-related genes. Our toolkit is available from our GitHub repository as Galaxy tool and dependency definitions and has been deployed using virtualization on multiple platforms including Docker.

Genetic heterogeneity in primary and relapsed mantle cell lymphomas: Impact of recurrent CARD11 mutations

The genetic mechanisms underlying disease progression, relapse and therapy resistance in mantle cell lymphoma (MCL) remain largely unknown. Whole-exome sequencing was performed in 27 MCL samples from 13 patients, representing the largest analyzed series of consecutive biopsies obtained at diagnosis and/or relapse for this type of lymphoma. Eighteen genes were found to be recurrently mutated in these samples, including known (ATM, MEF2B and MLL2) and novel mutation targets (S1PR1 and CARD11). CARD11, a scaffold protein required for B-cell receptor (BCR)-induced NF-κB activation, was subsequently screened in an additional 173 MCL samples and mutations were observed in 5.5% of cases. Based on in vitro cell line-based experiments, overexpression of CARD11 mutants were demonstrated to confer resistance to the BCR-inhibitor ibrutinib and NF-κB-inhibitor lenalidomide. Genetic alterations acquired in the relapse samples were found to be largely non-recurrent, in line with the branched evolutionary pattern of clonal evolution observed in most cases. In summary, this study highlights the genetic heterogeneity in MCL, in particular at relapse, and provides for the first time genetic evidence of BCR/NF-κB activation in a subset of MCL.

The genomic landscape of two Burkitt lymphoma cases and derived cell lines: comparison between primary and relapse samples

Relapse occurs in 10-40% of Burkitt lymphoma (BL) patients that have completed intensive chemotherapy regimens and is typically fatal. While treatment-naive BL has been characterized, the genomic landscape of BL at the time of relapse (rBL) has never been reported. Here, we present a genomic characterization of two rBL patients. The diagnostic samples had mutations common in BL, including MYC and CCND3. Additional mutations were detected at relapse, affecting important pathways such as NFκB (IKBKB) and MEK/ERK (NRAS) signaling, glutamine metabolism (SIRT4), and RNA processing (ZFP36L2). Genes implicated in drug resistance were also mutated at relapse (TP53, BAX, ALDH3A1, APAF1, FANCI). This concurrent genomic profiling of samples obtained at diagnosis and relapse has revealed mutations not previously reported in this disease. The patient-derived cell lines will be made available and, along with their detailed genetics, will be a valuable resource to examine the role of specific mutations in therapeutic resistance.

Targeted error-suppressed quantification of circulating tumor DNA using semi-degenerate barcoded adapters and biotinylated baits

Ultrasensitive methods for rare allele detection are critical to leverage the full potential offered by liquid biopsies. Here, we describe a novel molecular barcoding method for the precise detection and quantification of circulating tumor DNA (ctDNA). The major benefits of our design include straightforward and cost-effective production of barcoded adapters to tag individual DNA molecules before PCR and sequencing, and better control over cross-contamination between experiments. We validated our approach in a cohort of 24 patients with a broad spectrum of cancer diagnoses by targeting and quantifying single-nucleotide variants (SNVs), indels and genomic rearrangements in plasma samples. By using personalized panels targeting a priori known mutations, we demonstrate comprehensive error-suppression capabilities for SNVs and detection thresholds for ctDNA below 0.1%. We also show that our semi-degenerate barcoded adapters hold promise for noninvasive genotyping in the absence of tumor biopsies and monitoring of minimal residual disease in longitudinal plasma samples. The benefits demonstrated here include broad applicability, flexibility, affordability and reproducibility in the research and clinical settings.

Investigating the Genetic Causes of Sudden Unexpected Death in Children Through Targeted Next-Generation Sequencing Analysis

Background—

Inherited arrhythmia syndromes are responsible for a significant portion of autopsy-negative sudden unexpected death (SUD) cases, but molecular autopsy used to identify potentially causal variants is not routinely included in SUD investigations. We collaborated with a medical examiner's office to assist in finding a diagnosis for their autopsy-negative child SUD cases.

Methods and Results—

191 child SUD cases (<5 years of age) were selected for analyses. Our next generation sequencing panel incorporated 38 inherited arrhythmia syndrome candidate genes and another 33 genes not previously investigated for variants that may underlie SUDY pathophysiology. Overall, we identified 11 potentially causal disease-associated variants in 12 cases, for an overall yield of 6.3%. We also identified 31 variants of uncertain significance in 36 cases and 16 novel variants predicted to be pathogenic in silico in 15 cases. The disease-associated variants were reported to the medical examiner to notify surviving relatives and recommend clinical assessment.

Conclusions—

We have identified variants that may assist in the diagnosis of at least 6.3% of autopsy-negative child SUD cases and reduce risk of future SUD in surviving relatives. We recommend a cautious approach to variant interpretation. We also suggest inclusion of cardiomyopathy genes as well as other candidate SUD genes in molecular autopsy analyses.