Serres, MaKayla R.

Introduction: Molecular classifications have been developed to characterize the heterogeneity of diffuse large B cell lymphoma (DLBCL), with activation of BCL6 and Notch signaling pathways acting as driving pathogenic features in C1/BN2 subtype tumors. Among C1/BN2 defining variants are truncating/inactivating SPEN mutations. SPEN regulates the Notch signaling pathway by initiating formation of a repressive complex that facilitates transcriptional repression of Notch target genes. Our group has identified SPEN mutations enriched in DLBCL cases failing to achieve event-free status at 24 months (EFS24). This contrasted with generally favorable outcomes in C1/BN2 tumors, suggesting a potential high-risk molecular feature with unexplored significance. Presenting an opportunity to expand underlying mechanisms and risk stratification in C1/BN2 tumors, we aimed to use a functional genomic approach, combined with clinical data, to characterize the alternative biology stemming from SPEN truncating mutations. Methods: FFPE tumor biopsies from newly diagnosed DLBCL cases (N=444) enrolled in the Mayo Clinic/Iowa Lymphoma Molecular Epidemiology Resource (MER) program were analyzed by WES (n=404) and RNA-seq (n=321). Biologic and clinical variables were annotated for all cases. CRISPR/Cas9 introduced truncating mutations to the endogenous SPEN locus in OCI-LY3 & SU-DHL-2 ABC-type human cell lines to establish models of SPEN depletion. Functional in vitro analysis utilized standard methods for RNA-seq, cell proliferation, western blot, and flow cytometry. Results: WES identified 29/404 ndDLBCL (7.2%) harboring SPEN mutations; n=17 missense and n=12 nonsense or frameshift insertions/deletions (truncating variants - [trunc]), the later considered as the focal point of our study. Against all cases, SPEN trunc cases displayed enrichment for high-risk factors including non-GCB COO (p=0.06, OR=3.7, 95% CI [0.9, 22.2]) and primary refractory disease (PTR; p=0.02, OR=5.0, 95% CI [1.05, 19.9]), and demonstrated inferior OS (p=0.06, HR=2.2), where a majority of clinical progression events occurred within 12 months of diagnosis. Leveraging RNA-seq, we observed 11 of 11 SPEN trunc cases exhibiting a double-hit negative gene expression signature, 9/11 cases with inflammatory or depleted microenvironment (LME), and 8/11 with a signature associating with intermediate/high risk of EFS24 failure. Applying the LymphGen algorithm, 5/12 SPEN trunc cases received the BN2 classification, observing co-occurring mutations with CD70 (p<0.01, OR=9.68), BRCA2 (p=0.03, OR=5.84), TET2 (p=0.02, OR=4.91), IRF4 (p=0.02, OR=7.27), and PTPRD (p<0.01, OR=8.31) genes, and were mutually exclusive with ARID1A, BCL2, EZH2, FAS, HLA-B, IRF8, and TMEM30A mutations ( Figure 1A). Isolating SPEN trunc against BN2 cases without SPEN variants (n=14), despite similar clustering of non-GCB COO (p=0.7), PTR was exclusively observed in SPEN trunc cases, accompanied again by inferior OS (p=0.07, HR=4.2; Figure 1B). In the OCI-LY3 SPEN depletion model, GSEA analysis from RNA-seq revealed activation of E2F (NES=2.7; p=9.4E-20) and MYC target genes (NES=2.3; p=4.1E-6), while both OCI-LY3 and SU-DHL-2 models featured depletion of antigen presentation genes (NES= -1.8, p=5.7E-4; NES= -2.2, p=3.8E-7, respectively). A decrease of MHC class II proteins (HLA-DR, HLA-DP) was confirmed by flow cytometry in both models (p<0.01). Hypothesizing loss of SPEN enhances reliance on the BCL6-HDAC3 axis featured in BN2 tumors, we treated both SPEN depletion models with a class-I HDAC inhibitor highly selective for HDAC3, BRD-3308, and indeed observed a 1.4 and 2.4-fold increase in sensitivity compared to WT cells (p=0.07, p<0.01, respectively). Significance: We find SPEN trunc mutations in poor acting DLBCL cases with clinical and biologic characteristics that diverge from C1/BN2 biology. The enhanced rate of clinical progression highlights urgency to better characterize these tumors. Our data suggests loss of SPEN regulatory mechanisms may act cooperatively with C1/BN2 biology to facilitate tumor immune evasion and drive progression through Notch-mediated pathways. Importantly, in vitro evidence shows SPEN deficient lymphoma cells are susceptible to targeted therapeutics against dependent pathways, providing rationale for advancing individualized care strategies for affected patients.

Abstract Non‐follicular low‐grade B‐cell lymphomas (LGBCL) are biologically diverse entities that share clinical and histologic features that make definitive pathologic categorization challenging. While most patients with LGBCL have an indolent course, some experience aggressive disease, highlighting additional heterogeneity across these subtypes. To investigate the potential for shared biology across subtypes, we performed RNA sequencing and applied machine learning approaches that identified five clusters of patients that grouped independently of subtype. One cluster was characterized by inferior outcome, upregulation of cell cycle genes, and increased tumor immune cell content. Integration of whole exome sequencing identified novel LGBCL mutations and enrichment of TNFAIP3 and BCL2 alterations in the poor survival cluster. Building on this, we further refined a transcriptomic signature associated with early clinical failure in two independent cohorts. Taken together, this study identifies unique clusters of LGBCL defined by novel gene expression signatures and immune profiles associated with outcome across diagnostic subtypes.

The clinical and largely unpredictable heterogeneity of phenotypes in patients with mitochondrial disorders demonstrates the ongoing challenges in the understanding of this semi-autonomous organelle in biology and disease. Previously, we used the gene-breaking transposon to create 1200 transgenic zebrafish strains tagging protein-coding genes (Ichino et al., 2020), including the lrpprc locus. Here, we present and characterize a new genetic revertible animal model that recapitulates components of Leigh Syndrome French Canadian Type (LSFC), a mitochondrial disorder that includes diagnostic liver dysfunction. LSFC is caused by allelic variations in the LRPPRC gene, involved in mitochondrial mRNA polyadenylation and translation. lrpprc zebrafish homozygous mutants displayed biochemical and mitochondrial phenotypes similar to clinical manifestations observed in patients, including dysfunction in lipid homeostasis. We were able to rescue these phenotypes in the disease model using a liver-specific genetic model therapy, functionally demonstrating a previously under-recognized critical role for the liver in the pathophysiology of this disease.

Abstract Introduction: Low-grade B-cell lymphomas (LGBCL), aside from follicular lymphoma and chronic lymphocytic leukemia/small lymphocytic lymphoma, account for approximately 10% of B-cell non-Hodgkin lymphomas and consist of several subtypes. While a majority of LGBCL cases have an overall favorable prognosis, we have previously shown that cases who have an event (relapse or progression, transformation, or re-treatment) within 24 months of diagnosis (EFS24) have an inferior overall survival (OS) compared to those achieving EFS24 (Tracy et al., AJH 2019;94:658-66). However, the underlying biological characteristics associated with early failure and aggressive disease across LGBCL subtypes are unknown. In this study, we used matched transcriptomic, genomic, and immune profiling data from LGBCL cases, the largest cohort to date, and asked whether there were unique biological phenotypes across different LGBCL subtypes and whether we could identify signatures associated with aggressive LGBCL. Validation of the prognostic utility of this signature was performed on a previously published, independent cohort of 63 pre-treatment LGBCL cases. Methods: Tumors from 64 newly diagnosed LGBCL patients from the Molecular Epidemiology Resource of the University of Iowa/Mayo Clinic Lymphoma Specialized Program of Research Excellence were included in this study (SMZL (n = 48), NMZL (n = 6), LPL (n = 5), B-NOS (n = 3), EMZL (n = 2)). RNA sequencing (RNAseq) data from 61 LGBCL tumors and 5 benign CD19+CD27+ memory B samples was subjected to NMF clustering to define groups. Differential expression and pathway analysis were used to identify biological characteristics of each cluster. CIBERSORT was used to identify immune cells in the tumor microenvironment. Whole exome sequencing (WES) was performed on 61 tumor-normal pairs. Singscore was used to assign a single score per patient representing gene expression of the survival-associated transcriptomic signature identified in this study. Results: NMF analysis of RNAseq data identified 5 clusters of patients, denoted LGBCL1-5 (Fig 1A). Patients from the same diagnostic subtype did not exclusively cluster together, with all LGBCL clusters comprised of patients from multiple subtypes (Fig 1B). Exploring the association between patient cluster and outcome, we observed significantly inferior event-free survival (EFS) (HR 2.24; 95% CI 1.01-4.98) and overall survival (OS) (HR 5.59; 95% CI 2.00-15.63) in LGBCL5 patients compared to LGBCL1-4 (Fig 1C). In addition, 80% of the transformation cases in our cohort were classified as LGBCL5 (Fig 1D). Differential expression and pathway analysis showed distinct processes significantly upregulated in each cluster (FDR < 0.05), with LGBCL5 demonstrating enrichment of cell cycle and mitosis pathways. CIBERSORT identified increased immune cell content in LGBCL3 and LGBCL5 compared to other clusters, with high frequencies of mast cells in both (p = 0.0002), increased CD8 T cells in LGBCL3 (p < 0.0001), and increased T follicular helper cells in LGBCL5 (p = 0.004). WES identified previously reported alterations in NOTCH, NFkB, and chromatin remodeling pathways and novel variants in LGBCL, including mutations in HNRNPK, CLTC, HLA-A, HLA-B and HLA-C. Assessment of alterations by cluster showed significant enrichment of TNFAIP3 (OR 5.54; 95% CI 1.20-28.14) and BCL2 alterations (OR 5.49; 95% CI 1.07-32.02) in LGBCL5 cluster. Finally, we identified a cell cycle-related transcriptomic signature of 108 genes upregulated in LGBCL5 and EFS24 failure cases. Cases with high expression of this signature showed significantly inferior EFS (HR 14.25; 95% CI 4.90-41.38) and OS (HR 7.82; 95% CI 2.40-25.44) compared to cases with low expression in our discovery cohort. This observation was reproduced in an independent validation cohort, where patients with high expression of this signature demonstrated significantly inferior EFS (HR 5.70; 95% CI 1.49-21.79) and OS (HR 10.07; 95% CI 2.00-50.61). Conclusions: In this study, we are the first to define mechanisms of pathogenesis in LGBCL with shared transcriptomic, genomic, and immune profiles present across LGBCL subtypes. We then further defined a gene expression signature associated with inferior patient outcome, with application of this signature to an independent validation cohort demonstrating proof of concept and utility of this signature as a prognostic marker in LGBCL patients. Figure 1 Figure 1. Disclosures Maurer: Morphosys: Membership on an entity's Board of Directors or advisory committees, Research Funding; Genentech: Research Funding; Kite Pharma: Membership on an entity's Board of Directors or advisory committees; Pfizer: Membership on an entity's Board of Directors or advisory committees; BMS: Research Funding; Nanostring: Research Funding. Paludo: Karyopharm: Research Funding. Habermann: Tess Therapeutics: Other: Data Monitoring Committee; Morphosys: Other: Scientific Advisory Board; Incyte: Other: Scientific Advisory Board; Seagen: Other: Data Monitoring Committee; Loxo Oncology: Other: Scientific Advisory Board; Eli Lilly & Co.,: Other: Scientific Advisor. Link: MEI: Consultancy; Genentech/Roche: Consultancy, Research Funding; Novartis, Jannsen: Research Funding. Rimsza: NanoString Technologies: Other: Fee-for-service contract. Ansell: Bristol Myers Squibb, ADC Therapeutics, Seattle Genetics, Regeneron, Affimed, AI Therapeutics, Pfizer, Trillium and Takeda: Research Funding. Cerhan: Genentech: Research Funding; Regeneron Genetics Center: Other: Research Collaboration; Celgene/BMS: Other: Connect Lymphoma Scientific Steering Committee, Research Funding; NanoString: Research Funding. Novak: Celgene/BMS: Research Funding.

One key bottleneck in understanding the human genome is the relative under-characterization of 90% of protein coding regions. We report a collection of 1200 transgenic zebrafish strains made with the gene-break transposon (GBT) protein trap to simultaneously report and reversibly knockdown the tagged genes. Protein trap-associated mRFP expression shows previously undocumented expression of 35% and 90% of cloned genes at 2 and 4 days post-fertilization, respectively. Further, investigated alleles regularly show 99% gene-specific mRNA knockdown. Homozygous GBT animals in ryr1b , fras1 , tnnt2a , edar and hmcn1 phenocopied established mutants. 204 cloned lines trapped diverse proteins, including 64 orthologs of human disease-associated genes with 40 as potential new disease models. Severely reduced skeletal muscle Ca 2+ transients in GBT ryr1b homozygous animals validated the ability to explore molecular mechanisms of genetic diseases. This GBT system facilitates novel functional genome annotation towards understanding cellular and molecular underpinnings of vertebrate biology and human disease.

When vertebrates face acute stressors, their bodies rapidly undergo a repertoire of physiological and behavioral adaptations, which is termed the stress response. Rapid changes in heart rate and blood glucose levels occur via the interaction of glucocorticoids and their cognate receptors following hypothalamic‐pituitary‐adrenal axis activation. These physiological changes are observed within minutes of encountering a stressor and the rapid time domain rules out genomic responses that require gene expression changes. Although behavioral changes corresponding to physiological changes are commonly observed, it is not clearly understood to what extent hypothalamic‐pituitary‐adrenal axis activation dictates adaptive behavior. We hypothesized that rapid locomotor response to acute stressors in zebrafish requires hypothalamic‐pituitary‐interrenal (HPI) axis activation. In teleost fish, interrenal cells are functionally homologous to the adrenocortical layer. We derived eight frameshift mutants in genes involved in HPI axis function: two mutants in exon 2 ofmc2r(adrenocorticotropic hormone receptor), five in exon 2 or 5 ofnr3c1(glucocorticoid receptor [GR]) and two in exon 2 ofnr3c2(mineralocorticoid receptor [MR]). Exposing larval zebrafish to mild environmental stressors, acute changes in salinity or light illumination, results in a rapid locomotor response. We show that this locomotor response requires a functioning HPI axis via the action ofmc2rand the canonical GR encoded bynr3c1gene, but not MR (nr3c2). Our rapid behavioral assay paradigm based on HPI axis biology can be used to screen for genetic and environmental modifiers of the hypothalamic‐pituitary‐adrenal axis and to investigate the effects of corticosteroids and their cognate receptor interactions on behavior.