Artificial Intelligence Analysis of the Gene Expression of Follicular Lymphoma Predicted the Overall Survival and Correlated with the Immune Microenvironment Response Signatures

Follicular lymphoma (FL) is the second most common lymphoma in Western countries. FL is characterized by being incurable, usually having an indolent clinical course with frequent relapses, and an eventual patient’s death or transformation to Diffuse Large B-cell Lymphoma. The immune response and the tumoral immune microenvironment, including FOXP3+Tregs, PD-1+TFH cells, TNFRSF14 (HVEM), and BTLA play a role in the pathogenesis. We aimed to analyze the gene expression of FL by Artificial Intelligence (machine learning, deep learning), to identify genes associated with the prognosis of the patients and with the microenvironment in terms of overall survival (OS). A series of 184 cases of the GSE16131 dataset was analyzed by multilayer perceptron (MLP) and radial basis function (RBF) neural networks. In the analysis, MLP and RBF had a synergistic effect. From an initial set of 22,215 genes probes, a final set of 43 genes was highlighted. These 43 genes predicted the OS and correlated with the immune microenvironment: in a multivariate Cox analysis, 18 genes were associated with a poor prognosis (namely, MED8, KRT19, CDC40, SLC24A2, PRB1, KIAA0100, EVA1B, KLK10, TMEM70, BTN2A3P, TRPM4, MED6, FRYL, CBFA2T2, RANBP9, BNIP2, PTP4A2 and ALDH1L1) and 25 genes were associated with a good prognosis of the patients. Gene set enrichment analysis (GSEA) confirmed these findings and showed a typical sinusoidal-like shape. Some of the most relevant genes for poor OS were EVA1B, KRT19, BTN2A3P, KLK10, TRPM4, TMEM70, and SLC24A2 (hazard risk = from 1.7 to 4.3, p < 0.005) and for good OS, these were TDRD12 and ZNF230 (HR = 0.34 and 0.28, p < 0.001). EVA1B, KRT19, BTN2AP3, KLK10, and TRPM4 also associated with M2-like macrophage markers including CD163, MRC1 (CD206), and IL10 in the core enrichment for dead OS outcome by GSEA and to poor OS by Kaplan–Meier with Log rank test. The scientific literature showed that some of these genes also play a role in other types of cancer. In conclusion, by Artificial Intelligence, we have identified new biomarkers with prognostic relevance in FL

A Combination of Multilayer Perceptron, Radial Basis Function Artificial Neural Networks and Machine Learning Image Segmentation for the Dimension Reduction and the Prognosis Assessment of Diffuse Large B-Cell Lymphoma

The prognosis of diffuse large B-cell lymphoma (DLBCL) is heterogeneous. Therefore, we aimed to highlight predictive biomarkers. First, artificial intelligence was applied into a discovery series of gene expression of 414 patients (GSE10846). A dimension reduction algorithm aimed to correlate with the overall survival and other clinicopathological variables; and included a combination of Multilayer Perceptron (MLP) and Radial Basis Function (RBF) artificial neural networks, gene-set enrichment analysis (GSEA), Cox regression and other machine learning and predictive analytics modeling [C5.0 algorithm, logistic regression, Bayesian Network, discriminant analysis, random trees, tree-AS, Chi-squared Automatic Interaction Detection CHAID tree, Quest, classification and regression (C&R) tree and neural net)]. From an initial 54,613 gene-probes, a set of 488 genes and a final set of 16 genes were defined. Secondly, two identified markers of the immune checkpoint, PD-L1 (CD274) and IKAROS (IKZF4), were validated in an independent series from Tokai University, and the immunohistochemical expression was quantified, using a machine-learning-based Weka segmentation. High PD-L1 associated with poor overall and progression-free survival, non-GCB phenotype, Epstein–Barr virus infection (EBER+), high RGS1 expression and several clinicopathological variables, such as high IPI and absence of clinical response. Conversely, high expression of IKAROS was associated with a good overall and progression-free survival, GCB phenotype and a positive clinical response to treatment. Finally, the set of 16 genes (PAF1, USP28, SORT1, MAP7D3, FITM2, CENPO, PRCC, ALDH6A1, CSNK2A1, TOR1AIP1, NUP98, UBE2H, UBXN7, SLC44A2, NR2C2AP and LETM1), in combination with PD-L1, IKAROS, BCL2, MYC, CD163 and TNFAIP8, predicted the survival outcome of DLBCL with an overall accuracy of 82.1%. In conclusion, building predictive models of DLBCL is a feasible analytical strategy

The Use of the Random Number Generator and Artificial Intelligence Analysis for Dimensionality Reduction of Follicular Lymphoma Transcriptomic Data

Follicular lymphoma (FL) is one of the most frequent subtypes of non-Hodgkin lymphomas. This research predicted the prognosis of 184 untreated follicular lymphoma patients (LLMPP GSE16131 series), using gene expression data and artificial intelligence (AI) neural networks. A new strategy based on the random number generation was used to create 120 different and independent multilayer perceptron (MLP) solutions, and 22,215 gene probes were ranked according to their averaged normalized importance for predicting the overall survival. After dimensionality reduction, the final neural network architecture included (1) newly identified predictor genes related to cell adhesion and migration, cell signaling, and metabolism (EPB41L4B, MOCOS, SPIN2A, BTD, SRGAP3, CTNS, PRB1, L1CAM, and CEP57); (2) the international prognostic index (IPI); and (3) other relevant immuno-oncology, immune microenvironment, and checkpoint markers (CD163, CSF1R, FOXP3, PDCD1, TNFRSF14 (HVEM), and IL10). The performance of this neural network was good, with an area under the curve (AUC) of 0.89. A comparison with other machine learning techniques (C5 tree, logistic regression, Bayesian network, discriminant analysis, KNN algorithms, LSVM, random trees, SVM, tree-AS, XGBoost linear, XGBoost tree, CHAID, Quest, C&R tree, random forest, and neural network) was also made. In conclusion, the overall survival of follicular lymphoma was predicted with a neural network with high accuracy

Copy Number Alteration and Mutational Profile of High-Grade B-Cell Lymphoma with MYC and BCL2 and/or BCL6 Rearrangements, Diffuse Large B-Cell Lymphoma with MYC-Rearrangement, and Diffuse Large B-Cell Lymphoma with MYC-Cluster Amplification.

The authors thank the technicians and scientists of the Genomics Core Facility of the IDIBAPS for their assistance in performing the NGS analysis.Diffuse large B-cell lymphoma (DLBCL) with MYC alteration is classified as high-grade B-cell lymphoma with MYC and BCL2 and/or BCL6 rearrangements (double/triple-hit lymphoma; DHL/THL), DLBCL with MYC rearrangement (single-hit lymphoma; SHL) and DLBCL with MYC-cluster amplification (MCAD). To elucidate the genetic features of DHL/THL, SHL, and MCAD, 23 lymphoma cases from Tokai University Hospital were analyzed. The series included 10 cases of DHL/THL, 10 cases of SHL and 3 cases of MCAD. The analysis used whole-genome copy number microarray analysis (OncoScan) and a custom-made next-generation sequencing (NGS) panel of 115 genes associated with aggressive B-cell lymphomas. The copy number alteration (CNA) profiles were similar between DHL/THL and SHL. MCAD had fewer CNAs than those of DHL/THL and SHL, except for +8q24. The NGS profile characterized DHL/THL with a higher "mutation burden" than SHL (17 vs. 10, p = 0.010), and the most relevant genes for DHL/THL were BCL2 and SOCS1, and for SHL was DTX1. MCAD was characterized by mutations of DDX3X, TCF3, HLA-A, and TP53, whereas MYC was unmutated. In conclusion, DHL/THL, SHL, and MCAD have different profiles.S

High Expression of Caspase-8 Associated with Improved Survival in Diffuse Large B-Cell Lymphoma: Machine Learning and Artificial Neural Networks Analyses

High expression of the anti-apoptotic TNFAIP8 is associated with poor survival of the patients with diffuse large B-cell lymphoma (DLBCL), and one of the functions of TNFAIP8 is to inhibit the pro-apoptosis Caspase-8. We aimed to analyze the immunohistochemical expression of Caspase-8 (active subunit p18; CASP8) in a series of 97 cases of DLBCL from Tokai University Hospital, and to correlate with other Caspase-8 pathway-related markers, including cleaved Caspase-3, cleaved PARP, BCL2, TP53, MDM2, MYC, Ki67, E2F1, CDK6, MYB and LMO2. After digital image quantification, the correlation with several clinicopathological characteristics of the patients showed that high protein expression of Caspase-8 was associated with a favorable overall and progression-free survival (Hazard Risks = 0.3; p = 0.005 and 0.03, respectively). Caspase-8 also positively correlated with cCASP3, MDM2, E2F1, TNFAIP8, BCL2 and Ki67. Next, the Caspase-8 protein expression was modeled using predictive analytics, and a high overall predictive accuracy (>80%) was obtained with CHAID decision tree, Bayesian network, discriminant analysis, C5 tree, logistic regression, and Artificial Intelligence Neural Network methods (both Multilayer perceptron and Radial basis function); the most relevant markers were cCASP3, E2F1, TP53, cPARP, MDM2, BCL2 and TNFAIP8. Finally, the CASP8 gene expression was also successfully modeled in an independent DLBCL series of 414 cases from the Lymphoma/Leukemia Molecular Profiling Project (LLMPP). In conclusion, high protein expression of Caspase-8 is associated with a favorable prognosis of DLBCL. Predictive modeling is a feasible analytic strategy that results in a solution that can be understood (i.e., explainable artificial intelligence, “white-box” algorithms)

Integrative Statistics, Machine Learning and Artificial Intelligence Neural Network Analysis Correlated CSF1R with the Prognosis of Diffuse Large B-Cell Lymphoma

Tumor-associated macrophages (TAMs) of the immune microenvironment play an important role in the Diffuse Large B-cell Lymphoma (DLBCL) pathogenesis. This research aimed to characterize the expression of macrophage colony-stimulating factor 1 receptor (CSF1R) at the gene and protein level in correlation with survival. First, the immunohistochemical expression of CSF1R was analyzed in a series of 198 cases from Tokai University Hospital and two patterns of histological expression were found, a TAMs, and a diffuse B-lymphocytes pattern. The clinicopathological correlations showed that the CSF1R + TAMs pattern associated with a poor progression-free survival of the patients, disease progression, higher MYC proto-oncogene expression, lower MDM2 expression, BCL2 translocation, and a MYD88 L265P mutation. Conversely, a diffuse CSF1R + B-cells pattern was associated with a favorable progression-free survival. Second, the histological expression of CSF1R was also correlated with 10 CSF1R-related markers including CSF1, STAT3, NFKB1, Ki67, MYC, PD-L1, TNFAIP8, IKAROS, CD163, and CD68. CSF1R moderately correlated with STAT3, TNFAIP8, CD68, and CD163 in the cases with the CSF1R + TAMs pattern. In addition, machine learning modeling predicted the CSF1R immunohistochemical expression with high accuracy using regression, generalized linear, an artificial intelligence neural network (multilayer perceptron), and support vector machine (SVM) analyses. Finally, a multilayer perceptron analysis predicted the genes associated with the CSF1R gene expression using the GEO GSE10846 DLBCL series of the Lymphoma/Leukemia Molecular Profiling Project (LLMPP), with correlation to the whole set of 20,683 genes as well as with an immuno-oncology cancer panel of 1790 genes. In addition, CSF1R positively correlated with SIRPA and inversely with CD47. In conclusion, the CSF1R histological pattern correlated with the progression-free survival of the patients of the Tokai series, and predictive analytics is a feasible strategy in DLBCL

Artificial Neural Networks Predicted the Overall Survival and Molecular Subtypes of Diffuse Large B-Cell Lymphoma Using a Pancancer Immune-Oncology Panel

Diffuse large B-cell lymphoma (DLBCL) is one of the most frequent subtypes of non-Hodgkin lymphomas. We used artificial neural networks (multilayer perceptron and radial basis function), machine learning, and conventional bioinformatics to predict the overall survival and molecular subtypes of DLBCL. The series included 106 cases and 730 genes of a pancancer immune oncology panel (nCounter) as predictors. The multilayer perceptron predicted the outcome with high accuracy, with an area under the curve (AUC) of 0.98, and ranked all the genes according to their importance. In a multivariate analysis, ARG1, TNFSF12, REL, and NRP1 correlated with favorable survival (hazard risks: 0.3–0.5), and IFNA8, CASP1, and CTSG, with poor survival (hazard risks = 1.0–2.1). Gene set enrichment analysis (GSEA) showed enrichment toward poor prognosis. These high-risk genes were also associated with the gene expression of M2-like tumor-associated macrophages (CD163), and MYD88 expression. The prognostic relevance of this set of 7 genes was also confirmed within the IPI and MYC translocation strata, the EBER-negative cases, the DLBCL not-otherwise specified (NOS) (High-grade B-cell lymphoma with MYC and BCL2 and/or BCL6 rearrangements excluded), and an independent series of 414 cases of DLBCL in Europe and North America (GSE10846). The perceptron analysis also predicted molecular subtypes (based on the Lymph2Cx assay) with high accuracy (AUC = 1). STAT6, TREM2, and REL were associated with the germinal center B-cell (GCB) subtype, and CD37, GNLY, CD46, and IL17B were associated with the activated B-cell (ABC)/unspecified subtype. The GSEA had a sinusoidal-like plot with association to both molecular subtypes, and immunohistochemistry analysis confirmed the correlation of MAPK3 with the GCB subtype in another series of 96 cases (notably, MAPK3 also correlated with LMO2, but not with M2-like tumor-associated macrophage markers CD163, CSF1R, TNFAIP8, CASP8, PD-L1, PTX3, and IL-10). Finally, survival and molecular subtypes were successfully modeled using other machine learning techniques including logistic regression, discriminant analysis, SVM, CHAID, C5, C&R trees, KNN algorithm, and Bayesian network. In conclusion, prognoses and molecular subtypes were predicted with high accuracy using neural networks, and relevant genes were highlighted

High PTX3 expression is associated with a poor prognosis in diffuse large B-cell lymphoma

Tumor-associated macrophages (TAMs) are associated with a poor prognosis of diffuse large B-cell lymphoma (DLBCL). As macrophages are heterogeneous, the immune polarization and their pathological role warrant further study. We characterized the microenvironment of DLBCL by immunohistochemistry in a training set of 132 cases, which included 10 Epstein–Barr virus-encoded small RNA (EBER)-positive and five high-grade B-cell lymphomas, with gene expression profiling in a representative subset of 37 cases. Diffuse large B-cell lymphoma had a differential infiltration of TAMs. The high infiltration of CD68 (pan-macrophages), CD16 (M1-like), CD163, pentraxin 3 (PTX3), and interleukin (IL)-10-positive macrophages (M2c-like) and low infiltration of FOXP3-positive regulatory T lymphocytes (Tregs) correlated with poor survival. Activated B cell-like DLBCL was associated with high CD16, CD163, PTX3, and IL-10, and EBER-positive DLBCL with high CD163 and PTX3. Programmed cell death-ligand 1 positively correlated with CD16, CD163, IL-10, and RGS1. In a multivariate analysis of overall survival, PTX3 and International Prognostic Index were identified as the most relevant variables. The gene expression analysis showed upregulation of genes involved in innate and adaptive immune responses and macrophage and Toll-like receptor pathways in high PTX3 cases. The prognostic relevance of PTX3 was confirmed in a validation set of 159 cases. Finally, in a series from Europe and North America (GSE10846, R-CHOP-like treatment, n = 233) high gene expression of PTX3 correlated with poor survival, and moderately with CSF1R, CD16, MITF, CD163, MYC, and RGS1. Therefore, the high infiltration of M2c-like immune regulatory macrophages and low infiltration of FOXP3-positive Tregs is associated with a poor prognosis in DLBCL, for which PTX3 is a new prognostic biomarker

High TNFRSF14 and low BTLA are associated with poor prognosis in Follicular Lymphoma and in Diffuse Large B-cell Lymphoma transformation

The microenvironment influences the behavior of follicular lymphoma (FL) but the specific roles of the immunomodulatory BTLA and TNFRSF14 (HVEM) are unknown. Therefore, we examined their immunohistochemical expression in the intrafollicular, interfollicular and total histological compartments in 106 FL cases (57M/49F; median age 57-years), and in nine relapsed-FL with transformation to DLBCL (tFL). BTLA expression pattern was of follicular T-helper cells (TFH) in the intrafollicular and of T-cells in the interfollicular compartments. The mantle zones were BTLA+ in 35.6% of the cases with similar distribution of IgD. TNFRSF14 expression pattern was of neoplastic B lymphocytes (centroblasts) and "tingible body macrophages". At diagnosis, the averages of total BTLA and TNFRSF14-positive cells were 19.2%±12.4STD (range, 0.6%-58.2%) and 46.7 cells/HPF (1-286.5), respectively. No differences were seen between low-grade vs. high-grade FL but tFL was characterized by low BTLA and high TNFRSF14 expression. High BTLA correlated with good overall survival (OS) (total-BTLA, Hazard Risk=0.479, P=0.022) and with high PD-1 and FOXP3+Tregs. High TNFRSF14 correlated with poor OS and progression-free survival (PFS) (total-TNFRSF14, HR=3.9 and 3.2, respectively, P<0.0001), with unfavorable clinical variables and higher risk of transformation (OR=5.3). Multivariate analysis including BTLA, TNFRSF14 and FLIPI showed that TNFRSF14 and FLIPI maintained prognostic value for OS and TNFRSF14 for PFS. In the GSE16131 FL series, high TNFRSF14 gene expression correlated with worse prognosis and GSEA showed that NFkB pathway was associated with the "High-TNFRSF14/dead-phenotype". In conclusion, the BTLA-TNFRSF14 immune modulation pathway seems to play a role in the pathobiology and prognosis of FL