Transcriptome and cytogenetic profiling analysis of matched in situ/invasive cutaneous squamous cell carcinomas from immunocompetent patients

The Loss of PPARγ Expression and Signaling Is a Key Feature of Cutaneous Actinic Disease and Squamous Cell Carcinoma: Association with Tumor Stromal Inflammation

Given the importance of peroxisome proliferator-activated receptor (PPAR)-gamma in epidermal inflammation and carcinogenesis, we analyzed the transcriptomic changes observed in epidermal PPARγ-deficient mice (Pparg-/-epi). A gene set enrichment analysis revealed a close association with epithelial malignancy, inflammatory cell chemotaxis, and cell survival. Single-cell sequencing of Pparg-/-epi mice verified changes to the stromal compartment, including increased inflammatory cell infiltrates, particularly neutrophils, and an increase in fibroblasts expressing myofibroblast marker genes. A comparison of transcriptomic data from Pparg-/-epi and publicly available human and/or mouse actinic keratoses (AKs) and cutaneous squamous cell carcinomas (SCCs) revealed a strong correlation between the datasets. Importantly, PPAR signaling was the top common inhibited canonical pathway in AKs and SCCs. Both AKs and SCCs also had significantly reduced PPARG expression and PPARγ activity z-scores. Smaller reductions in PPARA expression and PPARα activity and increased PPARD expression but reduced PPARδ activation were also observed. Reduced PPAR activity was also associated with reduced PPARα/RXRα activity, while LPS/IL1-mediated inhibition of RXR activity was significantly activated in the tumor datasets. Notably, these changes were not observed in normal sun-exposed skin relative to non-exposed skin. Finally, Ppara and Pparg were heavily expressed in sebocytes, while Ppard was highly expressed in myofibroblasts, suggesting that PPARδ has a role in myofibroblast differentiation. In conclusion, these data provide strong evidence that PPARγ and possibly PPARα represent key tumor suppressors by acting as master inhibitors of the inflammatory changes found in AKs and SCCs.

Single-cell sequencing highlights heterogeneity and malignant progression in actinic keratosis and cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma (cSCC) is the second most frequent of the keratinocyte-derived malignancies with actinic keratosis (AK) as a precancerous lesion. To comprehensively delineate the underlying mechanisms for the whole progression from normal skin to AK to invasive cSCC, we performed single-cell RNA sequencing (scRNA-seq) to acquire the transcriptomes of 138,982 cells from 13 samples of six patients including AK, squamous cell carcinoma in situ (SCCIS), cSCC, and their matched normal tissues, covering comprehensive clinical courses of cSCC. We identified diverse cell types, including important subtypes with different gene expression profiles and functions in major keratinocytes. In SCCIS, we discovered the malignant subtypes of basal cells with differential proliferative and migration potential. Differentially expressed genes (DEGs) analysis screened out multiple key driver genes including transcription factors along AK to cSCC progression. Immunohistochemistry (IHC)/immunofluorescence (IF) experiments and single-cell ATAC sequencing (scATAC-seq) data verified the expression changes of these genes. The functional experiments confirmed the important roles of these genes in regulating cell proliferation, apoptosis, migration, and invasion in cSCC tumor. Furthermore, we comprehensively described the tumor microenvironment (TME) landscape and potential keratinocyte-TME crosstalk in cSCC providing theoretical basis for immunotherapy. Together, our findings provide a valuable resource for deciphering the progression from AK to cSCC and identifying potential targets for anticancer treatment of cSCC.

Genetic Studies of Actinic Keratosis Development: Where Are We Now?

Actinic keratosis (AK) is a common precancerous skin lesion that can develop into cutaneous squamous cell carcinoma (CSCC). AK is characterized by atypical keratinocytes in the skin's outer layer and is commonly found in sun-exposed areas. Like many precancerous lesions, the development of AK is closely associated with genetic mutations. The molecular biology and transcriptional mechanisms underlying AK development are not well understood. Ultraviolet (UV) light exposure, especially UVA and UVB radiation, is a significant risk factor for AK, causing DNA damage and mutagenic effects. Besides UV exposure, comorbidities like diabetes, rheumatoid arthritis, and psoriasis may also influence AK development. AK patients have shown associations with various internal malignancies, indicating potential vulnerability in cancer-associated genes. Treatment for AK includes cryosurgery, electrodesiccation and curettage, chemotherapeutic creams, photodynamic therapy, or topical immune-modulators. Genomic studies have identified genetic aberrations in AK, with common mutations found in genes like TP53, NOTCH1, and NOTCH2. The progression from AK to CSCC involves chromosomal aberrations and alterations in oncogenes and tumor-suppressor genes. The functional relationships among these genes are not fully understood, but network analysis provides insights into their potential mechanisms. Further research is needed to enhance our understanding of AK's pathogenesis and develop novel therapeutic approaches.

LINC01305 recruits basonuclin 1 to act on G-protein pathway suppressor 1 to promote esophageal squamous cell carcinoma

EsophageaL squamous cell carcinoma (ESCC) is one of the most common and lethal tumors, however, its underlying molecular mechanisms are not completely understood and new therapeutic targets are needed. Here, we found that the transcription factor basonuclin 1 (BNC1) was significantly upregulated and closely related to the differentiation and metastasis of ESCC. Furthermore, BNC1, LINC01305, and G-protein pathway suppressor 1 (GPS1) had significant oncogenic roles in ESCC. In addition, in vivo experiments showed that knockdown of BNC1 indeed significantly inhibited the proliferation and metastasis of ESCC. We also revealed the molecular mechanism by which LINC01305 recruits BNC1 to the promoter of GPS1, and then GPS1 could mediate the JNK signaling pathway to promote the proliferation and metastases of ESCC. Taken together, we discovered the novel molecular mechanism by which LINC01305/BNC1 upregulates GPS1 expression to promote the development of ESCC, providing a new therapeutic target for ESCC.

Dichotomous role of miR193b-3p in diabetic foot ulcers maintains inhibition of healing and suppression of tumor formation

Despite the hyperproliferative environment marked by activation of β-catenin and overexpression of c-myc, the epidermis surrounding chronic diabetic foot ulcers (DFUs) is clinically hypertrophic and nonmigratory yet does not undergo malignant transformation. We identified miR193b-3p as a master regulator that contributes to this unique cellular phenotype. We determined that induction of tumor suppressor miR193b-3p is a unique feature of DFUs that is not found in venous leg ulcers, acute wounds, or cutaneous squamous cell carcinoma (SCC). Genomic analyses of DFUs identified suppression of the miR193b-3p target gene network that orchestrates cell motility. Inhibition of migration and wound closure was further confirmed by overexpression of miR193b-3p in human organotypic and murine in vivo wound models, whereas miR193b-3p knockdown accelerated wound reepithelialization in human ex vivo and diabetic murine wounds in vivo. The dominant negative effect of miR193b-3p on keratinocyte migration was maintained in the presence of promigratory miR31-5p and miR15b-5p, which were also overexpressed in DFUs. miR193b-3p mediated antimigratory activity by disrupting stress fiber formation and by decreasing activity of GTPase RhoA. Conversely, miR193b-3p targets that typically participate in malignant transformation were found to be differentially regulated between DFUs and SCC, including the proto-oncogenes KRAS (Kirsten rat sarcoma viral proto-oncogene) and KIT (KIT proto-oncogene). Although miR193b-3p acts as a tumor suppressor contributing to low tumor incidence in DFUs, it also acts as a master inhibitor of cellular migration and epithelialization in DFUs. Thus, miR193b-3p may represent a target for wound healing induction, cancer therapeutics, and diagnostics.

Cancer Progression Gene Expression Profiling Identifies the Urokinase Plasminogen Activator Receptor as a Biomarker of Metastasis in Cutaneous Squamous Cell Carcinoma

Cutaneous squamous cell carcinoma (cSCC) of the head and neck region is the second most prevalent skin cancer, with metastases to regional lymph nodes occurring in 2%–5% of cases. To further our understanding of the molecular events characterizing cSCC invasion and metastasis, we conducted targeted cancer progression gene expression and pathway analysis in non-metastasizing (PRI-) and metastasizing primary (PRI+) cSCC tumors of the head and neck region, cognate lymph node metastases (MET), and matched sun-exposed skin (SES). The highest differentially expressed genes in metastatic (MET and PRI+) versus non-metastatic tumors (PRI-) and SES included PLAU, PLAUR, MMP1, MMP10, MMP13, ITGA5, VEGFA, and various inflammatory cytokine genes. Pathway enrichment analyses implicated these genes in cellular pathways and functions promoting matrix remodeling, cell survival and migration, and epithelial to mesenchymal transition, which were all significantly activated in metastatic compared to non-metastatic tumors (PRI-) and SES. We validated the overexpression of urokinase plasminogen activator receptor (uPAR, encoded by PLAUR) in an extended patient cohort by demonstrating higher uPAR staining intensity in metastasizing tumors. As pathway analyses identified epidermal growth factor (EGF) as a potential upstream regulator of PLAUR, the effect of EGF on uPAR expression levels and cell motility was functionally validated in human metastatic cSCC cells. In conclusion, we propose that uPAR is an important driver of metastasis in cSCC and represents a potential therapeutic target in this disease.

Emerging precision diagnostics in advanced cutaneous squamous cell carcinoma

Advanced cutaneous squamous cell carcinoma (cSCC) encompasses unresectable and metastatic disease. Although immune checkpoint inhibition has been approved for this entity recently, a considerable proportion of cases is associated with significant morbidity and mortality. Clinical, histopathological, and radiological criteria are used for current diagnostics, classification, and therapeutic decision-making. The identification of complex molecular biomarkers to accurately stratify patients is a not yet accomplished requirement to further shift current diagnostics and care to a personalized precision medicine. This article highlights new insights into the mutational profile of cSCC, summarizes current diagnostic and therapeutic standards, and discusses emerging diagnostic approaches with emphasis on liquid biopsy and tumor tissue-based analyses.

Molecular Mechanisms of Cutaneous Squamous Cell Carcinoma

Non-melanoma skin cancers are cutaneous malignancies representing the most common form of cancer in the United States. They are comprised predominantly of basal cell carcinomas and squamous cell carcinomas (cSCC). The incidence of cSCC is increasing, resulting in substantial morbidity and ever higher treatment costs; currently in excess of one billion dollars, per annum. Here, we review research defining the molecular basis and development of cSCC that aims to provide new insights into pathogenesis and drive the development of novel, cost and morbidity saving therapies.

Clinicopathological significance and underlying molecular mechanism of downregulation of basonuclin 1 expression in ovarian carcinoma

In this study, we aim to identify the clinical significance of basonuclin 1 (BNC1) expression in ovarian carcinoma (OV) and to explore its latent mechanisms. Via integrating in-house tissue microarrays, gene chips, and RNA-sequencing data, we explored the expression and clinical value of BNC1 in OV. Immunohistochemical staining was utilized to confirm the protein expression status of BNC1. A combined SMD of -2.339 (95% CI: -3.649 to -1.028, P < 0.001) identified that BNC1 was downregulated based on 1346 samples, and the sROC (AUC = 0.93) showed a favorable discriminatory ability of BNC1 in OV patients. We used univariate and multivariate Cox regulation to evaluate the prognostic role of BNC1 for OV patients, and a combined hazard ratio of 0.717 (95% CI: 0.445-0.989, P < 0.001) revealed that BNC1 was a protective factor for OV. Furthermore, the fraction of infiltrating naive B cells, memory B cells, and other immune cells showed statistical differences between the high- and low-BNC1 expression groups through cell-type identification by estimating relative subsets of RNA transcripts (CIBERSORT) algorithm. Enrichment analysis showed that BNC1 may have a relationship with immune-related items in OV. By predicting the potential regulatory transcription factors (TFs) of BNC1, friend leukemia virus integration 1 (FLI1) may be a potential upstream TF of BNC1. Corporately, a decreasing trend of BNC1 may serve as a tumor suppressor and prognostic biomarker in OV patients. Moreover, BNC1 may take part in immune-related pathways and influence the fraction of tumor-infiltrating immune cells.

Blood Differential Gene Expression in Patients with Chronic Heart Failure and Systemic Iron Deficiency: Pathways Involved in Pathophysiology and Impact on Clinical Outcomes

Background: Iron deficiency is a common disorder in patients with heart failure and is related with adverse outcomes and poor quality of life. Previous experimental studies have shown biological connections between iron homeostasis, mitochondrial metabolism, and myocardial function. However, the mechanisms involved in this crosstalk are yet to be unfolded. Methods: The present research attempts to investigate the intrinsic biological mechanisms between heart failure and iron deficiency and to identify potential prognostic biomarkers by determining the gene expression pattern in the blood of heart failure patients, using whole transcriptome and targeted TaqMan^® low-density array analyses. Results: We performed a stepwise cross-sectional longitudinal study in a cohort of chronic heart failure patients with and without systemic iron deficiency. First, the full transcriptome was performed in a nested case-control exploratory cohort of 7 paired patients and underscored 1128 differentially expressed transcripts according to iron status (cohort1#). Later, we analyzed the messenger RNA levels of 22 genes selected by their statistical significance and pathophysiological relevance, in a validation cohort of 71 patients (cohort 2#). Patients with systemic iron deficiency presented lower mRNA levels of mitochondrial ferritin, sirtuin-7, small integral membrane protein 20, adrenomedullin and endothelin converting enzyme-1. An intermediate mitochondrial ferritin gene expression and an intermediate or low sirtuin7 and small integral membrane protein 20 mRNA levels were associated with an increased risk of all-cause mortality and heart failure admission ((HR 2.40, 95% CI 1.04–5.50, p-value = 0.039), (HR 5.49, 95% CI 1.78–16.92, p-value = 0.003), (HR 9.51, 95% CI 2.69–33.53, p-value < 0.001), respectively). Conclusions: Patients with chronic heart failure present different patterns of blood gene expression depending on systemic iron status that affect pivotal genes involved in iron regulation, mitochondrial metabolism, endothelial function and cardiovascular physiology, and correlate with adverse clinical outcomes.

The Genomic Landscape of Actinic Keratosis

Actinic keratoses (AKs) are lesions of epidermal keratinocyte dysplasia and are precursors for invasive cutaneous squamous cell carcinoma (cSCC). Identifying the specific genomic alterations driving the progression from normal skin to skin with AK to skin with invasive cSCC is challenging because of the massive UVR-induced mutational burden characteristic at all stages of this progression. In this study, we report the largest AK whole-exome sequencing study to date and perform a mutational signature and candidate driver gene analysis on these lesions. We demonstrate in 37 AKs from both immunosuppressed and immunocompetent patients that there are significant similarities between AKs and cSCC in terms of mutational burden, copy number alterations, mutational signatures, and patterns of driver gene mutations. We identify 44 significantly mutated AK driver genes and confirm that these genes are similarly altered in cSCC. We identify azathioprine mutational signature in all AKs from patients exposed to the drug, providing further evidence for its role in keratinocyte carcinogenesis. cSCCs differ from AKs in having higher levels of intrasample heterogeneity. Alterations in signaling pathways also differ, with immune-related signaling and TGFβ signaling significantly more mutated in cSCC. Integrating our findings with independent gene expression datasets confirms that dysregulated TGFβ signaling may represent an important event in AK‒cSCC progression.

Assessing SPP1/Osteopontin (OPN) Splice Variants and Their Association to Nonmelanoma Skin Cancer by Absolute Quantification: Identification of OPN-5 Subvariants and Their Protein Coding Potential

The study evaluated whether SPP1/osteopontin (OPN) splice variants are differentially expressed in nonmelanoma skin cancer compared to normal skin. The absolute number of mRNA molecules of OPN-a predominated in normal skin and nonmelanoma skin cancer compared to OPN-b, OPN-c, and OPN-5. However, mRNAs of OPN-a, OPN-b, and OPN-c were expressed in higher levels in cutaneous squamous cell carcinomas (cSCCs) and basal cell carcinomas relative to normal skin. Additionally, OPN-5 expression was higher than OPN-b and OPN-c, and OPN-c, in normal skin and nonmelanoma skin cancer, respectively. Furthermore, we identified four OPN-5 splice variants, which were cloned and analyzed for protein expression.

Loss of retinoic acid receptor-related receptor alpha (Rorα) promotes the progression of UV-induced cSCC

Cutaneous squamous cell carcinoma (cSCC) is prevalent in the world, accounting for a huge part of non-melanoma skin cancer. Most cSCCs are associated with a distinct pre-cancerous lesion, the actinic keratosis (AK). However, the progression trajectory from normal skin to AK and cSCC has not been fully demonstrated yet. To identify genes involved in this progression trajectory and possible therapeutic targets for cSCC, here we constructed a UV-induced cSCC mouse model covering the progression from normal skin to AK to cSCC, which mimicked the solar UV radiation perfectly using the solar-like ratio of UVA and UVB, firstly. Then, transcriptome analysis and a series of bioinformatics analyses and cell experiments proved that Rorα is a key transcript factor during cSCC progression. Rorα could downregulate the expressions of S100a9 and Sprr2f in cSCC cells, which can inhibit the proliferation and migration in cSCC cells, but not the normal keratinocyte. Finally, further animal experiments confirmed the inhibitory effect of cSCC growth by Rorα in vivo. Our findings showed that Rorα would serve as a potential novel target for cSCC, which will facilitate the treatment of cSCC in the future.

Invasive squamous cell carcinomas and precursor lesions on UV-exposed epithelia demonstrate concordant genomic complexity in driver genes

Although squamous cell carcinomas (SCC) are the most frequent human solid tumor at many anatomic sites, the driving molecular alterations underlying their progression from precursor lesions are poorly understood, especially in the context of photodamage. Therefore, we used high-depth, targeted next-generation sequencing (NGS) of RNA and DNA from routine tissue samples to characterize the progression of both well- (cutaneous) and poorly (ocular) studied SCCs. We assessed 56 formalin-fixed paraffin-embedded (FFPE) cutaneous lesions (n = 8 actinic keratosis, n = 30 carcinoma in situ [CIS], n = 18 invasive) and 43 FFPE ocular surface lesions (n = 2 conjunctival/corneal intraepithelial neoplasia, n = 20 CIS, n = 21 invasive), from institutions in the US and Brazil. An additional seven cases of advanced cutaneous SCC were profiled by hybrid capture-based NGS of >1500 genes. The cutaneous and ocular squamous neoplasms displayed a predominance of UV-signature mutations. Precursor lesions had highly similar somatic genomic landscapes to SCCs, including chromosomal gains of 3q involving SOX2, and highly recurrent mutations and/or loss of heterozygosity events affecting tumor suppressors TP53 and CDKN2A. Additionally, we identify a novel molecular subclass of CIS with RB1 mutations. Among TP53 wild-type tumors, human papillomavirus transcript was detected in one matched pair of cutaneous CIS and SCC. Amplicon-based whole-transcriptome sequencing of select 20 cutaneous lesions demonstrated significant upregulation of pro-invasion genes in cutaneous SCCs relative to precursors, including MMP1, MMP3, MMP9, LAMC2, LGALS1, and TNFRSF12A. Together, ocular and cutaneous squamous neoplasms demonstrate similar alterations, supporting a common model for neoplasia in UV-exposed epithelia. Treatment modalities useful for cutaneous SCC may also be effective in ocular SCC given the genetic similarity between these tumor types. Importantly, in both systems, precursor lesions possess the full complement of major genetic changes seen in SCC, supporting non-genetic drivers of invasiveness.

Mutation of NEKL-4/NEK10 and TTLL genes suppress neuronal ciliary degeneration caused by loss of CCPP-1 deglutamylase function

Ciliary microtubules are subject to post-translational modifications that act as a "Tubulin Code" to regulate motor traffic, binding proteins and stability. In humans, loss of CCP1, a cytosolic carboxypeptidase and tubulin deglutamylating enzyme, causes infantile-onset neurodegeneration. In C. elegans, mutations in ccpp-1, the homolog of CCP1, result in progressive degeneration of neuronal cilia and loss of neuronal function. To identify genes that regulate microtubule glutamylation and ciliary integrity, we performed a forward genetic screen for suppressors of ciliary degeneration in ccpp-1 mutants. We isolated the ttll-5(my38) suppressor, a mutation in a tubulin tyrosine ligase-like glutamylase gene. We show that mutation in the ttll-4, ttll-5, or ttll-11 gene suppressed the hyperglutamylation-induced loss of ciliary dye filling and kinesin-2 mislocalization in ccpp-1 cilia. We also identified the nekl-4(my31) suppressor, an allele affecting the NIMA (Never in Mitosis A)-related kinase NEKL-4/NEK10. In humans, NEK10 mutation causes bronchiectasis, an airway and mucociliary transport disorder caused by defective motile cilia. C. elegans NEKL-4 localizes to the ciliary base but does not localize to cilia, suggesting an indirect role in ciliary processes. This work defines a pathway in which glutamylation, a component of the Tubulin Code, is written by TTLL-4, TTLL-5, and TTLL-11; is erased by CCPP-1; is read by ciliary kinesins; and its downstream effects are modulated by NEKL-4 activity. Identification of regulators of microtubule glutamylation in diverse cellular contexts is important to the development of effective therapies for disorders characterized by changes in microtubule glutamylation. By identifying C. elegans genes important for neuronal and ciliary stability, our work may inform research into the roles of the tubulin code in human ciliopathies and neurodegenerative diseases.

Molecular prediction of metastasis in cutaneous squamous cell carcinoma

PURPOSE OF REVIEW

Cutaneous squamous cell carcinoma (cSCC) is a highly prevalent malignancy frequently occurring on body surfaces chronically exposed to ultraviolet radiation. While a large majority of tumors remain localized to the skin and immediate subcutaneous tissue and are cured with surgical excision, a small subset of patients with cSCC will develop metastatic disease. Risk stratification for cSCC is performed using clinical staging systems, but given a high mutational burden and advances in targeted and immunotherapy, there is growing interest in molecular predictors of high-risk disease.

RECENT FINDINGS

Recent literature on the risk for metastasis in cSCC includes notable findings in genes involved in cell-cycle regulation, tumor suppression, tissue invasion and microenvironment, interactions with the host-immune system, and epigenetic regulation.

SUMMARY

cSCC is a highly mutated tumor with complex carcinogenesis. Regulators of tumor growth and local invasion are numerous and increasingly well-understood but drivers of metastasis are less established. Areas of importance include central system regulators (NOTCH, miRNAs), proteins involved in tissue invasion (podoplanin, E-cadherin), and targets of existing and emerging therapeutics (PD-1, epidermal growth factor receptor). Given the complexity of cSCC carcinogenesis, the use of machine learning algorithms and computational genomics may provide ultimate insight and prospective studies are needed to verify clinical relevance.