Follicular Lymphoma Tumor-Cell Transcriptional Programs Associate with Distinct Somatic Alterations and Tumor-Immune Microenvironments
Authors: Krull Jordan E. , Wenzl Kerstin , Manske Michelle K. , Hopper Melissa A. , Larson Melissa C. , Sarangi Vivekananda , Maurer Matthew J. , Yang Zhi-Zhang , Rimsza Lisa M. , Link Brian K. , Habermann Thomas M. , Ansell Stephen M. , King Rebecca L. , Cerhan James R. , Novak Anne J. 
Published: November 2021

Abstract
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&lt;0.001). VIPER master regulator activity inferencing revealed that these pathways were likely being controlled by differential activity in NF-kB, IRFs, STAT1, BCL6, and FOXO1. Each program significantly enriched for, but were not defined by, portions of specific germinal center programs, such as pre-memory (G1), light-zone-to-dark-zone transition (G2), and a pre-light-zone intermediate (G3). We next assessed the connection between B cell programs and the tumor microenvironment (TME) using available paired CyTOF data on 67 cases, which revealed an active TME in G1, with an abundance of CD8 T cell and NK cell populations, a wide variety of immune content in G2 that consisted mostly of Tfh and myeloid cells, and a poorly populated immune compartment in G3 compared to G1 and G2. Finally, somatic driver mutations and copy number alterations from WES were identified and associated with the three clusters. The three groups distinguished themselves by significant enrichment of copy number alterations (TNFAIP3-loss , 1q23-gain, 1q32-gain) in G2, while 10q-loss and mutations in BCL2 and chromatin modifiers (KMT2D and CREBBP) enriched in G3. G1, overall, had lower alteration burden and had weak associations with any specific alterations, suggesting an alternative mechanism for driving the G1 program.
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.