Identification of candidate B-lymphoma genes by cross-species gene expression profiling

miR-28-based combination therapy impairs aggressive B cell lymphoma growth by rewiring DNA replication

Diffuse large B cell lymphoma (DLBCL) is the most common aggressive B cell lymphoma and accounts for nearly 40% of cases of B cell non-Hodgkin lymphoma. DLBCL is generally treated with R-CHOP chemotherapy, but many patients do not respond or relapse after treatment. Here, we analyzed the therapeutic potential of the tumor suppressor microRNA-28 (miR-28) for DLBCL, alone and in combination with the Bruton's tyrosine kinase inhibitor ibrutinib. Combination therapy with miR-28 plus ibrutinib potentiated the anti-tumor effects of monotherapy with either agent by inducing a specific transcriptional cell-cycle arrest program that impairs DNA replication. The molecular actions of miR-28 and ibrutinib synergistically impair DNA replication by simultaneous inhibition of origin activation and fork progression. Moreover, we found that downregulation of the miR-28-plus-ibrutinib gene signature correlates with better survival of ABC-DLBCL patients. These results provide evidence for the effectiveness of a new miRNA-based ibrutinib combination therapy for DLBCL and unveil the miR-28-plus-ibrutinib gene signature as a new predictor of outcome in ABC-DLBCL patients.

Bioinformatics Analysis of miRNAs Targeting TRAF5 in DLBCL Involving in NF-κB Signaling Pathway and Affecting the Apoptosis and Signal Transduction

Background

Diffuse large B-cell lymphoma (DLBCL) is an aggressive B-cell lymphoma with high heterogeneity. There is an unmet need to investigate valid indicators for the diagnosis and therapy of DLBCL.

Methods

GEO database was utilized to screen for differentially expressed genes (DEGs) and differential miRNAs in DLBCL tissues. The Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) were applied to analyse DEGs. Then multiple databases were searched for related miRNAs within DLBCL, TNF receptor-associated factor 5 (TRAF5) and NF-kappa B (NF-κB) signaling pathways. The KOBAS database was used to assist in the screening of miRNAs of interest and construct the regulatory network of miRNA-mRNA. Finally, the expression level and diagnostic performance of miRNAs were analyzed with GEO datasets, and DEGs were identified from the GEPIA database.

Results

DEGs were significantly concentrated in the NF-κB signaling pathway and cytokine-cytokine receptor interaction, and involved in the process of immune response and protein binding. MiR-15a-5p, miR-147a, miR-192-5p, miR-197-3p, miR-532-5p, and miR-650 were revealed to be targeting TRAF5 and participating in NF-κB signaling pathway and might impact the apoptosis and signal transduction of DLBCL. In the GEPIA database, TRAF5 was significantly overexpressed in DLBCL. The expression of miR-197-3p was upregulated within GEO datasets, while the rest of the miRNAs were downregulated in DLBCL.

Conclusions

Subsets of miRNAs may participate in the NF-κB signaling pathway by co-targeting TRAF5 and could be prospective biomarkers exploring the pathogenesis of DLBCL.

Targetable vulnerability of deregulated FOXM1/PLK1 signaling axis in diffuse large B cell lymphoma

FOXM1 is a transcription factor that controls cell cycle regulation, cell proliferation, and differentiation. Overexpression of FOXM1 has been implicated in various cancer types. However, the activation status and functional significance of FOXM1 in diffuse large B cell lymphoma (DLBCL) have not been well investigated. Using proteomic approaches, we discovered that the protein expression levels of FOXM1 and PLK1 were positively correlated in DLBCL cell lines and primary DLBCL. Expression levels of FOXM1 and PLK1 mRNAs were also significantly higher in DLBCL than in normal human B cells and could predict poor prognosis of DLBCL, particularly in patients with germinal center B cell-like (GCB) DLBCL. Furthermore, proteomic studies defined a FOXM1-PLK1 signature that consisted of proteins upstream and downstream of that axis involved in the p38-MAPK-AKT pathway, cell cycle, and DNA damage/repair. Further studies demonstrated a mechanistic function of the FOXM1/PLK1 axis in connection with the DNA damage response pathways regulating the S/G2 checkpoint of the cell cycle. Therapeutic targeting of FOXM1/PLK1 using a FOXM1 or PLK1 inhibitor, as well as other clinically relevant small-molecule inhibitors targeting ATR-CHK1, was highly effective in DLBCL in vitro models. These findings are instrumental for lymphoma drug discovery aiming at the FOXM1/PLK1/ATR/CHK1 axis.

The Tumor Suppressor Protein TRAF3 Modulates GSK3 Activity and Susceptibility of B Lymphoma Cells to GSK3 Inhibition

TNF receptor-associated factor 3 (TRAF3) is an adapter protein that inhibits many signals that promote B cell survival and activation. Mice with a B cell-specific TRAF3 deficiency and humans with a rare haploinsufficiency in TRAF3 have enhanced development of BCLs as they age. Loss-of-function mutations in TRAF3 are common in B cell malignancies. Recent studies show that pharmacological inhibition of the enzyme glycogen synthase kinase 3 (GSK3), which regulates cellular growth, survival, and metabolism, inhibits growth and survival of BCL-derived B cells. In this study, we found that TRAF3 and GSK3 associated in B cells. The relative levels of TRAF3 in BCL cell lines correlated positively with the ratio of inactive to total GSK3β, and negatively correlated with susceptibility to GSK3 inhibition by the GSK3 inhibitory drug 9-ING-41, currently in clinical trials. Uniquely in BCLs with low TRAF3, GSK3 inhibition caused increased loss of the TRAF3-regulated, anti-apoptotic protein Mcl-1. GSK3 inhibition also blocked hyperresponsiveness to IL-6 receptor signaling in TRAF3-deficient BCL cells. Together, these results support the utility of 9-ING-41 as a treatment for BCL, and suggest that a decrease or loss of TRAF3 in BCLs could act as a biomarker for increased susceptibility to GSK3 inhibitor treatment.

TNF receptor-associated factor 3 restrains B-cell receptor signaling in normal and malignant B cells

TRAF3 has diverse signaling functions, which vary by cell type. Uniquely in B lymphocytes, TRAF3 inhibits homeostatic survival. Highlighting the role of TRAF3 as a tumor suppressor, loss-of-function TRAF3 mutations are associated with human B-cell malignancies, while B-cell-specific deletion of TRAF3 in mice leads to autoimmunity and lymphoma development. The role of TRAF3 in inhibiting noncanonical NF-κB activation, CD40 and BAFF-R signaling to B cells is well documented. In contrast, TRAF3 enhances many T-cell effector functions, through associating with and enhancing signaling by the T-cell receptor (TCR)-CD28 complex. The present study was designed to determine the role of TRAF3 in signaling via the B-cell antigen receptor (BCR). The BCR is crucial for antigen recognition, survival, proliferation, and antibody production, and defects in BCR signaling can promote abnormal survival of malignant B cells. Here, we show that TRAF3 is associated with both CD79B and the BCR-activated kinases Syk and Btk following BCR stimulation. BCR-induced phosphorylation of Syk and additional downstream kinases was increased in TRAF3^−/− B cells, with regulation observed in both follicular and marginal zone B-cell subsets. BCR stimulation of TRAF3^−/− B cells resulted in increased surface expression of MHC-II, CD80, and CD86 molecules. Interestingly, increased survival of TRAF3^−/− primary B cells was resistant to inhibition of Btk, while TRAF3-deficient malignant B-cell lines showed enhanced sensitivity. TRAF3 serves to restrain normal and malignant BCR signaling, with important implications for its role in normal B-cell biology and abnormal survival of malignant B cells.

Identification of hub genes associated with the pathogenesis of diffuse large B-cell lymphoma subtype one characterized by host response via integrated bioinformatic analyses

Background

Host response diffuse large B-cell lymphoma (HR DLBCL) shares features of histologically defined T-cell/histiocyte-rich B-cell lymphoma, including fewer genetic abnormalities, frequent splenic and bone marrow involvement, and younger age at presentation. HR DLBCL is inherently less responsive to the standard treatment for DLBCL. Moreover, the mechanism of infiltration of HR DLBCL with preexisting abundant T-cells and dendritic cells is unknown, and their associated underlying immune responses incompletely defined. Here, hub genes and pathogenesis associated with HR DLBCL were explored to reveal molecular mechanisms and treatment targets.

Methods

Differentially expressed genes were identified in three datasets (GSE25638, GSE44337, GSE56315). The expression profile of the genes in the GSE53786 dataset was used to constructed a co-expression network. Protein-protein interactions analysis in the modules of interest identified candidate hub genes. Then screening of real hub genes was carried out by survival analysis within the GSE53786 and GSE10846 datasets. Expression of hub genes was validated in the Gene expression profiling interactive analysis, Oncomine databases and human tissue specimens. Functional enrichment analysis and Gene set enrichment analysis were utilized to investigate the potential mechanisms. Tumor Immune Estimation Resource and The Cancer Genome Atlas were used to mine the association of the hub gene with tumor immunity, potential upstream regulators were predicted using bioinformatics tools.

Results

A total of 274 common differentially expressed genes were identified. Within the key module, we identified CXCL10 as a real hub gene. The validation of upregulated expression level of CXCL10 was consistent with our study. CXCL10 might have a regulatory effect on tumor immunity. The predicted miRNA (hsa-mir-6849-3p) and transcription factor (IRF9) might regulate gene expression in the hub module.

Differential transcriptome response to proton versus X-ray radiation reveals novel candidate targets for combinatorial PT therapy in lymphoma

BACKGROUND AND PURPOSE

Knowledge of biological responses to proton therapy (PT) in comparison to X-ray remains in its infancy. Identification of PT specific molecular signals is an important opportunity for the discovery of biomarkers and synergistic drugs to advance clinical application. Since PT is used for the treatment of lymphoma, we report here transcriptomic responses of lymphoma cell lines to PT vs X-ray and identify potential therapeutic targets.

MATERIALS AND METHODS

Two lymphoma cell lines of human (BL41) and murine (J3D) origin were irradiated by X-ray and PT. Differential transcriptome regulation was quantified by RNA sequencing for each radiation type at 12 hours post irradiation. Gene-set enrichment analysis revealed deregulated molecular pathways and putative targets for lymphoma cell sensitization to PT.

RESULTS

Transcriptomic gene set enrichment analyses uncovered pathways that contribute to the unfolded protein response (UPR) and mitochondrial transport. Functional validation at multiple time points demonstrated increased UPR activation and decreased protein translation, perhaps due to increased oxidative stress and oxidative protein damage after PT. PPARgamma was identified as a potential regulator of the PT transcriptomic response. Inhibition of PPARgamma by two compounds, T0070907 and SR2595, sensitized lymphoma cells to PT.

CONCLUSIONS

Proton vs X-ray radiation leads to the transcriptional regulation of a specific subset of genes in line with diminished protein translation and UPR activation that may be due to oxidative stress. This study demonstrates that different radiation qualities trigger distinct cellular responses in lymphoma cells, and identifies PPARgamma inhibition as a potential strategy for the sensitization of lymphoma to PT.

Sulfiredoxin as a Potential Therapeutic Target for Advanced and Metastatic Prostate Cancer

The incidence of prostate cancer (PCa) is increasing, and it is currently the second most frequent cause of death by cancer in men. Despite advancements in cancer therapies, new therapeutic approaches are still needed for treatment-refractory advanced metastatic PCa. Cross-species analysis presents a robust strategy for the discovery of new potential therapeutic targets. This strategy involves the integration of genomic data from genetically engineered mouse models (GEMMs) and human PCa datasets. Considering the role of antioxidant pathways in tumor initiation and progression, we searched oxidative stress-related genes for a potential therapeutic target for PCa. First, we analyzed RNA-sequencing data from Pb-Cre4; Pten^f/f mice and discovered an increase in sulfiredoxin (Srxn1) mRNA expression in high-grade prostatic intraepithelial neoplasia (PIN), well-differentiated adenocarcinoma (medium-stage tumors), and poor-differentiated adenocarcinoma (advanced-stage prostate tumors). The increase of SRXN1 protein expression was confirmed by immunohistochemistry in mouse prostate tumor paraffin samples. Analyses of human databases and prostate tissue microarrays demonstrated that SRXN1 is overexpressed in a subset of high-grade prostate tumors and correlates with aggressive PCa with worse prognosis and decreased survival. Analyses in vitro showed that SRXN1 expression is also higher in most PCa cell lines compared to normal cell lines. Furthermore, siRNA-mediated downregulation of SRXN1 led to decreased viability of PCa cells LNCaP. In conclusion, we identified the antioxidant enzyme SRXN1 as a potential therapeutic target for PCa. Our results suggest that the use of specific SRXN1 inhibitors may be an effective strategy for the adjuvant treatment of castration-resistant PCa with SRXN1 overexpression.

Identification of potential therapeutic targets in prostate cancer through a cross-species approach

Genetically engineered mouse models of cancer can be used to filter genome-wide expression datasets generated from human tumours and to identify gene expression alterations that are functionally important to cancer development and progression. In this study, we have generated RNAseq data from tumours arising in two established mouse models of prostate cancer, PB-Cre/PtenloxP/loxP and p53loxP/loxPRbloxP/loxP, and integrated this with published human prostate cancer expression data to pinpoint cancer-associated gene expression changes that are conserved between the two species. To identify potential therapeutic targets, we then filtered this information for genes that are either known or predicted to be druggable. Using this approach, we revealed a functional role for the kinase MELK as a driver and potential therapeutic target in prostate cancer. We found that MELK expression was required for cell survival, affected the expression of genes associated with prostate cancer progression and was associated with biochemical recurrence.

FOXM1 is a therapeutic target for high-risk multiple myeloma

The transcription factor forkhead box M1 (FOXM1) is a validated oncoprotein in solid cancers, but its role in malignant plasma cell tumors such as multiple myeloma (MM) is unknown. We analyzed publicly available MM data sets and found that overexpression of FOXM1 prognosticates inferior outcome in a subset (~15%) of newly diagnosed cases, particularly patients with high-risk disease based on global gene expression changes. Follow-up studies using human myeloma cell lines (HMCLs) as the principal experimental model system demonstrated that enforced expression of FOXM1 increased growth, survival and clonogenicity of myeloma cells, whereas knockdown of FOXM1 abolished these features. In agreement with that, constitutive upregulation of FOXM1 promoted HMCL xenografts in laboratory mice, whereas inducible knockdown of FOXM1 led to growth inhibition. Expression of cyclin-dependent kinase 6 (CDK6) and NIMA-related kinase 2 (NEK2) was coregulated with FOXM1 in both HMCLs and myeloma patient samples, suggesting interaction of these three genes in a genetic network that may lend itself to targeting with small-drug inhibitors for new approaches to myeloma therapy and prevention. These results establish FOXM1 as high-risk myeloma gene and provide support for the design and testing of FOXM1-targeted therapies specifically for the FOXM1(High) subset of myeloma.

Control of asymmetric Hopfield networks and application to cancer attractors

The asymmetric Hopfield model is used to simulate signaling dynamics in gene regulatory networks. The model allows for a direct mapping of a gene expression pattern into attractor states. We analyze different control strategies aimed at disrupting attractor patterns using selective local fields representing therapeutic interventions. The control strategies are based on the identification of signaling bottlenecks, which are single nodes or strongly connected clusters of nodes that have a large impact on the signaling. We provide a theorem with bounds on the minimum number of nodes that guarantee control of bottlenecks consisting of strongly connected components. The control strategies are applied to the identification of sets of proteins that, when inhibited, selectively disrupt the signaling of cancer cells while preserving the signaling of normal cells. We use an experimentally validated non-specific and an algorithmically-assembled specific B cell gene regulatory network reconstructed from gene expression data to model cancer signaling in lung and B cells, respectively. Among the potential targets identified here are TP53, FOXM1, BCL6 and SRC. This model could help in the rational design of novel robust therapeutic interventions based on our increasing knowledge of complex gene signaling networks.

The NF-κB genomic landscape in lymphoblastoid B cells

The nuclear factor κB (NF-κΒ) subunits RelA, RelB, cRel, p50, and p52 are each critical for B cell development and function. To systematically characterize their responses to canonical and noncanonical NF-κB pathway activity, we performed chromatin immunoprecipitation followed by high-throughput DNA sequencing (ChIP-seq) analysis in lymphoblastoid B cell lines (LCLs). We found a complex NF-κB-binding landscape, which did not readily reflect the two NF-κB pathway paradigms. Instead, 10 subunit-binding patterns were observed at promoters and 11 at enhancers. Nearly one-third of NF-κB-binding sites lacked κB motifs and were instead enriched for alternative motifs. The oncogenic forkhead box protein FOXM1 co-occupied nearly half of NF-κB-binding sites and was identified in protein complexes with NF-κB on DNA. FOXM1 knockdown decreased NF-κB target gene expression and ultimately induced apoptosis, highlighting FOXM1 as a synthetic lethal target in B cell malignancy. These studies provide a resource for understanding mechanisms that underlie NF-κB nuclear activity and highlight opportunities for selective NF-κB blockade.