Melanoblast transcriptome analysis reveals pathways promoting melanoma metastasis

A multi-omics analysis reveals KDELR1 promotes malignant progression and correlates with tumor microenvironment in head and neck squamous cell carcinoma

KDELR1, a constituent of the KDEL endoplasmic reticulum protein retention receptors family, is implicated in immune responses and cancers progression. In this study, we delineate the clinicopathological significance and oncogenic role of KDELR1 in head and neck squamous cell carcinoma (HNSCC) through a comprehensive multi-omics approach. KDELR1 expression is correlated with tumor grade, tumor stage, lymph node metastasis, clinical stage and poor prognosis in HNSCC. Moreover, our results indicate a marked upregulation of KDELR1 expression in primary tumor tissues compared to normal and dysplasia tissues. Furthermore, increased KDELR1 expression is associated with tumor progression and indicates an unfavorable prognosis in HNSCC. Functional enrichment analysis highlights the involvement of KDELR1 in HNSCC malignant progression. Depletion of KDELR1 inhibits proliferation, stemness, migration and invasion in HNSCC cells. Bioinformatics analysis results indicate that KDELR1 promotes HNSCC progression with regulating wnt signaling pathway and CTNNB1 expression. Of note, KDELR1 is associated with HNSCC tumor microenvironment compositions, especially positively correlating with cancer associated fibroblasts and negatively correlating with CD8+ T cells and B cells infiltration. Collectively, these findings indicate that KDELR1 as an oncogene in driving progression and correlates with tumor microenvironment, suggesting its potential as a promising biomarker and a therapeutic target in HNSCC.

MCAM Expression Facilitates Melanoma-Endothelial Interactions and Promotes Metastatic Disease Progression

Invasive growth and metastatic dissemination represent the primary cause of death in cancer patients. In order to successfully detach from the primary tumour and establish metastases in distant tissues, cancer cells need to dynamically rewire their cell adhesion machinery. Here we revisit the potential association of MCAM, a member of the immunoglobulin superfamily that was initially identified as a melanoma antigen, with disease progression. Using immunohistochemical stainings and bioinformatic analyses of published datasets, we find abundant MCAM expression both in primary and metastatic human melanomas. In additional bioinformatic analyses, we show that MCAM is highly expressed in foetal melanocytes and subsequently downregulated during melanocyte maturation. Bioinformatic inference of cellular communication networks reveals that melanoma cells with high MCAM expression more actively engage in signalling crosstalk with endothelial cells. Experimental investigations demonstrate that disruption of MCAM in melanoma cells inhibits their migration on endothelial cell surfaces in vitro and decreases their ability to develop spontaneous lung metastases in vivo. Taken together, our results could not confirm the notion that MCAM expression represents a useful biomarker for disease progression but provide evidence that MCAM expression might represent part of a reactivated embryonal transcriptional program that facilitates melanoma-endothelial cell interactions during metastatic progression.

KDELR3 overexpression as a novel prognostic and diagnostic biomarker in glioma: comprehensive bioinformatic analysis insights

Glioma, the most prevalent primary intracranial malignancy among adults, is distinguished by its high morbidity and recurrence rates, posing a considerable threat to patients' quality of life and survival prospects. Consequently, the pursuit of efficacious molecular prognostic markers holds paramount importance. The exploration of the role of the KDELR3 kinase family in various neoplastic conditions constitutes an emerging area of research. However, the biological functions of KDELR3 and its prognostic implications in brain tumors remain largely undocumented. This study endeavored to ascertain the potential of KDELR3 as a prognostic indicator for glioma. We integrated a comprehensive dataset encompassing 1127 glioma samples, sourced from our cohort, The Cancer Genome Atlas (TCGA), and the Chinese Glioma Genome Atlas (CGGA), to delve into the expression patterns of KDELR3 in glioma and their associated implications. Notably, KDELR3 was markedly overexpressed in glioma and demonstrated a positive correlation with clinical progression. By utilizing Kaplan-Meier survival analysis and the Cox proportional hazards regression model, we evaluated the prognostic significance of KDELR3, revealing it as an independent predictor of adverse outcomes in glioma patients. Furthermore, gene set enrichment analysis unveiled potential signaling pathways associated with KDELR3 expression in glioma, primarily encompassing Cytokine-cytokine receptor interaction, extracellular matrix (ECM)-receptor interaction, and complement and coagulation cascades. In summation, our findings provide profound insights into the potential role of KDELR3 and its application as a novel and promising prognostic biomarker for glioma.

Melanocyte lineage dynamics in development, growth and disease

Melanocytes evolved to produce the melanin that gives colour to our hair, eyes and skin. The melanocyte lineage also gives rise to melanoma, the most lethal form of skin cancer. The melanocyte lineage differentiates from neural crest cells during development, and most melanocytes reside in the skin and hair, where they are replenished by melanocyte stem cells. Because the molecular mechanisms necessary for melanocyte specification, migration, proliferation and differentiation are co-opted during melanoma initiation and progression, studying melanocyte development is directly relevant to human disease. Here, through the lens of advances in cellular omic and genomic technologies, we review the latest findings in melanocyte development and differentiation, and how these developmental pathways become dysregulated in disease.

Exploring the potential biological significance of KDELR family genes in lung adenocarcinoma

The Lys-Asp-Glu-Leu receptor (KDELR) family genes play critical roles in a variety of biological processes in different tumors. Our study aimed to provide a comprehensive analysis of the potential roles of KDELRs in lung adenocarcinoma (LUAD). Utilizing data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database, as well as clinical samples, we conducted a series of analyses and validations using R software tools and various online resources. The results showed that KDELR family genes and proteins were highly expressed and associated with a poor prognosis of LUAD. Promoter hypomethylation and the competing endogenous RNA (ceRNA) network of PCAT6/hsa-miR-326/KDELR1 might be potential causes of aberrant KDELR1 overexpression in LUAD. Three key Transcription factors (TFs) (SPI1, EP300, and MAZ) and a TFs-miRNAs-KDELRs network (involving 11 TFs) might be involved in modulating KDELRs expression abnormalities. Gene Set Enrichment Analysis (GSEA) indicated enrichment of genes highly expressing KDELR1, KDELR2, and KDELR3 in MTORC1_SIGNALING, P53_PATHWAY, and ANGIOGENESIS. Negative correlations between KDELRs expression and CD8 + T cell infiltration, as well as CTLA-4 expression. Our multiple analyses suggested that the KDELRs are important signaling molecules in LUAD. These results provided novel insights for developing prognostic markers and novel therapies of LUAD.

Molecular Markers in Melanoma Progression: A Study on the Expression of miRNA Gene Subtypes in Tumoral vs. Benign Nevi

This study investigates the differential expression of miRNA gene subtypes in tumoral versus benign nevi in individuals with melanoma, aiming to identify clinically significant correlations that could serve as reliable markers for assessing tumor stage and progression. Conducted between 2019 and 2022, this descriptive, quantitative observational research analyzed 90 formalin-fixed paraffin-embedded (FFPE) samples from the Pius Brinzeu County Emergency Clinical Hospital, Timisoara, including 45 samples of advanced-stage melanoma and 45 samples of pigmented nevi. miRNA purification and analysis were performed using the miRNeasy Kit and the Human Cancer PathwayFinder miScript miRNA PCR Array, with statistical analysis (including logistic regression) to determine associations with cancer staging, such as high Breslow index risk, number of mitoses, and vascular invasion. After the analysis and comparison of 180 miRNA gene subtypes, we selected 10 of the most upregulated and 10 most downregulated genes. The results revealed that hsa-miR-133b, hsa-miR-335-5p, hsa-miR-200a-3p, and hsa-miR-885-5p were significantly upregulated in melanoma samples, with fold changes ranging from 1.09 to 1.12. Conversely, hsa-miR-451a and hsa-miR-29b-3p showed notable downregulation in melanoma, with fold changes of 0.90 and 0.92, respectively. Additionally, logistic regression analysis identified hsa-miR-29b-3p (OR = 2.51) and hsa-miR-200a-3p (OR = 2.10) as significantly associated with an increased risk of a high Breslow index, while hsa-miR-127-3p and hsa-miR-451a were associated with a reduced risk. Conclusively, this study underscores the significant alterations in miRNA expression in melanoma compared to benign nevi and highlights the potential of specific miRNAs as biomarkers for melanoma progression. The identification of miRNAs with significant associations to melanoma characteristics suggests their utility in developing non-invasive, cost-effective diagnostic tools and in guiding therapeutic decisions, potentially improving patient outcomes in melanoma management.

Nonepithelial cancer dissemination: specificities and challenges

Epithelial cancers have served as a paradigm to study tumor dissemination but recent data have highlighted significant differences with nonepithelial cancers. Here, we review the current knowledge on nonepithelial tumor dissemination, drawing examples from the latest developments in melanoma, glioma, and sarcoma research. We underscore the importance of the reactivation of developmental processes during cancer progression and describe the nongenetic mechanisms driving nonepithelial tumor spread. We also outline therapeutic opportunities and ongoing clinical approaches to fight disseminating cancers. Finally, we discuss remaining challenges and emerging questions in the field. Defining the core principles underlying nonepithelial cancer dissemination may uncover actionable vulnerabilities of metastatic tumors and help improve the prognosis of patients with cancer.

Endoplasmic reticulum stress regulators exhibit different prognostic, therapeutic and immune landscapes in pancreatic adenocarcinoma

Endoplasmic reticulum stress (ERS) and unfolded protein response are the critical processes of tumour biology. However, the roles of ERS regulatory genes in pancreatic adenocarcinoma (PAAD) remain elusive. A novel ERS-related risk signature was constructed using the Lasso regression analysis. Its prognostic value, immune effect, metabolic influence, mutational feature and therapeutic correlation were comprehensively analysed through multiple bioinformatic approaches. The biofunctions of KDELR3 and YWHAZ in pancreatic cancer (PC) cells were also investigated through colony formation, Transwell assays, flow cytometric detection and a xenograft model. The upstream miRNA regulatory mechanism of KDELR3 was predicted and validated. ERS risk score was identified as an independent prognostic factor and could improve traditional prognostic model. Meanwhile, it was closely associated with metabolic reprogramming and tumour immune. High ERS risk enhanced glycolysis process and nucleotide metabolism, but was unfavourable for anti-tumour immune response. Moreover, ERS risk score could act as a potential biomarker for predicting the efficacy of ICBs. Overexpression of KDELR3 and YWHAZ stimulated the proliferation, migration and invasion of SW1990 and BxPC-3 cells. Silencing KDELR3 suppressed tumour growth in a xenograft model. miR-137 could weaken the malignant potentials of PC cells through inhibiting KDELR3 (5'-AGCAAUAA-3'). ERS risk score greatly contributed to PAAD clinical assessment. KDELR3 and YWHAZ possessed cancer-promoting capacities, showing promise as a novel treatment target.

Transcriptome-Based Treatment for Melanoma With Brain Metastasis: A Case Report

Malignant melanoma with brain metastasis has a high mortality rate. New approaches for diagnosis and treatment are urgently required to improve prognosis. Here we present a 60-year-old male with metastatic melanoma to the brain. Using a transcriptomics pipeline, we analyzed whole blood and resected tumor tissue, identifying enriched gene expression biomarkers and pathways - including seven upregulated ( BRAF, CDK4, EIF1AX, IK, NRAS, PIK3R2, and TP53) and 11 downregulated (CASP8, CDK10, CDKN2A, CTLA4, GNA11, HERC2, IRF4, MC1R, PLA2G6, RREB1, and TPCN2) genes in the blood (across 15 pathways), showing 14% enrichment, and 16 upregulated (CCND1, CDK4, CTLA4, EIF1AX, IK, IRF4, MITF, NRAS, PIK3CB, PIK3R2, PMEL, RREB1, SLC45A2, SOX10, TYR, and TYRP1) and three downregulated ( GNA11, KITLG, and PLA2G6) genes in tissue (across 17 pathways), showing 33% enrichment, with five shared markers and 12 shared pathways. The model connected CDK4 pathway overactivity observed in both samples to inhibitors like ribociclib, abemaciclib, and palbociclib as putative treatments. By enabling objective personalized therapy selection, this approach shows great promise for advancing patient outcomes.

LncRNAs in melanoma phenotypic plasticity: emerging targets for promising therapies

Long non-coding RNAs (lncRNAs) have received growing attention due to their diverse regulatory roles in cancer, including in melanoma, an aggressive type of skin cancer. The plasticity and phenotypic adaptability of melanoma cells are crucial factors contributing to therapeutic resistance. The identification of molecules playing key roles in melanoma cell plasticity could unravel novel and more effective therapeutic targets. This review presents current concepts of melanoma cell plasticity, illustrating its fluidity and dismissing the outdated notion of epithelial-mesenchymal-like transition as a simplistic binary process. Emphasis is placed on the pivotal role of lncRNAs in orchestrating cell plasticity, employing various mechanisms recently elucidated and unveiling their potential as promising targets for novel therapeutic strategies. Insights into the molecular mechanisms coordinated by lncRNAs in melanoma pave the way for the development of RNA-based therapies, holding great promise for enhancing treatment outcomes and offering a glimpse into a more effective approach to melanoma treatment.

The crucial roles of m6A RNA modifications in cutaneous cancers: Implications in pathogenesis, metastasis, drug resistance, and targeted therapies

N6-methyladenosine (m6A) is the most abundant internal modification on RNA. It is a dynamical and reversible process, which is regulated by m6A methyltransferase and m6A demethylase. The m6A modified RNA can be specifically recognized by the m6A reader, leading to RNA splicing, maturation, degradation or translation. The abnormality of m6A RNA modification is closely related to a variety of biological processes, especially the occurrence and development of tumors. Recent studies have shown that m6A RNA modification is involved in the pathogenesis of skin cancers. However, the precise molecular mechanisms of m6A-mediated cutaneous tumorigenesis have not been fully elucidated. Therefore, this review will summarize the biological characteristics of m6A modification, its regulatory role and mechanism in skin cancers, and the recent research progress of m6A-related molecular drugs, aiming to provide new ideas for clinical diagnosis and targeted therapy of cutaneous cancers.

A deep learning approach based on multi-omics data integration to construct a risk stratification prediction model for skin cutaneous melanoma

PURPOSE

Skin cutaneous melanoma (SKCM) is a highly aggressive melanocytic carcinoma whose high heterogeneity and complex etiology make its prognosis difficult to predict. This study aimed to construct a risk subtype typing model for SKCM.

METHODS

The study proposes a deep learning framework combining early fusion feature autoencoder (AE) and late fusion feature AE for risk subtype prediction of SKCM. The deep learning framework integrates mRNA, miRNA, and DNA methylation data of SKCM patients from The Cancer Genome Atlas (TCGA), and clusters the screened multi-omics features associated with survival prognosis to identify risk subtypes. Differential expression analysis and functional enrichment analysis were performed between risk subtypes, while SVM classifiers were constructed between differentially expressed genes (DEGs) obtained by Least Absolute Shrinkage and Selection Operator (LASSO) logistic regression screening and risk subtype labels inferred from multi-omics data, and the predictive robustness of risk subtypes inferred from the risk subtype classification prediction model was validated using two independent datasets.

RESULTS

The deep learning framework that combined early fusion feature AE with late fusion feature AE distinguished the two best risk subtypes compared to the multi-omics integration approach with single strategy AE or PCA. A promising C-index (C-index = 0.748) and a significant difference in survival (log-rank P value = 4.61 × 10-9) were found between the identified risk subtypes. The DEGs with the top significance values together with differentially expressed miRNAs provided the biological interpretation of risk subtypes on SKCM. Finally, the framework was applied to predict risk subtypes in two independent test datasets of SKCM patients, all of which showed good predictive power (C-index > 0.680) and significant survival differences (log-rank P value < 0.01).

CONCLUSION

The SKCM risk subtypes identified by integrating multi-omics data based on deep learning can not only improve the understanding of the molecular mechanisms of SKCM, but also provide clinicians with assistance in treatment decisions.

Activation of autophagy by in situ Zn2+ chelation reaction for enhanced tumor chemoimmunotherapy

Chemotherapy can induce a robust T cell antitumor immune response by triggering immunogenic cell death (ICD), a process in which tumor cells convert from nonimmunogenic to immunogenic forms. However, the antitumor immune response of ICD remains limited due to the low immunogenicity of tumor cells and the immunosuppressive tumor microenvironment. Although autophagy is involved in activating tumor immunity, the synergistic role of autophagy in ICD remains elusive and challenging. Herein, we report an autophagy amplification strategy using an ion-chelation reaction to augment chemoimmunotherapy in cancer treatments based on zinc ion (Zn²⁺)-doped, disulfiram (DSF)-loaded mesoporous silica nanoparticles (DSF@Zn-DMSNs). Upon pH-sensitive biodegradation of DSF@Zn-DMSNs, Zn²⁺ and DSF are coreleased in the mildly acidic tumor microenvironment, leading to the formation of toxic Zn²⁺ chelate through an in situ chelation reaction. Consequently, this chelate not only significantly stimulates cellular apoptosis and generates damage-associated molecular patterns (DAMPs) but also activates autophagy, which mediates the amplified release of DAMPs to enhance ICD. In vivo results demonstrated that DSF@Zn-DMSNs exhibit strong therapeutic efficacy via in situ ion chelation and possess the ability to activate autophagy, thus enhancing immunotherapy by promoting the infiltration of T cells. This study provides a smart in situ chelation strategy with tumor microenvironment-responsive autophagy amplification to achieve high tumor chemoimmunotherapy efficacy and biosafety.

Reactivation of embryonic genetic programs in tissue regeneration and disease

Embryonic genetic programs are reactivated in response to various types of tissue damage, providing cell plasticity for tissue regeneration or disease progression. In acute conditions, these programs remedy the damage and then halt to allow a return to homeostasis. In chronic situations, including inflammatory diseases, fibrosis and cancer, prolonged activation of embryonic programs leads to disease progression and tissue deterioration. Induction of progenitor identity and cell plasticity, for example, epithelial-mesenchymal plasticity, are critical outcomes of reactivated embryonic programs. In this Review, we describe molecular players governing reactivated embryonic genetic programs, their role during disease progression, their similarities and differences and lineage reversion in pathology and discuss associated therapeutics and drug-resistance mechanisms across many organs. We also discuss the diversity of reactivated programs in different disease contexts. A comprehensive overview of commonalities between development and disease will provide better understanding of the biology and therapeutic strategies.

Noncoding RNA circuitry in melanoma onset, plasticity, and therapeutic response

Melanoma, the cancer of the melanocyte, is the deadliest form of skin cancer with an aggressive nature, propensity to metastasize and tendency to resist therapeutic intervention. Studies have identified that the re-emergence of developmental pathways in melanoma contributes to melanoma onset, plasticity, and therapeutic response. Notably, it is well known that noncoding RNAs play a critical role in the development and stress response of tissues. In this review, we focus on the noncoding RNAs, including microRNAs, long non-coding RNAs, circular RNAs, and other small RNAs, for their functions in developmental mechanisms and plasticity, which drive onset, progression, therapeutic response and resistance in melanoma. Going forward, elucidation of noncoding RNA-mediated mechanisms may provide insights that accelerate development of novel melanoma therapies.

The Genomic Landscape of Melanoma and Its Therapeutic Implications

Melanoma is one of the most aggressive malignancies of the skin. The genetic composition of melanoma is complex and varies among different subtypes. With the aid of recent technologies such as next generation sequencing and single-cell sequencing, our understanding of the genomic landscape of melanoma and its tumor microenvironment has become increasingly clear. These advances may provide explanation to the heterogenic treatment outcomes of melanoma patients under current therapeutic guidelines and provide further insights to the development of potential new therapeutic targets. Here, we provide a comprehensive review on the genetics related to melanoma tumorigenesis, metastasis, and prognosis. We also review the genetics affecting the melanoma tumor microenvironment and its relation to tumor progression and treatment.

The journey from melanocytes to melanoma

Over the past decade, melanoma has led the field in new cancer treatments, with impressive gains in on-treatment survival but more modest improvements in overall survival. Melanoma presents heterogeneity and transcriptional plasticity that recapitulates distinct melanocyte developmental states and phenotypes, allowing it to adapt to and eventually escape even the most advanced treatments. Despite remarkable advances in our understanding of melanoma biology and genetics, the melanoma cell of origin is still fiercely debated because both melanocyte stem cells and mature melanocytes can be transformed. Animal models and high-throughput single-cell sequencing approaches have opened new opportunities to address this question. Here, we discuss the melanocytic journey from the neural crest, where they emerge as melanoblasts, to the fully mature pigmented melanocytes resident in several tissues. We describe a new understanding of melanocyte biology and the different melanocyte subpopulations and microenvironments they inhabit, and how this provides unique insights into melanoma initiation and progression. We highlight recent findings on melanoma heterogeneity and transcriptional plasticity and their implications for exciting new research areas and treatment opportunities. The lessons from melanocyte biology reveal how cells that are present to protect us from the damaging effects of ultraviolet radiation reach back to their origins to become a potentially deadly cancer.

An epigenetic switch controls an alternative NR2F2 isoform that unleashes a metastatic program in melanoma

Metastatic melanoma develops once transformed melanocytic cells begin to de-differentiate into migratory and invasive melanoma cells with neural crest cell (NCC)-like and epithelial-to-mesenchymal transition (EMT)-like features. However, it is still unclear how transformed melanocytes assume a metastatic melanoma cell state. Here, we define DNA methylation changes that accompany metastatic progression in melanoma patients and discover Nuclear Receptor Subfamily 2 Group F, Member 2 - isoform 2 (NR2F2-Iso2) as an epigenetically regulated metastasis driver. NR2F2-Iso2 is transcribed from an alternative transcriptional start site (TSS) and it is truncated at the N-terminal end which encodes the NR2F2 DNA-binding domain. We find that NR2F2-Iso2 expression is turned off by DNA methylation when NCCs differentiate into melanocytes. Conversely, this process is reversed during metastatic melanoma progression, when NR2F2-Iso2 becomes increasingly hypomethylated and re-expressed. Our functional and molecular studies suggest that NR2F2-Iso2 drives metastatic melanoma progression by modulating the activity of full-length NR2F2 (Isoform 1) over EMT- and NCC-associated target genes. Our findings indicate that DNA methylation changes play a crucial role during metastatic melanoma progression, and their control of NR2F2 activity allows transformed melanocytes to acquire NCC-like and EMT-like features. This epigenetically regulated transcriptional plasticity facilitates cell state transitions and metastatic spread.

Broad misappropriation of developmental splicing profile by cancer in multiple organs

Oncogenesis mimics key aspects of embryonic development. However, the underlying mechanisms are incompletely understood. Here, we demonstrate that the splicing events specifically active during human organogenesis, are broadly reactivated in the organ-specific tumor. Such events are associated with key oncogenic processes and predict proliferation rates in cancer cell lines as well as patient survival. Such events preferentially target nitrosylation and transmembrane-region domains, whose coordinated splicing in multiple genes respectively affect intracellular transport and N-linked glycosylation. We infer critical splicing factors potentially regulating embryonic splicing events and show that such factors are potential oncogenic drivers and are upregulated specifically in malignant cells. Multiple complementary analyses point to MYC and FOXM1 as potential transcriptional regulators of critical splicing factors in brain and liver. Our study provides a comprehensive demonstration of a splicing-mediated link between development and cancer, and suggest anti-cancer targets including splicing events, and their upstream splicing and transcriptional regulators.

Fingerprinting of skin cells by live cell Raman spectroscopy reveals melanoma cell heterogeneity and cell-type-specific responses to UVR

Raman spectroscopy is an emerging dermatological technique with the potential to discriminate biochemically between cell types in a label-free and non-invasive manner. Here, we use live single-cell Raman spectroscopy and principal component analysis (PCA) to fingerprint mouse melanoblasts, melanocytes, keratinocytes and melanoma cells. We show the differences in their spectra are attributable to biomarkers in the melanin biosynthesis pathway and that melanoma cells are a heterogeneous population that sit on a trajectory between undifferentiated melanoblasts and differentiated melanocytes. We demonstrate the utility of Raman spectroscopy as a highly sensitive tool to probe the melanin biosynthesis pathway and its immediate response to ultraviolet (UV) irradiation revealing previously undescribed opposing responses to UVA and UVB irradiation in melanocytes. Finally, we identify melanocyte-specific accumulation of β-carotene correlated with a stabilisation of the UVR response in lipids and proteins consistent with a β-carotene-mediated photoprotective mechanism. In summary, our data show that Raman spectroscopy can be used to determine the differentiation status of cells of the melanocyte lineage and describe the immediate and temporal biochemical changes associated with UV exposure which differ depending on cell type, differentiation status and competence to synthesise melanin. Our work uniquely applies Raman spectroscopy to discriminate between cell types by biological function and differentiation status while they are growing in culture. In doing so, we demonstrate for the first time its utility as a tool with which to probe the melanin biosynthesis pathway.

KDELR3 Is a Prognostic Biomarker Related to the Immune Infiltration and Chemoresistance of Anticancer Drugs in Uveal Melanoma

Uveal melanoma (UM) is an intraocular malignancy in adults in which approximately 50% of patients develop metastatic diseases and have a poor clinical outcome. Immunotherapies are quickly becoming a need, and recent research has produced some amazing achievements in this area. In the current investigation, an attempt was made to evaluate the prognostic usefulness of KDELR3 in UM, particularly its connection with tumor-infiltrating lymphocytes (TILs). The expression patterns of mRNAs and related clinical data of 80 UM patients were obtained from The Cancer Genome Atlas (TCGA). By using RT-PCR, we were able to investigate whether or not UM cells and D78 cells expressed KDELR3. The Kaplan-Meier approach, as well as univariate and multivariate tests, was utilized in order to investigate the potential predictive significance of KDELR3 expression. The associations between KDELR3 and TILs and immunological checkpoints were analyzed in order to evaluate the effect that KDELR3 may have on UM immunotherapy. On the basis of the differential expression of KDELR3, a distribution of the half-maximal inhibitory concentration (IC50) of various targeted medicines was observed. In this study, we found that the expression of KDELR3 was distinctly increased in most types of tumors. In addition, KDELR3 was highly expressed in UM cells. Moreover, patients with high KDELR3 expression exhibited a shorter overall survival and disease-free survival than those with low KDELR3 expression. Multivariate analyses confirmed that KDELR3 expression was an independent prognostic factor for overall survival and disease-free survival in patients with UM. Furthermore, KDELR3 expression was demonstrated to be positively correlated with macrophage M1, T cell CD8, T cell follicular helper, dendritic cell resting, and T cell CD4 memory activated. Meanwhile, the expression of KDELR3 was related to several immune checkpoints. The IC50 of AP-24534, BHG712, bleomycin, camptothecin, cisplatin, cytarabine, GSK1070916, and tipifarnib was higher in the KDELR3 high-expression group. In conclusion, KDELR3 may be applied as a potential diagnostic and prognostic biomarker for UM patients.

TRPM4 and TRPV2 are two novel prognostic biomarkers and promising targeted therapy in UVM

Uveal melanoma (UVM) is the most common primary intraocular malignancy tumor in adults. Almost 50% of UVM patients develop metastatic disease, and is usually fatal within 1 year. However, the mechanism of etiology remains unclear. The lack of prognostic, diagnostic and therapeutic biomarkers is a main limitation for clinical diagnosis and treatment. The transient receptor potential (TRP) channels play important roles in the occurrence and development of tumors, which may have the potential as a therapeutic target for UVM. This current study aimed to identify the potential effect and function of the TRPs that could provide survival prediction and new insight into therapy for UVM. Based on the transcriptome data and potential key genes of UVM were screened using the Cancer Genome Atlas (TCGA) databases, Gene expression analysis showed the expression of TRPM4, TRPV2 and other TRPs was high levels in UVM. Using survival analysis, we screened out that the high expression of TRPM4 and TRPV2 was negatively correlated with the prognosis of UVM patients. Cox regression analysis and functional enrichment analysis further indicated that TRPM4 and TRPV2 were the most convincing therapeutic targets of UVM, and the majority of genes involved in ferroptosis pathways in UVM showed positively correlated with the expression levels of TRPM4 and TRPV2. In conclusion, TRPM4 and TRPV2 were considered as two novel prognostic biomarkers and a promising targeted therapy in UVM.

Synergistic effects of nab-PTX and anti-PD-1 antibody combination against lung cancer by regulating the Pi3K/AKT pathway through the Serpinc1 gene

Lung cancer is a type of cancer with higher morbidity and mortality. In spite of the impressive response rates of nab-paclitaxel (nab-PTX) or programmed cell death-1 (PD-1) and its ligand inhibitors, the effective treatment remains limited. Currently, alternative strategies aim at drug combination of nab-PTX and PD-1/PD-L1 inhibitors. Even as the clinical impact of the combined agents continues to increase, basic research studies are still limited and the mechanisms underlying this synergy are not well studied. In this study, we evaluated the antitumor efficacy and the molecular mechanisms of action of nab-PTX in combination with anti-PD-1 antibody, using Lewis lung carcinoma (LLC) cell and subcutaneously transplanted tumor models. The combination of nab-PTX and anti-PD-1 antibody displayed stronger antitumor effects, manifested at tumor volume, proliferation and apoptosis through Ki67 and TUNEL staining. In-vivo experiments showed significant increases in CD4+ T cells, CD8+ T cells, IFN-γ, TNF-α, IL-2, PF, and Gzms-B, exerting antitumor effects with reductions in MDSCs and IL-10 after the treatments. Furthermore, transcriptomic analysis indicated 20 overlapped differentially expressed genes, and Serpin peptidase inhibitor clade C Member 1 (Serpinc1) was downregulated during treatment in vivo, whose expression level was markedly related to metastasis and overall survival of lung cancer patients. Functional enrichment analysis of the target gene revealed primary GO terms related to tumor, which warrants further investigation. We also found that Serpinc1 overexpression promoted cell proliferation, migration, and invasion and inhibited cell apoptosis of LLC cells in vitro, possibly regulating the associated factors via the Pi3K/AKT pathway. In summary, our results reveal the synergistic antitumor responses of nab-PTX combined with anti-PD-1 antibody, in which Serpinc1 may play an important role, providing a target gene for combination treatment strategy.

RIG-I acts as a tumor suppressor in melanoma via regulating the activation of the MKK/p38MAPK signaling pathway

Studies have indicated that RIG-I may act as a tumor suppressor and participate in the tumorigenesis of some malignant diseases. However, RIG-I induces distinct cellular responses via different downstream signaling pathways depending on the cell type. To investigate the biological function and underlying molecular mechanism of RIG-I in the tumorigenesis of melanoma, we constructed RIG-I knockout, RIG-I-overexpressing B16-F10 and RIG-I knockdown A375 melanoma cell lines, and analyzed the RIG-I-mediated change in the biological behavior of tumor cells in spontaneous and poly (I:C)-induced RIG-I activation. Cell proliferation, cell cycling, apoptosis and migration were detected by CCK-8 assay, BrdU incorporation assay, Annexin V-PI staining assay and Transwell assay, respectively. In vivo tumorigenicity was evaluated by tumor xenograft growth in nude mice and subsequently by Ki67 staining and TUNEL assays. Furthermore, Western blotting was utilized to explore the underlying mechanism of RIG-I in melanoma cells. Our data showed that RIG-I promotes apoptosis and inhibits proliferation by G1 phase cell cycle arrest in the melanoma cell lines. Mechanistically, RIG-I induced the phosphorylation of p38 MAPK and MAPK kinases MKK3 and MKK4. In conclusion, the current study demonstrated that RIG-I suppressed the development of melanoma by regulating the activity of the MKK/p38 MAPK signaling pathway, which is relevant to research on novel therapeutic targets for this malignant disease.

KDELR1 Is an Independent Prognostic Predictor and Correlates With Immunity in Glioma

Background

Gliomas are the most malignant central nervous system tumors. With the development of sequencing technology, more potential biomarkers related to the treatment, prognosis, and molecular classification of glioma have been identified. Here, we intend to investigate the potential biological function and clinical value of a new biomarker in glioma.

Methods

KDELR1 expression data and the corresponding clinical information were downloaded from public databases and then preprocessed using R language. Correlation, Kaplan–Meier survival, and Cox regression analyses were performed to explore the clinical significance of KDELR1 in glioma patients. Furthermore, the immune infiltration and microenvironment parameters were evaluated via TIMER and CIBERSORT. Immunohistochemistry was conducted to confirm the KDELR1 expression and its correlation with immunity infiltration and prognosis.

Results

KDELR1 was upregulated in glioma samples compared with normal brain tissues, and its expression was significantly correlated with age, the World Health Organization (WHO) grade, recurrence, necrosis, microvascular proliferation, molecular classification, isocitrate dehydrogenase (IDH) mutation, and 1p/19q codeletion status. In addition, survival analysis showed that glioma patients with KDELR1 overexpression had shorter overall survival (OS) and disease-free survival times, and Cox regression analysis revealed that KDELR1 acted as an independent prognostic factor of OS in glioma patients. Gene set enrichment analysis indicated a significant enrichment of metabolism-associated pathways. KDELR1 expression was positively associated with immune infiltration (including infiltration by CD8+ T cells, CD4+ T cells, macrophages, and so on) and microenvironment parameters (including stromal, immune, and ESTIMATE scores) in gliomas. The expression of KDELR1 and its correlation with the tumor grade and prognosis were confirmed by immunohistochemistry in clinical samples (n = 119, P < 0.05).

Conclusions

Taken together, these findings suggest that KDELR1 is correlated with the tumor grade, molecular classifications, and immune infiltration; highlighting that KDELR1 is a novel and promising biomarker for molecular classification, treatment, and prognostic assessment may further indicate the treating effect of immune therapy.

KDEL Receptors: Pathophysiological Functions, Therapeutic Options, and Biotechnological Opportunities

KDEL receptors (KDELRs) are ubiquitous seven-transmembrane domain proteins encoded by three mammalian genes. They bind to and retro-transport endoplasmic reticulum (ER)-resident proteins with a C-terminal Lys-Asp-Glu-Leu (KDEL) sequence or variants thereof. In doing this, KDELR participates in the ER quality control of newly synthesized proteins and the unfolded protein response. The binding of KDEL proteins to KDELR initiates signaling cascades involving three alpha subunits of heterotrimeric G proteins, Src family kinases, protein kinases A (PKAs), and mitogen-activated protein kinases (MAPKs). These signaling pathways coordinate membrane trafficking flows between secretory compartments and control the degradation of the extracellular matrix (ECM), an important step in cancer progression. Considering the basic cellular functions performed by KDELRs, their association with various diseases is not surprising. KDELR mutants unable to bind the collagen-specific chaperon heat-shock protein 47 (HSP47) cause the osteogenesis imperfecta. Moreover, the overexpression of KDELRs appears to be linked to neurodegenerative diseases that share pathological ER-stress and activation of the unfolded protein response (UPR). Even immune function requires a functional KDELR1, as its mutants reduce the number of T lymphocytes and impair antiviral immunity. Several studies have also brought to light the exploitation of the shuttle activity of KDELR during the intoxication and maturation/exit of viral particles. Based on the above, KDELRs can be considered potential targets for the development of novel therapeutic strategies for a variety of diseases involving proteostasis disruption, cancer progression, and infectious disease. However, no drugs targeting KDELR functions are available to date; rather, KDELR has been leveraged to deliver drugs efficiently into cells or improve antigen presentation.

Loss of 'Epidermal Melanin Unit' Integrity in Human Skin During Melanoma-Genesis

Cutaneous melanoma can be a most challenging neoplasm of high lethality, in part due to its extreme heterogeneity and characteristic aggressive and invasive nature. Indeed, its moniker 'the great masquerader' reflects that not all melanomas are created equal in terms of their originating cellular contexts, but also that melanoma cells in the malignant tumor can adopt a wide range of different cell states and variable organotropism. In this review, we focus on the early phases of melanomagenesis by discussing how the originating pigment cell of the melanocyte lineage can be influenced to embark on a wide range of tumor fates with distinctive microanatomical pathways. In particular, we assess how cells of the melanocyte lineage can differ by maturation status (stem cell; melanoblast; transiently amplifying cell; differentiated; post-mitotic; terminally-differentiated) as well as by micro-environmental niche (in the stratum basale of the epidermis; within skin appendages like hair follicle, eccrine gland, etc). We discuss how the above variable contexts may influence the susceptibility of the epidermal-melanin unit (EMU) to become unstable, which may presage cutaneous melanoma development. We also assess how unique features of follicular-melanin unit(s) (FMUs) can, by contrast, protect melanocytes from melanomagenesis. Lastly, we postulate how variable melanocyte fates in vitiligo, albinism, psoriasis, and alopecia areata may provide new insights into immune-/non immune-mediated outcomes for melanocytes in cutaneous melanin units.

Aldh2 is a lineage-specific metabolic gatekeeper in melanocyte stem cells

Melanocyte stem cells (McSCs) in zebrafish serve as an on-demand source of melanocytes during growth and regeneration, but metabolic programs associated with their activation and regenerative processes are not well known. Here, using live imaging coupled with scRNA-sequencing, we discovered that, during regeneration, quiescent McSCs activate a dormant embryonic neural crest transcriptional program followed by an aldehyde dehydrogenase (Aldh) 2 metabolic switch to generate progeny. Unexpectedly, although ALDH2 is well known for its aldehyde-clearing mechanisms, we find that, in regenerating McSCs, Aldh2 activity is required to generate formate – the one-carbon (1C) building block for nucleotide biosynthesis – through formaldehyde metabolism. Consequently, we find that disrupting the 1C cycle with low doses of methotrexate causes melanocyte regeneration defects. In the absence of Aldh2, we find that purines are the metabolic end product sufficient for activated McSCs to generate progeny. Together, our work reveals McSCs undergo a two-step cell state transition during regeneration, and that the reaction products of Aldh2 enzymes have tissue-specific stem cell functions that meet metabolic demands in regeneration.

Summary: In zebrafish melanocyte regeneration, quiescent McSCs respond by re-expressing a neural crest identity, followed by an Aldh2-dependent metabolic switch to generate progeny.

Endogenous and Exogenous Regulatory Signaling in the Secretory Pathway: Role of Golgi Signaling Molecules in Cancer

The biosynthetic transport route that constitutes the secretory pathway plays a fundamental role in the cell, providing to the synthesis and transport of around one third of human proteins and most lipids. Signaling molecules within autoregulatory circuits on the intracellular membranes of the secretory pathway regulate these processes, especially at the level of the Golgi complex. Indeed, cancer cells can hijack several of these signaling molecules, and therefore also the underlying regulated processes, to bolster their growth or gain more aggressive phenotypes. Here, we review the most important autoregulatory circuits acting on the Golgi, emphasizing the role of specific signaling molecules in cancer. In fact, we propose to draw awareness to highlight the Golgi-localized regulatory systems as potential targets in cancer therapy.

Revisiting miRNA Association with Melanoma Recurrence and Metastasis from a Machine Learning Point of View

The diagnostic and prognostic value of miRNAs in cutaneous melanoma (CM) has been broadly studied and supported by advanced bioinformatics tools. From early studies using miRNA arrays with several limitations, to the recent NGS-derived miRNA expression profiles, an accurate diagnostic panel of a comprehensive pre-specified set of miRNAs that could aid timely identification of specific cancer stages is still elusive, mainly because of the heterogeneity of the approaches and the samples. Herein, we summarize the existing studies that report several miRNAs as important diagnostic and prognostic biomarkers in CM. Using publicly available NGS data, we analyzed the correlation of specific miRNA expression profiles with the expression signatures of known gene targets. Combining network analytics with machine learning, we developed specific non-linear classification models that could successfully predict CM recurrence and metastasis, based on two newly identified miRNA signatures. Subsequent unbiased analyses and independent test sets (i.e., a dataset not used for training, as a validation cohort) using our prediction models resulted in 73.85% and 82.09% accuracy in predicting CM recurrence and metastasis, respectively. Overall, our approach combines detailed analysis of miRNA profiles with heuristic optimization and machine learning, which facilitates dimensionality reduction and optimization of the prediction models. Our approach provides an improved prediction strategy that could serve as an auxiliary tool towards precision treatment.

Tfap2b specifies an embryonic melanocyte stem cell that retains adult multifate potential

Melanocytes, the pigment-producing cells, are replenished from multiple stem cell niches in adult tissue. Although pigmentation traits are known risk factors for melanoma, we know little about melanocyte stem cell (McSC) populations other than hair follicle McSCs and lack key lineage markers with which to identify McSCs and study their function. Here we find that Tfap2b and a select set of target genes specify an McSC population at the dorsal root ganglia in zebrafish. Functionally, Tfap2b is required for only a few late-stage embryonic melanocytes, and is essential for McSC-dependent melanocyte regeneration. Fate mapping data reveal that tfap2b⁺ McSCs have multifate potential, and are the cells of origin for large patches of adult melanocytes, two other pigment cell types (iridophores and xanthophores), and nerve-associated cells. Hence, Tfap2b confers McSC identity in early development, distinguishing McSCs from other neural crest and pigment cell lineages, and retains multifate potential in the adult zebrafish.

Transcriptional profile and chromatin accessibility in zebrafish melanocytes and melanoma tumors

Transcriptional and epigenetic characterization of melanocytes and melanoma cells isolated from their in vivo context promises to unveil key differences between these developmentally related normal and cancer cell populations. We therefore engineered an enhanced Danio rerio (zebrafish) melanoma model with fluorescently labeled melanocytes to allow for isolation of normal (wild type) and premalignant (BRAFV600E-mutant) populations for comparison to fully transformed BRAFV600E-mutant, p53 loss-of-function melanoma cells. Using fluorescence-activated cell sorting to isolate these populations, we performed high-quality RNA- and ATAC-seq on sorted zebrafish melanocytes vs. melanoma cells, which we provide as a resource here. Melanocytes had consistent transcriptional and accessibility profiles, as did melanoma cells. Comparing melanocytes and melanoma, we note 4128 differentially expressed genes and 56,936 differentially accessible regions with overall gene expression profiles analogous to human melanocytes and the pigmentation melanoma subtype. Combining the RNA- and ATAC-seq data surprisingly revealed that increased chromatin accessibility did not always correspond with increased gene expression, suggesting that though there is widespread dysregulation in chromatin accessibility in melanoma, there is a potentially more refined gene expression program driving cancerous melanoma. These data serve as a resource to identify candidate regulators of the normal vs. diseased states in a genetically controlled in vivo context.

Genomic and Transcriptomic Underpinnings of Melanoma Genesis, Progression, and Metastasis

Melanoma is a deadly skin cancer with rapidly increasing incidence worldwide. The discovery of the genetic drivers of melanomagenesis in the last decade has led the World Health Organization to reclassify melanoma subtypes by their molecular pathways rather than traditional clinical and histopathologic features. Despite this significant advance, the genomic and transcriptomic drivers of metastatic progression are less well characterized. This review describes the known molecular pathways of cutaneous and uveal melanoma progression, highlights recently identified pathways and mediators of metastasis, and touches on the influence of the tumor microenvironment on metastatic progression and treatment resistance. While targeted therapies and immune checkpoint blockade have significantly aided in the treatment of advanced disease, acquired drug resistance remains an unfortunately common problem, and there is still a great need to identify potential prognostic markers and novel therapeutic targets to aid in such cases.

Identification of the Key Genes Associated with the Yak Hair Follicle Cycle

The development of hair follicles in yak shows significant seasonal cycles. In our previous research, transcriptome data including mRNAs and lncRNAs in five stages during the yak hair follicles (HFs) cycle were detected, but their regulation network and the hub genes in different periods are yet to be explored. This study aimed to screen and identify the hub genes during yak HFs cycle by constructing a mRNA-lncRNA co-expression network. A total of 5000 differently expressed mRNA (DEMs) and 729 differently expressed long noncoding RNA (DELs) were used to construct the co-expression network, based on weighted genes co-expression network analysis (WGCNA). Four temporally specific modules were considered to be significantly associated with the HFs cycle of yak. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that the modules are enriched into Wnt, EMC-receptor interaction, PI3K-Akt, focal adhesion pathways, and so on. The hub genes, such as FER, ELMO1, PCOLCE, and HOXC13, were screened in different modules. Five hub genes (WNT5A, HOXC13, DLX3, FOXN1, and OVOL1) and part of key lncRNAs were identified for specific expression in skin tissue. Furthermore, immunofluorescence staining and Western blotting results showed that the expression location and abundance of DLX3 and OVOL1 are changed following the process of the HFs cycle, which further demonstrated that these two hub genes may play important roles in HFs development.

Identification of a novel signature based on unfolded protein response-related gene for predicting prognosis in bladder cancer

BACKGROUND

The unfolded protein response (UPR) served as a vital role in the progression of tumors, but the molecule mechanisms of UPR in bladder cancer (BLCA) have been not fully investigated.

METHODS

We identified differentially expressed unfolded protein response-related genes (UPRRGs) between BLCA samples and normal bladder samples in the Cancer Genome Atlas (TCGA) database. Univariate Cox analysis and the least absolute shrinkage and selection operator penalized Cox regression analysis were used to construct a prognostic signature in the TCGA set. We implemented the validation of the prognostic signature in GSE13507 from the Gene Expression Omnibus database. The ESTIMATE, CIBERSORT, and ssGSEA algorithms were used to explore the correlation between the prognostic signature and immune cells infiltration as well as key immune checkpoints (PD-1, PD-L1, CTLA-4, and HAVCR2). GDSC database analyses were conducted to investigate the chemotherapy sensitivity among different groups. GSEA analysis was used to explore the potential mechanisms of UPR-based signature.

RESULTS

A prognostic signature comprising of seven genes (CALR, CRYAB, DNAJB4, KDELR3, CREB3L3, HSPB6, and FBXO6) was constructed to predict the outcome of BLCA. Based on the UPRRGs signature, the patients with BLCA could be classified into low-risk groups and high-risk groups. Patients with BLCA in the low-risk groups showed the more favorable outcomes than those in the high-risk groups, which was verified in GSE13507 set. This signature could serve as an autocephalous prognostic factor in BLCA. A nomogram based on risk score and clinical characteristics was established to predict the over survival of BLCA patients. Furthermore, the signature was closely related to immune checkpoints (PD-L1, CTLA-4, and HAVCR2) and immune cells infiltration including CD8⁺ T cells, follicular helper T cells, activated dendritic cells, and M2 macrophages. GSEA analysis indicated that immune and carcinogenic pathways were enriched in high-risk group.

CONCLUSIONS

We identified a novel unfolded protein response-related gene signature which could predict the over survival, immune microenvironment, and chemotherapy response of patients with bladder cancer.

Transcriptomic Profiling Revealed Plexin A2 Downregulation With Migration and Invasion Alteration in Dacarbazine-Treated Primary Melanoma Cells

Melanoma is highly heterogeneous type of malignant neoplasm that is responsible for the majority of deaths among other types of skin cancer. In the present study, we screened a list of differentially expressed genes in two primary, drug-naïve melanoma cell lines derived from patients with melanoma following treatment of the cells with the chemotherapeutic agent dacarbazine. The aim was to determine the transcriptomic profiles and associated alterations in the cell phenotype. We found the vascular endothelial growth factor A/vascular endothelial growth factor receptor 2, phosphoinositide 3-kinase/protein kinase B and focal adhesion signaling pathways to be top altered after dacarbazine treatment. In addition, we observed the expression levels of genes associated with tumor dissemination, integrin β8 and matrix metalloproteinase-1, to be diminished in both cell lines studied, the results of which were confirmed by reverse transcription-quantitative polymerase chain reaction. By contrast, plexin A2 expression was found to be upregulated in K2303 cells, where reduced migration and invasion were also observed, following dacarbazine treatment. Plexin A2 downregulation was associated with the promotion of migrative and invasive capacities in B0404 melanoma cells. Since plexin A2 is semaphorin co-receptor that is involved in focal adhesion and cell migration regulation, the present study suggested that plexin A2 may be implicated in the dacarbazine-mediated phenotypic shift of melanoma cells. We propose that the signature of cancer cell invasiveness can be revealed by using a combination of transcriptomic and functional approaches, which should be applied in the development of personalized therapeutic strategies for each patient with melanoma.

Neural Crest-Like Stem Cell Transcriptome Analysis Identifies LPAR1 in Melanoma Progression and Therapy Resistance

Metastatic melanoma is challenging to clinically address. Although standard-of-care targeted therapy has high response rates in patients with BRAF-mutant melanoma, therapy relapse occurs in most cases. Intrinsically resistant melanoma cells drive therapy resistance and display molecular and biologic properties akin to neural crest-like stem cells (NCLSC) including high invasiveness, plasticity, and self-renewal capacity. The shared transcriptional programs and vulnerabilities between NCLSCs and cancer cells remains poorly understood. Here, we identify a developmental LPAR1-axis critical for NCLSC viability and melanoma cell survival. LPAR1 activity increased during progression and following acquisition of therapeutic resistance. Notably, genetic inhibition of LPAR1 potentiated BRAFi ± MEKi efficacy and ablated melanoma migration and invasion. Our data define LPAR1 as a new therapeutic target in melanoma and highlights the promise of dissecting stem cell-like pathways hijacked by tumor cells. SIGNIFICANCE: This study identifies an LPAR1-axis critical for melanoma invasion and intrinsic/acquired therapy resistance.

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.

Human melanocyte development and melanoma dedifferentiation at single-cell resolution

In humans, epidermal melanocytes are responsible for skin pigmentation, defence against ultraviolet radiation and the deadliest common skin cancer, melanoma. Although there is substantial overlap in melanocyte development pathways between different model organisms, species-dependent differences are frequent and the conservation of these processes in human skin remains unresolved. Here, we used a single-cell enrichment and RNA-sequencing pipeline to study human epidermal melanocytes directly from the skin, capturing transcriptomes across different anatomical sites, developmental age, sexes and multiple skin tones. We uncovered subpopulations of melanocytes that exhibit anatomical site-specific enrichment that occurs during gestation and persists through adulthood. The transcriptional signature of the volar-enriched subpopulation is retained in acral melanomas. Furthermore, we identified human melanocyte differentiation transcriptional programs that are distinct from gene signatures generated from model systems. Finally, we used these programs to define patterns of dedifferentiation that are predictive of melanoma prognosis and response to immune checkpoint inhibitor therapy.

Dosage-Dependent Gynoecium Development and Gene Expression in Brassica napus-Orychophragmus violaceus Addition Lines

Distant hybridization usually leads to female sterility of the hybrid but the mechanism behind this is poorly understood. Complete pistil abortion but normal male fertility was shown by one Brassica napus-Orychophragmus violaceus monosomic alien addition line (MA, AACC + 1 IO, 2n = 39) produced previously. To study the effect of a single O. violaceus chromosome addition on pistil development in different genetic backgrounds, hybrids between the MA and B. carinata (BBCC), B. juncea (AABB), and two synthetic hexaploids (AABBCC) were firstly produced in this study which show complete female sterility. A microspore culture was further performed to produce the haploid monosomic alien addition line (HMA, AC + 1 IO, 2n = 20) and disomic addition line (DA, AACC + 2 IO, 2n = 40) together with haploid (H, AC, 2n = 19) and double haploid (DH, AACC, 2n = 38) plants of B. napus from MA to investigate the dosage effect of the alien O. violaceus chromosome on pistil development and gene expression. Compared to MA, the development of the pistils of DA and HMA was completely or partially recovered, in which the pistils could swell and elongate to a normal shape after open pollination, although no seeds were produced. Comparative RNA-seq analyses revealed that the numbers of the differentially expressed genes (DEGs) were significantly different, dosage-dependent, and consistent with the phenotypic difference in pairwise comparisons of HMA vs. H, DA vs. DH, MA vs. DH, MA vs. DA, and MA vs. HMA. The gene ontology (GO) enrichment analysis of DEGs showed that a number of genes involved in the development of the gynoecium, embryo sac, ovule, and integuments. Particularly, several common DEGs for pistil development shared in HMA vs. H and DA vs. DH showed functions in genotoxic stress response, auxin transport, and signaling and adaxial/abaxial axis specification. The results provided updated information for the molecular mechanisms behind the gynoecium development of B. napus responding to the dosage of alien O. violaceus chromosomes.

Development and validation of an intra-tumor heterogeneity-related signature to predict prognosis of bladder cancer: a study based on single-cell RNA-seq

Intra-tumor heterogeneity (ITH) was a potential mechanism of progression and drug resistance in bladder cancer (BCa). However, the understanding of ITH in BCa remains insufficient. Single-cell RNA sequencing (scRNA-seq) profiles of 2075 cells were analyzed, and 2940 cell markers were screened. The ITH of 396 cases was evaluated, and 96 ITH-related genes were identified. Based on the gene-pair strategy, 96 genes were cyclically paired, and an 8-gene-pair model was successfully established to evaluate the overall survival of BCa through Lasso and multivariate Cox regressions. The risk model showed high predictive value in the training dataset (n = 396, p = 0) and external validation datasets (n = 165, p = 2.497e-02; n = 224, p = 3.423e-02). The model was also valuable for the prediction of clinical treatment outcomes. Totally, a prognostic model based on scRNA-seq and ITH was successfully constructed and validated in large cohorts, providing novel clues for ITH study of BCa.

Targeting the SHP2 phosphatase promotes vascular damage and inhibition of tumor growth

The tyrosine phosphatase SHP2 is oncogenic in cancers driven by receptor-tyrosine-kinases, and SHP2 inhibition reduces tumor growth. Here, we report that SHP2 is an essential promoter of endothelial cell survival and growth in the remodeling tumor vasculature. Using genetic and chemical approaches to inhibit SHP2 activity in endothelial cells, we show that SHP2 inhibits pro-apoptotic STAT3 and stimulates proliferative ERK1/2 signaling. Systemic SHP2 inhibition in mice bearing tumor types selected for SHP2-independent tumor cell growth promotes degeneration of the tumor vasculature and blood extravasation; reduces tumor vascularity and blood perfusion; and increases tumor necrosis. Reduction of tumor growth ensues, independent of SHP2 targeting in the tumor cells, blocking immune checkpoints, or recruiting macrophages. We also show that inhibiting the Angiopoietin/TIE2/AKT cascade magnifies the vascular and anti-tumor effects of SHP2 inhibition by blocking tumor endothelial AKT signaling, not a target of SHP2. Since the SHP2 and Ang2/TIE2 pathways are active in vascular endothelial cells of human melanoma and colon carcinoma, SHP2 inhibitors alone or with Ang2/TIE2 inhibitors hold promise to effectively target the tumor endothelium.

Prioritizing and characterizing functionally relevant genes across human tissues

Knowledge of genes that are critical to a tissue's function remains difficult to ascertain and presents a major bottleneck toward a mechanistic understanding of genotype-phenotype links. Here, we present the first machine learning model-FUGUE-combining transcriptional and network features, to predict tissue-relevant genes across 30 human tissues. FUGUE achieves an average cross-validation auROC of 0.86 and auPRC of 0.50 (expected 0.09). In independent datasets, FUGUE accurately distinguishes tissue or cell type-specific genes, significantly outperforming the conventional metric based on tissue-specific expression alone. Comparison of tissue-relevant transcription factors across tissue recapitulate their developmental relationships. Interestingly, the tissue-relevant genes cluster on the genome within topologically associated domains and furthermore, are highly enriched for differentially expressed genes in the corresponding cancer type. We provide the prioritized gene lists in 30 human tissues and an open-source software to prioritize genes in a novel context given multi-sample transcriptomic data.

KDELR2 knockdown synergizes with temozolomide to induce glioma cell apoptosis through the CHOP and JNK/p38 pathways

BACKGROUND

The C-terminal tetrapeptide Lys-Asp-Glu-Leu receptors (KDELRs) are transmembrane proteins that regulate ER stress (ERS) response, growth, differentiation, and immune responses. There is an association between KDELR2and promotion of glioblastoma tumorigenesis. The aim of the present study was to explore the functional mechanism of KDELR2 in glioma and during response to chemotherapy to temozolomide (TMZ).

METHODS

The expression of KDELR2 in glioma tissues and cells was evaluated by immunohistochemistry, western blot and RT-qPCR assay. Then role of KDELR2 was demonstrated by CCK8, colony formation, flow cytometry and Hochest 33258 assays. The expression of genes (ATF4, ATF6, PERK, eIF2-α, GRP78 and CHOP) in U373 cells was evaluated by RT-qPCR. The protein expression of genes (cleaved caspase 3, caspase 3, cleaved PARP, PARP, Bax, Bcl-2, JNK, p-JNK, p38, p-p38, ATF4, ATF6, XBP-1s, PERK, p-PERK, GRP78 and CHOP) was measured by western blot assay.

RESULTS

The expression of KDELR2 was upregulated in high-grade gliomas tissues. KDELR2 knockdown suppressed cell proliferation but increased cell apoptosis. Further, Knockdown of KDELR2 also activated the ER stress (ERS)-dependent CHOP pathway, and resulted in increased levels of phosphorylated c-Jun N-terminal kinase (JNK) and p38. Moreover, the combination of KDELR2 knockdown and TMZ application showed a synergistic cytotoxic effect in U373 cells through the ERS-dependent CHOP and JNK/p38 pathways.

CONCLUSIONS

KDELR2 knockdown induces apoptosis and sensitizes glioma cells to TMZ, which is mediated by the ERS-dependent CHOP and JNK/p38 pathways.

Circulating Tumor DNA as a Marker for Treatment Response in Metastatic Melanoma Patients Using Next-Generation Sequencing-A Prospective Feasibility Study

We prospectively performed a longitudinal analysis of circulating tumor DNA (ctDNA) from 149 plasma samples and CT scans in Stage III and IV metastatic melanoma patients (n = 20) treated with targeted agents or immunotherapy using two custom next-generation sequencing (NGS) Ion AmpliSeq™ HD panels including 60 and 81 amplicons in 18 genes, respectively. Concordance of matching cancer-associated mutations in tissue and plasma was 73.3%. Mutant allele frequency (MAF) levels showed a range from 0.04% to 28.7%, well detectable with NGS technologies utilizing single molecule tagging like the AmpliSeq™ HD workflow. Median followup time of the tissue and/or plasma positive cohort (n = 15) was 24.6 months and median progression-free survival (PFS) was 7.8 months. Higher MAF ≥ 1% at baseline was not significantly associated with a risk of progression (Odds Ratio = 0.15; p = 0.155). Although a trend could be seen, MAF levels did not differ significantly over time between patients with and without a PFS event (p = 0.745). Depending on the cell-free DNA amount, NGS achieved a sensitivity down to 0.1% MAF and allowed for parallel analysis of multiple mutations and previously unknown mutations. Our study indicates that NGS gene panels could be useful for monitoring disease burden during therapy with ctDNA in melanoma patients.

Rethinking the biology of metastatic melanoma: a holistic approach

Over the past decades, melanoma-related mortality has remained nearly stable. The main reason is treatment failure of metastatic disease and the inherently linked knowledge gap regarding metastasis formation. In order to elicit invasion, melanoma cells manipulate the tumor microenvironment, gain motility, and adhere to the extracellular matrix and cancer-associated fibroblasts. Melanoma cells thereby express different cell adhesion molecules like laminins, integrins, N-cadherin, and others. Epithelial-mesenchymal transition (EMT) is physiological during embryologic development, but reactivated during malignancy. Despite not being truly epithelial, neural crest-derived malignancies like melanoma share similar biological programs that enable tumorigenesis, invasion, and metastasis. This complex phenomenon is termed phenotype switching and is intertwined with oncometabolism as well as dormancy escape. Additionally, it has been shown that primary melanoma shed exosomes that create a favorable premetastatic niche in the microenvironment of secondary organs and lymph nodes. Although the growing body of literature describes the aforementioned concepts separately, an integrative holistic approach is missing. Using melanoma as a tumor model, this review will shed light on these complex biological principles in an attempt to clarify the mechanistic metastatic pathways that dictate tumor and patient fate.

Development and Validation of a Hypoxia-Related Signature for Predicting Survival Outcomes in Patients With Bladder Cancer

Objectives

This study aimed to develop and validate a hypoxia signature for predicting survival outcomes in patients with bladder cancer.

Methods

We downloaded the RNA sequence and the clinicopathologic data of the patients with bladder cancer from The Cancer Genome Atlas (TCGA) (https://portal.gdc.cancer.gov/repository?facetTab=files) and the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) databases. Hypoxia genes were retrieved from the Molecular Signatures Database (https://www.gsea-msigdb.org/gsea/msigdb/index.jsp). Differentially expressed hypoxia-related genes were screened by univariate Cox regression analysis and Lasso regression analysis. Then, the selected genes constituted the hypoxia signature and were included in multivariate Cox regression to generate the risk scores. After that, we evaluate the predictive performance of this signature by multiple receiver operating characteristic (ROC) curves. The CIBERSORT tool was applied to investigate the relationship between the hypoxia signature and the immune cell infiltration, and the maftool was used to summarize and analyze the mutational data. Gene-set enrichment analysis (GSEA) was used to investigate the related signaling pathways of differentially expressed genes in both risk groups. Furthermore, we developed a model and presented it with a nomogram to predict survival outcomes in patients with bladder cancer.

Results

Eight genes (AKAP12, ALDOB, CASP6, DTNA, HS3ST1, JUN, KDELR3, and STC1) were included in the hypoxia signature. The patients with higher risk scores showed worse overall survival time than the ones with lower risk scores in the training set (TCGA) and two external validation sets (GSE13507 and GSE32548). Immune infiltration analysis showed that two types of immune cells (M0 and M1 macrophages) had a significant infiltration in the high-risk group. Tumor mutation burden (TMB) analysis showed that the risk scores between the wild types and the mutation types of TP53, MUC16, RB1, and FGFR3 were significantly different. Gene-Set Enrichment Analysis (GSEA) showed that immune or cancer-associated pathways belonged to the high-risk groups and metabolism-related signal pathways were enriched into the low-risk group. Finally, we constructed a predictive model with risk score, age, and stage and validated its performance in GEO datasets.

Conclusion

We successfully constructed and validated a novel hypoxia signature in bladder cancer, which could accurately predict patients’ prognosis.