In silico analyses reveal common cellular pathways affected by loss of heterozygosity (LOH) events in the lymphomagenesis of Non-Hodgkin’s lymphoma (NHL)
BMC Genomics volume 15, Article number: 390 (2014)
The analysis of cellular networks and pathways involved in oncogenesis has increased our knowledge about the pathogenic mechanisms that underlie tumour biology and has unmasked new molecular targets that may lead to the design of better anti-cancer therapies. Recently, using a high resolution loss of heterozygosity (LOH) analysis, we identified a number of potential tumour suppressor genes (TSGs) within common LOH regions across cases suffering from two of the most common forms of Non-Hodgkin’s lymphoma (NHL), Follicular Lymphoma (FL) and Diffuse Large B-cell Lymphoma (DLBCL). From these studies LOH of the protein tyrosine phosphatase receptor type J (PTPRJ) gene was identified as a common event in the lymphomagenesis of these B-cell lymphomas. The present study aimed to determine the cellular pathways affected by the inactivation of these TSGs including PTPRJ in FL and DLBCL tumourigenesis.
Pathway analytical approaches identified that candidate TSGs located within common LOH regions participate within cellular pathways, which may play a crucial role in FL and DLBCL lymphomagenesis (i.e., metabolic pathways). These analyses also identified genes within the interactome of PTPRJ (i.e. PTPN11 and B2M) that when inactivated in NHL may play an important role in tumourigenesis. We also detected genes that are differentially expressed in cases with and without LOH of PTPRJ, such as NFATC3 (nuclear factor of activated T-cells, cytoplasmic, calcineurin-dependent 3). Moreover, upregulation of the VEGF, MAPK and ERBB signalling pathways was also observed in NHL cases with LOH of PTPRJ, indicating that LOH-driving events causing inactivation of PTPRJ, apart from possibly inducing a constitutive activation of these pathways by reduction or abrogation of its dephosphorylation activity, may also induce upregulation of these pathways when inactivated. This finding implicates these pathways in the lymphomagenesis and progression of FL and DLBCL.
The evidence obtained in this research supports findings suggesting that FL and DLBCL share common pathogenic mechanisms. Also, it indicates that PTPRJ can play a crucial role in the pathogenesis of these B-cell tumours and suggests that activation of PTPRJ might be an interesting novel chemotherapeutic target for the treatment of these B-cell tumours.
The identification of altered pathways in tumour cells has provided a more holistic understanding of the pathogenic mechanisms that underlie the genesis, progression and chemoresponse of cancer. As a consequence, malfunction of a gene or a group of genes must be analysed as part of a complex network of components that are highly related to each other. However, this analysis is limited by factors such as the lack of analytic tools to convert and integrate, in a feasible and reliable manner, large amounts of data derived from different high-throughput genomic approaches into outputs with more biological meaning.
Non-Hodgkin’s lymphoma (NHL) represents a highly biological and clinical heterogeneous group of blood cancers [1, 2]. Two of the most common NHL subtypes are Diffuse Large B-cell Lymphoma (DLBCL), an aggressive lymphoma, and Follicular lymphoma (FL), a slow-growing type of lymphoma, which account for around 50% of NHL cases [1, 3]. The poor treatment outcomes obtained between 20-40% of cases with these NHL lymphomas has prompted studies aimed at the discovery of genes and pathways that could act as novel targets for a therapy that increase the survival rates of patients suffering from NHL [4–6]. However, the high genetic variability observed within NHL subtypes, has limited the understanding of the pathology of these NHL subtypes, as well as the discovery of new molecular targets for therapeutic development.
In a previous study, the integration of copy number (CNV) and gene expression profiling (GEP) data from DLBCL and FL cases, allowed us to identify common and disease-specific genetic alterations targeting known oncogenic pathways, such as the mitogen activated protein kinase (MAPK) and apoptosis signalling pathways, unmasking common pathogenic mechanisms underlying the malignant phenotype of these biologically and genetically distinct NHL subtypes . Likewise, in a recent study, using a high resolution loss of heterozygosity (LOH) analysis in FLs and DLBCLs, we also identified candidate tumour suppressor genes (TSGs) within common LOH regions across these NHL subtypes and implicated them in the lymphomagenesis of these B-cell lymphomas . In this study, we have performed pathway analysis of these candidate TSGs, in order to identify common cellular networks that might be altered by the inactivation of one or more TSGs in the lymphomagenesis of these B-cell lymphomas.
As part of our recent LOH studies, we also implicated PTPRJ (protein tyrosine phosphatase receptor type J) as a novel TSG in the tumourigenesis of FL and DLBCL, with LOH of PTPRJ identified as a common event in FL and DLBCL. LOH of PTPRJ was also confirmed by a decrease of heterozygosity of a microsatellite targeting PTPRJ loci in these NHL cases. In addition, FL cases with LOH exhibited a significant downregulation of PTPRJ . Several lines of evidence support the tumour suppressive role of PTPRJ, as this protein tyrosine phosphatase has been implicated in the oncogenesis of breast, lung, colorectal, thyroid and meningioma cancers [9–12]. Furthermore, PTPRJ regulates signalling pathways involved in cell growth, proliferation and angiogenesis, such as MAPK (ERK1/2), PLCG1, PI3K (p85), FLT3, B-cell receptor (BCR), PDGFRB and VEGFR2 signalling [13–19]. Nonetheless, the natural ligands of PTPRJ, Syndecan-2 (SDC2), a transmembrane heparan sulfate proteoglycan, and Thrombospondin-1 (THBS1), a homotrimeric glycoprotein, have been identified to induce cell adhesion and inhibit cell growth and angiogenesis, respectively [20, 21]. This anti-tumour activity of PTPRJ has also been demonstrated in In vitro experiments, using agonist peptides of PTPRJ and the oncogenic silencing of PTPRJ expression by microRNA-328 expression [22–24].
Since the role of PTPRJ in normal and malignant B-cell differentiation is poorly understood, we used pathway and genomic analyses to identify cellular pathways that may be altered by the inactivation of PTPRJ in order to provide a better understanding of the role of PTPRJ in the lymphomagenesis of FL and DLBCL. These analyses identified metabolic pathways as one of the most enriched and affected pathways resulting from the inactivation of candidate TSGs, indicating that these cellular pathways might play an important role in FL and DLBCL tumourigenesis. Additionally, inactivation of PTPRJ was shown to affect the expression of a number of genes and pathways that are regulated by PTPRJ through protein-protein interactions.
Candidate tumour suppressor genes (TSGs) targeted by LOH events participate within common cellular networks that may orchestrate NHL lymphomagenesis
In order to determine whether the candidate tumour suppressor genes (TSGs) located within common LOH regions across DLBCL and FL cases interact and participate within common cellular networks, a global interactome of a total of 262 genes affected by LOH events , was created using the VisANT (v. 4.06) platform (Figure 1). As a result of this analysis, the METABOLIC pathway (KEGG hsa-01100) was identified as the most enriched pathway by these candidate TSGs. This approach also revealed a high level of interaction amongst these candidate TSGs targeted by LOH events, indicating that NHL tumourigenesis might be orchestrated by the possible inactivation of these candidate TSGs, which could lead to the constitutive activation of oncogenic pathways. In this context, it is important to highlight that the genes MAPK6 (mitogen-activated protein kinase 6), PTPN11 (protein tyrosine phosphatase, non-receptor type 11) and ANXA7 (annexin A7) showed the highest number of links with other genes (170, 122 and 106 respectively). The high number of protein links of these proteins suggests that a possible inactivation of these genes may cause a major deregulation of patwhays by altering the function of a high number of genes. Thus, these genes are shown as the most interesting candidate TSGs in NHL lymphomagenesis.
In addition to the global interactome study, a gene set enrichment analysis (GSEA) was performed to classify the 262 candidate TSGs targeted by LOH events, thenceforth also referred in this study as LOH genes, into curated gene families and also to determine the pathways that are mostly enriched with these genes. This approach categorized these LOH genes into gene families (Additional file 1). For instance, FASLG was classified into the cytokines and growth factor family and together with PTPRJ into the cell differentiation markers family. EP300 was found in the translocated cancer gene, transcription factor and tumour suppressor gene families. Additionally, MAPK6, PTPN11 and NFATC3 were identified as members of the protein kinases, oncogenes and transcription factors families, respectively.
Furthermore, 215 out of 262 genes within LOH regions were identified and contrasted with collections from the molecular signature database (MsigDB, v3.1). Several cellular pathways were identified to be enriched with genes located within LOH regions (Table 1). The most enriched pathways for LOH genes were the, KEGG_ARGININE_AND_PROLINE METABOLISM, REACTOME_DOUBLE STRAND_BREAK_REPAIR and MICROTUBULE_ASSOCIATED_COMPLEX GO (gene ontology) gene sets. Interestingly, the TSGs TP53BP1 and B2M were identified as members of the REACTOME_DOUBLE_STRAND_BREAK_REPAIR and BIOCARTA_CTL PATHWAY pathways, respectively.
Significant enrichment of LOH genes was also observed in a collection of gene sets with chemical and genetic alterations (CGP), indicating that these candidate TSGs have been previously implicated in multiple forms of cancer (Additional file 2). In this collection, MARSON_BOUND_BY_FOXP3_STIMULATED, DIAZ_CHRONIC MEYLOGENOUS_LEUKEMIA_UP and GRAESSMANN_RESPONSE_TO_MC AND _DOXORUBICIN_UP were the one of most significant gene sets. To highlight, PTPRJ was included in the ROYLANCE_BREAST_CANCER_16Q_COPY_NUMBER_DN and TP53BP1 was listed in the PUJANA_BREAST_CANCER_LIT_INT_NETWORK gene sets.
Interactome of PTPRJ overlaps with the interactome of PTPN11 (protein tyrosine phosphatase, non-receptor type 11), a candidate TSG also found to be targeted by common LOH events in NHL
In a recent study, the protein tyrosine phosphatase receptor type J (PTPRJ) gene was implicated as a novel TSG in the lymphomagenesis of DLBCL and FL . In this previous study, using publically available gene expression profiling data from our cohort of NHL cases and subsequent validation by qualitative PCR (qPCR), a significant downregulation of PTPRJ expression in FL cases with LOH, indicating that the LOH-driving events targeting PTPRJ might have an effect on mechanisms regulating the expression of PTPRJ [7, 8]. This significant reduction in PTPRJ transcript abundance may result in haploinsufficiency of this TSG and thus, play a crucial role in the FL and DLBCL tumourigenesis. Based on this, we performed in silico pathway and gene expression analyses in order to provide more information about the unclear role and function of PTPRJ in normal and malignant B-cell development and to understand the effect of LOH of PTPRJ in the lymphomagenesis of FL and DLBCL from a more cellular perspective.
Initially, we studied the interactome of PTPRJ to determine whether candidate TSGs, affected by LOH events, interact directly or indirectly with PTPRJ and play a role together with PTPRJ, within common cellular networks. The analysis of the gene-gene interactions of PTPRJ revealed that none of the 41 genes, which are known to interact with PTPRJ, were previously identified among the genes located within the common LOH regions across FLs and DLBCLs (Figure 2) . However, 22 genes that interact indirectly (one level) with PTPRJ were found to be commonly inactivated in NHL cases. Among these candidate TSGs, B2M (Beta-2_microglobulin) and PTPN11 (protein tyrosine phosphatase, non-receptor type 11) were highlighted. B2M was found to interact indirectly with PTPRJ, through the PTPRJ-GRB2-B2M interaction. PTPN11, another protein tyrosine phosphatase (PTP), was identified to share 24 interacting genes with PTPRJ. Additionally, using previously published high resolution LOH data, LOH of PTPN11 was identified in 74% of NHL, 71% of FL and 76% of DLBCL cases .
The genes showing interaction with both PTPRJ and PTPN11 are: MAPK3 (mitogen-activated protein kinase 3), GRB2 (growth factor receptor-bound protein 2), ERBB2 (v-erb-b2 erythroblastic leukemia viral oncogene homolog 2), MET (met proto-oncogene (hepatocyte growth factor receptor), PDGFRB (platelet-derived growth factor receptor, beta polypeptide), CTNNB1 (catenin (cadherin-associated protein), beta 1, 88 kDa), PLCG1 (phospholipase C, gamma 1), CBL (Cas-Br-M (murine) ecotropic retroviral transforming sequence), KIT (v-kit Hardy-Zuckerman 4 feline sarcoma viral oncogene homolog), NTRK1 (neurotrophic tyrosine kinase, receptor, type 1), KDR (kinase insert domain receptor (a type III receptor tyrosine kinase)), PXN (paxillin), GAB1 (GRB2-associated binding protein 1), LEPR (leptin receptor), EPOR (erythropoietin receptor), GHR (growth hormone receptor), CXCR4 (chemokine (C-X-C motif) receptor 4), CSK (c-src tyrosine kinase), INSR (insulin receptor), JAK2 (Janus kinase 2), LCK (lymphocyte-specific protein tyrosine kinase), GRIN2B (glutamate receptor, ionotropic, N-methyl D-aspartate 2B), TEK (TEK tyrosine kinase, endothelial) and TIE1 (tyrosine kinase with immunoglobulin-like and EGF-like domains 1). Interestingly, these 24 genes were found to significantly enrich important signalling pathways in lymphocyte biology, which suggests that a double inactivation of PTPRJ and PTPN11 might have an aberrant effect on the function of these pathways (Table 2).
LOH of PTPRJinduces deregulation of genes and signalling pathways in NHL tumours
We also used gene expression profiling (GEP) data from NHL patients with known LOH status for PTPRJ , to detect genes that are significantly deregulated by LOH of PTPRJ. Using this approach, we identified differentially expressed genes between cases with and without LOH of PTPRJ. This result indicates that LOH of PTPRJ could have an effect on the expression pattern of some genes (Figure 3A), suggesting that some pathways might be differentially enriched between the two categories. In order to prove this hypothesis, a Gene Set Enrichment Analysis (GSEA) identified some cellular pathways with significant upregulation in cases with retention and some others, in cases with LOH (Figure 3B and 3C). The genes that were upregulated in the most significant gene sets are described in Table 3. The lists of the gene sets enriched in each LOH status (RET or LOH) are detailed in the Additional file 3 and Additional file 4. Nonetheless, it is important to highlight that despite all enrichment analyses exhibiting FDR values equal to 1, the results from these gene-set enrichments, described in Figure 3 and Table 3, were analysed based on the nominal significance (P value lower than 0.05) instead of the FDR values, due to the low population size, which could be driving the high FDR scores, and the high biological relevance of the gene-sets that were significantly upregulated in cases with LOH of PTPRJ.
NFATC3 is a commonly inactivated TSG in NHL cases and is significantly upregulated in cases with LOH of PTPRJ
In order to identify genes, previously found to be inactivated in our high resolution LOH study , and that are differentially expressed in the context of LOH of PTPRJ, we compared 262 candidate TSGs located within LOH regions against the 400 genes with differential expression between cases with retention and LOH of PTPRJ. This approach identified only the NFATC3 (nuclear factor of activated T-cells, cytoplasmic, calcineurin-dependent 3) gene, located at 16q22, to be commonly inactivated and differentially expressed in NHL cases with LOH of PTPRJ (Figure 4). Comparison of the mean fluorescence intensity values for NFATC3 between NHL cases with retention and LOH of PTPRJ revealed that NFATC3 was significantly upregulated in cases with LOH of PTPRJ (325 vs. 375; P = 0.042). Furthermore, It is important to mention that LOH of NFATC3 was found in 58% of NHL, 48% of FL and 67% of DLBCL cases .
The implementation of pathway analyses on candidate tumour suppressor genes (TSGs) found to be targeted by loss of heterozygosity in NHL, provided a holistic perspective of the shared cellular mechanisms underlying the lymphomagenesis of the commonest forms of NHL, FL and DLBCL. This analytical approach unmasked the pathways that are mostly altered by the genetic inactivation of TSGs potentially caused by LOH events in NHL and revealed common pathogenic mechanisms between FL and DLBCL. The analysis of candidate tumour suppressor genes, previously identified to be located within common LOH regions across FL and DLBCL cases , via interactome and pathway analytical approaches, identified the KEGG METABOLIC pathway (hsa-01100) and the ARGININE AND PROLINE METABOLISM pathway (hsa-00330) as the most enriched pathways in relation to these candidate TSGs. In addition, this finding suggests that the regulation of these pathways might be commonly altered in NHL. These candidate TSGs encode proteins (i.e., enzymes) that may directly or indirectly (i.e., transcriptional regulators) disrupt the activity of these pathways ,either by reduced or null expression of these genes with tumour suppressor activity or by the encoding of dysfunctional proteins, resulting from the LOH-induced inactivation of these TSGs.
Disruption of metabolic pathways such as energy production or biosynthesis of amino acids, nucleosides, etc. are one of the hallmarks in cancer, as the malignant cells demand a higher activity in energy metabolism and biosynthesis in order to sustain their constitutive and accelerated proliferation and their malignant performance [25, 26]. The purpose of this metabolic reprogramming of tumour cells is not only to increase the energy output in these cells but it is also to maximize the biosynthetic pathways of the malignant cells to turn them into more efficient and self-sustaining organisms [25–27]. For instance, reprogramming of the glutamine-proline-arginine metabolic circuit has been associated with cancer and found to be regulated by TP53 and MYC, in which proline acts as the regulatory axis of this circuit [27, 28]. Therefore, the significant enrichment of candidate TSGs in the metabolic pathways related to arginine and proline metabolism in our samples is evidence that this reprogramming of energetic and biosynthetic pathways is important in NHL tumour cells. However, further functional studies are required to corroborate the effect of LOH events targeting these TSGs in the function of this metabolic pathway and its role in the tumourigenesis of FL and DLBCL.
Using an annotated network database enriched with protein-protein interactions and associations, we were able to visualize the interactome of PTPRJ in relation to candidate TSGs located within LOH regions. This analysis determined that PTPRJ does not interact directly with any other possible LOH genes; however, this TSG indirectly interacts with 22 possible inactivated TSGs by LOH in NHL cases. This finding suggests that the malignant phenotype of tumoural B-cells is orchestrated by an uncontrolled activation of oncogenic signaling pathways, which may be caused by the inactivation of the dephosphorylation activity of PTPRJ together with the inactivation of other TSGs that might regulate direct targets of PTPRJ within these pathways. For instance, PTPRJ interacts directly with GRB2, a pivotal protein in signal transduction that activates the RAS-MAPK pathway, which is inactivated by the dephosphorylation activity of PTPRJ. GRB2 was also found to interact directly with B2M in an unclear manner . Interestingly, B2M was found to be inactivated in NHLs and its inactivation has been associated as a mechanism to evade the immune surveillance in DLBCL and other types of cancer [30, 31]. In addition, this finding also indicates that lymphomagenesis could be orchestrated by a network of TSGs that are inactivated in malignant B-cells to block tumour suppressor signalling and constitutively maintain pathogenic signals.
PTPN11, highly targeted by LOH events in NHL, was also identified in the PTPRJ interactome. PTPN11 encodes a protein tyrosine phosphatase (PTP) with two SH2 domains which acts as an intracellular signalling transducer of growth factors and cytokines receptors by regulation of the MAPK pathway . Furthermore, it has been demonstrated that this PTP plays a crucial role in hematopoiesis and that mutations in PTPN11 block the T-cell and B-cell development . Mutations in the PTPN11 gene have been detected in low frequencies in leukemia cases; however, alterations in this gene predispose patients with Noolan syndrome to several types of leukemia (mainly juvenile myelomonocytic leukemia JMML) and other types of cancer. [34–36]. Interestingly, PTPN11 and PTPRJ share around 24 gene interactions, which enriched important pathways in cancer such as ERBB signalling pathway, cytokine-cytokine interactions pathway, JAK/STAT signalling pathway and others. Thus, this finding indicates that the inactivation of these two PTPs may be a key factor in the malignant B-cell development of these NHL subtypes, as double inactivation of PTPRJ and PTPN11 may induce a constitutive activation of oncogenic pathways regulated by these PTPs.
Additionally, the gene NFATC3 was identified as being significantly upregulated in NHL cases with LOH of PTPRJ. Interestingly, LOH of NFATC3 was identified to be very common across FL and DLBCL cases and copy-neutral events were found to be the driving cause of LOH of NFATC3 in these cases, based on the previously described high resolution LOH approach in these NHL cases . NFATC3 is a Ca++ -dependent protein that regulates T-cell activation and migration, proliferation and angiogenesis by regulation of VEGF . Moreover, the abrogation of NFATC3 expression in a murine model suggested the implication of this TSG in T-cell lymphomagenesis; whereas, its implication in B-cell lymphomas has not been studied . In addition, NFATC3 was hypothesized as being part of a mechanism whereby intratumoural CD4+CD25+ T-cells (Treg cells) interact with activated CD4+ T-cells to suppress the anti-tumour activity of infiltrated CD4+ T-cells in B-cell NHL tumours and thus, induce immune tolerance to these tumours . Interestingly, these Treg cells also suppress the cytotoxic activity of CD8+ T-cells . These findings suggest that B-cell tumours may escape the immune surveillance through inactivation of B2M, the malignant B-cell-mediated recruitment of Treg cells and regulation of infiltrating CD4+, possibly through NFATC3, and CD8+ T-cells. Further studies are required to determine the role of NFATC3 and its inactivation in B-cell lymphomagenesis, the copy-neutral events driving LOH of NFATC3 and its relation with the inactivation of PTPRJ, which is known to inhibit the TCR-mediated T-cell activation by dephosphorylation of LAT and PLCG1 pathways [14, 41].
Although gene expression data was used to infer pathways of a protein tyrosine phosphatase, which lacks a DNA-binding regulatory motif, we were able to identify genes with different patterns of expression between cases with retention and LOH of PTPRJ, which suggests that the inactivation of PTPRJ might affect pathways whose ultimate goals are to switch on/off transcription factors and as a result of this, to induce an aberrant expression of these genes. In addition, using a Gene Set Enrichment Analysis (GSEA) of the differentially expressed genes between cases with LOH and retention; several pathways related to cancer and lymphocyte function were identified to be significantly upregulated in cases with LOH. This upregulation may be due to an indirect effect of the LOH of PTPRJ, which may cause a constitutive activation of pathways and may result in an abnormal expression of genes downstream of these pathways. Nonetheless, further functional studies (i.e., immunoblot) are required to validate these results and confirm the effect of LOH of PTPRJ in the expression of these oncogenic pathways.
For instance, the VEGF signalling pathway was found as the most significant upregulated pathway in cases with LOH of PTPRJ. This finding was expected as PTPRJ dephosphorylates VEGFR-2 and inhibits the VEGF-mediated cell proliferation, migration, angiogenesis and anti-apoptosis signalling . Malignant activation of the VEGF signalling pathway has been widely implicated in the tumour growth and lymphangiogenesis of NHLs and the expression of VEGF has been considered as a poor prognostic factor in these hematological malignancies and has been implicated in transformation of FL to DLBCL [42–44]. Furthermore, the ERBB and MAPK signalling pathway were also found to be upregulated in cases with LOH. These oncogenic pathways have been previously described to be dephosphorylated by PTPRJ [17, 45, 46]. Likewise, as PTPRJ also affects changes in Ca++ concentration [19, 47], it was expected that inactivation of PTPRJ upregulates Ca++ signalling pathways. Unexpectedly, signalling mediated by the Fc epsilon receptor I was found to be also upregulated in LOH cases; however, despite that this receptor for IgE is expressed in early stages of B-cells , its relation with PTPRJ is unclear.
On the other hand, the most enriched pathways in cases with retention of PTPRJ were DNA polymerase, TP53 signaling and proteasome pathways. It is likely that the upregulation of these pathways in cells with functional PTPRJ are the consequence of the activation of tumour suppressor activities, such as controlling cell proliferation and proapoptotic signalling through dephosphorylation of VEGFR, MAPK and PI3K pathways [16–18]. The upregulation of proteasome pathway in cases with retention of PTPRJ correlates with a previous study, which found that a rat homolog of PTPRJ controlled the proteasome-mediated degradation rate of its regulator (p27Ki) by activation of MAPK pathway . Furthermore, an analysis using freely available gene expression profiling data from FL and DLBCL cases expressing high and low levels of PTPRJ [50, 51], determined a low correlation of genes affected by PTPRJ transcript abundance between FL and DLBCL cases. This discrepancy may indicate differences in the LOH-driving events targeting this TSG in FLs and DLBCLs, which is supported by the significant downregulation of PTPRJ in FL cases, but not in DLBCL with LOH .
In conclusion, the use of pathway analytical approaches has provided more evidence supporting the presence of common pathogenic mechanisms underlying the lymphomagenesis of FL and DLBCL. Specifically, this work identified genes and pathways affected by LOH of PTPRJ. Furthermore, these findings suggest that PTPRJ plays a crucial role in the lymphomagenesis of FL and DLBCL as this TSG was found to induce aberrant expression of genes and PTPRJ-regulated pathways in NHL cases with LOH. Finally, it is important to highlight that the role of the most promising TSGs and pathways affected by LOH in FL and DLBCL tumourigenesis will need to be validated in further in vitro modeling experiments.
LOH pathway analyses
A total of 42 NHL cases, 21 DLBCL and 21 FL, were analysed using a high resolution LOH approach, which identified 46 common LOH regions across FL and DLBCL cases, harbouring 262 candidate tumour suppressor genes . The global interactome of genes affected by LOH events in NHL cases was built inserting these 262 candidate TSGs in the VisANT v. 4.06 (Integrative Visual Analysis Tool for Biological Networks and Pathways) platform to determine the level interaction amongst them. Genes without any interaction were eliminated from the analysis. In addition, a gene set enrichment analysis (GSEA) was performed on these 262 candidate TSGs, in order to investigate candidate TSGs with collections of curated gene sets, such as canonical pathways, Biocarta, KEGG, and Reactome, using the molecular signature database (MsigDB, v3.1). The significance level for this GSEA analysis was set at α < 0.05. Both analyses were based on the assumption that all the genes located within the inferred LOH regions, previously described , were being targeted by LOH events in the cohort of NHL cases.
Interactome of PTPRJ
In order to extend on recently published findings implicating PTPRJ as a novel TSG in the lymphomagenesis of FL and DLBCL  and to understand the unclear role of PTPRJ in B-cell tumourigenesis, we used bioinformatics tools to investigate the effect of LOH of PTPRJ on genes and pathways. Thus, to identify candidate TSGs affected by LOH events that interact directly or indirectly with PTPRJ and play a role together with PTPRJ, within common cellular networks, we extracted and analyzed the interactome of PTPRJ from the global interactome enriched with the 262 candidate TSGs , using the VisANT v. 4.06 platform (Integrative Visual Analysis Tool for Biological Networks and Pathways). In addition, the interactome of PTPRJ was enriched with log2-transformed data from an Illumina Sentrix Human-6 (v2.0) Expression Beadchip from NHL cases (19 FLs and 20 DLBCLs) with a known LOH status (15 RET and 14 LOH) [7, 8].
Selection of differentially expressed genes
To determine differentially expressed genes between cases with retention and LOH of PTPRJ, we used our previously described linear gene expression profiling (GEP) data from cases (19 FLs and 20 DLBCLs) with a known LOH status (15 RET and 14 LOH) [7, 8]. The ComparativeMarkerSelection module of the Genepattern platform  was employed for this selection and the analysis was carried out using the default options (2-sided T-Test and 10000 permutations). Results from this analysis were visualized using the ComparativeMarkerSelectionViewer module. A total of 400 (200 per status) differentially expressed genes were selected by the ExtractorComparativeMarkerResults module based on their score. Expression patterns were visualized using the HeatMapViewer module. A Venn diagram was used to identify common genes found to be differentially expressed from the comparison of cases with RET and LOH of PTPRJ and the 262 candidate TSGs identified within common LOH regions in NHL cases . Comparison of the expression levels of NFATC3 (nuclear factor of activated T-cells, cytoplasmic 3) between cases with retention and LOH of PTPRJ was carried out using an independent t-test (α < 0.05).
Gene set enrichment analysis (GSEA)
The identification of gene sets that were upregulated in cases with RET and LOH of PTPRJ was performed using the above mentioned GEP data from NHL cases with a known LOH status (15 RET and 14 LOH) for PTPRJ in the GSEA v 2.0 software. The collection of curated pathways from the KEGG database and the default settings (without collapsing GEP data) were used for this analysis. Significance level was set at an α < 0.05.
Swerdlow SH, Campo E, Harris NL, Jaffe ES, Pileri SA, Stein H, Thiele J, Vardiman JW: WHO Classification of Tumours of Haematopoietic and Lymphoid Tissues. 2008, Lyon: International Agency for Research on Cancer (IARC)
Turner JJ, Morton LM, Linet MS, Clarke CA, Kadin ME, Vajdic CM, Monnereau A, Maynadie M, Chiu BC, Marcos-Gragera R, Costantini AS, Cerhan JR, Weisenburger DD: InterLymph hierarchical classification of lymphoid neoplasms for epidemiologic research based on the WHO classification (2008): update and future directions. Blood. 2010, 116 (20): e90-98. 10.1182/blood-2010-06-289561.
Jayasekara H, Karahalios A, Juneja S, Thursfield V, Farrugia H, English DR, Giles GG: Incidence and survival of lymphohematopoietic neoplasms according to the World Health Organization classification: a population-based study from the Victorian Cancer Registry in Australia. Leuk Lymphoma. 2010, 51 (3): 456-468. 10.3109/10428190903552104.
Federico M, Luminari S, Dondi A, Tucci A, Vitolo U, Rigacci L, Di Raimondo F, Carella AM, Pulsoni A, Merli F, Arcaini L, Angrilli F, Stelitano C, Gaidano G, Dell'Olio M, Marcheselli L, Franco V, Galimberti S, Sacchi S, Brugiatelli M: R-CVP Versus R-CHOP Versus R-FM for the initial treatment of patients with advanced-stage follicular lymphoma: results of the FOLL05 trial conducted by the Fondazione Italiana Linfomi. J Clin Oncol. 2013, 31: 1506-1513. 10.1200/JCO.2012.45.0866.
Friedberg JW, Fisher RI: Diffuse large B-cell lymphoma. Hematol Oncol Clin North Am. 2008, 22 (5): 941-952. 10.1016/j.hoc.2008.07.002. ix
Feugier P, Van Hoof A, Sebban C, Solal-Celigny P, Bouabdallah R, Ferme C, Christian B, Lepage E, Tilly H, Morschhauser F, Gaulard P, Salles G, Bosly A, Gisselbrecht C, Reyes F, Coiffier B: Long-term results of the R-CHOP study in the treatment of elderly patients with diffuse large B-cell lymphoma: a study by the Groupe d’Etude des Lymphomes de l’Adulte. J Clin Oncol. 2005, 23 (18): 4117-4126. 10.1200/JCO.2005.09.131.
Green MR, Aya-Bonilla C, Gandhi MK, Lea RA, Wellwood J, Wood P, Marlton P, Griffiths LR: Integrative genomic profiling reveals conserved genetic mechanisms for tumorigenesis in common entities of non-Hodgkin’s lymphoma. Genes Chromosomes Cancer. 2011, 50 (5): 313-326. 10.1002/gcc.20856.
Aya-Bonilla C, Green MR, Camilleri E, Benton M, Keane C, Marlton P, Lea R, Gandhi MK, Griffiths LR: High-resolution loss of heterozygosity screening implicates PTPRJ as a potential tumor suppressor gene that affects susceptibility to non-hodgkin’s lymphoma. Genes Chromosom Cancer. 2013, 52 (5): 467-479. 10.1002/gcc.22044.
Iuliano R, Le Pera I, Cristofaro C, Baudi F, Arturi F, Pallante P, Martelli ML, Trapasso F, Chiariotti L, Fusco A: The tyrosine phosphatase PTPRJ/DEP-1 genotype affects thyroid carcinogenesis. Oncogene. 2004, 23 (52): 8432-8438. 10.1038/sj.onc.1207766.
Petermann A, Haase D, Wetzel A, Balavenkatraman KK, Tenev T, Guhrs KH, Friedrich S, Nakamura M, Mawrin C, Bohmer FD: Loss of the protein-tyrosine phosphatase DEP-1/PTPRJ drives meningioma cell motility. Brain Pathol. 2011, 21 (4): 405-418. 10.1111/j.1750-3639.2010.00464.x.
Ruivenkamp C, Hermsen M, Postma C, Klous A, Baak J, Meijer G, Demant P: LOH of PTPRJ occurs early in colorectal cancer and is associated with chromosomal loss of 18q12-21. Oncogene. 2003, 22 (22): 3472-3474. 10.1038/sj.onc.1206246.
Ruivenkamp CA, van Wezel T, Zanon C, Stassen AP, Vlcek C, Csikos T, Klous AM, Tripodis N, Perrakis A, Boerrigter L, Groot PC, Lindeman J, Mooi WJ, Meijjer GA, Scholten G, Dauwerse H, Paces V, van Zandwijk N, van Ommen GJ, Demant P: Ptprj is a candidate for the mouse colon-cancer susceptibility locus Scc1 and is frequently deleted in human cancers. Nat Genet. 2002, 31 (3): 295-300. 10.1038/ng903.
Arora D, Stopp S, Bohmer SA, Schons J, Godfrey R, Masson K, Razumovskaya E, Ronnstrand L, Tanzer S, Bauer R, Bohmer FD, Muller JP: Protein-tyrosine phosphatase DEP-1 controls receptor tyrosine kinase FLT3 signaling. J Biol Chem. 2011, 286 (13): 10918-10929. 10.1074/jbc.M110.205021.
Baker JE, Majeti R, Tangye SG, Weiss A: Protein tyrosine phosphatase CD148-mediated inhibition of T-cell receptor signal transduction is associated with reduced LAT and phospholipase Cgamma1 phosphorylation. Mol Cell Biol. 2001, 21 (7): 2393-2403. 10.1128/MCB.21.7.2393-2403.2001.
Kovalenko M, Denner K, Sandstrom J, Persson C, Gross S, Jandt E, Vilella R, Bohmer F, Ostman A: Site-selective dephosphorylation of the platelet-derived growth factor beta-receptor by the receptor-like protein-tyrosine phosphatase DEP-1. J Biol Chem. 2000, 275 (21): 16219-16226. 10.1074/jbc.275.21.16219.
Lampugnani MG, Zanetti A, Corada M, Takahashi T, Balconi G, Breviario F, Orsenigo F, Cattelino A, Kemler R, Daniel TO, Dejana E: Contact inhibition of VEGF-induced proliferation requires vascular endothelial cadherin, beta-catenin, and the phosphatase DEP-1/CD148. J Cell Biol. 2003, 161 (4): 793-804. 10.1083/jcb.200209019.
Sacco F, Tinti M, Palma A, Ferrari E, Nardozza AP, van Huijsduijnen Hooft R, Takahashi T, Castagnoli L, Cesareni G: Tumor suppressor density-enhanced phosphatase-1 (DEP-1) inhibits the RAS pathway by direct dephosphorylation of ERK1/2 kinases. J Biol Chem. 2009, 284 (33): 22048-22058. 10.1074/jbc.M109.002758.
Tsuboi N, Utsunomiya T, Roberts RL, Ito H, Takahashi K, Noda M, Takahashi T: The tyrosine phosphatase CD148 interacts with the p85 regulatory subunit of phosphoinositide 3-kinase. Biochem J. 2008, 413 (1): 193-200. 10.1042/BJ20071317.
Zhu JW, Brdicka T, Katsumoto TR, Lin J, Weiss A: Structurally distinct phosphatases CD45 and CD148 both regulate B cell and macrophage immunoreceptor signaling. Immunity. 2008, 28 (2): 183-196. 10.1016/j.immuni.2007.11.024.
Takahashi K, Mernaugh RL, Friedman DB, Weller R, Tsuboi N, Yamashita H, Quaranta V, Takahashi T: Thrombospondin-1 acts as a ligand for CD148 tyrosine phosphatase. Proc Natl Acad Sci. 2012, 109 (6): 1985-1990. 10.1073/pnas.1106171109.
Whiteford JR, Xian X, Chaussade C, Vanhaesebroeck B, Nourshargh S, Couchman JR: Syndecan-2 is a novel ligand for the protein tyrosine phosphatase receptor CD148. Mol Biol Cell. 2011, 22 (19): 3609-3624. 10.1091/mbc.E11-02-0099.
Ortuso F, Paduano F, Carotenuto A, Gomez-Monterrey I, Bilotta A, Gaudio E, Sala M, Artese A, Vernieri E, Dattilo V, Iuliano R, Brancaccio D, Bertamino A, Musella S, Alcaro S, Grieco P, Perrotti N, Croce CM, Novellino E, Fusco A, Campiglia P, Trapasso F: Discovery of PTPRJ agonist peptides that effectively inhibit in vitro cancer cell proliferation and tube formation. ACS Chem Biol. 2013, 8 (7): 1497-1506. 10.1021/cb3007192.
Paduano F, Dattilo V, Narciso D, Bilotta A, Gaudio E, Menniti M, Agosti V, Palmieri C, Perrotti N, Fusco A, Trapasso F, Iuliano R: Protein tyrosine phosphatase PTPRJ is negatively regulated by microRNA-328. FEBS J. 2013, 280 (2): 401-412. 10.1111/j.1742-4658.2012.08624.x.
Paduano F, Ortuso F, Campiglia P, Raso C, Iaccino E, Gaspari M, Gaudio E, Mangone G, Carotenuto A, Bilotta A, Narciso D, Palmieri C, Agosti V, Artese A, Gomez-Monterrey I, Sala M, Cuda G, Iuliano R, Perrotti N, Scala G, Viglietto G, Alcaro S, Croce CM, Novellino E, Fusco A, Trapasso F: Isolation and functional characterization of peptide agonists of PTPRJ, a tyrosine phosphatase receptor endowed with tumor suppressor activity. ACS Chem Biol. 2012, 7 (10): 1666-1676. 10.1021/cb300281t.
Hanahan D, Weinberg RA: Hallmarks of cancer: the next generation. Cell. 2011, 144 (5): 646-674. 10.1016/j.cell.2011.02.013.
Vander Heiden MG, Cantley LC, Thompson CB: Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science. 2009, 324 (5930): 1029-1033. 10.1126/science.1160809.
Phang JM, Liu W, Hancock C, Christian KJ: The proline regulatory axis and cancer. Front Oncol. 2012, 2: 60-
Catchpole G, Platzer A, Weikert C, Kempkensteffen C, Johannsen M, Krause H, Jung K, Miller K, Willmitzer L, Selbig J, Weikert S: Metabolic profiling reveals key metabolic features of renal cell carcinoma. J Cell Mol Med. 2011, 15 (1): 109-118. 10.1111/j.1582-4934.2009.00939.x.
Wang J, Yuan Y, Zhou Y, Guo L, Zhang L, Kuai X, Deng B, Pan Z, Li D, He F: Protein interaction data set highlighted with human Ras-MAPK/PI3K signaling pathways. J Proteome Res. 2008, 7 (9): 3879-3889. 10.1021/pr8001645.
Challa-Malladi M, Lieu YK, Califano O, Holmes AB, Bhagat G, Murty VV, Dominguez-Sola D, Pasqualucci L, Dalla-Favera R: Combined genetic inactivation of beta2-Microglobulin and CD58 reveals frequent escape from immune recognition in diffuse large B cell lymphoma. Cancer Cell. 2011, 20 (6): 728-740. 10.1016/j.ccr.2011.11.006.
Hill DM, Kasliwal T, Schwarz E, Hebert AM, Chen T, Gubina E, Zhang L, Kozlowski S: A dominant negative mutant beta 2-microglobulin blocks the extracellular folding of a major histocompatibility complex class I heavy chain. J Biol Chem. 2003, 278 (8): 5630-5638. 10.1074/jbc.M208381200.
Cunnick JM, Meng S, Ren Y, Desponts C, Wang HG, Djeu JY, Wu J: Regulation of the mitogen-activated protein kinase signaling pathway by SHP2. J Biol Chem. 2002, 277 (11): 9498-9504. 10.1074/jbc.M110547200.
Qu CK, Nguyen S, Chen J, Feng GS: Requirement of Shp-2 tyrosine phosphatase in lymphoid and hematopoietic cell development. Blood. 2001, 97 (4): 911-914. 10.1182/blood.V97.4.911.
Hugues L, Cave H, Philippe N, Pereira S, Fenaux P, Preudhomme C: Mutations of PTPN11 are rare in adult myeloid malignancies. Haematologica. 2005, 90 (6): 853-854.
Jongmans MC, van der Burgt I, Hoogerbrugge PM, Noordam K, Yntema HG, Nillesen WM, Kuiper RP, Ligtenberg MJ, van Kessel AG, van Krieken JH, Kiemeney LA, Hoogerbrugge N: Cancer risk in patients with Noonan syndrome carrying a PTPN11 mutation. Eur J Hum Genet. 2011, 19 (8): 870-874. 10.1038/ejhg.2011.37.
Tartaglia M, Martinelli S, Cazzaniga G, Cordeddu V, Iavarone I, Spinelli M, Palmi C, Carta C, Pession A, Arico M, Masera G, Basso G, Sorcini M, Gelb BD, Biondi A: Genetic evidence for lineage-related and differentiation stage-related contribution of somatic PTPN11 mutations to leukemogenesis in childhood acute leukemia. Blood. 2004, 104 (2): 307-313. 10.1182/blood-2003-11-3876.
Urso K, Alfranca A, Martínez-Martínez S, Escolano A, Ortega I, Rodríguez A, Redondo JM: NFATc3 regulates the transcription of genes involved in T-cell activation and angiogenesis. Blood. 2011, 118 (3): 795-803. 10.1182/blood-2010-12-322701.
Glud SZ, Sorensen AB, Andrulis M, Wang B, Kondo E, Jessen R, Krenacs L, Stelkovics E, Wabl M, Serfling E, Palmetshofer A, Pedersen FS: A tumor-suppressor function for NFATc3 in T-cell lymphomagenesis by murine leukemia virus. Blood. 2005, 106 (10): 3546-3552. 10.1182/blood-2005-02-0493.
Yang ZZ, Novak AJ, Stenson MJ, Witzig TE, Ansell SM: Intratumoral CD4+CD25+ regulatory T-cell-mediated suppression of infiltrating CD4+ T cells in B-cell non-Hodgkin lymphoma. Blood. 2006, 107 (9): 3639-3646. 10.1182/blood-2005-08-3376.
Yang ZZ, Novak AJ, Ziesmer SC, Witzig TE, Ansell SM: Attenuation of CD8(+) T-cell function by CD4(+)CD25(+) regulatory T cells in B-cell non-Hodgkin’s lymphoma. Cancer Res. 2006, 66 (20): 10145-10152. 10.1158/0008-5472.CAN-06-1822.
Tangye SG, Phillips JH, Lanier LL, de Vries JE, Aversa G: CD148: a receptor-type protein tyrosine phosphatase involved in the regulation of human T cell activation. J Immunol. 1998, 161 (7): 3249-3255.
Paydas S, Seydaoglu G, Ergin M, Erdogan S, Yavuz S: The prognostic significance of VEGF-C and VEGF-A in non-Hodgkin lymphomas. Leuk Lymphoma. 2009, 50 (3): 366-373. 10.1080/10428190802706665.
Ruan J, Hajjar K, Rafii S, Leonard JP: Angiogenesis and antiangiogenic therapy in non-Hodgkin’s lymphoma. Ann Oncol. 2009, 20 (3): 413-424.
Shipp MA, Ross KN, Tamayo P, Weng AP, Kutok JL, Aguiar RC, Gaasenbeek M, Angelo M, Reich M, Pinkus GS, Ray TS, Koval MA, Last KW, Norton A, Lister TA, Mesirov J, Neuberg DS, Lander ES, Aster JC, Golub TR: Diffuse large B-cell lymphoma outcome prediction by gene-expression profiling and supervised machine learning. Nat Med. 2002, 8 (1): 68-74. 10.1038/nm0102-68.
Takahashi T, Takahashi K, Mernaugh RL, Tsuboi N, Liu H, Daniel TO: A monoclonal antibody against CD148, a receptor-like tyrosine phosphatase, inhibits endothelial-cell growth and angiogenesis. Blood. 2006, 108 (4): 1234-1242. 10.1182/blood-2005-10-4296.
Tarcic G, Boguslavsky SK, Wakim J, Kiuchi T, Liu A, Reinitz F, Nathanson D, Takahashi T, Mischel PS, Ng T, Yarden Y: An unbiased screen identifies DEP-1 tumor suppressor as a phosphatase controlling EGFR endocytosis. Curr Biol. 2009, 19 (21): 1788-1798. 10.1016/j.cub.2009.09.048.
de la Fuente-Garcia MA, Nicolas JM, Freed JH, Palou E, Thomas AP, Vilella R, Vives J, Gaya A: CD148 is a membrane protein tyrosine phosphatase present in all hematopoietic lineages and is involved in signal transduction on lymphocytes. Blood. 1998, 91 (8): 2800-2809.
Toba K, Hanawa H, Sakaue M, Yoshida K, Itoh H, Tsuchiyama J, Maruyama S, Narita M, Takahashi M, Watanabe K, Aizawa Y: Fc epsilon RI and CD22 mRNA are expressed in early B-lineage and myeloid leukemia cell lines. Leuk Res. 2003, 27 (2): 173-182. 10.1016/S0145-2126(02)00086-3.
Trapasso F, Iuliano R, Boccia A, Stella A, Visconti R, Bruni P, Baldassarre G, Santoro M, Viglietto G, Fusco A: Rat protein tyrosine phosphatase eta suppresses the neoplastic phenotype of retrovirally transformed thyroid cells through the stabilization of p27(Kip1). Mol Cell Biol. 2000, 20 (24): 9236-9246. 10.1128/MCB.20.24.9236-9246.2000.
Dave SS, Wright G, Tan B, Rosenwald A, Gascoyne RD, Chan WC, Fisher RI, Braziel RM, Rimsza LM, Grogan TM, Miller TP, LeBlanc M, Greiner TC, Weisenburger DD, Lynch JC, Vose J, Armitage JO, Smeland EB, Kvaloy S, Holte H, Delabie J, Connors JM, Lansdorp PM, Ouyang Q, Lister TA, Davies AJ, Norton AJ, Muller-Hermelink HK, Ott G, Campo E, et al: Prediction of survival in follicular lymphoma based on molecular features of tumor-infiltrating immune cells. N Engl J Med. 2004, 351 (21): 2159-2169. 10.1056/NEJMoa041869.
Lenz G, Wright G, Dave SS, Xiao W, Powell J, Zhao H, Xu W, Tan B, Goldschmidt N, Iqbal J, Vose J, Bast M, Fu K, Weisenburger DD, Greiner TC, Armitage JO, Kyle A, May L, Gascoyne RD, Connors JM, Troen G, Holte H, Kvaloy S, Dierickx D, Verhoef G, Delabie J, Smeland EB, Jares P, Martinez A, Lopez-Guillermo A, et al: Stromal gene signatures in large-B-cell lymphomas. N Engl J Med. 2008, 359 (22): 2313-2323. 10.1056/NEJMoa0802885.
Gould J, Getz G, Monti S, Reich M, Mesirov JP: Comparative gene marker selection suite. Bioinformatics. 2006, 22 (15): 1924-1925. 10.1093/bioinformatics/btl196.
We thank all participating institutions and patients for supporting our research. This study was supported by funding from the Cancer Council of Queensland, the Genomics Lymphoma Project and the Herbert family. The Australasian Leukaemia and Lymphoma Group (ALLG) Tissue Bank is supported by the Leukaemia Foundation of Australia and the National Health and Medical Research Council.
The authors declare that they have no competing interests.
CAB designed, performed experiments, data analyses and drafted the manuscript. EC analyzed data. RL and LMH analyzed data. MKG designed experiments and analyzed data. LRG coordinated the research and designed analyses, analyzed data and wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: List of gene sets with chemical and genetic perturbations (CGP) enriched with genes within LOH regions.(DOC 46 KB)
About this article
Cite this article
Aya-Bonilla, C., Camilleri, E., Haupt, L.M. et al. In silico analyses reveal common cellular pathways affected by loss of heterozygosity (LOH) events in the lymphomagenesis of Non-Hodgkin’s lymphoma (NHL). BMC Genomics 15, 390 (2014). https://doi.org/10.1186/1471-2164-15-390