Non-IG aberrations of FOXP1 in B-cell malignancies lead to an aberrant expression of N-truncated isoforms of FOXP1

Comparison of B-Cell Lupus and Lymphoma Using a Novel Immune Imbalance Transcriptomics Algorithm Reveals Potential Therapeutic Targets

BACKGROUND/OBJECTIVES

Systemic lupus erythematosus (lupus) and B-cell lymphoma (lymphoma) co-occur at higher-than-expected rates and primarily depend on B cells for their pathology. These observations implicate shared inflammation-related B cell molecular mechanisms as a potential cause of co-occurrence.

METHODS

We consequently implemented a novel Immune Imbalance Transcriptomics (IIT) algorithm and applied IIT to lupus, lymphoma, and healthy B cell RNA-sequencing (RNA-seq) data to find shared and contrasting mechanisms that are potential therapeutic targets.

RESULTS

We observed 7143 significantly dysregulated genes in both lupus and lymphoma. Of those genes, we found 5137 to have a significant immune imbalance, defined as a significant dysregulation by both diseases, as analyzed by IIT. Gene Ontology (GO) term and pathway enrichment of the IIT genes yielded immune-related "Neutrophil Degranulation" and "Adaptive Immune System", which validates that the IIT algorithm isolates biologically relevant genes in immunity and inflammation. We found that 344 IIT gene products are known targets for established and/or repurposed drugs. Among our results, we found 48 known and 296 novel lupus targets, along with 151 known and 193 novel lymphoma targets. Known disease drug targets in our IIT results further validate that IIT isolates genes with disease-relevant mechanisms.

CONCLUSIONS

We anticipate the IIT algorithm, together with the shared and contrasting gene mechanisms uncovered here, will contribute to the development of immune-related therapeutic options for lupus and lymphoma patients.

Joint Secondary Transcriptomic Analysis of Non-Hodgkin's B-Cell Lymphomas Predicts Reliance on Pathways Associated with the Extracellular Matrix and Robust Diagnostic Biomarkers

Approximately 450,000 cases of Non-Hodgkin's lymphoma are annually diagnosed worldwide, resulting in ~240,000 deaths. An augmented understanding of the common mechanisms of pathology among larger numbers of B-cell Non-Hodgkin's Lymphoma (BCNHL) patients is sorely needed. We consequently performed a large joint secondary transcriptomic analysis of the available BCNHL RNA-sequencing projects from GEO, consisting of 322 relevant samples across ten distinct public studies, to find common underlying mechanisms and biomarkers across multiple BCNHL subtypes and patient subpopulations; limitations may include lack of diversity in certain ethnicities and age groups and limited clinical subtype diversity due to sample availability. We found ~10,400 significant differentially expressed genes (FDR-adjusted p-value < 0.05) and 33 significantly modulated pathways (Bonferroni-adjusted p-value < 0.05) when comparing BCNHL samples to non-diseased B-cell samples. Our findings included a significant class of proteoglycans not previously associated with lymphomas as well as significant modulation of genes that code for extracellular matrix-associated proteins. Our drug repurposing analysis predicted new candidates for repurposed drugs including ocriplasmin and collagenase. We also used a machine learning approach to identify robust BCNHL biomarkers that include YES1, FERMT2, and FAM98B, which have not previously been associated with BCNHL in the literature, but together provide ~99.9% combined specificity and sensitivity for differentiating lymphoma cells from healthy B-cells based on measurement of transcript expression levels in B-cells. This analysis supports past findings and validates existing knowledge while providing novel insights into the inner workings and mechanisms of transformed B-cell lymphomas that could give rise to improved diagnostics and/or therapeutics.

Preprocessing of Public RNA-Sequencing Datasets to Facilitate Downstream Analyses of Human Diseases

Publicly available RNA-sequencing (RNA-seq) data are a rich resource for elucidating the mechanisms of human disease; however, preprocessing these data requires considerable bioinformatic expertise and computational infrastructure. Analyzing multiple datasets with a consistent computational workflow increases the accuracy of downstream meta-analyses. This collection of datasets represents the human intracellular transcriptional response to disorders and diseases such as acute lymphoblastic leukemia (ALL), B-cell lymphomas, chronic obstructive pulmonary disease (COPD), colorectal cancer, lupus erythematosus; as well as infection with pathogens including Borrelia burgdorferi, hantavirus, influenza A virus, Middle East respiratory syndrome coronavirus (MERS-CoV), Streptococcus pneumoniae, respiratory syncytial virus (RSV), severe acute respiratory syndrome coronavirus (SARS-CoV), and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We calculated the statistically significant differentially expressed genes and Gene Ontology terms for all datasets. In addition, a subset of the datasets also includes results from splice variant analyses, intracellular signaling pathway enrichments as well as read mapping and quantification. All analyses were performed using well-established algorithms and are provided to facilitate future data mining activities, wet lab studies, and to accelerate collaboration and discovery.

Biological Difference Between Epstein-Barr Virus Positive and Negative Post-transplant Lymphoproliferative Disorders and Their Clinical Impact

Epstein–Barr virus (EBV) infection is correlated with several lymphoproliferative disorders, including Hodgkin disease, Burkitt lymphoma, diffuse large B-cell lymphoma (DLBCL), and post-transplant lymphoproliferative disorder (PTLD). The oncogenic EBV is present in 80% of PTLD. EBV infection influences immune response and has a causative role in the oncogenic transformation of lymphocytes. The development of PTLD is the consequence of an imbalance between immunosurveillance and immunosuppression. Different approaches have been proposed to treat this disorder, including suppression of the EBV viral load, reduction of immune suppression, and malignant clone destruction. In some cases, upfront chemotherapy offers better and durable clinical responses. In this work, we elucidate the clinicopathological and molecular-genetic characteristics of PTLD to clarify the biological differences of EBV(+) and EBV(–) PTLD. Gene expression profiling, next-generation sequencing, and microRNA profiles have recently provided many data that explore PTLD pathogenic mechanisms and identify potential therapeutic targets. This article aims to explore new insights into clinical behavior and pathogenesis of EBV(–)/(+) PTLD with the hope to support future therapeutic studies.

Novel markers in pediatric-type follicular lymphoma

The aim of this study was to review the histopathological, phenotypic, and molecular characteristics of pediatric-type follicular lymphoma (PTFL) and to assess the diagnostic value of novel immunohistochemical markers in distinguishing PTFL from follicular hyperplasia (FH). A total of 13 nodal PTFLs were investigated using immunohistochemistry, fluorescence in situ hybridization (FISH), and PCR and were compared with a further 20 reactive lymph nodes showing FH. Morphologically, PTFL cases exhibited a follicular growth pattern with irregular lymphoid follicles in which the germinal centers were composed of numerous blastoid cells showing a starry-sky appearance. Immunohistochemistry highlighted preserved CD10 (13/13) and BCL6 (13/13) staining, CD20 (13/13) positivity, a K light chain predominance (7/13), and partial BCL2 expression in 6/13 cases (using antibodies 124, E17, and SP66). The germinal center (GC)–associated markers stathmin and LLT-1 were positive in most of the cases (12/13 and 12/13, respectively). Interestingly, FOXP-1 was uniformly positive in PTFL (12/13 cases) in contrast to reactive GCs in FH, where only a few isolated positive cells were observed. FISH revealed no evidence of BCL2, BCL6, or MYC rearrangements in the examined cases. By PCR, clonal immunoglobulin gene rearrangements were detected in 100% of the tested PTFL cases. Our study confirmed the unique morphological and immunophenotypic features of PTFL and suggests that FOXP-1 can represent a novel useful diagnostic marker in the differential diagnosis between PTFL and FH.

FOXP1 expression is a prognostic biomarker in follicular lymphoma treated with rituximab and chemotherapy

Follicular lymphoma (FL) is a clinically and molecularly highly heterogeneous disease, yet prognostication relies predominantly on clinical tools. We recently demonstrated that integration of mutation status of 7 genes, including EZH2 and MEF2B, improves risk stratification. We mined gene expression data to uncover genes that are differentially expressed in EZH2- and MEF2B-mutated cases. We focused on FOXP1 and assessed its protein expression by immunohistochemistry (IHC) in 763 tissue biopsies. For outcome correlation, a population-based training cohort of 142 patients with FL treated with rituximab, cyclophosphamide, vincristine, and prednisone, and a clinical trial validation cohort comprising 395 patients treated with cyclophosphamide, doxorubicin, vincristine, and prednisone (CHOP) ± rituximab were used. We found FOXP1 to be significantly downregulated in both EZH2- and MEF2B-mutated cases. By IHC, 76 specimens in the training cohort (54%) had high FOXP1 expression (>10%), which was associated with reduced 5-year failure-free survival (FFS) rates (55% vs 70%). In the validation cohort, high FOXP1 expression status was observed in 248 patients (63%) and correlated with significantly shorter FFS in patients treated with R-CHOP (hazard ratio [HR], 1.95; P = .017) but not in patients treated with CHOP (HR, 1.15; P = .44). The impact of high FOXP1 expression on FFS in immunochemotherapy-treated patients was additional to the Follicular Lymphoma International Prognostic Index. High FOXP1 expression was associated with distinct molecular features such as TP53 mutations, expression of IRF4, and gene expression signatures reminiscent of dark zone germinal center or activated B cells. In summary, FOXP1 is a downstream phenotypic commonality of gene mutations and predicts outcome following rituximab-containing regimens.

Novel developments in the pathogenesis and diagnosis of extranodal marginal zone lymphoma

Extranodal marginal zone lymphoma (EMZL), mostly represented by mucosa-associated lymphoid tissue (MALT) type, also referred to as MALT lymphoma, is a clinically heterogeneous entity within the group of low-grade B cell lymphomas that arises in a wide range of different extranodal sites, including the stomach, lung, ocular adnexa, and skin. It represents the third most common non-Hodgkin lymphoma in the Western world, and the median age of occurrence is around 60 years. One characteristic aspect in a subset of EMZL detectable in about 25% of the cases is the presence of specific chromosomal translocations involving the genes MALT1 and BCL10, which lead to activation of the NF-κB signaling pathway. Another unique aspect is that several infectious agents, such as Helicobacter pylori in the case of gastric EMZL, and autoimmune disorders, like Sjögren syndrome, have been implicated in the pathogenesis of this cancer. Recent findings as summarized in this review have further improved our understanding of the complex pathobiology of this disease and have been essential to better define novel treatment strategies. In addition, many of these specific features are currently being implemented for the diagnosis of EMZL.

PUMILIO/FOXP1 signaling drives expansion of hematopoietic stem/progenitor and leukemia cells

RNA-binding proteins (RBPs) have emerged as important regulators of invertebrate adult stem cells, but their activities remain poorly appreciated in mammals. Using a short hairpin RNA strategy, we demonstrate here that the 2 mammalian RBPs, PUMILIO (PUM)1 and PUM2, members of the PUF family of posttranscriptional regulators, are essential for hematopoietic stem/progenitor cell (HSPC) proliferation and survival in vitro and in vivo upon reconstitution assays. Moreover, we found that PUM1/2 sustain myeloid leukemic cell growth. Through a proteomic approach, we identified the FOXP1 transcription factor as a new target of PUM1/2. Contrary to its canonical repressive activity, PUM1/2 rather promote FOXP1 expression by a direct binding to 2 canonical PUM responsive elements present in the FOXP1-3' untranslated region (UTR). Expression of FOXP1 strongly correlates with PUM1 and PUM2 levels in primary HSPCs and myeloid leukemia cells. We demonstrate that FOXP1 by itself supports HSPC and leukemic cell growth, thus mimicking PUM activities. Mechanistically, FOXP1 represses the expression of the p21^-CIP1 and p27^-KIP1 cell cycle inhibitors. Enforced FOXP1 expression reverses shPUM antiproliferative and proapoptotic activities. Altogether, our results reveal a novel regulatory pathway, underscoring a previously unknown and interconnected key role of PUM1/2 and FOXP1 in regulating normal HSPC and leukemic cell growth.

The small FOXP1 isoform predominantly expressed in activated B cell-like diffuse large B-cell lymphoma and full-length FOXP1 exert similar oncogenic and transcriptional activity in human B cells

The forkhead transcription factor FOXP1 is generally regarded as an oncogene in activated B cell-like diffuse large B-cell lymphoma. Previous studies have suggested that a small isoform of FOXP1 rather than full-length FOXP1, may possess this oncogenic activity. Corroborating those studies, we herein show that activated B cell-like diffuse large B-cell lymphoma cell lines and primary activated B cell-like diffuse large B-cell lymphoma cells predominantly express a small FOXP1 isoform, and that the 5′-end of the Foxp1 gene is a common insertion site in murine lymphomas in leukemia virus- and transposon-mediated insertional mutagenesis screens. By combined mass spectrometry, (quantative) reverse transcription polymerase chain reaction/sequencing, and small interfering ribonucleic acid-mediated gene silencing, we determined that the small FOXP1 isoform predominantly expressed in activated B cell-like diffuse large B-cell lymphoma lacks the N-terminal 100 amino acids of full-length FOXP1. Aberrant overexpression of this FOXP1 isoform (ΔN100) in primary human B cells revealed its oncogenic capacity; it repressed apoptosis and plasma cell differentiation. However, no difference in potency was found between this small FOXP1 isoform and full-length FOXP1. Furthermore, overexpression of full-length FOXP1 or this small FOXP1 isoform in primary B cells and diffuse large B-cell lymphoma cell lines resulted in similar gene regulation. Taken together, our data indicate that this small FOXP1 isoform and full-length FOXP1 have comparable oncogenic and transcriptional activity in human B cells, suggesting that aberrant expression or overexpression of FOXP1, irrespective of the specific isoform, contributes to lymphomagenesis. These novel insights further enhance the value of FOXP1 for the diagnostics, prognostics, and treatment of diffuse large B-cell lymphoma patients.

The significance of FOXP1 in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common subtype of mature B-cell lymphoma. While the majority of patients are cured with immunochemotherapy incorporating the anti-CD20 monoclonal antibody rituximab (R-CHOP), relapsed and refractory patients still have a dismal prognosis. DLBCL subtypes including an aggressive activated B-cell-like (ABC) and a more favorable prognosis germinal center-like (GCB) DLBCL have been identified by gene expression profiling and are characterized by distinct genetic abnormalities and oncogenic pathways. This identification of novel molecular targets is now enabling clinical trials to evaluate more effective personalized approaches to DLBCL therapy. The forkhead transcription factor FOXP1 is highly expressed in the ABC-DLBCL gene signature and has been extensively studied within the context of DLBCL for more than a decade. Here, we review the significance of FOXP1 in the pathogenesis of DLBCL, summarizing data supporting its utility as a prognostic and subtyping marker, its targeting by genetic aberrations, the importance of specific isoforms, and emerging data demonstrating a functional role in lymphoma biology. FOXP1 is one of the critical transcription factors whose deregulated expression makes important contributions to DLBCL pathogenesis. Thus, FOXP1 warrants further study as a potential theranostic in ABC-DLBCL.

N-terminally truncated FOXP1 protein expression and alternate internal FOXP1 promoter usage in normal and malignant B cells

Strong FOXP1 protein expression is a poor risk factor in diffuse large B-cell lymphoma and has been linked to an activated B-cell-like subtype, which preferentially expresses short FOXP1 (FOXP1S) proteins. However, both short isoform generation and function are incompletely understood. Here we prove by mass spectrometry and N-terminal antibody staining that FOXP1S proteins in activated B-cell-like diffuse large B-cell lymphoma are N-terminally truncated. Furthermore, a rare strongly FOXP1-expressing population of normal germinal center B cells lacking the N-terminus of the regular long protein (FOXP1L) was identified. Exon-targeted silencing and transcript analyses identified three alternate 5' non-coding exons [FOXP1-Ex6b(s), FOXP1-Ex7b and FOXP1-Ex7c], downstream of at least two predicted promoters, giving rise to FOXP1S proteins. These were differentially controlled by B-cell activation and methylation, conserved in murine lymphoma cells, and significantly correlated with FOXP1S protein expression in primary diffuse large B-cell lymphoma samples. Alternatively spliced isoforms lacking exon 9 (e.g. isoform 3) did not encode FOXP1S, and an alternate long human FOXP1 protein (FOXP1AL) likely generated from a FOXP1-Ex6b(L) transcript was detected. The ratio of FOXP1L:FOXP1S isoforms correlated with differential expression of plasmacytic differentiation markers in U-2932 subpopulations, and altering this ratio was sufficient to modulate CD19 expression in diffuse large B-cell lymphoma cell lines. Thus, the activity of multiple alternate FOXP1 promoters to produce multiple protein isoforms is likely to regulate B-cell maturation.

EBV-Positive and EBV-Negative Posttransplant Diffuse Large B Cell Lymphomas Have Distinct Genomic and Transcriptomic Features

The molecular pathogenesis of posttransplant diffuse large B cell lymphoma (PT-DLBCL) is largely unknown. We have recently shown that Epstein-Barr virus-positive (EBV(+)) and -negative (EBV(-)) PT-DLBCL have distinct gene expression profiles, and the transcriptomic profile of EBV(-) PT-DLBCL is similar to that of DLBCL in immunocompetent individuals (IC-DLBCL). To validate these observations at the genomic level, we performed array-comparative genome hybridization (aCGH) analysis of 21 EBV(+) PT-DLBCL, 6 EBV(-) PT-DLBCL, and 11 control IC-DLBCL, and subsequently combined genomic and transcriptomic data. The analysis showed that EBV(+) and EBV(-) PT-DLBCL have distinct aCGH profiles and shared only one recurrent imbalance. EBV(-) PT-DLBCL, however, displayed at least 10 aberrations recurrent in IC-DLBCL, among which characteristic gain of 3/3q and 18q, and loss of 6q23/TNFAIP3 as well as 9p21/CDKN2A. The most prevalent aberration in EBV(+) PT-DLBCL was gain/amplification of 9p24.1 targeting PDCD1LG2/PDL2. Our data indicate that the FOXP1 oncogene and the tumor suppressor CDKNA2 implicated in EBV(-) DLBCL, do not play a critical role in the pathogenesis of EBV(+) PT-DLBCL. Altogether, genomic profiling of PT-/IC-DLBCL confirms that EBV(-) and EBV(+) PT-DLBCL are distinct entities, while EBV(-) PT-DLBCL has features in common with IC-DLBCL. These findings support the hypothesis that EBV(-) PT-DLBCL are de novo lymphomas in transplant recipients.

The hematopoietic oncoprotein FOXP1 promotes tumor cell survival in diffuse large B-cell lymphoma by repressing S1PR2 signaling

Aberrant expression of the oncogenic transcription factor forkhead box protein 1 (FOXP1) is a common feature of diffuse large B-cell lymphoma (DLBCL). We have combined chromatin immunoprecipitation and gene expression profiling after FOXP1 depletion with functional screening to identify targets of FOXP1 contributing to tumor cell survival. We find that the sphingosine-1-phosphate receptor 2 (S1PR2) is repressed by FOXP1 in activated B-cell (ABC) and germinal center B-cell (GCB) DLBCL cell lines with aberrantly high FOXP1 levels; S1PR2 expression is further inversely correlated with FOXP1 expression in 3 patient cohorts. Ectopic expression of wild-type S1PR2, but not a point mutant incapable of activating downstream signaling pathways, induces apoptosis in DLBCL cells and restricts tumor growth in subcutaneous and orthotopic models of the disease. The proapoptotic effects of S1PR2 are phenocopied by ectopic expression of the small G protein Gα13 but are independent of AKT signaling. We further show that low S1PR2 expression is a strong negative prognosticator of patient survival, alone and especially in combination with high FOXP1 expression. The S1PR2 locus has previously been demonstrated to be recurrently mutated in GCB DLBCL; the transcriptional silencing of S1PR2 by FOXP1 represents an alternative mechanism leading to inactivation of this important hematopoietic tumor suppressor.

FOXP1 suppresses immune response signatures and MHC class II expression in activated B-cell-like diffuse large B-cell lymphomas

The FOXP1 (forkhead box P1) transcription factor is a marker of poor prognosis in diffuse large B-cell lymphoma (DLBCL). Here microarray analysis of FOXP1-silenced DLBCL cell lines identified differential regulation of immune response signatures and major histocompatibility complex class II (MHC II) genes as some of the most significant differences between germinal center B-cell (GCB)-like DLBCL with full-length FOXP1 protein expression versus activated B-cell (ABC)-like DLBCL expressing predominantly short FOXP1 isoforms. In an independent primary DLBCL microarray data set, multiple MHC II genes, including human leukocyte antigen DR alpha chain (HLA-DRA), were inversely correlated with FOXP1 transcript expression (P<0.05). FOXP1 knockdown in ABC-DLBCL cells led to increased cell-surface expression of HLA-DRA and CD74. In R-CHOP (rituximab, cyclophosphamide, doxorubicin, vincristine and prednisone)-treated DLBCL patients (n=150), reduced HLA-DRA (<90% frequency) expression correlated with inferior overall survival (P=0.0003) and progression-free survival (P=0.0012) and with non-GCB subtype stratified by the Hans, Choi or Visco-Young algorithms (all P<0.01). In non-GCB DLBCL cases with <90% HLA-DRA, there was an inverse correlation with the frequency (P=0.0456) and intensity (P=0.0349) of FOXP1 expression. We propose that FOXP1 represents a novel regulator of genes targeted by the class II MHC transactivator CIITA (MHC II and CD74) and therapeutically targeting the FOXP1 pathway may improve antigen presentation and immune surveillance in high-risk DLBCL patients.

FOXP1 potentiates Wnt/β-catenin signaling in diffuse large B cell lymphoma

The transcription factor FOXP1 (forkhead box protein P1) is a master regulator of stem and progenitor cell biology. In diffuse large B cell lymphoma (DLBCL), copy number amplifications and chromosomal translocations result in overexpression of FOXP1. Increased abundance of FOXP1 in DLBCL is a predictor of poor prognosis and resistance to therapy. We developed a genome-wide, mass spectrometry-coupled, gain-of-function genetic screen, which revealed that FOXP1 potentiates β-catenin-dependent, Wnt-dependent gene expression. Gain- and loss-of-function studies in cell models and zebrafish confirmed that FOXP1 was a general and conserved enhancer of Wnt signaling. In a Wnt-dependent fashion, FOXP1 formed a complex with β-catenin, TCF7L2 (transcription factor 7-like 2), and the acetyltransferase CBP [CREB (adenosine 3',5'-monophosphate response element-binding protein)-binding protein], and this complex bound the promoters of Wnt target genes. FOXP1 promoted the acetylation of β-catenin by CBP, and acetylation was required for FOXP1-mediated potentiation of β-catenin-dependent transcription. In DLBCL, we found that FOXP1 promoted sensitivity to Wnt pathway inhibitors, and knockdown of FOXP1 or blocking β-catenin transcriptional activity slowed xenograft tumor growth. These data connect excessive FOXP1 with β-catenin-dependent signal transduction and provide a molecular rationale for Wnt-directed therapy in DLBCL.