Identifying and targeting determinants of melanoma cellular invasion

Dangerous liaisons: Loss of keratinocyte control over melanocytes in melanomagenesis

Melanomas arise from transformed melanocytes, positioned at the dermal-epidermal junction in the basal layer of the epidermis. Melanocytes are completely surrounded by keratinocyte neighbors, with which they communicate through direct contact and paracrine signaling to maintain normal growth control and homeostasis. UV radiation from sunlight reshapes this communication network to drive a protective tanning response. However, repeated rounds of sun exposure result in accumulation of mutations in melanocytes that have been considered as primary drivers of melanoma initiation and progression. It is now clear that mutations in melanocytes are not sufficient to drive tumor formation-the tumor environment plays a critical role. This review focuses on changes in melanocyte-keratinocyte communication that contribute to melanoma initiation and progression, with a particular focus on recent mechanistic insights that lay a foundation for developing new ways to intercept melanoma development.

Melanoma redox biology and the emergence of drug resistance

Melanoma is the deadliest form of skin cancer, with the loss of approximately 60,000 lives world-wide each year. Despite the development of targeted therapeutics, including compounds that have selectivity for mutant oncoproteins expressed only in cancer cells, many patients are either unresponsive to initial therapy or their tumors acquire resistance. This results in five-year survival rates of below 25%. New strategies that either kill drug-resistant melanoma cells or prevent their emergence would be extremely valuable. Melanoma, like other cancers, has long been described as being under increased oxidative stress, resulting in an increased reliance on antioxidant defense systems. Changes in redox homeostasis are most apparent during metastasis and during the metabolic reprogramming associated with the development of treatment resistance. This review discusses oxidative stress in melanoma, with a particular focus on targeting antioxidant pathways to limit the emergence of drug resistant cells.

Integrative analysis of bulk and single-cell RNA-seq reveals the molecular characterization of the immune microenvironment and oxidative stress signature in melanoma

Background

The immune microenvironment and oxidative stress of melanoma show significant heterogeneity, which affects tumor growth, invasion and treatment response. Single-cell and bulk RNA-seq data were used to explore the heterogeneity of the immune microenvironment and oxidative stress of melanoma.

Methods

The R package Seurat facilitated the analysis of the single-cell dataset, while Harmony, another R package, was employed for batch effect correction. Cell types were classified using Uniform Manifold Approximation and Projection (UMAP). The Secreted Signaling algorithm from CellChatDB.human was applied to elucidate cell-to-cell communication patterns within the single-cell data. Consensus clustering analysis for the skin cutaneous melanoma (SKCM) samples was executed with the R package ConsensusClusterPlus. To quantify immune infiltrating cells, we utilized CIBERSORT, ESTIMATE, and TIMERxCell algorithms provided by the R package Immuno-Oncology Biological Research (IOBR). Single nucleotide variant (SNV) analysis was conducted using Maftools, an R package specifically designed for this purpose. Subsequently, the expression levels of PXDN and PAPSS2 genes were assessed in melanoma tissues compared to adjacent normal tissues. Furthermore, in vitro experiments were conducted to evaluate the proliferation and reactive oxygen species expression in melanoma cells following transfection with siRNA targeting PXDN and PAPSS2.

Results

Malignant tumor cell populations were reclassified based on a comprehensive single-cell dataset analysis, which yielded six distinct tumor subsets. The specific marker genes identified for these subgroups were then used to interrogate the Cancer Genome Atlas Skin Cutaneous Melanoma (TCGA-SKCM) cohort, derived from bulk RNA sequencing data, resulting in the delineation of two immune molecular subtypes. Notably, patients within the cluster2 (C2) subtype exhibited a significantly more favorable prognosis compared to those in the cluster1 (C1) subtype. An alignment of immune characteristics was observed between the C2 subtype and unique immune functional tumor cell subsets. Genes differentially expressed across these subtypes were subsequently leveraged to construct a predictive risk model. In vitro investigations further revealed elevated expression levels of PXDN and PAPSS2 in melanoma tissue samples. Functional assays indicated that modulation of PXDN and PAPSS2 expression could influence the production of reactive oxygen species (ROS) and the proliferative capacity of melanoma cells.

Conclusion

The constructed six-gene signature can be used as an immune response and an oxidative stress marker to guide the clinical diagnosis and treatment of melanoma.

Methylation of the epigenetic JMJD2D protein by SET7/9 promotes prostate tumorigenesis

How the function of the JMJD2D epigenetic regulator is regulated or whether it plays a role in prostate cancer has remained elusive. We found that JMJD2D was overexpressed in prostate tumors, stimulated prostate cancer cell growth and became methylated by SET7/9 on K427. Mutation of this lysine residue in JMJD2D reduced the ability of DU145 prostate cancer cells to grow, invade and form tumors and elicited extensive transcriptomic changes. This included downregulation of CBLC, a ubiquitin ligase gene with hitherto unknown functions in prostate cancer, and upregulation of PLAGL1, a transcription factor with reported tumor suppressive characteristics in the prostate. Bioinformatic analyses indicated that CBLC expression was elevated in prostate tumors. Further, downregulation of CBLC largely phenocopied the effects of the K427 mutation on DU145 cells. In sum, these data have unveiled a novel mode of regulation of JMJD2D through lysine methylation, illustrated how this can affect oncogenic properties by influencing expression of the CBLC gene, and established a pro-tumorigenic role for CBLC in the prostate. A corollary is that JMJD2D and CBLC inhibitors could have therapeutic benefits in the treatment of prostate and possibly other cancers.

The role of peroxidasin in solid cancer progression

Peroxidasin is a heme-containing peroxidase enzyme that plays a vital role in the cross-linking of collagen IV molecules in basement membranes. Collagen IV cross-links are essential for providing structure and mechanical stability throughout tissue development, homeostasis, and wound healing. During cancer progression, the basement membrane is degraded, and proteins typically found in the basement membrane, including peroxidasin and collagen IV, can be found spread throughout the tumour microenvironment where they interact with cancer cells and alter cell behaviour. Whilst peroxidasin is reported to be up-regulated in a number of different cancers, the role that it plays in disease progression and metastasis has only recently begun to be studied. This review highlights the current literature exploring the known roles of peroxidasin in normal tissues and cancer progression, regulators of peroxidasin expression, and the reported relationships between peroxidasin expression and patient outcome in cancer.

GLIS Family Zinc Finger 3 Promotes Triple-Negative Breast Cancer Progression by Inducing Cell Proliferation, Migration and Invasion, and Activating the NF-κB Signaling Pathway

Triple-negative breast cancer (TNBC) puts a great threat to women's health. GLIS family zinc finger 3 (GLIS3) belongs to the GLI transcription factor family and acts as a critical factor in cancer progression. Nevertheless, the part of GLIS3 played in TNBC is not known. Immunohistochemical (IHC) staining analysis displayed that GLIS3 was highly expressed in TNBC tissues. The effect of GLIS3 on the malignant phenotype of TNBC was tested in two different cell lines according to GLIS3 regulation. Upregulation of GLIS3 promoted the proliferation, migration, and invasion of TNBC cell lines, whereas the knockdown of GLIS3 suppressed these tumor activities. Inhibition of GLIS3 induced TNBC cell apoptosis. Furthermore, study as immunofluorescence and electrophoretic mobility shift assay confirmed that the nuclear factor-κB (NF-κB) signaling pathway activated by GLIS3 played an important role in TNBC cells' malignant phenotype. In conclusion, the present work demonstrated that GLIS3 acts as a crucial element in TNBC progression via activating the NF-κB signaling pathway. Accordingly, above mentioned findings indicated that modulation of GLIS3 expression is a potential tactic to interfere with the progression of TNBC.

A Composite Bioinformatic Analysis to Explore Endoplasmic Reticulum Stress-Related Prognostic Marker and Potential Pathogenic Mechanisms in Glioma by Integrating Multiomics Data

In recent years, abnormal endoplasmic reticulum stress (ERS) response, as an important regulator of immunity, may play a vital role in the occurrence, development, and treatment of glioma. Weighted correlation network analysis (WGCNA) based on six glioma datasets was used to screen eight prognostic-related differentially expressed ERS-related genes (PR-DE-ERSGs) and to construct a prognostic model. BMP2 and HEY2 were identified as protective factors (HR < 1), and NUP107, DRAM1, F2R, PXDN, RNF19A, and SCG5 were identified as risk factors for glioma (HR > 1). QRT-PCR further supported significantly higher DRAM1 and lower SCG5 relative mRNA expression in gliomas. Our model has demonstrated excellent performance in predicting the prognosis of glioma patients from numerous datasets. In addition, the model shows good stability in multiple tests. Our model also shows broad clinical promise in predicting drug treatment effects. More immune cells/processes in the high-risk population with poor prognosis illustrate the importance of the tumor immunosuppressive environment in glioma. The potential role of the HEY2-based competitive endogenous RNA (ceRNA) regulatory network in glioma was validated and revealed the possible important role of glycolysis in glioma ERS. IDH1 and TP53 mutations with better prognosis were strongly associated with the risk score and PR-DE-ERSGs expression in the model. mDNAsi was also closely related to the risk score and clinical characteristics.

A Composite Bioinformatic Analysis to Explore Endoplasmic Reticulum Stress-Related Prognostic Marker and Potential Pathogenic Mechanisms in Glioma by Integrating Multiomics Data

In recent years, abnormal endoplasmic reticulum stress (ERS) response, as an important regulator of immunity, may play a vital role in the occurrence, development, and treatment of glioma. Weighted correlation network analysis (WGCNA) based on six glioma datasets was used to screen eight prognostic-related differentially expressed ERS-related genes (PR-DE-ERSGs) and to construct a prognostic model. BMP2 and HEY2 were identified as protective factors (HR < 1), and NUP107, DRAM1, F2R, PXDN, RNF19A, and SCG5 were identified as risk factors for glioma (HR > 1). QRT-PCR further supported significantly higher DRAM1 and lower SCG5 relative mRNA expression in gliomas. Our model has demonstrated excellent performance in predicting the prognosis of glioma patients from numerous datasets. In addition, the model shows good stability in multiple tests. Our model also shows broad clinical promise in predicting drug treatment effects. More immune cells/processes in the high-risk population with poor prognosis illustrate the importance of the tumor immunosuppressive environment in glioma. The potential role of the HEY2-based competitive endogenous RNA (ceRNA) regulatory network in glioma was validated and revealed the possible important role of glycolysis in glioma ERS. IDH1 and TP53 mutations with better prognosis were strongly associated with the risk score and PR-DE-ERSGs expression in the model. mDNAsi was also closely related to the risk score and clinical characteristics.

A Composite Bioinformatic Analysis to Explore Endoplasmic Reticulum Stress-Related Prognostic Marker and Potential Pathogenic Mechanisms in Glioma by Integrating Multiomics Data

In recent years, abnormal endoplasmic reticulum stress (ERS) response, as an important regulator of immunity, may play a vital role in the occurrence, development, and treatment of glioma. Weighted correlation network analysis (WGCNA) based on six glioma datasets was used to screen eight prognostic-related differentially expressed ERS-related genes (PR-DE-ERSGs) and to construct a prognostic model. BMP2 and HEY2 were identified as protective factors (HR < 1), and NUP107, DRAM1, F2R, PXDN, RNF19A, and SCG5 were identified as risk factors for glioma (HR > 1). QRT-PCR further supported significantly higher DRAM1 and lower SCG5 relative mRNA expression in gliomas. Our model has demonstrated excellent performance in predicting the prognosis of glioma patients from numerous datasets. In addition, the model shows good stability in multiple tests. Our model also shows broad clinical promise in predicting drug treatment effects. More immune cells/processes in the high-risk population with poor prognosis illustrate the importance of the tumor immunosuppressive environment in glioma. The potential role of the HEY2-based competitive endogenous RNA (ceRNA) regulatory network in glioma was validated and revealed the possible important role of glycolysis in glioma ERS. IDH1 and TP53 mutations with better prognosis were strongly associated with the risk score and PR-DE-ERSGs expression in the model. mDNAsi was also closely related to the risk score and clinical characteristics.

Pre-diagnostic DNA methylation in blood leucocytes in cutaneous melanoma; a nested case-control study within the Norwegian Women and Cancer cohort

The prognosis of cutaneous melanoma depends on early detection, and good biomarkers for melanoma risk may provide a valuable tool to detect melanoma development at a pre-clinical stage. By studying the epigenetic profile in pre-diagnostic blood samples of melanoma cases and cancer free controls, we aimed to identify DNA methylation sites conferring melanoma risk. DNA methylation was measured at 775,528 CpG sites using the Illumina EPIC array in whole blood in incident melanoma cases (n = 183) and matched cancer-free controls (n = 183) in the Norwegian Women and Cancer cohort. Phenotypic information and ultraviolet radiation exposure were obtained from questionnaires. Epigenome wide association (EWAS) was analyzed in future melanoma cases and controls with conditional logistic regression, with correction for multiple testing using the false discovery rate (FDR). We extended the analysis by including a public data set on melanoma (GSE120878), and combining these different data sets using a version of covariate modulated FDR (AdaPT). The analysis on future melanoma cases and controls did not identify any genome wide significant CpG sites (0.85 ≤ padj ≤ 0.99). In the restricted AdaPT analysis, 7 CpG sites were suggestive at the FDR level of 0.15. These CpG sites may potentially be used as pre-diagnostic biomarkers of melanoma risk.

Measuring peroxidasin activity in live cells using bromide addition for signal amplification

Peroxidasin (PXDN) is involved in the crosslinking of collagen IV, a major constituent of basement membranes. Disruption of basement membrane integrity as observed in genetic alterations of collagen IV or PXDN can result in developmental defects and diverse pathologies. Hence, the study of PXDN activity in (patho)physiological contexts is highly relevant. So far, measurements of PXDN activity have been reported from purified proteins, cell lysates and de-cellularized extracellular matrix. Here, for the first time we report the measurement of PXDN activity in live cells using the Amplex Red assay with a signal amplifying modification. We observe that bromide addition enhances the obtained signal, most likely due to formation of HOBr. Abrogation of signal amplification by the HOBr scavenger carnosine supports this hypothesis. Both, pharmacological inhibition as well as complementary genetic approaches confirm that the obtained signal is indeed related to PXDN activity. We validate the modified assay by investigating the effect of Brefeldin A, to inhibit the secretory pathway and thus the access of PXDN to the extracellular Amplex Red dye. Our method opens up new possibilities to investigate the activity of PXDN in (patho)physiological contexts.

Uric Acid Reacts with Peroxidasin, Decreases Collagen IV Crosslink, Impairs Human Endothelial Cell Migration and Adhesion

Uric acid is considered the main substrate for peroxidases in plasma. The oxidation of uric acid by human peroxidases generates urate free radical and urate hydroperoxide, which might affect endothelial function and explain, at least in part, the harmful effects of uric acid on the vascular system. Peroxidasin (PXDN), the most recent heme-peroxidase described in humans, catalyzes the formation of hypobromous acid, which mediates collagen IV crosslinks in the extracellular matrix. This enzyme has gained increasing scientific interest since it is associated with cardiovascular disease, cancer, and renal fibrosis. The main objective here was to investigate whether uric acid would react with PXDN and compromise the function of the enzyme in human endothelial cells. Urate decreased Amplex Red oxidation and brominating activity in the extracellular matrix (ECM) from HEK293/PXDN overexpressing cells and in the secretome of HUVECs. Parallelly, urate was oxidized to 5-hydroxyisourate. It also decreased collagen IV crosslink in isolated ECM from PFHR9 cells. Urate, the PXDN inhibitor phloroglucinol, and the PXDN knockdown impaired migration and adhesion of HUVECs. These results demonstrated that uric acid can affect extracellular matrix formation by competing for PXDN. The oxidation of uric acid by PXDN is likely a relevant mechanism in the endothelial dysfunction related to this metabolite.

Differential DNA methylation analysis of SUMF2, ADAMTS5, and PXDN provides novel insights into colorectal cancer prognosis prediction in Taiwan

BACKGROUND

Patients with colorectal cancer (CRC) undergo surgery, as well as perioperative chemoradiation or adjuvant chemotherapy primarily based on the tumor–node– metastasis (TNM) cancer staging system. However, treatment responses and prognostic outcomes of patients within the same stage vary markedly. The potential use of novel biomarkers can improve prognostication and shared decision making before implementation into certain therapies.

AIM

To investigate whether SUMF2, ADAMTS5, and PXDN methylation status could be associated with CRC prognosis.

METHODS

We conducted a Taiwan region cohort study involving 208 patients with CRC recruited from Tri-Service General Hospital and applied the candidate gene approach to identify three genes involved in oncogenesis pathways. A methylation-specific polymerase chain reaction (MS-PCR) and EpiTYPER DNA methylation analysis were employed to detect methylation status and to quantify the methylation level of candidate genes in tumor tissue and adjacent normal tissue from participants. We evaluated SUMF2, ADAMTS5, and PXDN methylation as predictors of prognosis, including recurrence-free survival (RFS), progression-free survival (PFS), and overall survival (OS), using a Cox regression model and Kaplan–Meier analysis.

RESULTS

We revealed various outcomes related to methylation and prognosis. Significantly shorter PFS and OS were associated with the CpG_3+CpG_7 hypermethylation of SUMF2 from tumor tissue compared with CpG_3+CpG_7 hypomethylation [hazard ratio (HR) = 2.24, 95% confidence interval (CI) = 1.03-4.85 for PFS, HR = 2.56 and 95%CI = 1.08-6.04 for OS]. By contrast, a significantly longer RFS was associated with CpG_2 and CpG_13 hypermethylation of ADAMTS5 from normal tissue compared with CpG_2 and CpG_13 hypomethylation [HR (95%CI) = 0.15 (0.03-0.71) for CpG_2 and 0.20 (0.04-0.97) for CpG_13]. The relationship between the methylation status of PXDN and the prognosis of CRC did not reach statistical significance.

CONCLUSION

Our study found that CpG_3+CpG_7 hypermethylation of SUMF2 from tumor tissue was associated with significantly shorter PFS and OS compared with CpG_3+CpG_7 hypomethylation. CpG_2 and CpG_13 hypermethylation of ADAMTS5 from normal tissue was associated with a significantly longer RFS compared with CpG_2 and CpG_13 hypomethylation. These methylation-related biomarkers which have implications for CRC prognosis prediction may aid physicians in clinical decision-making.

Redox-Related Proteins in Melanoma Progression

Melanoma is the most aggressive type of skin cancer. Despite the available therapies, the minimum residual disease is still refractory. Reactive oxygen and nitrogen species (ROS and RNS) play a dual role in melanoma, where redox imbalance is involved from initiation to metastasis and resistance. Redox proteins modulate the disease by controlling ROS/RNS levels in immune response, proliferation, invasion, and relapse. Chemotherapeutics such as BRAF and MEK inhibitors promote oxidative stress, but high ROS/RNS amounts with a robust antioxidant system allow cells to be adaptive and cooperate to non-toxic levels. These proteins could act as biomarkers and possible targets. By understanding the complex mechanisms involved in adaptation and searching for new targets to make cells more susceptible to treatment, the disease might be overcome. Therefore, exploring the role of redox-sensitive proteins and the modulation of redox homeostasis may provide clues to new therapies. This study analyzes information obtained from a public cohort of melanoma patients about the expression of redox-generating and detoxifying proteins in melanoma during the disease stages, genetic alterations, and overall patient survival status. According to our analysis, 66% of the isoforms presented differential expression on melanoma progression: NOS2, SOD1, NOX4, PRX3, PXDN and GPX1 are increased during melanoma progression, while CAT, GPX3, TXNIP, and PRX2 are decreased. Besides, the stage of the disease could influence the result as well. The levels of PRX1, PRX5 and PRX6 can be increased or decreased depending on the stage. We showed that all analyzed isoforms presented some genetic alteration on the gene, most of them (78%) for increased mRNA expression. Interestingly, 34% of all melanoma patients showed genetic alterations on TRX1, most for decreased mRNA expression. Additionally, 15% of the isoforms showed a significant reduction in overall patient survival status for an altered group (PRX3, PRX5, TR2, and GR) and the unaltered group (NOX4). Although no such specific antioxidant therapy is approved for melanoma yet, inhibitors or mimetics of these redox-sensitive proteins have achieved very promising results. We foresee that forthcoming investigations on the modulation of these proteins will bring significant advances for cancer therapy.

Peroxidasin Enhances Basal Phenotype and Inhibits Branching Morphogenesis in Breast Epithelial Progenitor Cell Line D492

The human breast is composed of terminal duct lobular units (TDLUs) that are surrounded by stroma. In the TDLUs, basement membrane separates the stroma from the epithelial compartment, which is divided into an inner layer of luminal epithelial cells and an outer layer of myoepithelial cells. Stem cells and progenitor cells also reside within the epithelium and drive a continuous cycle of gland remodelling that occurs throughout the reproductive period. D492 is an epithelial cell line originally isolated from the stem cell population of the breast and generates both luminal and myoepithelial cells in culture. When D492 cells are embedded into 3D reconstituted basement membrane matrix (3D-rBM) they form branching colonies mimicking the TDLUs of the breast, thereby providing a well-suited in vitro model for studies on branching morphogenesis and breast development. Peroxidasin (PXDN) is a heme-containing peroxidase that crosslinks collagen IV with the formation of sulfilimine bonds. Previous studies indicate that PXDN plays an integral role in basement membrane stabilisation by crosslinking collagen IV and as such contributes to epithelial integrity. Although PXDN has been linked to fibrosis and cancer in some organs there is limited information on its role in development, including in the breast. In this study, we demonstrate expression of PXDN in breast epithelium and stroma and apply the D492 cell line to investigate the role of PXDN in cell differentiation and branching morphogenesis in the human breast. Overexpression of PXDN induced basal phenotype in D492 cells, loss of plasticity and inhibition of epithelial-to-mesenchymal transition as is displayed by complete inhibition of branching morphogenesis in 3D culture. This is supported by results from RNA-sequencing which show significant enrichment in genes involved in epithelial differentiation along with significant negative enrichment of EMT factors. Taken together, we provide evidence for a novel role of PXDN in breast epithelial differentiation and mammary gland development.

Peroxidasin protein expression and enzymatic activity in metastatic melanoma cell lines are associated with invasive potential

Peroxidasin, a heme peroxidase, has been shown to play a role in cancer progression. mRNA expression has been reported to be upregulated in metastatic melanoma cell lines and connected to the invasive phenotype, but little is known about how peroxidasin acts in cancer cells. We have analyzed peroxidasin protein expression and activity in eight metastatic melanoma cell lines using an ELISA developed with an in-house peroxidasin binding protein. RNAseq data analysis confirmed high peroxidasin mRNA expression in the five cell lines classified as invasive and low expression in the three non-invasive cell lines. Protein levels of peroxidasin were higher in the cell lines with an invasive phenotype. Active peroxidasin was secreted to the cell culture medium, where it accumulated over time, and peroxidasin protein levels in the medium were also much higher in invasive than non-invasive cell lines. The only well-established physiological role of peroxidasin is in the formation of a sulfilimine bond, which cross-links collagen IV in basement membranes via catalyzed oxidation of bromide to hypobromous acid. We found that peroxidasin secreted from melanoma cells formed sulfilimine bonds in uncross-linked collagen IV, confirming peroxidasin activity and hypobromous acid formation. Moreover, 3-bromotyrosine, a stable product of hypobromous acid reacting with tyrosine residues, was detected in invasive melanoma cells, substantiating that their expression of peroxidasin generates hypobromous acid, and showing that it does not exclusively react with collagen IV, but also with other biomolecules.

Genetic and Genomic Pathways of Melanoma Development, Invasion and Metastasis

Melanoma is a serious form of skin cancer that accounts for 80% of skin cancer deaths. Recent studies have suggested that melanoma invasiveness is attributed to phenotype switching, which is a reversible type of cell behaviour with similarities to epithelial to mesenchymal transition. Phenotype switching in melanoma is reported to be independent of genetic alterations, whereas changes in gene transcription, and epigenetic alterations have been associated with invasiveness in melanoma cell lines. Here, we review mutational, transcriptional, and epigenomic alterations that contribute to tumour heterogeneity in melanoma, and their potential to drive melanoma invasion and metastasis. We also discuss three models that are hypothesized to contribute towards aspects of tumour heterogeneity and tumour progression in melanoma, namely the clonal evolution model, the cancer stem cell model, and the phenotype switching model. We discuss the merits and disadvantages of each model in explaining tumour heterogeneity in melanoma, as a precursor to invasion and metastasis.

Basement membrane stiffness determines metastases formation

The basement membrane (BM) is a special type of extracellular matrix and presents the major barrier cancer cells have to overcome multiple times to form metastases. Here we show that BM stiffness is a major determinant of metastases formation in several tissues and identify netrin-4 (Net4) as a key regulator of BM stiffness. Mechanistically, our biophysical and functional analyses in combination with mathematical simulations show that Net4 softens the mechanical properties of native BMs by opening laminin node complexes, decreasing cancer cell potential to transmigrate this barrier despite creating bigger pores. Our results therefore reveal that BM stiffness is dominant over pore size, and that the mechanical properties of 'normal' BMs determine metastases formation and patient survival independent of cancer-mediated alterations. Thus, identifying individual Net4 protein levels within native BMs in major metastatic organs may have the potential to define patient survival even before tumour formation. The ratio of Net4 to laminin molecules determines BM stiffness, such that the more Net4, the softer the BM, thereby decreasing cancer cell invasion activity.

Identification of novel drug resistance mechanisms by genomic and transcriptomic profiling of glioblastoma cells with mutation-activated EGFR

AIMS

Epidermal growth factor receptor (EGFR) is not only involved in carcinogenesis, but also in chemoresistance. We characterized U87.MGΔEGFR glioblastoma cells with constitutively active EGFR due to deletion at the ligand binding domain in terms of gene expression profiling and chromosomal aberrations. Wild-type U87.MG cells served as control.

MATERIALS AND METHODS

RNA sequencing and network analyses (Ingenuity Pathway Analysis) were performed to identify novel drug resistance mechanisms related to expression of mutation activated EGFR. Chromosomal aberrations were characterized by multicolor fluorescence in situ hybridization (mFISH) and array comparative genomic hybridization (aCGH).

KEY FINDINGS

U87.MGΔEGFR cells presented much more chromosomal aberrations, amplifications and deletions than wild-type U87.MG cells. Still, both cell lines were near-triploid. Numerous genes were overexpressed in U87.MGΔEGFR cells, some of which have been already linked to drug resistance. PXDN, which is associated with epithelial mesenchymal transition, was the most upregulated gene (901.8-fold). TENM1 was 331.6-fold upregulated, and it was previously reported to modulate neural development. EGFR-AS1 (161.2-fold upregulated) has been reported to increase the EGFR mRNA stability and its expression - in accordance with that of EGFR - was upregulated (85.5-fold). In addition to well-known resistance genes, numerous novel genes and genomic aberrations were identified. ANGPT2 upregulation and CPM downregulation were validated by Western blotting.

SIGNIFICANCE

Transcriptomics and genomics analyses in U87.MGΔEGFR cells unraveled a range of novel drug resistance mechanisms including apoptosis, DNA repair, ferroptosis, glutathione related gene activities, heat shock, oxidative stress, transcription factor activities, which may have important implications for future treatment strategies.

Genome-wide DNA methylation and RNA expression differences correlate with invasiveness in melanoma cell lines

Aims & objectives: The aim of this study was to investigate the role of DNA methylation in invasiveness in melanoma cells. Materials & methods: The authors carried out genome-wide transcriptome (RNA sequencing) and reduced representation bisulfite sequencing methylome profiling between noninvasive (n = 4) and invasive melanoma cell lines (n = 5). Results: The integration of differentially expressed genes and differentially methylated fragments (DMFs) identified 12 DMFs (two in AVPI1, one in HMG20B, two in BCL3, one in NTSR1, one in SYNJ2, one in ROBO2 and four in HORMAD2) that overlapped with either differentially expressed genes (eight DMFs and six genes) or cis-targets of lncRNAs (five DMFs associated with cis-targets and four differentially expressed lncRNAs). Conclusions: DNA methylation changes are associated with a number of transcriptional differences observed in noninvasive and invasive phenotypes in melanoma.

Integrated analysis of circRNA-miRNA-mRNA network reveals potential prognostic biomarkers for radiotherapies with X-rays and carbon ions in non-small cell lung cancer

Background

This work was aimed at exploring the regulatory network of non-coding RNA (ncRNA) especially circular RNA (circRNA) and microRNA (miRNA), in the sensitivity of non-small cell lung cancer (NSCLC) cells to low linear energy transfer (LET) X-ray and high-LET carbon ion irradiations.

Methods

The radioresistant NSCLC cell line A549-R11 was obtained from its parental cell line A549 through irradiation with X-rays of 2.0 Gy per fraction for 30 times. The sensitivities of A549, A549-R11 and H1299 cells exposed to X-rays and carbon ions were verified using the colony formation assay. A comprehensive circRNA-miRNA-mRNA network was constructed through the sequencing data in parental A549, acquired radioresistant A549-R11 and intrinsic radioresistant H1299 cells, and the network was further optimized according to the prognostic results from the TCGA and GEO databases.

Results

Based on high-throughput sequencing of circRNAs, we found that 40 circRNAs were up-regulated while 184 circRNAs were down-regulated in the intersection of the sets of A549-R11 and H1299 cells. Subsequently, a circRNA- miRNA-mRNA network, including 14 interactive pairs and 8 circRNAs, 4 overall survival-associated miRNAs, and 4 mRNAs, was constructed through the high-throughput data screening and bioinformatics methods.

Conclusions

Our results provide a complete understanding to the regulatory mechanism of the sensitivities to low-LET X-ray and high-LET carbon ion irradiations, and might be helpful to screen potential biomarkers for predicting the Carbon-ion radiotherapy (CIRT) and X-ray radiotherapy responses in NSCLC.

eQTL Colocalization Analyses Identify NTN4 as a Candidate Breast Cancer Risk Gene

Breast cancer genome-wide association studies (GWASs) have identified 150 genomic risk regions containing more than 13,000 credible causal variants (CCVs). The CCVs are predominantly noncoding and enriched in regulatory elements. However, the genes underlying breast cancer risk associations are largely unknown. Here, we used genetic colocalization analysis to identify loci at which gene expression could potentially explain breast cancer risk phenotypes. Using data from the Breast Cancer Association Consortium (BCAC) and quantitative trait loci (QTL) from the Genotype-Tissue Expression (GTEx) project and The Cancer Genome Project (TCGA), we identify shared genetic relationships and reveal novel associations between cancer phenotypes and effector genes. Seventeen genes, including NTN4, were identified as potential mediators of breast cancer risk. For NTN4, we showed the rs61938093 CCV at this region was located within an enhancer element that physically interacts with the NTN4 promoter, and the risk allele reduced NTN4 promoter activity. Furthermore, knockdown of NTN4 in breast cells increased cell proliferation in vitro and tumor growth in vivo. These data provide evidence linking risk-associated variation to genes that may contribute to breast cancer predisposition.

Aberrant pregnancy-associated plasma protein-A expression in breast cancers prognosticates clinical outcomes

Elevated levels of pregnancy-associated plasma protein-A (PAPP-A) have been implicated in the pathogenesis of various malignancies, including breast cancers. Breast cancer is one of the most frequent carcinomas and is the second most common cancer type detected in women of child-bearing age. Throughout pregnancy PAPP-A is produced and secreted by the placental syncytiotrophoblast cells; co-incidentally pregnancy-associated breast cancers often have an aggressive clinical course. The components of the PAPP-A/IGF axis was assessed in a panel of breast cancer cell lines. Using neutralising antibodies the impact of PAPP-A/IGF axis on cell motility was evaluated. PAPP-A was expressed in four of the twelve breast cancer cell lines tested. Blocking PAPP-A and IGFBP4 with neutralising antibodies significantly decreased motiliy of MDA-MB-231 cells. Upregulation of PAPP-A expression in breast tumours resulted in a trend towards worse overall survival. Notably, PAPP-A expression also positively correlated with epithelial-to-mesenchymal transition markers. In conclusion, these results indicate that PAPP-A plays an important role in breast cancer progression and it may be a promising therapeutic target in breast cancer.

Peroxidan Plays a Tumor-Promoting Role in Oral Squamous Cell Carcinoma

Despite dramatic progress in cancer diagnosis and treatment, the five-year survival rate of oral squamous cell carcinoma (OSCC) is still only about 50%. Thus, the need for elucidating the molecular mechanisms underlying OSCC is urgent. We previously identified the peroxidasin gene (PXDN) as one of several novel genes associated with OSCC. Although the PXDN protein is known to act as a tumor-promoting factor associated with the Warburg effect, its function and role in OSCC are poorly understood. In this study, we investigated the expression, function, and relationship with the Warburg effect of PXDN in OSCC. In immunohistochemical analysis of OSCC specimens, we observed that elevated PXDN expression correlated with lymph node metastasis and a diffuse invasion pattern. High PXDN expression was confirmed as an independent predictor of poor prognosis by multivariate analysis. The PXDN expression level correlated positively with that of pyruvate kinase (PKM2) and heme oxygenase-1 (HMOX1) and with lactate and ATP production. No relationship between PXDN expression and mitochondrial activation was observed, and PXDN expression correlated inversely with reactive oxygen species (ROS) production. These results suggest that PXDN might be a tumor progression factor causing a Warburg-like effect in OSCC.

Peroxidasin mediates bromination of tyrosine residues in the extracellular matrix

Peroxidasin is a heme peroxidase that oxidizes bromide to hypobromous acid (HOBr), a powerful oxidant that promotes the formation of the sulfilimine crosslink in collagen IV in basement membranes. We investigated whether HOBr released by peroxidasin leads to other oxidative modifications of proteins, particularly bromination of tyrosine residues, in peroxidasin-expressing PFHR9 cells. Using stable isotope dilution LC-MS/MS, we detected the formation of 3-bromotyrosine, a specific biomarker of HOBr-mediated protein modification. The level of 3-bromotyrosine in extracellular matrix proteins from normally cultured cells was 1.1 mmol/mol tyrosine and decreased significantly in the presence of the peroxidasin inhibitor, phloroglucinol. A negligible amount of 3-bromotyrosine was detected in peroxidasin-knockout cells. 3-Bromotyrosine formed both during cell growth in culture and in the isolated decellularized extracellular matrix when embedded peroxidasin was supplied with hydrogen peroxide and bromide. The level of 3-bromotyrosine was significantly higher in extracellular matrix than intracellular proteins, although a low amount was detected intracellularly. 3-Bromotyrosine levels increased with higher bromide concentrations and decreased in the presence of physiological concentrations of thiocyanate and urate. However, these peroxidase substrates showed moderate to minimal inhibition of collagen IV crosslinking. Our findings provide evidence that peroxidasin promotes the formation of 3-bromotyrosine in proteins. They show that HOBr produced by peroxidasin is selective for, but not limited to, the crosslinking of collagen IV. Based on our findings, the use of 3-bromotyrosine as a specific biomarker of oxidative damage by HOBr warrants further investigation in clinical conditions linked to high peroxidasin expression.

Proteomic analysis enables distinction of early- versus advanced-stage lung adenocarcinomas

Background

A gel‐free proteomic approach was utilized to perform in‐depth tissue protein profiling of lung adenocarcinoma (ADC) and normal lung tissues from early and advanced stages of the disease. The long‐term goal of this study is to generate a large‐scale, label‐free proteomics dataset from histologically well‐classified lung ADC that can be used to increase further our understanding of disease progression and aid in identifying novel biomarkers.

Methods and results

Cases of early‐stage (I‐II) and advanced‐stage (III‐IV) lung ADCs were selected and paired with normal lung tissues from 22 patients. The histologically and clinically stratified human primary lung ADCs were analyzed by liquid chromatography‐tandem mass spectrometry. From the analysis of ADC and normal specimens, 4863 protein groups were identified. To examine the protein expression profile of ADC, a peak area‐based quantitation method was used. In early‐ and advanced‐stage ADC, 365 and 366 proteins were differentially expressed, respectively, between normal and tumor tissues (adjusted P‐value < .01, fold change ≥ 4). A total of 155 proteins were dysregulated between early‐ and advanced‐stage ADCs and 18 were suggested as early‐specific stage ADC. In silico functional analysis of the upregulated proteins in both tumor groups revealed that most of the enriched pathways are involved in mRNA metabolism. Furthermore, the most overrepresented pathways in the proteins that were unique to ADC are related to mRNA metabolic processes.

Conclusions

Further analysis of these data may provide an insight into the molecular pathways involved in disease etiology and may lead to the identification of biomarker candidates and potential targets for therapy. Our study provides potential diagnostic biomarkers for lung ADC and novel stage‐specific drug targets for rational intervention.

Human peroxidasin 1 promotes angiogenesis through ERK1/2, Akt, and FAK pathways

Aims

The term angiogenesis refers to sprouting of new blood vessels from pre-existing ones. The angiogenic process involves cell migration and tubulogenesis requiring interaction between endothelial cells and the extracellular matrix. Human peroxidasin 1 (hsPxd01) is a multidomain heme peroxidase found embedded in the basement membranes. As it promotes the stabilization of extracellular matrix, we investigated its possible role in angiogenesis both in vitro and in vivo.

Methods and results

We analysed the effects of peroxidasin 1 gene silencing and supplementation by recombinant hsPxd01 in TeloHAEC endothelial cells on cell migration, tubulogenesis in matrigel, and intracellular signal transduction as assessed by kinase phosphorylation and expression of pro-angiogenic genes as measured by qRT-PCR. We further evaluated the angiogenic potential of recombinant peroxidasin in a chicken chorioallantoic membrane model. RNA silencing of endogenous hsPxd01 significantly reduced tube formation and cell migration, whereas supplementation by the recombinant peroxidase promoted tube formation in vitro and stimulated vascularization in vivo through its catalytic activity. Moreover, recombinant hsPxd01 promoted phosphorylation of Extracellular signal-Regulated Kinases (ERK1/2), Protein kinase B (Akt), and Focal Adhesion Kinase (FAK), and induced the expression of pro-angiogenic downstream genes: Platelet Derived Growth Factor Subunit B (PDGFB), endothelial-derived Heparin Binding EGF-like growth factor (HB-EGF), CXCL-1, Hairy-Related Transcription Factor 1 (HEY-1), DNA-binding protein inhibitor (ID-2), Snail Family Zinc Finger 1 (SNAI-1), as well as endogenous hsPxd01. However, peroxidasin silencing significantly reduced Akt and FAK phosphorylation but induced ERK1/2 activation after supplementation by recombinant hsPxd01. While hsPxd01 silencing significantly reduced expression of HEY-1, ID-2, and PDGFB, it did not affect expression of SNAI-1, HB-EGF, and CXCL-1 after supplementation by recombinant hsPxd01.

Conclusion

Our findings suggest a role of enzymatically active peroxidasin 1 as a pro-angiogenic peroxidase and a modulator of ERK1/2, Akt and FAK signalling.

Germline Mutations in Familial Papillary Thyroid Cancer

Thyroid cancer, predominantly of papillary histology (PTC), is a common cancer mostly diagnosed sporadically. Hereditary PTC is encountered in ~ 5% of cases and may present at an earlier age, with greater risks of metastasis and recurrence, compared with sporadic cases. The molecular basis of hereditary PTC is unknown in most cases. In this study, the genetic basis of hereditary PTC in three Brazilian families was investigated. Whole exome sequencing (WES) was carried out for probands in each family, and validated, pathogenic/likely pathogenic sequence variants (P/LPSVs) were genotyped in additional family members to establish their putative pathogenic role. Overall, seven P/LPSVs in seven novel genes were detected: p.D283N*ANXA3, p.Y157S*NTN4, p.G172W*SERPINA1, p.G188S*FKBP10, p.R937C*PLEKHG5, p.L32Q*P2RX5, and p.Q76*SAPCD1. These results indicate that these novel genes are seemingly associated with hereditary PTC, but extension and validation in other PTC families are required.

Emerging Roles of GLI-Similar Krüppel-like Zinc Finger Transcription Factors in Leukemia and Other Cancers

GLI-similar 1-3 (GLIS1-3), a subfamily of Krüppel-like zinc finger transcription factors, function as key regulators of several biological processes important to oncogenesis, including control of cell proliferation, differentiation, self-renewal, and epithelial-mesenchymal transition. This review provides a short overview of the critical roles genetic changes in GLIS1-3 play in the development of several malignancies. This includes intrachromosomal translocations involving GLIS2 and ETO2/CBFA2T3 in the development of pediatric non-Down's syndrome (DS), acute megakaryoblastic leukemia (AMKL), a malignancy with poor prognosis, and an association of interchromosomal translocations between GLIS3, GLIS1, and PAX8, and between GLIS3 and CLPTM1L with hyalinizing trabecular tumors (HTTs) and fibrolamellar hepatocellular carcinoma (FHCC), respectively. Targeting upstream signaling pathways that regulate GLIS signaling may offer new therapeutic strategies in the management of cancer.

Proteomics-Metabolomics Combined Approach Identifies Peroxidasin as a Protector against Metabolic and Oxidative Stress in Prostate Cancer

Peroxidasin (PXDN), a human homolog of Drosophila PXDN, belongs to the family of heme peroxidases and has been found to promote oxidative stress in cardiovascular tissue, however, its role in prostate cancer has not been previously elucidated. We hypothesized that PXDN promotes prostate cancer progression via regulation of metabolic and oxidative stress pathways. We analyzed PXDN expression in prostate tissue by immunohistochemistry and found increased PXDN expression with prostate cancer progression as compared to normal tissue or cells. PXDN knockdown followed by proteomic analysis revealed an increase in oxidative stress, mitochondrial dysfunction and gluconeogenesis pathways. Additionally, Liquid Chromatography with tandem mass spectrometry (LC-MS/MS)-based metabolomics confirmed that PXDN knockdown induced global reprogramming associated with increased oxidative stress and decreased nucleotide biosynthesis. We further demonstrated that PXDN knockdown led to an increase in reactive oxygen species (ROS) associated with decreased cell viability and increased apoptosis. Finally, PXDN knockdown decreased colony formation on soft agar. Overall, the data suggest that PXDN promotes progression of prostate cancer by regulating the metabolome, more specifically, by inhibiting oxidative stress leading to decreased apoptosis. Therefore, PXDN may be a biomarker associated with prostate cancer and a potential therapeutic target.

Silencing of carbonic anhydrase I enhances the malignant potential of exosomes secreted by prostatic tumour cells

We report results showing that the silencing of carbonic anhydrase I (siCA1) in prostatic (PC3) tumour cells has a significant impact on exosome formation. An increased diameter, concentration and diversity of the produced exosomes were noticed as a consequence of this knock‐down. The protein composition of the exosomes' cargo was also altered. Liquid chromatography and mass spectrometry analyses identified 42 proteins significantly altered in PC3 siCA1 exosomes compared with controls. The affected proteins are mainly involved in metabolic processes, biogenesis, cell component organization and defense/immunity. Interestingly, almost all of them have been described as ‘enhancers' of tumour development through the promotion of cell proliferation, migration and invasion. Thus, our results indicate that the reduced expression of the CA1 protein enhances the malignant potential of PC3 cells.

Netrin Family: Role for Protein Isoforms in Cancer

Netrins form a family of secreted and membrane-associated proteins. Netrins are involved in processes for axonal guidance, morphogenesis, and angiogenesis by regulating cell migration and survival. These processes are of special interest in tumor biology. From the netrin genes various isoforms are translated and regulated by alternative splicing. We review here the diversity of isoforms of the netrin family members and their known and potential roles in cancer.

MiR-203a is differentially expressed during branching morphogenesis and EMT in breast progenitor cells and is a repressor of peroxidasin

MicroRNAs regulate developmental events such as branching morphogenesis, epithelial to mesenchymal transition (EMT) and its reverse process mesenchymal to epithelial transition (MET). In this study, we performed small RNA sequencing of a breast epithelial progenitor cell line (D492), and its mesenchymal derivative (D492M) cultured in three-dimensional microenvironment. Among the most downregulated miRNAs in D492M was miR-203a, a miRNA that plays an important role in epithelial differentiation. Increased expression of miR-203a was seen in D492, concomitant with increased complexity of branching. When miR-203a was overexpressed in D492M, a partial reversion towards epithelial phenotype was seen. Gene expression analysis of D492M and D492MmiR-203a revealed peroxidasin, a collagen IV cross-linker, as the most significantly downregulated gene in D492MmiR-203a. Collectively, we demonstrate that miR-203a expression temporally correlates with branching morphogenesis and is suppressed in D492M. Overexpression of miR-203a in D492M induces a partial MET and reduces the expression of peroxidasin. Furthermore, we demonstrate that miR-203a is a novel repressor of peroxidasin. MiR-203-peroxidasin axis may be an important regulator in branching morphogenesis, EMT/MET and basement membrane remodeling.

Peroxidasin is a novel target of the redox-sensitive transcription factor Nrf2

Peroxidasin (PXDN) facilitates peroxidative reactions via utilisation of hydrogen peroxide (H2O2) and has been shown to crosslink collagen IV through sulfilimine bond formation in the presence of hypohalous acids. Aberrant PXDN expression has been associated with kidney fibrosis, cancer, congenital eye defects and various cardiovascular disorders. Since PXDN expression is modified by H2O2, we hypothesized that a major antioxidant response transcription factor, nuclear factor erythroid 2-related factor 2 (Nrf2), may regulate PXDN expression. PXDN expression in response to H2O2 and the Nrf2-specific inducers, tert-butylhydroquinone (tBHQ) and sulforaphane (SFN), was determined by western blotting and immunofluorescence microscopy, in HeLa and HEK293 cells. Chromatin immunoprecipitation and luciferase reporter assays were used to investigate the regulation of PXDN by Nrf2. We observed elevated Nrf2 nuclear translocation and increased PXDN protein expression in response to H2O2, tBHQ and SFN, in both cell lines. We found that Nrf2 binds to and increases luciferase reporter gene expression from the PXDN promoter via a putative Nrf2-binding site. In summary, we show that PXDN is a novel target of the redox sensitive transcription factor Nrf2. This finding further highlights the role of PXDN in redox-related processes and compliments the currently understood pathophysiological functions of PXDN.

Emerging roles of GLIS3 in neonatal diabetes, type 1 and type 2 diabetes

GLI-similar 3 (GLIS3), a member of the Krüppel-like zinc finger protein subfamily, is predominantly expressed in the pancreas, thyroid and kidney. Glis3 mRNA can be initially detected in mouse pancreas at embryonic day 11.5 and is largely restricted to β cells, pancreatic polypeptide-expressing cells, as well as ductal cells at later stage of pancreas development. Mutations in GLIS3 cause a neonatal diabetes syndrome, characterized by neonatal diabetes, congenital hypothyroidism and polycystic kidney. Importantly, genome-wide association studies showed that variations of GLIS3 are strongly associated with both type 1 diabetes (T1D) and type 2 diabetes (T2D) in multiple populations. GLIS3 cooperates with pancreatic and duodenal homeobox 1 (PDX1), v-maf musculoaponeurotic fibrosarcoma oncogene family, protein A (MAFA), as well as neurogenic differentiation 1 (NEUROD1) and potently controls insulin gene transcription. GLIS3 also plays a role in β cell survival and likely in insulin secretion. Any perturbation of these functions may underlie all three forms of diabetes. GLIS3, synergistically with hepatocyte nuclear factor 6 (HNF6) and forkhead box A2 (FOXA2), controls fetal islet differentiation via transactivating neurogenin 3 (NGN3) and impairment of this function leads to neonatal diabetes. In addition, GLIS3 is also required for the compensatory β cell proliferation and mass expansion in response to insulin resistance, which if disrupted may predispose to T2D. The increasing understanding of the mechanisms of GLIS3 in β cell development, survival and function maintenance will provide new insights into disease pathogenesis and potential therapeutic target identification to combat diabetes.