Deletion of Ptprd and Cdkn2a cooperate to accelerate tumorigenesis

Protein tyrosine phosphatase delta is a STAT3-phosphatase and suppressor of metabolic liver disease

ABSTRACT

Objective

Impaired hepatic expression of protein tyrosine phosphatase delta (PTPRD) is associated with increased STAT3 transcriptional activity and reduced survival from hepatocellular carcinoma in patients with chronic hepatitis C virus infection. However, the PTPRD-expressing hepatic cell types, signalling pathways responsive to PTPRD and their role in non-viral liver disease are largely unknown.

Methods

We studied PTPRD expression in single-cell and bulk liver transcriptomic data from mice and humans, and established a Ptprd-deficient mouse model for metabolic dysfunction-associated steatohepatitis (MASH). Identified pathways were validated by perturbation studies in human hepatocytes and PTPRD substrates by pull-down assays. The clinical relevance was further explored in a cohort with metabolic disease by ranking patients according to PTPRD expression and analysing its association with metabolic disease markers.

Results

The analysis of individuals ranked according to PTPRD expression and Ptprd-deficient mice, showed that PTPRD levels were associated with hepatic glucose/lipid signalling and peroxisome function. Hepatic PTPRD expression is impaired in aetiologies of chronic liver diseases that are associated with metabolic disease. We further validated PTPRD as a STAT3 phosphatase in the liver, acting as a regulator of peroxisomal fatty acid metabolism. During MASH, low PTPRD led to increased liver steatosis in Ptprd+/- mice and a pronounced unfolded protein response, which impacts insulin signalling. Accordingly, silencing of PTPRD blunted insulin-induced AKT phosphorylation. Patients with obesity and low hepatic PTPRD expression exhibit increased levels of metabolic risk factors.

Conclusion

Our data revealed an important regulatory role of the hepatic PTPRD-STAT3 axis in maintaining glucose/lipid homeostasis, which is recapitulated in clinical manifestations of metabolic liver disease.

Baseline circulating tumour DNA and interim PET predict response in relapsed/refractory classical Hodgkin lymphoma

Reliable biomarkers for early identification of treatment failure in relapsed/refractory (r/r) classical Hodgkin lymphoma (cHL) are lacking. Circulating tumour DNA (ctDNA) profiling has emerged as a powerful predictive and prognostic tool in several haemopoietic and non-haemopoietic malignancies and may guide rational treatment choices in r/r cHL. To assess the predictive and prognostic value of ctDNA, we performed a retrospective analysis on 55 r/r cHL patients treated with the bendamustine, gemcitabine and vinorelbine (BEGEV) regimen and additionally evaluated the potential utility of integrating ctDNA with interim [18 F]-FDG positron emission tomography (iPET). Baseline ctDNA genotyping in r/r cHL mirrored gene mutations and pathways involved in newly diagnosed cHL. We found that baseline ctDNA quantification and serial ctDNA monitoring have prognostic value in r/r cHL receiving salvage chemotherapy. Lastly, integrating ctDNA quantification with iPET evaluation may improve the early identification of patients at high risk of failing standard salvage therapy, who may benefit from an early switch to immunotherapeutic agents. Collectively, our results support the implementation of non-invasive methods to detect minimal residual disease in recurrent cHL and justify its prospective evaluation in appropriately designed clinical trials.

The genomic and transcriptomic landscape of advanced renal cell cancer for individualized treatment strategies

Differences in the clinical course and treatment responses in individual patients with advanced renal cell carcinoma (RCC) can largely be explained by the different genomics of this disease. To improve the personalized treatment strategy and survival outcomes for patients with advanced RCC, the genomic make-up in patients with advanced RCC was investigated to identify putative actionable variants and signatures. In this prospective multicenter study (NCT01855477), whole-genome sequencing (WGS) data of locally advanced and metastatic tissue biopsies and matched whole-blood samples were collected from 91 patients with histopathologically confirmed RCC. WGS data were analyzed for small somatic variants, copy-number alterations and structural variants. For a subgroup of patients, RNA sequencing (RNA-Seq) data could be analyzed. RNA-Seq data were clustered on immunogenic and angiogenic gene expression patterns according to a previously developed angio-immunogenic gene signature. In all patients with papillary and clear cell RCC, putative actionable drug targets were detected by WGS, of which 94% were on-label available. RNA-Seq data of clear cell and papillary RCC were clustered using a previously developed angio-immunogenic gene signature. Analyses of driver mutations and RNA-Seq data revealed clear differences among different RCC subtypes, showing the added value of WGS and RNA-Seq over clinicopathological data. By improving both histological subtyping and the selection of treatment according to actionable targets and immune signatures, WGS and RNA-Seq may improve therapeutic decision making for most patients with advanced RCC, including patients with non-clear cell RCC for whom no standard treatment is available to data. Prospective clinical trials are needed to evaluate the impact of genomic and transcriptomic diagnostics on survival outcome for advanced RCC patients.

Whole genome analysis reveals the genomic complexity in metastatic cutaneous squamous cell carcinoma

Metastatic cutaneous squamous cell carcinoma (CSCC) is a highly morbid disease requiring radical surgery and adjuvant therapy, which is associated with a poor prognosis. Yet, compared to other advanced malignancies, relatively little is known of the genomic landscape of metastatic CSCC. We have previously reported the mutational signatures and mutational patterns of CCCTC-binding factor (CTCF) regions in metastatic CSCC. However, many other genomic components (indel signatures, non-coding drivers, and structural variants) of metastatic CSCC have not been reported. To this end, we performed whole genome sequencing on lymph node metastases and blood DNA from 25 CSCC patients with regional metastases of the head and neck. We designed a multifaceted computational analysis at the whole genome level to provide a more comprehensive perspective of the genomic landscape of metastatic CSCC. In the non-coding genome, 3′ untranslated region (3′UTR) regions of EVC (48% of specimens), PPP1R1A (48% of specimens), and ABCA4 (20% of specimens) along with the tumor-suppressing long non-coding RNA (lncRNA) LINC01003 (64% of specimens) were significantly functionally altered (Q-value < 0.05) and represent potential non-coding biomarkers of CSCC. Recurrent copy number loss in the tumor suppressor gene PTPRD was observed. Gene amplification was much less frequent, and few genes were recurrently amplified. Single nucleotide variants driver analyses from three tools confirmed TP53 and CDKN2A as recurrently mutated genes but also identified C9 as a potential novel driver in this disease. Furthermore, indel signature analysis highlighted the dominance of ID signature 13 (ID13) followed by ID8 and ID9. ID9 has previously been shown to have no association with skin melanoma, unlike ID13 and ID8, suggesting a novel pattern of indel variation in metastatic CSCC. The enrichment analysis of various genetically altered candidates shows enrichment of “TGF-beta regulation of extracellular matrix” and “cell cycle G1 to S check points.” These enriched terms are associated with genetic instability, cell proliferation, and migration as mechanisms of genomic drivers of metastatic CSCC.

Association of PTPRD/PTPRT Mutation With Better Clinical Outcomes in NSCLC Patients Treated With Immune Checkpoint Blockades

The common gamma receptor–dependent cytokines and their JAK-STAT pathways play important roles in T cell immunity and have been demonstrated to be related with response to immune checkpoint blockades (ICBs). PTPRD and PTPRT are phosphatases involved in JAK-STAT pathway. However, their clinical significance for non-small cell lung cancer (NSCLC) treated with ICBs is still unclear. Genomic and survival data of NSCLC patients administrated with anti–PD-1/PD-L1 or anti–CTLA-4 antibodies (Rizvi2015; Hellmann2018; Rizvi2018 Samstein2019) were retrieved from publicly accessible data. Genomic, survival and mRNA data of 1007 patients with NSCLC were obtained from The Cancer Genome Atlas (TCGA). PTPRD/PTPRT mutation was significantly associated with better progression-free survival (PFS) in three independent Rizvi2015, Hellmann2018 and Rizvi2018 cohorts. The median PFS for PTPRD/PTPRT mutant-type vs. wild-type NSCLC patients were not reached vs. 6.3 months (Rizvi2015, HR = 0.16; 95% CI, 0.02-1.17; P=0.03), 24.0 vs. 5.4 months (Hellmann2018, HR, 0.49; 95% CI, 0.26-0.94; P=0.03), 5.6 vs. 3.0 months (Rizvi2018, HR = 0.64; 95% CI, 0.44-0.92; P=0.01) and 6.8 vs. 3.5 months (Pooled cohort, HR, 0.54; 95% CI, 0.39-0.73; P<0.0001) respectively. PTPRD/PTPRT mutation was an independent predictive factor for PFS in pooled cohort (P = 0.01). Additionally, PTPRD/PTPRT mutation associated with better overall survival (OS) in Samstein2019 cohort (19 vs. 10 months, P=0.03). While similar clinical benefits were not observed in patients without ICBs treatment (TCGA cohort, P=0.78). In the further exploratory analysis, PTPRD/PTPRT mutation was significantly associated with increased tumor mutation burden and higher mRNA expression of JAK1 and STAT1. Gene Set Enrichment Analysis revealed prominent enrichment of signatures related to antigen processing and presentation in patients with PTPRD/PTPRT mutation. This work suggested that PTPRD/PTPRT mutation might be a potential positive predictor for ICBs in NSCLC. These results need to be further confirmed in future.

Decreased HLF Expression Predicts Poor Survival in Lung Adenocarcinoma

Background

Lung adenocarcinoma (LUAD) is a type of non-small cell carcinoma. Its pathogenesis is being explored and there is no cure for the disease.

Material/Methods

The Gene Expression Omnibus (GEO) was searched to obtain data on expression of messenger RNA. GEO2R, an interactive web tool, was used to calculate the differentially expressed genes (DEGs) in LUAD. All the DEGs from different datasets were imported into VENNY 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html) to identify the intersection of the DEGs. An online analysis tool, the Database for Annotation, Visualization, and Integrated Discovery (DAVID), was used to help understand the biological meaning of DEG enrichment in LUAD. Cytoscape 3.7.2 was used to perform centrality analysis and visualize hub genes and related networks. Furthermore, the prognostic value of the hub genes was evaluated with the Kaplan-Meier plotter survival analysis tool.

Results

The GEO database was used to obtain RNA sequencing information for LUAD and normal tissue from the GSE118370, GSE136043, and GSE140797 datasets. A total of 376 DEGs were identified from GSE118370, 248 were identified from GSE136403, and 718 DEGs were identified from GSE140797. The 10 genes with the highest degrees of expression – the hub genes – were CAV1, TEK, SLIT2, RHOJ, DGSX, HLF, MEIS1, PTPRD, FOXF1, and ADRB2. In addition, Kaplan-Meier survival evaluation showed that CAV1, TEK, SLIT2, HLF, MEIS1, PTPRD, FOXF1, and ADRB2 were associated with favorable outcomes for LUAD.

Conclusions

CAV1, TEK, SLIT2, HLF, MEIS1, PTPRD, FOXF1, and ADRB2 are hub genes in the DEG interaction network for LUAD and are involved in the development of and prognosis for the disease. The mechanisms underlying these genes should be the subject of further studies.

Biomarker-guided targeted and immunotherapies in malignant pleural mesothelioma

Malignant pleural mesothelioma (MPM) is a lethal thoracic malignancy whose incidence is still increasing worldwide. MPM is characterized by frequent inactivation of tumor-suppressor genes (TSGs), e.g., the homozygous deletion of CDKN2A/2B and various genetic alterations that inactivate BAP1, NF2, LATS1/2, and TP53. The leading cause for the poor prognosis of patients with MPM is the lack of effective treatment options, with conventional chemotherapy being the standard of care in the clinic, which has remained unchanged for almost 20 years. Precision oncology, a burgeoning effort to provide precise cancer treatment tailored to unique molecular changes in individual patients, has made tremendous progress in the last decade in several cancers, but not in MPM. Recent studies indicate a high degree of tumor heterogeneity in MPM and the importance to optimize histological and molecular classifications for improved treatment. In this review, we provide an up-to-date overview of recent advances in MPM by focusing on new stratifications of tumor subgroups, specific vulnerabilities associated with functional loss of TSGs and other biomarkers, and potential clinical implications. The molecularly based subdivisions not only deepen our understanding of MPM pathobiology, but more importantly, they may raise unprecedented new hopes for personalized treatment of MPM patients with biomarker-guided targeted and immunotherapies.

IL-6 induces tumor suppressor protein tyrosine phosphatase receptor type D by inhibiting miR-34a to prevent IL-6 signaling overactivation

Protein tyrosine phosphatase receptor type D (PTPRD) is a tumor suppressor gene that is epigenetically silenced and mutated in several cancers, including breast cancer. Since IL-6/STAT3 signaling is often hyperactivated in breast cancer and STAT3 is a direct PTPRD substrate, we investigated the role of PTPRD in breast cancer and the association between PTPRD and IL-6/STAT3 signaling. We found that PTPRD acts as a tumor suppressor in breast cancer tissues and that high PTPRD expression is positively associated with tumor size, lymph node metastasis, PCNA expression, and patient survival. Moreover, breast cancers with high PTPRD expression tend to exhibit high IL-6 and low phosphorylated-STAT3 expression. IL-6 was found to inhibit miR-34a transcription and induce PTPRD expression in breast cancer and breast epithelial cells, whereas PTPRD was shown to mediate activated STAT3 dephosphorylation and to be a conserved, direct target of miR-34a. IL-6-induced PTPRD upregulation was blocked by miR-34a mimics, whereas experimental PTPRD overexpression suppressed MDA-MB-231 cell migration, invasion, and epithelial to mesenchymal transition, decreased STAT3 phosphorylation, and increased miR-34a transcription. Our findings suggest that PTPRD mediates activated STAT3 dephosphorylation and is induced by the IL-6/STAT3-mediated transcriptional inhibition of miR-34a, thereby establishing a negative feedback loop that inhibits IL-6/STAT3 signaling overactivation.

Integrated analysis of relapsed B-cell precursor Acute Lymphoblastic Leukemia identifies subtype-specific cytokine and metabolic signatures

Recent efforts reclassified B-Cell Precursor Acute Lymphoblastic Leukemia (BCP-ALL) into more refined subtypes. Nevertheless, outcomes of relapsed BCP-ALL remain unsatisfactory, particularly in adult patients where the molecular basis of relapse is still poorly understood. To elucidate the evolution of relapse in BCP-ALL, we established a comprehensive multi-omics dataset including DNA-sequencing, RNA-sequencing, DNA methylation array and proteome MASS-spec data from matched diagnosis and relapse samples of BCP-ALL patients (n = 50) including the subtypes DUX4, Ph-like and two aneuploid subtypes. Relapse-specific alterations were enriched for chromatin modifiers, nucleotide and steroid metabolism including the novel candidates FPGS, AGBL and ZNF483. The proteome expression analysis unraveled deregulation of metabolic pathways at relapse including the key proteins G6PD, TKT, GPI and PGD. Moreover, we identified a novel relapse-specific gene signature specific for DUX4 BCP-ALL patients highlighting chemotaxis and cytokine environment as a possible driver event at relapse. This study presents novel insights at distinct molecular levels of relapsed BCP-ALL based on a comprehensive multi-omics integrated data set including a valuable proteomics data set. The relapse specific aberrations reveal metabolic signatures on genomic and proteomic levels in BCP-ALL relapse. Furthermore, the chemokine expression signature in DUX4 relapse underscores the distinct status of DUX4-fusion BCP-ALL.

p16 overexpression and 9p21 deletion are linked to unfavorable tumor phenotype in breast cancer

Overexpression of the p16 tumor suppressor, but also deletion of its gene locus 9p21, is linked to unfavorable tumor phenotype and poor prognosis in breast cancer. To better understand these contradictory observations, and to clarify the prognostic impact of p16 expression and 9p21 deletion, a tissue microarray (TMA) with 2,197 breast cancers was analyzed by fluorescence in-situ hybridization and immunohistochemistry (FISH) for 9p21 deletion and p16 expression. p16 immunostaining was weak in 25.6%, moderate in 7.1%, and strong in 12.7% of 1,684 evaluable cancers. Strong p16 staining was linked to advanced tumor stage (p = 0.0003), high-grade (p < 0.0001), high tumor cell proliferation (p < 0.0001), negative hormone receptor (ER/PR) status (p < 0.0001 each), and shorter overall survival (p = 0.0038). 9p21 deletion was found in 15.3% of 1,089 analyzable breast cancers, including 1.7% homozygous and 13.6% heterozygous deletions. 9p21 deletion was linked to adverse tumor features, including high-grade (p < 0.0001) and nodal positive cancers (p = 0.0063), high cell proliferation (p < 0.0001), negative hormone receptor (ER/PR) status (p ≤ 0.0006), and HER2 amplification (p = 0.0078). Patient outcome was worse in 9p21 deleted than in undeleted cancers (p = 0.0720). p16 expression was absent in cancers harboring homozygous 9p21 deletions, but no difference in p16 expression was found between cancers with (59.2% p16 positive) and without heterozygous 9p21 deletion (51.3% p16 positive, p = 0.0256). In summary, p16 expression is unrelated to partial 9p21 deletion, but both alterations are linked to aggressive breast cancer phenotype. High-level p16 expression is a strong predictor of unfavorable disease course in breast cancer.

Deletion and downregulation of MTAP contribute to the motility of esophageal squamous carcinoma cells

Esophageal squamous cell carcinoma (ESCC) is among the most common malignancies, with a low 5-year overall survival rate. In previous studies, we and others have found that 9p21.3 was the most frequently deleted region in ESCC. The MTAP gene, which is located close to CDKN2A/B in 9p21.3, encodes methylthioadenosine phosphorylase. This enzyme plays an important role during the process of adenosine transfer. In the present study, we found that MTAP is deleted at the genomic level in 19.1% (64/341) of primary ESCC tumors, and decreased mRNA and protein expression were present in 31.1% (28/90) and 33.3% (6/18) of ESCCs, respectively. Further statistical analysis showed a positive correlation between deletion and decreased mRNA expression of MTAP in the ESCC tissues tested (coefficient: 0.826; P=1.17×10⁻²³). Knockdown of MTAP expression using small interfering RNA-mediated silencing promoted the invasion and migration of ESCC cells. Also, overexpression of MATP using pcDNA3.1-MTAP plasmid decreased the cell invasion and migration. At the molecular level, MTAP knockdown downregulated E-cadherin and p-GSK3β but upregulated Slug expression. Our results indicated that MTAP deletion results in the decreased expression in ESCCs and that it plays a role in promoting the mobility and inducing the epithelial-to-mesenchymal transition of ESCC cells via the GSK3β/Slug/E-cadherin axis. The data suggest that MTAP might function as a tumor suppressor gene in ESCC.

Genetic defects of the IRF1-mediated major histocompatibility complex class I antigen presentation pathway occur prevalently in the JAK2 gene in non-small cell lung cancer

Recognition of major histocompatibility complex (MHC) class I antigens on tumor cells by cytotoxic T cells is involved in T cell-mediated tumor immune surveillance and immune checkpoint therapy. The interferon-γ (IFNγ)-IRF1 signaling pathway regulates MHC class I antigen presentation. To examine genetic defects of the IFNγ-IRF1 pathway in non-small cell lung cancer (NSCLC), we analyzed The Cancer Genome Atlas (TCGA) lung adenocarcinoma (LuAd) and lung squamous cell carcinoma (LuSc) data. Loss-of-function (LOF) genetic alterations of the IFNγ-IRF1 pathway genes (IFNGR1, IFNGR2, JAK1, JAK2, STAT1, IRF1) were found in 64 (6.3%) of 1,016 patients. These genetic defects occur prevalently in JAK2 (33 cases) and often through deletions (29 cases) of chromosome 9p24.1. JAK2 deletions were frequently, but not always, associated with deletions of PD-L1 gene (CD274), PD-L2 gene (PDCD1LG2), PTPRD, and CDKN2A/CDKN2B at the chromosome 9p24.1-9p21.3 region. IRF1 expression was correlated with immune cytolytic activity markers GZMA and PRF1 in NSCLC. IFNγ induced IRF1 expression and cell surface HLA-A/HLA-B/HLA-C (HLA-ABC) in A549, H661, H292, and H2172 cells that contained the wildtype JAK2, but not in H1573 and H1623 cells that were JAK2 defective. Deletion of JAK2 or inhibition of the JAK2 kinase activity resulted in loss of IFNγ-induced IRF1 and cell surface HLA-ABC in JAK2 wildtype NSCLC cells, whereas expression of exogenous JAK2 in H1573 cells restored the IFNγ responses. These findings show that JAK2 deficiency is the major mechanism of genetic defects of the IFNγ-IRF1 pathway in NSCLC and reveal a previously unrecognized significance of chromosome 9p deletion in NSCLC.

Recurrent epigenetic silencing of the PTPRD tumor suppressor in laryngeal squamous cell carcinoma

Cellular processes like differentiation, mitotic cycle, and cell growth are regulated by tyrosine kinases with known oncogenic potential and tyrosine phosphatases that downmodulate the first. Therefore, tyrosine phosphatases are recurrent targets of gene alterations in human carcinomas. We and others suggested recently a tumor suppressor function of the PTPRD tyrosine phosphatase and reported homozygous deletions of the PTPRD locus in laryngeal squamous cell carcinoma. In this study, we investigated other gene-inactivating mechanisms potentially targeting PTPRD, including loss-of-function mutations and also epigenetic alterations like promoter DNA hypermethylation. We sequenced the PTPRD gene in eight laryngeal squamous cell carcinoma cell lines but did not identify any inactivating mutations. In contrast, by bisulfite pyrosequencing of the gene promoter region, we identified significantly higher levels of methylation (p = 0.001 and p = 0.0002, respectively) in 9/14 (64%) laryngeal squamous cell carcinoma cell lines and 37/79 (47%) of primary laryngeal squamous cell carcinoma tumors as compared to normal epithelium of the upper aerodigestive tract. There was also a strong correlation (p = 0.0001) between methylation and transcriptional silencing for the PTPRD gene observed in a cohort of 497 head and neck tumors from The Cancer Genome Atlas dataset suggesting that DNA methylation is the main mechanism of PTPRD silencing in these tumors. In summary, our data provide further evidence of the high incidence of PTPRD inactivation in laryngeal squamous cell carcinoma. We suggest that deletions and loss-of-function mutations are responsible for PTPRD loss only in a fraction of cases, whereas DNA methylation is the dominating mechanism of PTPRD inactivation.

Molecular pathogenesis of splenic and nodal marginal zone lymphoma

Genomic studies have improved our understanding of the biological basis of splenic (SMZL) and nodal (NMZL) marginal zone lymphoma by providing a comprehensive and unbiased view of the genes/pathways that are deregulated in these diseases. Consistent with the physiological involvement of NOTCH, NF-κB, B-cell receptor and toll-like receptor signaling in mature B-cells differentiation into the marginal zone B-cells, many oncogenic mutations of genes involved in these pathways have been identified in SMZL and NMZL. Beside genetic lesions, also epigenetic and post-transcriptional modifications contribute to the deregulation of marginal zone B-cell differentiation pathways in SMZL and NMZL. This review describes the progress in understanding the molecular mechanism underlying SMZL and NMZL, including molecular and post-transcriptional modifications, and discusses how information gained from these efforts has provided new insights on potential targets of diagnostic, prognostic and therapeutic relevance in SMZL and NMZL.

The genetics of nodal marginal zone lymphoma

Nodal marginal zone lymphoma (NMZL) is a rare, indolent B-cell tumor that is distinguished from splenic marginal zone lymphoma (SMZL) by the different pattern of dissemination. NMZL still lacks distinct markers and remains orphan of specific cancer gene lesions. By combining whole-exome sequencing, targeted sequencing of tumor-related genes, whole-transcriptome sequencing, and high-resolution single nucleotide polymorphism array analysis, we aimed at disclosing the pathways that are molecularly deregulated in NMZL and we compare the molecular profile of NMZL with that of SMZL. These analyses identified a distinctive pattern of nonsilent somatic lesions in NMZL. In 35 NMZL patients, 41 genes were found recurrently affected in ≥3 (9%) cases, including highly prevalent molecular lesions of MLL2 (also known as KMT2D; 34%), PTPRD (20%), NOTCH2 (20%), and KLF2 (17%). Mutations of PTPRD, a receptor-type protein tyrosine phosphatase regulating cell growth, were enriched in NMZL across mature B-cell tumors, functionally caused the loss of the phosphatase activity of PTPRD, and were associated with cell-cycle transcriptional program deregulation and increased proliferation index in NMZL. Although NMZL shared with SMZL a common mutation profile, NMZL harbored PTPRD lesions that were otherwise absent in SMZL. Collectively, these findings provide new insights into the genetics of NMZL, identify PTPRD lesions as a novel marker for this lymphoma across mature B-cell tumors, and support the distinction of NMZL as an independent clinicopathologic entity within the current lymphoma classification.

Protein tyrosine phosphatases receptor type D is a potential tumour suppressor gene inactivated by deoxyribonucleic acid methylation in paediatric acute myeloid leukaemia

AIM

Protein tyrosine phosphatases receptor type D (PTPRD) is a tumour suppressor gene, and its epigenetic silencing is frequently found in glioblastoma. As aberrant deoxyribonucleic acid (DNA) methylation patterning has been shown to play a role in leukaemogenesis, we studied the promoter methylation, expression profiles and molecular functions of PTPRD in paediatric patients with acute myeloid leukaemia (AML).

METHODS

Bone marrow specimens were obtained from 32 Chinese patients with a mean age of 7.2 years (range 1.1-16.5). PTPRD and methylation status were evaluated by real-time polymerase chain reaction (PCR) and methylation-specific PCR. Western blot and flow cytometry techniques were also used.

RESULTS

PTPRD expression was decreased by promoter region methylation in six AML cells and methylated in 21 (65.6%) of the 32 samples. In addition, PTPRD expression could be induced by the DNA demethylating agent 5-aza-2'-deoxycytidine. Furthermore, functional studies showed that overexpression of PTPRD in AML cells inhibited cell proliferation and clonogenicity as well as inducing apoptosis. However, PTPRD knockdown increased cell proliferation. These effects were associated with downregulation of cyclin D1, c-myc and upregulation of Bax.

CONCLUSION

The results of this study demonstrated that PTPRD was a potential tumour suppressor gene inactivated by DNA methylation in paediatric AML.

PTPRD is homozygously deleted and epigenetically downregulated in human hepatocellular carcinomas

PTPRD (protein tyrosine phosphatase, receptor type, D) is a tumor suppressor gene, frequently inactivated through deletions or epigenetic mechanisms in several cancers with importance for global health. In this study, we provide new and functionally integrated evidence on genetic and epigenetic alterations of PTPRD gene in hepatocellular carcinomas (HCCs). Importantly, HCC is the sixth most common malignancy and the third most common cause of cancer-related mortality worldwide. We used a high throughput single nucleotide polymorphism (SNP) microarray assay (Affymetrix, 10K2.0 Assay) covering the whole genome to screen an extensive panel of HCC cell lines (N=14 in total) to detect DNA copy number changes. PTPRD expression was determined in human HCCs by Q-RT-PCR and immunohistochemistry. Promoter hypermethylation was assessed by combined bisulfite restriction analysis (COBRA). DNA methyl transferase inhibitor 5-azacytidine (5-AzaC) and/or histone deacetylase inhibitor Trichostain A (TSA) were used to restore the expression. We identified homozygous deletions in Mahlavu and SNU475 cells, in the 5'UTR and coding regions, respectively. PTPRD mRNA expression was downregulated in 78.5% of cell lines and 82.6% of primary HCCs. PTPRD protein expression was also found to be lost or reduced in HCC tumor tissues. We found promoter hypermethylation in 22.2% of the paired HCC samples and restored PTPRD expression by 5-AzaC and/or TSA treatments. In conclusion, PTPRD is homozygously deleted and epigenetically downregulated in HCCs. We hypothesize PTPRD as a tumor suppressor candidate and potential cancer biomarker in human HCCs. This hypothesis is consistent with compelling evidences in other organ systems, as discussed in this article. Further functional assays in larger samples may ascertain the contribution of PTPRD to hepatocarcinogenesis in greater detail, not to forget its broader importance for diagnostic medicine and the emerging field of personalized medicine in oncology.

Identification of Copy Number Aberrations in Breast Cancer Subtypes Using Persistence Topology

DNA copy number aberrations (CNAs) are of biological and medical interest because they help identify regulatory mechanisms underlying tumor initiation and evolution. Identification of tumor-driving CNAs (driver CNAs) however remains a challenging task, because they are frequently hidden by CNAs that are the product of random events that take place during tumor evolution. Experimental detection of CNAs is commonly accomplished through array comparative genomic hybridization (aCGH) assays followed by supervised and/or unsupervised statistical methods that combine the segmented profiles of all patients to identify driver CNAs. Here, we extend a previously-presented supervised algorithm for the identification of CNAs that is based on a topological representation of the data. Our method associates a two-dimensional (2D) point cloud with each aCGH profile and generates a sequence of simplicial complexes, mathematical objects that generalize the concept of a graph. This representation of the data permits segmenting the data at different resolutions and identifying CNAs by interrogating the topological properties of these simplicial complexes. We tested our approach on a published dataset with the goal of identifying specific breast cancer CNAs associated with specific molecular subtypes. Identification of CNAs associated with each subtype was performed by analyzing each subtype separately from the others and by taking the rest of the subtypes as the control. Our results found a new amplification in 11q at the location of the progesterone receptor in the Luminal A subtype. Aberrations in the Luminal B subtype were found only upon removal of the basal-like subtype from the control set. Under those conditions, all regions found in the original publication, except for 17q, were confirmed; all aberrations, except those in chromosome arms 8q and 12q were confirmed in the basal-like subtype. These two chromosome arms, however, were detected only upon removal of three patients with exceedingly large copy number values. More importantly, we detected 10 and 21 additional regions in the Luminal B and basal-like subtypes, respectively. Most of the additional regions were either validated on an independent dataset and/or using GISTIC. Furthermore, we found three new CNAs in the basal-like subtype: a combination of gains and losses in 1p, a gain in 2p and a loss in 14q. Based on these results, we suggest that topological approaches that incorporate multiresolution analyses and that interrogate topological properties of the data can help in the identification of copy number changes in cancer.

Loss-of-Function PTPRD Mutations Lead to Increased STAT3 Activation and Sensitivity to STAT3 Inhibition in Head and Neck Cancer

BACKGROUND

Protein tyrosine phosphatase receptor type D (PTPRD) is a putative tumor suppressor in several cancers including head and neck squamous cell carcinoma (HNSCC). STAT3 is a frequently hyperactivated oncogene in HNSCC. As STAT3 is a direct substrate of PTPRD, we sought to determine the genetic or epigenetic alterations of PTPRD that contribute to overactive STAT3 in HNSCC.

METHODS

We analyzed data from The Cancer Genome Atlas (TCGA) and our previous whole-exome sequencing study and summarized the mutation, methylation, and copy number status of PTPRD in HNSCC and other cancers. In vitro studies involved standard transfection and MTT protocols, as well as methylation-specific PCR.

RESULTS

Our findings indicate that PTPRD mutation, rather than methylation or copy number alteration, is the primary mechanism by which PTPRD function is lost in HNSCC. We demonstrate that overexpression of wild-type PTPRD in HNSCC cells significantly inhibits growth and STAT3 activation while PTPRD mutants do not, suggesting that mutation may lead to loss of function and subsequent hyper-phosphorylation of PTPRD substrates, especially STAT3. Importantly, we determined that HNSCC cells harboring an endogenous PTPRD mutation are more sensitive to STAT3 blockade than PTPRD wild-type cells. We additionally found that PTPRD mRNA expression does not correlate with pSTAT3 expression, suggesting that alterations that manifest through altered mRNA expression, including hypermethylation and gene copy number alterations, do not significantly contribute to STAT3 overactivation in HNSCC.

CONCLUSION

PTPRD mutation, but not methylation or copy number loss, may serve as a predictive biomarker of sensitivity to STAT3 inhibitors in HNSCC.

Canine spontaneous head and neck squamous cell carcinomas represent their human counterparts at the molecular level

Spontaneous canine head and neck squamous cell carcinoma (HNSCC) represents an excellent model of human HNSCC but is greatly understudied. To better understand and utilize this valuable resource, we performed a pilot study that represents its first genome-wide characterization by investigating 12 canine HNSCC cases, of which 9 are oral, via high density array comparative genomic hybridization and RNA-seq. The analyses reveal that these canine cancers recapitulate many molecular features of human HNSCC. These include analogous genomic copy number abnormality landscapes and sequence mutation patterns, recurrent alteration of known HNSCC genes and pathways (e.g., cell cycle, PI3K/AKT signaling), and comparably extensive heterogeneity. Amplification or overexpression of protein kinase genes, matrix metalloproteinase genes, and epithelial-mesenchymal transition genes TWIST1 and SNAI1 are also prominent in these canine tumors. This pilot study, along with a rapidly growing body of literature on canine cancer, reemphasizes the potential value of spontaneous canine cancers in HNSCC basic and translational research.