Search

Supplementary Figure 6 Boxplots showing the expression of indicated metabolites in C1 (red), C2 (blue) and C3d (green). Differences between groups were calculated with the Mann-Whitney U test. ns, non-significant.

Supplemental Fig. 2. Gene expression of non-tumor cells by cluster.

Purpose:IgM monoclonal gammopathy of undetermined significance (MGUS) and Waldenström macroglobulinemia (WM) represent a disease spectrum with highly varied therapeutic management, ranging from observation to chemoimmunotherapy. The current classification relies solely on clinical features and does not explain the heterogeneity that exists within each of these conditions. Further investigation is warranted to shed light on the biology that may account for the clinical differences.Experimental Design:We used bone marrow (BM) clonal CD19+ and/or CD138+ sorted cells, matched BM supernatant, and peripheral blood serum from 32 patients (7 MGUS, 25 WM) to perform the first multi-omics approach including whole-exome sequencing, RNA sequencing, proteomics, metabolomics, and mass cytometry.Results:We identified three clusters with distinct pathway activation, immune content, metabolomic, and clinical features. Cluster 1 included only patients with WM and was characterized by transcriptional silencing of genes involved in cell cycle and immune response, enrichment of mitochondrial metabolism, infiltration of senescent T effector memory cells, and aggressive clinical behavior. Genetic/structural alterations of TNFAIP3 were distinct events of this cluster. Cluster 2 comprised both MGUS and WM patients with upregulation of inflammatory response, senescence and glycolysis signatures, increased activated T follicular helper and T regulatory cells, and indolent clinical behavior. Cluster 3 also included both MGUS and WM patients and exhibited intermediate features, including proliferative and inflammatory signaling, as well as glycolysis and mitochondrial metabolism.Conclusions:We have identified three distinct molecular clusters, suggesting a potential biologic classification that may have therapeutic implications.

Supplementary Figure 4. Schematic representation of the mechanism associated with the three clusters.

Supplementary Figure 8 Boxplots showing the expression of CD69 (A), BTLA (B), and CD25+CTLA-4+ (C) in CD4+ T cells.

Supplementary Material

Supplementary Figure 5 A. Heat map of Spearman's correlation coefficients among paired proteins in the bone marrow (BM) supernatant and peripheral blood (PB) serum of 16 patients with WM (n=11) and IgM MGUS (n=5). Positive correlations are shown in red, whereas negative correlations are indicated in blue. B. Heat maps showing expression of (i) 12 proteins significantly different between C1 vs C2 plus C3, and (i) 74 proteins significantly different between C2 vs C1 plus C3 in the PB serum. The left bar shows correlation coefficients between BM and PB proteins. Gray squares indicate non significance.

Supplementary Figure 3 Boxplot showing that 6q deletion resulted in decreased expression of TNFAIP3.

Supplementary Figure 9 SPADE maps of CD4+ CTLA-4+ T cells from concatenated files (C1, n=7; C2, n=4). The dots in red and blue indicate high and low expression, respectively. The dot size is representative of cell events, with the big size containing more cell events.

Supplementary Figure 7 A. Heat map of Spearman's correlation coefficients among paired metabolites in the bone marrow (BM) supernatant and peripheral blood (PB) serum of 12 patients with WM. Positive correlations are shown in red, whereas negative correlations are indicated in blue. B. Heat map showing different expression of 44 metabolites in the BM supernatant among the three molecular clusters. The left bar shows correlation coefficients between BM and PB metabolites. Gray squares indicate non significance.

Purpose:IgM monoclonal gammopathy of undetermined significance (MGUS) and Waldenström macroglobulinemia (WM) represent a disease spectrum with highly varied therapeutic management, ranging from observation to chemoimmunotherapy. The current classification relies solely on clinical features and does not explain the heterogeneity that exists within each of these conditions. Further investigation is warranted to shed light on the biology that may account for the clinical differences.Experimental Design:We used bone marrow (BM) clonal CD19+ and/or CD138+ sorted cells, matched BM supernatant, and peripheral blood serum from 32 patients (7 MGUS, 25 WM) to perform the first multi-omics approach including whole-exome sequencing, RNA sequencing, proteomics, metabolomics, and mass cytometry.Results:We identified three clusters with distinct pathway activation, immune content, metabolomic, and clinical features. Cluster 1 included only patients with WM and was characterized by transcriptional silencing of genes involved in cell cycle and immune response, enrichment of mitochondrial metabolism, infiltration of senescent T effector memory cells, and aggressive clinical behavior. Genetic/structural alterations of TNFAIP3 were distinct events of this cluster. Cluster 2 comprised both MGUS and WM patients with upregulation of inflammatory response, senescence and glycolysis signatures, increased activated T follicular helper and T regulatory cells, and indolent clinical behavior. Cluster 3 also included both MGUS and WM patients and exhibited intermediate features, including proliferative and inflammatory signaling, as well as glycolysis and mitochondrial metabolism.Conclusions:We have identified three distinct molecular clusters, suggesting a potential biologic classification that may have therapeutic implications.

SLCO1B1 (solute carrier organic anion transporter family member 1B1) is an important transmembrane hepatic uptake transporter. Genetic variants in the SLCO1B1 gene have been associated with altered protein folding, resulting in protein degradation and decreased transporter activity. Next-generation sequencing (NGS) of pharmacogenes is being applied increasingly to associate variation in drug response with genetic sequence variants. However, it is difficult to link variants of unknown significance with functional phenotypes using “one-at-a-time” functional systems. Deep mutational scanning (DMS) using a “landing pad cell–based system” is a high-throughput technique designed to analyze hundreds of gene open reading frame (ORF) missense variants in a parallel and scalable fashion. We have applied DMS to analyze 137 missense variants in the SLCO1B1 ORF obtained from the Exome Aggregation Consortium project. ORFs containing these variants were fused to green fluorescent protein and were integrated into “landing pad” cells. Florescence-activated cell sorting was performed to separate the cells into four groups based on fluorescence readout indicating protein expression at the single cell level. NGS was then performed and SLCO1B1 variant frequencies were used to determine protein abundance. We found that six variants not previously characterized functionally displayed less than 25% and another 12 displayed approximately 50% of wild-type protein expression. These results were then functionally validated by transporter studies. Severely damaging variants identified by DMS may have clinical relevance for SLCO1B1-dependent drug transport, but we need to exercise caution since the relatively small number of severely damaging variants identified raise questions with regard to the application of DMS to intrinsic membrane proteins such as organic anion transporter protein 1B1. Significance Statement The functional implications of a large numbers of open reading frame (ORF) “variants of unknown significance” (VUS) in transporter genes have not been characterized. This study applied deep mutational scanning to determine the functional effects of VUS that have been observed in the ORF of SLCO1B1(solute carrier organic anion transporter family member 1B1). Several severely damaging variants were identified, studied, and validated. These observations have implications for both the application of deep mutational scanning to intrinsic membrane proteins and for the clinical effect of drugs and endogenous compounds transported by SLCO1B1.

Introduction Genome-wide association studies (GWAS) have confirmed the leucine-rich repeat kinase 2 (LRRK2) gene as a susceptibility locus for idiopathic Parkinson's disease (PD) in Caucasians. Though the rs1491942 and rs76904798 variants have shown the strongest associations, the causal variant(s) remains unresolved. Therefore, the aim of this study was to identify variants that may be driving the LRRK2 GWAS signal by sequencing the entire LRRK2 gene in Caucasian PD patients and controls. Methods A discovery series (287 PD patients, 294 controls) and replication series (362 PD patients, 168 controls) were included. The entire LRRK2 gene as well as 10 Kb upstream/downstream was sequenced. Candidate potential causal variants were considered to be those that (a) were in at least weak linkage disequilibrium with the two GWAS-nominated variants (rs1491942 and rs76904798), and (b) displayed an association odds ratio (OR) that is stronger than the two GWAS variants. Results Thirty-four candidate variants (all intronic/intergenic) that may drive the LRRK2 PD GWAS signal were identified in the discovery series. However, examination of the replication series for these variants did not reveal any with a consistently stronger OR than both PD GWAS variants. Evaluation of public databases to determine which candidate variants are most likely to have a direct functional effect on LRRK2 expression was inconclusive. Conclusion Though our findings provide novel insights into the LRRK2 GWAS association, a clear causal variant was not identified. The identified candidate variants can form the basis for future experiments and functional studies that can more definitively assess causal LRRK2 variants.

The human angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2) proteins play key roles in the cellular internalization of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the coronavirus responsible for the coronavirus disease of 2019 (COVID-19) pandemic. We set out to functionally characterize the ACE2 and TMPRSS2 protein abundance for variant alleles encoding these proteins that contained non-synonymous single-nucleotide polymorphisms (nsSNPs) in their open reading frames (ORFs). Specifically, a high-throughput assay, deep mutational scanning (DMS), was employed to test the functional implications of nsSNPs, which are variants of uncertain significance in these two genes. Specifically, we used a ‘landing pad’ system designed to quantify the protein expression for 433 nsSNPs that have been observed in the ACE2 and TMPRSS2 ORFs and found that 8 of 127 ACE2, 19 of 157 TMPRSS2 isoform 1 and 13 of 149 TMPRSS2 isoform 2 variant proteins displayed less than ~25% of the wild-type protein expression, whereas 4 ACE2 variants displayed 25% or greater increases in protein expression. As a result, we concluded that nsSNPs in genes encoding ACE2 and TMPRSS2 might potentially influence SARS-CoV-2 infectivity. These results can now be applied to DNA sequence data for patients infected with SARS-CoV-2 to determine the possible impact of patient-based DNA sequence variation on the clinical course of SARS-CoV-2 infection.

Single nucleotide variants in the open reading frames (ORFs) of pharmacogenes are important causes of interindividual variability in drug response. The functional characterization of variants of unknown significance within ORFs remains a major challenge for pharmacogenomics. Deep mutational scanning (DMS) is a high‐throughput technique that makes it possible to analyze the functional effect of hundreds of variants in a parallel and scalable fashion. We adapted a “landing pad” DMS system to study the function of missense variants in the ORFs of cytochrome P450 family 2 subfamily C member 9 (CYP2C9) and cytochrome P450 family 2 subfamily C member 19 (CYP2C19). We studied 230 observed missense variants in the CYP2C9 and CYP2C19 ORFs and found that 19 of 109 CYP2C9 and 36 of 121 CYP2C19 variants displayed less than ~ 25% of the wild‐type protein expression, a level that may have clinical relevance. Our results support DMS as an efficient method for the identification of damaging ORF variants that might have potential clinical pharmacogenomic application.

Accurate discovery of somatic mutations in a cell is a challenge that partially lays in immaturity of dedicated analytical approaches. Approaches comparing a cell’s genome to a control bulk sample miss common mutations, while approaches to find such mutations from bulk suffer from low sensitivity. We developed a tool, All2, which enables accurate filtering of mutations in a cell without the need for data from bulk(s). It is based on pair-wise comparisons of all cells to each other where every call for base pair substitution and indel is classified as either a germline variant, mosaic mutation, or false positive. As All2 allows for considering dropped-out regions, it is applicable to whole genome and exome analysis of cloned and amplified cells. By applying the approach to a variety of available data, we showed that its application reduces false positives, enables sensitive discovery of high frequency mutations, and is indispensable for conducting high resolution cell lineage tracing.

We previously identified five single nucleotide polymorphisms (SNPs) at four susceptibility loci for diffuse large B-cell lymphoma (DLBCL) in individuals of European ancestry through a large genome-wide association study (GWAS). To further elucidate genetic susceptibility to DLBCL, we sought to validate two loci at 3q13.33 and 3p24.1 that were suggestive in the original GWAS with additional genotyping. In the meta-analysis (5662 cases and 9237 controls) of the four original GWAS discovery scans and three replication studies, the 3q13.33 locus (rs9831894; minor allele frequency [MAF] = 0.40) was associated with DLBCL risk [odds ratio (OR) = 0.83, P = 3.62 × 10−13]. rs9831894 is in linkage disequilibrium (LD) with additional variants that are part of a super-enhancer that physically interacts with promoters of CD86 and ILDR1. In the meta-analysis (5510 cases and 12 817 controls) of the four GWAS discovery scans and four replication studies, the 3p24.1 locus (rs6773363; MAF = 0.45) was also associated with DLBCL risk (OR = 1.20, P = 2.31 × 10−12). This SNP is 29 426-bp upstream of the nearest gene EOMES and in LD with additional SNPs that are part of a highly lineage-specific and tumor-acquired super-enhancer that shows long-range interaction with AZI2 promoter. These loci provide additional evidence for the role of immune function in the etiology of DLBCL, the most common lymphoma subtype.

Background: DNA polymerase epsilon (POLE) is encoded by the POLE gene, and POLE-driven tumors are characterized by high mutational rates. POLE-driven tumors are relatively common in endometrial and colorectal cancer, and their presence is increasingly recognized in ovarian cancer (OC) of endometrioid type. POLE-driven cases possess an abundance of TCT>TAT and TCG>TTG somatic mutations characterized by mutational signature 10 from the Catalog of Somatic Mutations in Cancer (COSMIC). By quantifying the contribution of COSMIC mutational signature 10 in RNA sequencing (RNA-seq) we set out to identify POLE-driven tumors in a set of unselected Mayo Clinic OC. Methods: Mutational profiles were calculated using expressed single-nucleotide variants (eSNV) in the Mayo Clinic OC tumors (n=195), The Cancer Genome Atlas (TCGA) OC tumors (n=419), and the Genotype-Tissue Expression (GTEx) normal ovarian tissues (n=84). Non-negative Matrix Factorization (NMF) of the mutational profiles inferred the contribution per sample of four distinct mutational signatures, one of which corresponds to COSMIC mutational signature 10. Results: In the Mayo Clinic OC cohort we identified six tumors with a predicted contribution from COSMIC mutational signature 10 of over five mutations per megabase. These six cases harbored known POLE hotspot mutations (P286R, S297F, V411L, and A456P) and were of endometrioid histotype (P=5e-04). These six tumors were hypermutated with a higher tumor mutation load (mean, 54.02 mutations per megabase) compared to non-POLE endometrioid OC cases (mean, 7.69 mutations per megabase; P=5e-04), and had an early onset (average age of patients at onset, 48.33 years) when compared to non-POLE endometrioid OC cohort (average age at onset, 60.13 years; P=.008). Samples from TCGA and GTEx had a low COSMIC signature 10 contribution (median 0.16 mutations per megabase; maximum 1.78 mutations per megabase) and carried no POLE hotspot mutations.Conclusions: From the largest cohort of RNA-seq from endometrioid OC to date (n=53), we identified six hypermutated samples likely driven by POLE (frequency, 11%). Our result suggests the clinical need to screen for POLE driver mutations in endometrioid OC, which can guide enrollment in immunotherapy clinical trials.

By analyzing whole-exome data from the Alzheimer’s disease sequencing project (ADSP), we identify a set of 4 genes that show highly significant association with Alzheimer’s disease (AD). These genes were identified within a human TREM2 co-expression network using a novel approach wherein prioritized polygenic score analyses were performed sequentially to identify significant polygenic components. Two of the 4 genes (TREM2, RIN3) have previously been linked to AD and two (ATP8B4, IL17RA) are novel. Like TREM2, the 2 novel AD genes are selectively expressed in human microglial cells. The most significant variants in ATP8B4 and IL17RA are non-synonymous variants with strong effects comparable to the APOE ε4 and ε2 alleles. These protein-altering variants will provide unique opportunities to further explore the biological role of microglial cells in AD and help inform future immune modulatory therapeutic development for AD.

Post-zygotic mosaic mutations can be used to track cell lineages in humans. By using cell cloning and induced pluripotent cell lines, we analyzed early cell lineages in two living individuals (a patient and a control), and a postmortem human specimen. Of ten reconstructed post-zygotic divisions, none resulted in balanced contributions of daughter lineages to tissues. In both living individuals one of two lineages from the first cleavage was dominant across tissues, with 90% frequency in blood. We propose that the efficiency of DNA repair contributes to lineage imbalance. Allocation of lineages in postmortem brain correlated with anterior-posterior axis, associating lineage history with cell fate choices in embryos. Recurrence of germline variants as mosaic suggested that certain loci may be particularly susceptible to mutagenesis. We establish a minimally invasive framework for defining cell lineages in any living individual, which paves the way for studying their relevance in health and disease.

Background The study of mosaic mutation is important since it has been linked to cancer and various disorders. Single cell sequencing has become a powerful tool to study the genome of individual cells for the detection of mosaic mutations. The amount of DNA in a single cell needs to be amplified before sequencing and multiple displacement amplification (MDA) is widely used owing to its low error rate and long fragment length of amplified DNA. However, the phi29 polymerase used in MDA is sensitive to template fragmentation and presence of sites with DNA damage that can lead to biases such as allelic imbalance, uneven coverage and over representation of C to T mutations. It is therefore important to select cells with uniform amplification to decrease false positives and increase sensitivity for mosaic mutation detection. Results We propose a method, Scellector (single cell selector), which uses haplotype information to detect amplification quality in shallow coverage sequencing data. We tested Scellector on single human neuronal cells, obtained in vitro and amplified by MDA. Qualities were estimated from shallow sequencing with coverage as low as 0.3× per cell and then confirmed using 30× deep coverage sequencing. The high concordance between shallow and high coverage data validated the method. Conclusion Scellector can potentially be used to rank amplifications obtained from single cell platforms relying on a MDA-like amplification step, such as Chromium Single Cell profiling solution.

Double/triple hit lymphoma (DH/TH), known as high-grade B-cell lymphoma (HGBL), is an aggressive diffuse large B cell lymphoma (DLBCL), defined as having concurrent MYC, BCL2, and/or BCL6 gene rearrangements. While gene rearrangements represent significant genetic events in cancer, copy number alterations (CNAs) also play an important role, and their contributions to rearrangements have yet to be fully elucidated. Using FISH and high-resolution CNA data, we defined the landscape of concurrent gene rearrangements and copy gains in MYC, BCL2, and BCL6, in a cohort of 479 newly diagnosed DLBCL. We also show that concurrent translocations and copy number alterations, in combinations similar to DH/TH, identify a unique subset of DLBCL, alternative DH/TH, that have survival outcomes similar to DH/TH DLBCL patients.

Tumors acquire numerous mutations during development and progression. When translated into proteins, these mutations give rise to neoantigens that can be recognized by T cells and generate antibodies, representing an exciting direction of cancer immunotherapy. While neoantigens have been reported in many cancer types, the profiling of neoantigens often focused on the class-I subtype that are presented to CD8 + T cells, and the relationship between neoantigen load and clinical outcomes was often inconsistent among cancer types. In this study, we described an informatics workflow, REAL-neo, for identification, quality control (QC), and prioritization of both class-I and class-II human leukocyte antigen (HLA) bound neoantigens that arise from somatic single nucleotide mutations (SNM), small insertions and deletions (INDEL), and gene fusions. We applied REAL-neo to 835 primary breast tumors in the Cancer Genome Atlas (TCGA) and performed comprehensive profiling and characterization of the detected neoantigens. We found recurrent HLA class-I and class-II restricted neoantigens across breast cancer cases, and uncovered associations between neoantigen load and clinical traits. Both class-I and class-II neoantigen loads from SNM and INDEL were found to predict overall survival independent of tumor mutational burden (TMB), breast cancer subtypes, tumor-infiltrating lymphocyte (TIL) levels, tumor stage, and age at diagnosis. Our study highlighted the importance of accurate and comprehensive neoantigen profiling and QC, and is the first to report the predictive value of neoantigen load for overall survival in breast cancer.

Additional file 1. File containing supplementary figures.

Background: Follicular lymphoma (FL) exhibits significant clinical, cellular, molecular, and genetic heterogeneity. Our understanding of FL biology and molecular classifications of FL, to date, has largely been driven by pathologic classification (Grade 1-3b), genetic classification (m7-FLIPI), or gene expression profiling (IR-1/2; Huet-23), along with limited studies on small cohorts or targeted panels. In order to further understand the biological underpinnings and complexity of FL, large-scale and integrated whole exome sequencing (WES) and RNA sequencing (RNAseq) studies are needed. Using a highly-annotated cohort of 93 FL tumors with matched RNAseq, WES, and CyTOF data, we have explored the transcriptomic signature of purified FL B cells and identified unique molecular subsets that are defined by distinct pathway activation, immune content, and genomic signatures. Methods: Frozen cell suspensions from 93 untreated FL (Grade 1-3b) patients' tumor biopsies, enrolled in the University of Iowa/Mayo Clinic Lymphoma SPORE, were used for the study. DNA was isolated from whole tumor cell suspensions, and RNA was isolated from both whole tumor and B cell-enriched cell suspensions. RNAseq and WES were performed in the Mayo Clinic Genome Analysis Core. RNAseq and WES data were processed using a standard pipeline and novel driver genes were identified using Chasm+ driver analysis. Copy number variants were identified from WES data using GISTIC 2.0. NMF clustering and single sample gene set testing for B cell lineage and tumor microenvironment (TME) signatures were performed in R using the NMF and singscore packages. Results: Unsupervised clustering of RNAseq data identified three distinct expression programs which separated patient B cell samples into 3 groups: Group 1 (G1, n=20), Group 2 (G2, n=24), Group 3 (G3, n=43). While no clinical attributes were defined by any single group, G1 consisted of cases that had less aggressive characteristics (63% Stage I-II, 79% FLIPI 0-1). To identify unique transcriptional pathways driving the three expression programs, we scored gene level contributions to NMF expression programs and employed gene set enrichment analysis. This revealed significant pathway association with type-I IFN signaling (G1), DNA repair and stress response (G2), and epigenetic modulation (G3) as differentiating programs between the 3 groups (FDR q Conclusion: In this study, we have identified three unique FL tumor B cell groups, defined by specific transcriptional programs. Pathways such as inflammation, DNA damage response, and chromatin modification contribute most to the differences between B cell samples and group membership. Additionally, each program associated with specific genetic events as well as TME composition, highlighting potential drivers of these B cell states. This study improves the understanding of the mechanisms driving FL tumors and motivates further investigation into transcriptional consequences of genetic events as well as potential tumor intrinsic factors that may influence the TME. Figure 1 Figure 1. Disclosures Maurer: BMS: Research Funding; Genentech: Research Funding; Morphosys: Membership on an entity's Board of Directors or advisory committees, 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; Nanostring: Research Funding. Rimsza: NanoString Technologies: Other: Fee-for-service contract. Link: MEI: Consultancy; Genentech/Roche: Consultancy, Research Funding; Novartis, Jannsen: Research Funding. Habermann: Tess Therapeutics: Other: Data Monitoring Committee; Seagen: Other: Data Monitoring Committee; Incyte: Other: Scientific Advisory Board; Morphosys: Other: Scientific Advisory Board; Loxo Oncology: Other: Scientific Advisory Board; Eli Lilly & Co.,: Other: Scientific Advisor. Ansell: Bristol Myers Squibb, ADC Therapeutics, Seattle Genetics, Regeneron, Affimed, AI Therapeutics, Pfizer, Trillium and Takeda: Research Funding. King: Celgene/BMS: 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.

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.

Background Diffuse large B-cell lymphoma (DLBCL) is the most common aggressive non-Hodgkin lymphoma subtype with a 5-year survival of ~64%. While DLBCL is treated using immunochemotherapy (IC) with curative intent, 20%-40% of patients do not reach remission or relapse post IC. Clonal somatic mutations have been associated with aging, hematologic malignancies (predominately myeloid), and reduced OS in the general population. The objective of this study was to evaluate the association of clonal somatic mutations with event-free survival (EFS) and overall survival (OS). Methods We studied newly diagnosed DLBCL cases treated with IC who were enrolled in the Molecular Epidemiology Resource (MER) of the University of Iowa/Mayo Clinic Lymphoma Specialized Program of Research Excellence (SPORE). Clinical and treatment data were abstracted from medical records and all patients were systematically followed for disease outcomes. Pre-treatment DNA was extracted from matched peripheral blood and paraffin-embedded tumor tissue, and whole exome sequencing was conducted at 100x coverage. From the peripheral blood, allele counts from GVCF files produced by HaplotypeCaller in GATK were extracted for 42 genes commonly associated with clonal hematopoiesis (e.g. DNMT3A, TET2, and ASXL1). Mutations were deemed clonal somatic if the population minor allele frequency was Results The study consisted of 261 DLBCL patients treated with IC. The median age at diagnosis was 65 years (range 20-90) and 56% were male. With a median follow-up time of 5.1 years (range 0.1-15.0), there were 100 events and 80 deaths. A total of 17 (6.5%) patients had clonal somatic mutations, and16 patients were over the age of 60. A total of 8 (of 42) genes had clonal somatic mutations, with SF3B1, ASXL1, and TET2, being the most frequent (4 individuals per gene). VAFs ranged from 0.10 to 0.28 and none of the patients had multiple mutations. Additionally, the clonal somatic variants found in the peripheral blood were abscent in the tumor sample. Of the 17 patients with clonal somatic mutations, 12 had an event while 88 patients without a mutation had an event. In a univariate analysis, clonal somatic mutations were associated with inferior EFS (HR=2.55, 95% CI 1.39-4.68, p=0.002; Figure 1A). After adjusting for age, sex, and IPI, clonal somatic mutations remained associated with inferior EFS (HR=2.02, 95% CI 1.09-3.74, p=0.026). Clonal somatic mutations were also associated with inferior OS in the univariate analysis (HR=2.06, 95% CI 0.99-4.29, p=0.053), which attenuated after multivariate adjustment (HR=1.59, 95% CI: 0.76-3.34, p=0.22, Figure 1B). Although based on small numbers, mutations in SF3B1 were associated with inferior EFS (HR=3.25, 95% CI 1.16-9.12, p=0.025), but did not reach significance for OS (HR=2.56, 95% CI 0.78-8.38, p=0.120). Conclusions In this novel study of newly diagnosed DLBCL patients, clonal somatic mutations were identified in 6.5% of patients and were associated with inferior outcomes. Additional research is required at deeper sequencing to validate these findings and integration with tumor genomics is required to understand the prognosis of DLBCL patients with smaller clonal populations. Disclosures Shah: Dren Bio: Consultancy. Maurer:Celgene / BMS: Research Funding; Pfizer: Membership on an entity's Board of Directors or advisory committees; Morphosys: Membership on an entity's Board of Directors or advisory committees; Kite: Membership on an entity's Board of Directors or advisory committees; Nanostring: Research Funding. Martin:Celgene: Consultancy; Karyopharm: Consultancy, Research Funding; Morphosys: Consultancy; Regeneron: Consultancy; Incyte: Consultancy; Kite: Consultancy; Cellectar: Consultancy; Bayer: Consultancy; Beigene: Consultancy; Sandoz: Consultancy; I-MAB: Consultancy; Janssen: Consultancy; Teneobio: Consultancy. Witzig:AbbVie: Consultancy; Incyte: Consultancy; Acerta: Research Funding; Karyopharm Therapeutics: Research Funding; Immune Design: Research Funding; Spectrum: Consultancy; Celgene: Consultancy, Research Funding; MorphSys: Consultancy. Nowakowski:NanoString: Research Funding; Seattle Genetics: Consultancy; Curis: Consultancy; Kymera: Consultancy; Kite: Consultancy; Ryvu: Consultancy, Membership on an entity's Board of Directors or advisory committees; MorphoSys: Consultancy, Research Funding; Celgene/BMS: Consultancy, Research Funding. Novak:Celgene/BMS: Research Funding. Cerhan:NanoString: Research Funding; BMS/Celgene: Research Funding.

Background: Follicular Lymphoma (FL) is the second most common non-Hodgkin lymphoma and presents with significant clinical, cellular, molecular, and genetic heterogeneity. Despite the name and defining microanatomic location, the genetic and molecular identity and pathogenesis of the FL tumor cell is largely uncharacterized. Prior clinical and molecular classifications of FL have been primarily driven by pathologic classification (Grade 1-3b), genetic classification (M7-FLIPI), or gene expression profiling (IR-1/2). Using a unique cohort of 93 FL tumors, we have explored the transcriptomic signature of purified FL B cells, along with their matched whole tumor, and identified unique molecular subsets which are defined by distinct pathway activation, immune content, and genomic signatures identified through whole exome sequencing (WES). Methods: Frozen tumor biopsies from 93 untreated FL (Grade 1-3b) patients enrolled in the University of Iowa/Mayo Clinic Lymphoma SPORE were used for the study. DNA was isolated from whole tumor cell suspensions and RNA was isolated from both whole tumor and B cell enriched cell suspensions. RNA sequencing (RNAseq) and WES was performed in the Mayo Clinic Genome Analysis Core. RNAseq and WES data were processed using the Mayo Clinic standard pipeline and novel driver genes were identified using 20/20+ driver analysis. Copy number variants were identified using GISTIC 2.0. NMF clustering and single sample gene set testing, for B cell lineage and tumor microenvironment (TME) signatures, was performed in R using the NMF and SingScore packages. Results: NMF consensus clustering of FL B cell RNAseq data identified two distinct subsets, C1 (n=32) and C2 (n=57). Clinically, C1 was associated with being FL grade 3 (p5% frequency) were examined. TNFAIP3, TP53, and BCL6 alterations were enriched in C1 samples, whereas C2 associated with alterations in BCL2, KMT2D, CREBBP, REL, and MYC. Finally, B cell clusters were analyzed for TME signatures. C1 samples displayed significant enrichment of macrophage, cytotoxic cell, gamma-delta-Tcell, and endothelial cell TME elements (p Conclusion: Our results suggest that B cells from FL patients display two distinct transcriptomic signatures. C1 identifies an immunologically active tumor, driven by TNFAIP3 alterations, with pre-PBL characteristics, DNA replication and repair, inflammatory cytokine secretion/signaling, and hyper-metabolic characteristics. C2 identifies an immunologically quiet tumor, driven by alterations in BCL2 and chromatin modifiers, with an intermediate GC phenotype, repressed cytokine signaling, and active cell cycle progression and cytoskeleton rearrangement. This study improves our understanding of the mechanisms driving FL tumors and motivates further investigation into the relationship between tumor intrinsic factors that may influence the TME. Disclosures Maurer: Morphosys: Membership on an entity's Board of Directors or advisory committees; Kite: Membership on an entity's Board of Directors or advisory committees; Pfizer: Membership on an entity's Board of Directors or advisory committees; Nanostring: Research Funding; Celgene / BMS: Research Funding. Ansell:ADC Therapeutics: Research Funding; Trillium: Research Funding; Affimed: Research Funding; Regeneron: Research Funding; AI Therapeutics: Research Funding; Takeda: Research Funding; Seattle Genetics: Research Funding; Bristol Myers Squibb: Research Funding. Cerhan:NanoString: Research Funding; BMS/Celgene: Research Funding. Novak:Celgene/BMS: Research Funding.

Amplification of a CAG trinucleotide motif (CTG18.1) within the TCF4 gene has been strongly associated with Fuchs Endothelial Corneal Dystrophy (FECD). Nevertheless, a small minority of clinically unaffected elderly patients who have expanded CTG18.1 sequences have been identified. To test the hypothesis that the CAG expansions in these patients are protected from FECD because they have interruptions within the CAG repeats, we utilized a combination of an amplification-free, long-read sequencing method and a new target-enrichment sequence analysis tool developed by Pacific Biosciences to interrogate the sequence structure of expanded repeats. The sequencing was successful in identifying a previously described interruption within an unexpanded allele and provided sequence data on expanded alleles greater than 2000 bases in length. The data revealed considerable heterogeneity in the size distribution of expanded repeats within each patient. Detailed analysis of the long sequence reads did not reveal any instances of interruptions to the expanded CAG repeats, but did reveal novel variants within the AGG repeats that flank the CAG repeats in two of the five samples from clinically unaffected patients with expansions. This first examination of the sequence structure of CAG repeats in CTG18.1 suggests that factors other than interruptions to the repeat structure account for the absence of disease in some elderly patients with repeat expansions in the TCF4 gene.