BACKGROUND
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  represents a highly biological and clinical heterogeneous group of blood cancers  <cit> .
two of the most common nhl subtypes are diffuse large b-cell lymphoma , an aggressive lymphoma, and follicular lymphoma , a slow-growing type of lymphoma, which account for around 50% of nhl cases  <cit> .
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 .
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  and gene expression profiling  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  and apoptosis signalling pathways, unmasking common pathogenic mechanisms underlying the malignant phenotype of these biologically and genetically distinct nhl subtypes  <cit> .
likewise, in a recent study, using a high resolution loss of heterozygosity  analysis in fls and dlbcls, we also identified candidate tumour suppressor genes  within common loh regions across these nhl subtypes and implicated them in the lymphomagenesis of these b-cell lymphomas  <cit> .
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  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  <cit> .
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 .
furthermore, ptprj regulates signalling pathways involved in cell growth, proliferation and angiogenesis, such as mapk , plcg <dig>  pi3k , flt <dig>  b-cell receptor , pdgfrb and vegfr <dig> signalling .
nonetheless, the natural ligands of ptprj, syndecan- <dig> , a transmembrane heparan sulfate proteoglycan, and thrombospondin- <dig> , a homotrimeric glycoprotein, have been identified to induce cell adhesion and inhibit cell growth and angiogenesis, respectively  <cit> .
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- <dig> expression .
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.
RESULTS
candidate tumour suppressor genes  targeted by loh events participate within common cellular networks that may orchestrate nhl lymphomagenesis
in order to determine whether the candidate tumour suppressor genes  located within common loh regions across dlbcl and fl cases interact and participate within common cellular networks, a global interactome of a total of  <dig> genes affected by loh events  <cit> , was created using the visant  platform .
as a result of this analysis, the metabolic pathway  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 mapk <dig> , ptpn <dig>  and anxa <dig>  showed the highest number of links with other genes .
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.figure  <dig>
global interactome of genes commonly affected by loh across nhl patients.
an initial list containing  <dig> genes was used; however, only those with more than one interaction are shown.
a total of  <dig> nodes belonging to  <dig> pathways were mapped, which indicates the high level of interaction among the genes targeted by loh-driving events in nhl tumours.
the genes with labels correspond to the questioned loh genes.
in addition to the global interactome study, a gene set enrichment analysis  was performed to classify the  <dig> 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 .
for instance, faslg was classified into the cytokines and growth factor family and together with ptprj into the cell differentiation markers family.
ep <dig> was found in the translocated cancer gene, transcription factor and tumour suppressor gene families.
additionally, mapk <dig>  ptpn <dig> and nfatc <dig> were identified as members of the protein kinases, oncogenes and transcription factors families, respectively.
furthermore,  <dig> out of  <dig> genes within loh regions were identified and contrasted with collections from the molecular signature database .
several cellular pathways were identified to be enriched with genes located within loh regions .
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 sets.
interestingly, the tsgs tp53bp <dig> and b2m were identified as members of the reactome_double_strand_break_repair and biocarta_ctl pathway pathways, respectively.table  <dig>
list of cellular pathways enriched with genes within loh regions
*k/k is the ratio between the number of genes in overlap  and the number of genes in gene set .
this table combines the results from the collection of gene sets from canonical pathways, biocarta, kegg, reactome, and go  gene sets.
significant enrichment of loh genes was also observed in a collection of gene sets with chemical and genetic alterations , indicating that these candidate tsgs have been previously implicated in multiple forms of cancer .
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 tp53bp <dig> was listed in the pujana_breast_cancer_lit_int_network gene sets.
interactome of ptprj overlaps with the interactome of ptpn <dig> , 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  gene was implicated as a novel tsg in the lymphomagenesis of dlbcl and fl  <cit> .
in this previous study, using publically available gene expression profiling data from our cohort of nhl cases and subsequent validation by qualitative pcr , 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  <cit> .
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  <dig> genes, which are known to interact with ptprj, were previously identified among the genes located within the common loh regions across fls and dlbcls   <cit> .
however,  <dig> genes that interact indirectly  with ptprj were found to be commonly inactivated in nhl cases.
among these candidate tsgs, b2m  and ptpn <dig>  were highlighted.
b2m was found to interact indirectly with ptprj, through the ptprj-grb2-b2m interaction.
ptpn <dig>  another protein tyrosine phosphatase , was identified to share  <dig> interacting genes with ptprj.
additionally, using previously published high resolution loh data, loh of ptpn <dig> was identified in 74% of nhl, 71% of fl and 76% of dlbcl cases  <cit> .figure  <dig>
interactome of ptprj indicates that the inactivation of
ptprj
may affect cellular networks and that these networks are not only altered by inactivation of
ptprj
.
ptprj gene network showing a) all the known genes that interact with ptprj and b) only those genes that were also found to be inactivated in nhl cases.
two levels of interaction were used to filter the interactions.
interestingly, we found that the gene ptpn <dig>  which was also found inactivated in nhl cases, shares some target genes with ptprj; moreover, b2m has an indirect interaction with ptprj .
the genes showing interaction with both ptprj and ptpn <dig> are: mapk <dig> , grb <dig> , erbb <dig> , met , pdgfrb , ctnnb <dig> , beta  <dig>  88 kda), plcg <dig> , cbl  ecotropic retroviral transforming sequence), kit , ntrk <dig> , kdr ), pxn , gab <dig> , lepr , epor , ghr , cxcr <dig>  receptor 4), csk , insr , jak <dig> , lck , grin2b , tek  and tie <dig> .
interestingly, these  <dig> genes were found to significantly enrich important signalling pathways in lymphocyte biology, which suggests that a double inactivation of ptprj and ptpn <dig> might have an aberrant effect on the function of these pathways .table  <dig>
list of gene sets commonly regulated by ptprj and ptpn <dig> based on gene expression
*k/k is the ratio between the number of genes in overlap  and the number of genes in gene set .
loh of ptprj induces deregulation of genes and signalling pathways in nhl tumours
we also used gene expression profiling  data from nhl patients with known loh status for ptprj  <cit> , 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 , suggesting that some pathways might be differentially enriched between the two categories.
in order to prove this hypothesis, a gene set enrichment analysis  identified some cellular pathways with significant upregulation in cases with retention and some others, in cases with loh .
the genes that were upregulated in the most significant gene sets are described in table  <dig>  the lists of the gene sets enriched in each loh status  are detailed in the additional file  <dig> and additional file  <dig>  nonetheless, it is important to highlight that despite all enrichment analyses exhibiting fdr values equal to  <dig>  the results from these gene-set enrichments, described in figure  <dig> and table  <dig>  were analysed based on the nominal significance  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.figure  <dig>
genes and pathways affected by loh of
ptprj
.
a) heatmap of genes differentially expressed between cases with retention  and loh calls.
b) dna polymerase , tp <dig> signaling pathway  and proteasome  were identified as the most enriched gene sets in cases with retention of ptprj.
c) in cases with loh, the gene sets vegf signaling pathway , melanoma , erbb signaling pathway , calcium signaling pathway , fc epsilon ri signaling pathway  and mapk signaling pathway  were the most enriched.
despite an fdr  of  <dig> for all the enrichments, the results from this gene-set enrichment analysis were analysed based on the p values instead of the fdr scores.
the high fdr values obtained in this analysis may be due to the low population size.
the lists of all gene sets from this analysis are provided in additional file  <dig> and additional file  <dig>
list of genes enriched within the gene sets differentially expressed between cases with retention and loh of ptprj
polq,prim <dig> pold <dig> pole <dig> polk,pole <dig> pole <dig> pols,rev <dig> pola <dig> rfc <dig> pold <dig> polb,poli,pole
cdc <dig> pten,tnfrsf10b,chek <dig> pmaip <dig> cdk <dig> cdk <dig> sfn,perp,p53aip <dig> fas,ccnb <dig> serpinb <dig> casp <dig> ccnb <dig> ei <dig> ccnd <dig> tp53i <dig> ccne <dig> rfwd <dig> ccnb <dig> mdm2
psmb <dig> psma <dig> psma <dig> psma <dig> psmb <dig> psmd <dig> psmd <dig> psmd <dig> psmc <dig> psma <dig> psma <dig> psmd <dig> psma <dig> psma6
mapkapk <dig> pla2g <dig> pxn,pla2g <dig> pla2g2f,map2k <dig> map2k <dig> raf <dig> prkcg,pik3r <dig> pla2g1b,mapk <dig> mapk <dig> sphk <dig> akt <dig> pla2g2a,bad,rac <dig> pla2g <dig> pik3cd,nfatc <dig> prkca,pik3cg
fgf <dig> ccnd <dig> met,fgf <dig> fgf <dig> map2k <dig> map2k <dig> igf <dig> raf <dig> pik3r <dig> mapk <dig> cdk <dig> akt <dig> fgf <dig> bad,fgf <dig> araf,fgf <dig> fgfr <dig> pik3cd,pdgfra,pdgfd,fgf <dig> pik3cg
camk2g,shc <dig> map2k <dig> map2k <dig> raf <dig> stat5b,prkcg,pak <dig> nrg <dig> rps6kb <dig> camk2b,pik3r <dig> pak <dig> nrg <dig> mapk <dig> erbb <dig> akt <dig> tgfa,bad,araf,cbl,erbb <dig> pik3cd,prkca,pik3cg, pak3
gna <dig> htr <dig> cckbr,chp,adra1d,itpr <dig> ryr <dig> slc25a <dig> phka <dig> nos <dig> phka <dig> slc8a <dig> cacna1c,atp2b <dig> mylk,agtr <dig> camk2g,chrna <dig> camk <dig> mylk <dig> prkcg,camk2b,grm <dig> bdkrb <dig> adcy <dig> ryr <dig> erbb <dig> sphk <dig> trhr,calm <dig> prkacg,adora2a,pln,gna <dig> tbxa2r,ptafr,erbb <dig> adcy <dig> p2rxl <dig> tacr <dig> cacna1d,pdgfra,hrh <dig> ntsr <dig> plcb <dig> htr2b,prkca, cacna1i,slc25a <dig> ptger <dig> itpkb
pla2g2f,vav <dig> map2k <dig> map2k <dig> raf <dig> pik3r <dig> pla2g1b,mapk <dig> mapk <dig> akt <dig> pla2g2a,rac <dig> il <dig> pla2g <dig> pik3cd,il <dig>  pdk <dig> prkca,pik3cg
dusp <dig> map2k <dig> fgf <dig> cacna2d <dig> mef2c,fgf <dig> rps6ka <dig> cacna2d <dig> tgfbr <dig> taok <dig> cacna2d <dig> rac <dig> nfkb <dig> chp,flna,srf,fgf <dig> daxx,rps6ka <dig> rras,mapkapk <dig> pla2g <dig> ntrk <dig> cacna1c,pla2g <dig> fgf <dig> fgf <dig> rap1a,flnc,map3k <dig> pla2g2f,cacnb <dig> map2k <dig> gadd45a,rps6ka <dig> nf <dig> map2k <dig> rps6ka <dig> raf <dig> prkcg,pak <dig> taok <dig> mos,mapk8ip <dig> pla2g1b, mapk <dig> mapk <dig> rras <dig> prkacg,gadd45g,akt <dig> ntf <dig> fgf <dig>  pla2g2a,flnb,rac <dig> fgf <dig> fgf <dig> rasgrf <dig> fgfr <dig> pla2g <dig> map3k8nlk ptpn <dig> cacna1d il1r <dig> pdgfra arrb <dig> map3k <dig> rps6ka <dig> ppp5c fgf <dig> gna <dig> cacng <dig> map3k <dig> dusp <dig> prkca map2k <dig> cacna1i mapk8ip2
vegfb,rictor,pgf,akt <dig> prkaa <dig> rps6ka <dig> rps6ka <dig> rps6ka <dig> igf <dig> rps6kb <dig> pik3r <dig> eif4b,mapk <dig> vegfc,akt <dig> figf, pik3cd,pik3cg
pim <dig> nfkb <dig> stat <dig> ccnd <dig> flt <dig> map2k <dig> cebpa,map2k <dig> raf <dig> rps6kb <dig> jup,mapk <dig> tcf <dig> akt <dig> bad,araf,pik3cd, pik3cg,kit
nfatc <dig> 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  <cit> , and that are differentially expressed in the context of loh of ptprj, we compared  <dig> candidate tsgs located within loh regions against the  <dig> genes with differential expression between cases with retention and loh of ptprj.
this approach identified only the nfatc <dig>  gene, located at 16q <dig>  to be commonly inactivated and differentially expressed in nhl cases with loh of ptprj .
comparison of the mean fluorescence intensity values for nfatc <dig> between nhl cases with retention and loh of ptprj revealed that nfatc <dig> was significantly upregulated in cases with loh of ptprj .
furthermore, it is important to mention that loh of nfatc <dig> was found in 58% of nhl, 48% of fl and 67% of dlbcl cases  <cit> .figure  <dig>
identification of genes located within the common loh regions in nhl cases whose expression is affected by the loh status of
ptprj.
a) comparison of the differentially expressed genes between nhl cases with retention and loh of ptprj  and genes located within the common loh regions in nhl cases  discovered that nfatc <dig> gene in addition to being inactivated in nhl subtypes, b) is upregulated in cases with loh of ptprj.
nfatc <dig> mean fluorescence intensity values were obtained from an illumina sentrix human- <dig>  expression beadchip  <cit> .
discussion
the implementation of pathway analyses on candidate tumour suppressor genes  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  <cit> , via interactome and pathway analytical approaches, identified the kegg metabolic pathway  and the arginine and proline metabolism pathway  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  that may directly or indirectly  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  <cit> .
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 .
for instance, reprogramming of the glutamine-proline-arginine metabolic circuit has been associated with cancer and found to be regulated by tp <dig> and myc, in which proline acts as the regulatory axis of this circuit  <cit> .
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  <dig> 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 grb <dig>  a pivotal protein in signal transduction that activates the ras-mapk pathway, which is inactivated by the dephosphorylation activity of ptprj.
grb <dig> was also found to interact directly with b2m in an unclear manner  <cit> .
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  <cit> .
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.
ptpn <dig>  highly targeted by loh events in nhl, was also identified in the ptprj interactome.
ptpn <dig> encodes a protein tyrosine phosphatase  with two sh <dig> domains which acts as an intracellular signalling transducer of growth factors and cytokines receptors by regulation of the mapk pathway  <cit> .
furthermore, it has been demonstrated that this ptp plays a crucial role in hematopoiesis and that mutations in ptpn <dig> block the t-cell and b-cell development  <cit> .
mutations in the ptpn <dig> 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  and other types of cancer.
.
interestingly, ptpn <dig> and ptprj share around  <dig> 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 ptpn <dig> may induce a constitutive activation of oncogenic pathways regulated by these ptps.
additionally, the gene nfatc <dig> was identified as being significantly upregulated in nhl cases with loh of ptprj.
interestingly, loh of nfatc <dig> 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 nfatc <dig> in these cases, based on the previously described high resolution loh approach in these nhl cases  <cit> .
nfatc <dig> is a ca++ -dependent protein that regulates t-cell activation and migration, proliferation and angiogenesis by regulation of vegf  <cit> .
moreover, the abrogation of nfatc <dig> 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  <cit> .
in addition, nfatc <dig> was hypothesized as being part of a mechanism whereby intratumoural cd4+cd25+ t-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  <cit> .
interestingly, these treg cells also suppress the cytotoxic activity of cd8+ t-cells  <cit> .
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 nfatc <dig>  and cd8+ t-cells.
further studies are required to determine the role of nfatc <dig> and its inactivation in b-cell lymphomagenesis, the copy-neutral events driving loh of nfatc <dig> and its relation with the inactivation of ptprj, which is known to inhibit the tcr-mediated t-cell activation by dephosphorylation of lat and plcg <dig> pathways  <cit> .
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  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  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- <dig> and inhibits the vegf-mediated cell proliferation, migration, angiogenesis and anti-apoptosis signalling  <cit> .
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 .
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  <cit> .
likewise, as ptprj also affects changes in ca++ concentration  <cit> , 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  <cit> , its relation with ptprj is unclear.
on the other hand, the most enriched pathways in cases with retention of ptprj were dna polymerase, tp <dig> 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 .
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  by activation of mapk pathway  <cit> .
furthermore, an analysis using freely available gene expression profiling data from fl and dlbcl cases expressing high and low levels of ptprj  <cit> , 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  <cit> .
CONCLUSIONS
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.

 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 represents a highly biological and clinical heterogeneous group of blood cancers. two of the most common nhl subtypes are diffuse large b-cell lymphoma, an aggressive lymphoma, and follicular lymphoma, a slow-growing type of lymphoma, which account for around 50% of nhl cases. 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. 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 and gene expression profiling 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 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 analysis in fls and dlbcls, we also identified candidate tumour suppressor genes 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 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. furthermore, ptprj regulates signalling pathways involved in cell growth, proliferation and angiogenesis, such as mapk, plcg pi3k, flt b-cell receptor, pdgfrb and vegfr signalling. nonetheless, the natural ligands of ptprj, syndecan-, a transmembrane heparan sulfate proteoglycan, and thrombospondin-, a homotrimeric glycoprotein, have been identified to induce cell adhesion and inhibit cell growth and angiogenesis, respectively. 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- expression. 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 targeted by loh events participate within common cellular networks that may orchestrate nhl lymphomagenesis in order to determine whether the candidate tumour suppressor genes located within common loh regions across dlbcl and fl cases interact and participate within common cellular networks, a global interactome of a total of genes affected by loh events, was created using the visant platform. as a result of this analysis, the metabolic pathway 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 mapk, ptpn and anxa showed the highest number of links with other genes. 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.figure global interactome of genes commonly affected by loh across nhl patients. an initial list containing genes was used; however, only those with more than one interaction are shown. a total of nodes belonging to pathways were mapped, which indicates the high level of interaction among the genes targeted by loh-driving events in nhl tumours. the genes with labels correspond to the questioned loh genes. in addition to the global interactome study, a gene set enrichment analysis was performed to classify the 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. for instance, faslg was classified into the cytokines and growth factor family and together with ptprj into the cell differentiation markers family. ep was found in the translocated cancer gene, transcription factor and tumour suppressor gene families. additionally, mapk ptpn and nfatc were identified as members of the protein kinases, oncogenes and transcription factors families, respectively. furthermore, out of genes within loh regions were identified and contrasted with collections from the molecular signature database. several cellular pathways were identified to be enriched with genes located within loh regions. 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 sets. interestingly, the tsgs tp53bp and b2m were identified as members of the reactome_double_strand_break_repair and biocarta_ctl pathway pathways, respectively.table list of cellular pathways enriched with genes within loh regions *k/k is the ratio between the number of genes in overlap and the number of genes in gene set. this table combines the results from the collection of gene sets from canonical pathways, biocarta, kegg, reactome, and go gene sets. significant enrichment of loh genes was also observed in a collection of gene sets with chemical and genetic alterations, indicating that these candidate tsgs have been previously implicated in multiple forms of cancer. 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 tp53bp was listed in the pujana_breast_cancer_lit_int_network gene sets. interactome of ptprj overlaps with the interactome of ptpn, 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 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, 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. 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 genes, which are known to interact with ptprj, were previously identified among the genes located within the common loh regions across fls and dlbcls. however, genes that interact indirectly with ptprj were found to be commonly inactivated in nhl cases. among these candidate tsgs, b2m and ptpn were highlighted. b2m was found to interact indirectly with ptprj, through the ptprj-grb2-b2m interaction. ptpn another protein tyrosine phosphatase, was identified to share interacting genes with ptprj. additionally, using previously published high resolution loh data, loh of ptpn was identified in 74% of nhl, 71% of fl and 76% of dlbcl cases.figure interactome of ptprj indicates that the inactivation of ptprj may affect cellular networks and that these networks are not only altered by inactivation of ptprj. ptprj gene network showing a) all the known genes that interact with ptprj and b) only those genes that were also found to be inactivated in nhl cases. two levels of interaction were used to filter the interactions. interestingly, we found that the gene ptpn which was also found inactivated in nhl cases, shares some target genes with ptprj; moreover, b2m has an indirect interaction with ptprj. the genes showing interaction with both ptprj and ptpn are: mapk, grb, erbb, met, pdgfrb, ctnnb, beta 88 kda), plcg, cbl ecotropic retroviral transforming sequence), kit, ntrk, kdr ), pxn, gab, lepr, epor, ghr, cxcr receptor 4), csk, insr, jak, lck, grin2b, tek and tie. interestingly, these genes were found to significantly enrich important signalling pathways in lymphocyte biology, which suggests that a double inactivation of ptprj and ptpn might have an aberrant effect on the function of these pathways.table list of gene sets commonly regulated by ptprj and ptpn based on gene expression *k/k is the ratio between the number of genes in overlap and the number of genes in gene set. loh of ptprj induces deregulation of genes and signalling pathways in nhl tumours we also used gene expression profiling 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, suggesting that some pathways might be differentially enriched between the two categories. in order to prove this hypothesis, a gene set enrichment analysis identified some cellular pathways with significant upregulation in cases with retention and some others, in cases with loh. the genes that were upregulated in the most significant gene sets are described in table  the lists of the gene sets enriched in each loh status are detailed in the additional file and additional file nonetheless, it is important to highlight that despite all enrichment analyses exhibiting fdr values equal to the results from these gene-set enrichments, described in figure  and table  were analysed based on the nominal significance 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.figure genes and pathways affected by loh of ptprj. a) heatmap of genes differentially expressed between cases with retention and loh calls. b) dna polymerase, tp signaling pathway and proteasome were identified as the most enriched gene sets in cases with retention of ptprj. c) in cases with loh, the gene sets vegf signaling pathway, melanoma, erbb signaling pathway, calcium signaling pathway, fc epsilon ri signaling pathway and mapk signaling pathway were the most enriched. despite an fdr of for all the enrichments, the results from this gene-set enrichment analysis were analysed based on the p values instead of the fdr scores. the high fdr values obtained in this analysis may be due to the low population size. the lists of all gene sets from this analysis are provided in additional file and additional file list of genes enriched within the gene sets differentially expressed between cases with retention and loh of ptprj polq,prim pold pole polk,pole pole pols,rev pola rfc pold polb,poli,pole cdc pten,tnfrsf10b,chek pmaip cdk cdk sfn,perp,p53aip fas,ccnb serpinb casp ccnb ei ccnd tp53i ccne rfwd ccnb mdm2 psmb psma psma psma psmb psmd psmd psmd psmc psma psma psmd psma psma6 mapkapk pla2g pxn,pla2g pla2g2f,map2k map2k raf prkcg,pik3r pla2g1b,mapk mapk sphk akt pla2g2a,bad,rac pla2g pik3cd,nfatc prkca,pik3cg fgf ccnd met,fgf fgf map2k map2k igf raf pik3r mapk cdk akt fgf bad,fgf araf,fgf fgfr pik3cd,pdgfra,pdgfd,fgf pik3cg camk2g,shc map2k map2k raf stat5b,prkcg,pak nrg rps6kb camk2b,pik3r pak nrg mapk erbb akt tgfa,bad,araf,cbl,erbb pik3cd,prkca,pik3cg, pak3 gna htr cckbr,chp,adra1d,itpr ryr slc25a phka nos phka slc8a cacna1c,atp2b mylk,agtr camk2g,chrna camk mylk prkcg,camk2b,grm bdkrb adcy ryr erbb sphk trhr,calm prkacg,adora2a,pln,gna tbxa2r,ptafr,erbb adcy p2rxl tacr cacna1d,pdgfra,hrh ntsr plcb htr2b,prkca, cacna1i,slc25a ptger itpkb pla2g2f,vav map2k map2k raf pik3r pla2g1b,mapk mapk akt pla2g2a,rac il pla2g pik3cd,il pdk prkca,pik3cg dusp map2k fgf cacna2d mef2c,fgf rps6ka cacna2d tgfbr taok cacna2d rac nfkb chp,flna,srf,fgf daxx,rps6ka rras,mapkapk pla2g ntrk cacna1c,pla2g fgf fgf rap1a,flnc,map3k pla2g2f,cacnb map2k gadd45a,rps6ka nf map2k rps6ka raf prkcg,pak taok mos,mapk8ip pla2g1b, mapk mapk rras prkacg,gadd45g,akt ntf fgf pla2g2a,flnb,rac fgf fgf rasgrf fgfr pla2g map3k8nlk ptpn cacna1d il1r pdgfra arrb map3k rps6ka ppp5c fgf gna cacng map3k dusp prkca map2k cacna1i mapk8ip2 vegfb,rictor,pgf,akt prkaa rps6ka rps6ka rps6ka igf rps6kb pik3r eif4b,mapk vegfc,akt figf, pik3cd,pik3cg pim nfkb stat ccnd flt map2k cebpa,map2k raf rps6kb jup,mapk tcf akt bad,araf,pik3cd, pik3cg,kit nfatc 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 candidate tsgs located within loh regions against the genes with differential expression between cases with retention and loh of ptprj. this approach identified only the nfatc gene, located at 16q to be commonly inactivated and differentially expressed in nhl cases with loh of ptprj. comparison of the mean fluorescence intensity values for nfatc between nhl cases with retention and loh of ptprj revealed that nfatc was significantly upregulated in cases with loh of ptprj. furthermore, it is important to mention that loh of nfatc was found in 58% of nhl, 48% of fl and 67% of dlbcl cases.figure identification of genes located within the common loh regions in nhl cases whose expression is affected by the loh status of ptprj. a) comparison of the differentially expressed genes between nhl cases with retention and loh of ptprj and genes located within the common loh regions in nhl cases discovered that nfatc gene in addition to being inactivated in nhl subtypes, b) is upregulated in cases with loh of ptprj. nfatc mean fluorescence intensity values were obtained from an illumina sentrix human- expression beadchip. discussion the implementation of pathway analyses on candidate tumour suppressor genes 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 and the arginine and proline metabolism pathway 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 that may directly or indirectly 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. 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. for instance, reprogramming of the glutamine-proline-arginine metabolic circuit has been associated with cancer and found to be regulated by tp and myc, in which proline acts as the regulatory axis of this circuit. 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 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 grb a pivotal protein in signal transduction that activates the ras-mapk pathway, which is inactivated by the dephosphorylation activity of ptprj. grb 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. 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. ptpn highly targeted by loh events in nhl, was also identified in the ptprj interactome. ptpn encodes a protein tyrosine phosphatase with two sh 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 ptpn block the t-cell and b-cell development. mutations in the ptpn 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 and other types of cancer.. interestingly, ptpn and ptprj share around 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 ptpn may induce a constitutive activation of oncogenic pathways regulated by these ptps. additionally, the gene nfatc was identified as being significantly upregulated in nhl cases with loh of ptprj. interestingly, loh of nfatc 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 nfatc in these cases, based on the previously described high resolution loh approach in these nhl cases. nfatc is a ca++ -dependent protein that regulates t-cell activation and migration, proliferation and angiogenesis by regulation of vegf. moreover, the abrogation of nfatc 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, nfatc was hypothesized as being part of a mechanism whereby intratumoural cd4+cd25+ t-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 nfatc and cd8+ t-cells. further studies are required to determine the role of nfatc and its inactivation in b-cell lymphomagenesis, the copy-neutral events driving loh of nfatc and its relation with the inactivation of ptprj, which is known to inhibit the tcr-mediated t-cell activation by dephosphorylation of lat and plcg pathways. 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 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 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- 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. 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. likewise, as ptprj also affects changes in ca++ concentration, 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, tp 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. 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 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, 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. 
these analyses also identified genes within the interactome of ptprj that when inactivated in nhl may play an important role in tumourigenesis. recently, using a high resolution loss of heterozygosity analysis, we identified a number of potential tumour suppressor genes within common loh regions across cases suffering from two of the most common forms of non-hodgkin’s lymphoma, follicular lymphoma and diffuse large b-cell lymphoma. 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. the evidence obtained in this research supports findings suggesting that fl and dlbcl share common pathogenic mechanisms. from these studies loh of the protein tyrosine phosphatase receptor type j 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. 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 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.
 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 analysis, we identified a number of potential tumour suppressor genes within common loh regions across cases suffering from two of the most common forms of non-hodgkin’s lymphoma, follicular lymphoma and diffuse large b-cell lymphoma. from these studies loh of the protein tyrosine phosphatase receptor type j 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. these analyses also identified genes within the interactome of ptprj 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 nfatc. 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. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords lohptprjinteractomepathway analysisnhlissue-copyright-statement© the author 2014 