A comprehensive analysis of coding and non-coding transcriptomic changes in cutaneous squamous cell carcinoma

Gene expression landscape of cutaneous squamous cell carcinoma progression

BACKGROUND

Cutaneous squamous cell carcinomas (cSCCs) are the second most common human cancer and have been characterized by RNA sequencing (RNA-Seq); however, the transferability of findings from individual studies may be limited by small sample sizes and diverse analysis protocols.

OBJECTIVES

To define the transcriptome landscape at different stages in the progression of normal skin to cSCC via a meta-analysis of publicly available RNA-Seq samples.

METHODS

Whole-transcriptome data from 73 clinically normal skin samples, 46 actinic keratoses (AK) samples, 16 in situ SCC samples, 13 keratoacanthoma (KA) samples and 147 cSCC samples [including 30 samples from immunocompromised patients and 8 from individuals with recessive dystrophic epidermolysis bullosa (RDEB)] were uniformly processed to harmonize gene expression. Differential expression, fusion detection and cell-type deconvolution analyses were performed.

RESULTS

Individual RNA-Seq studies of cSCC demonstrated study-specific clustering and varied widely in their differential gene expression detection. Following batch correction, we defined a consensus set of differentially expressed genes (DEGs), including those altered in the preinvasive stages of cSCC development, and used single-cell RNA-Seq data to demonstrate that DEGs are often - but not always - expressed by tumour-specific keratinocytes (TSKs). Analysis of the cellular composition of cSCC, KA and RDEB-cSCC identified an increase in differentiated keratinocytes in KA, while RDEB-cSCC contained the most TSKs. Compared with cSCC arising in immunocompetent individuals, cSCC samples from immunosuppressed patients demonstrated fewer memory B cells and CD8+ T cells. A comprehensive and unbiased search for fusion transcripts in cSCC and intermediate disease stages identified few candidates that recurred in >1% of all specimens, suggesting that most cSCC are not driven by oncogenic gene fusions. Finally, using Genotype-Tissue Expression (GTEx) data, we distilled a novel 300-gene signature of chronic sun exposure that affirms greater cumulative ultraviolet (UV) exposure in later stages of cSCC development.

CONCLUSIONS

Our results define the gene expression landscape of cSCC progression, characterize cell subpopulation heterogeneity in cSCC subtypes that contribute to their distinct clinical phenotypes, demonstrate that gene fusions are not a common cause of cSCC and identify UV-responsive genes associated with cSCC development.

Transcriptomic analysis of genes associated with vitamin D receptor signalling reveals differences between skin cancers

Vitamin D activates the vitamin D receptor (VDR), which dimerizes preferentially with the retinoid X receptor-α (RXRα). This heterodimer connects with genetic elements responsive to vitamin D, inhibiting or stimulating gene activity. We performed Nanostring® analysis of VDR/RXRα to compare the mRNA expression of this heterodimer and their correlated transcriptomes in non-melanoma skin cancer (basal cell carcinomas (BCC) and squamous cell carcinomas (SCC)) and melanocytic lesions (intradermal nevi (IN), and melanomas (MM)) with control skin. To evaluate VDR, RXRα and other 22 correlated genes in BCC, SCC, IN and MM, paraffin samples had their transcriptomes analysed using Nanostring®, a platform that allows multiple mRNA analyses. There were 46 samples, including 11 BCC, 10 SCC, 10 IN, 12 MM and 3 pools of control skins. Most mRNAs differed between the lesion groups and the control group. BCC and SCC NCOR2 were upregulated; in MM and IN, RXRγ was higher than in the control group. TP53, FOXO3 and MED1 showed a significant difference when we compared the BCC group to the SCC group. Melanoma and intradermal nevi differed only in AhR. VDR and RXRα were lower than the control in all groups. The panel shows a clear difference between the non-melanocytic cancers and, on the other hand, a slight difference between the melanocytic lesions. The study of vitamin D's influence through its receptor and RXRα is an exciting issue for understanding the importance of this pathway, and the present study can impact the prevention and treatment strategies, mainly in non-melanocytic tumours.

The Mutational and Microenvironmental Landscape of Cutaneous Squamous Cell Carcinoma: A Review

Cutaneous squamous cell carcinoma (cSCC) manifests through the complex interactions of UV-induced DNA damage, genetic mutations, and alterations in the tumor microenvironment. A high mutational burden is present in cSCC, as well as both cSCC precursors and normal skin, making driver genes difficult to differentiate. Despite this, several key driver genes have been identified, including TP53, the NOTCH family, CDKN2A, PIK3CA, and EGFR. In addition to mutations, the tumor microenvironment and the manipulation and evasion of the immune system play a critical role in cSCC progression. Novel therapeutic approaches, such as immunotherapy and EGFR inhibitors, have been used to target these dysregulations, and have shown promise in treating advanced cSCC cases, emphasizing the need for targeted interventions considering both genetic and microenvironmental factors for improved patient outcomes.

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.

Clustering of RNA co-expression network identifies novel long non-coding RNA biomarkers in squamous cell carcinoma

Long non-coding RNAs (lncRNAs) have emerged as important players in cancer progression. Cutaneous squamous cell carcinoma (cSCC) is the most common metastatic skin cancer with increasing incidence worldwide. The prognosis of the metastatic cSCC is poor, and currently there are no established biomarkers to predict metastasis risk or specific therapeutic targets for advanced or metastatic cSCC. To elucidate the role of lncRNAs in cSCC, RNA sequencing of patient derived cSCC cell lines and normal human epidermal keratinocytes was performed. The correlation analysis of differentially expressed lncRNAs and protein-coding genes revealed six distinct gene clusters with one of the upregulated clusters featuring genes associated with cell motility. Upregulation of the expression of lncRNAs linked to cSCC cell motility in cSCC and head and neck SCC (HNSCC) cells was confirmed using qRT-PCR. Elevated expression of HOTTIP and LINC00543 was also noted in SCC tumors in vivo and was associated with poorer prognosis in HNSCC and lung SCC cohorts within TCGA data, respectively. Altogether, these findings uncover a novel set of lncRNAs implicated in cSCC cell locomotion. These lncRNAs may serve as potential novel biomarkers and as putative therapeutic targets for locally advanced and metastatic cSCC.

Long noncoding RNA plasmacytoma variant translocation 1 is overexpressed in cutaneous squamous cell carcinoma and exon 2 is critical for its oncogenicity

BACKGROUND

Cutaneous squamous cell carcinoma (cSCC) is one of the most common and fastest increasing forms of cancer worldwide with metastatic potential. Long noncoding RNAs (lncRNAs) are a group of RNA molecules with essential regulatory functions in both physiological and pathological processes.

OBJECTIVES

To investigate the function and mode of action of lncRNA plasmacytoma variant translocation 1 (PVT1) in cSCC.

METHODS

Quantitative reverse transcriptase polymerase chain reaction and single-molecule in situ hybridization were used to quantify the expression level of PVT1 in normal skin, premalignant skin lesions, actinic keratosis (AK) and primary and metastatic cSCCs. The function of PVT1 in cSCC was investigated both in vivo (tumour xenografts) and in vitro (competitive cell growth assay, 5-ethynyl-2'-deoxyuridine incorporation assay, colony formation assay and tumour spheroid formation assay) upon CRISPR-Cas9-mediated knockout of the entire PVT1 locus, the knockout of exon 2 of PVT1, and locked nucleic acid (LNA) gapmer-mediated PVT1 knockdown. RNA sequencing analysis was conducted to identify genes and processes regulated by PVT1.

RESULTS

We identified PVT1 as a lncRNA upregulated in cSCC in situ and cSCC, associated with the malignant phenotype of cSCC. We showed that the expression of PVT1 in cSCC was regulated by MYC. Both CRISPR-Cas9 deletion of the entire PVT1 locus and LNA gapmer-mediated knockdown of PVT1 transcript impaired the malignant behaviour of cSCC cells, suggesting that PVT1 is an oncogenic transcript in cSCC. Furthermore, knockout of PVT1 exon 2 inhibited cSCC tumour growth both in vivo and in vitro, demonstrating that exon 2 is a critical element for the oncogenic role of PVT1. Mechanistically, we showed that PVT1 was localized in the cell nucleus and its deletion resulted in cellular senescence, increased cyclin-dependent kinase inhibitor 1 (p21/CDKN1A) expression and cell cycle arrest.

CONCLUSIONS

Our study revealed a previously unrecognized role for exon 2 of PVT1 in its oncogenic role and that PVT1 suppresses cellular senescence in cSCC. PVT1 may be a potential biomarker and therapeutic target in 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.

MicroRNA-23b Plays a Tumor-Suppressive Role in Cutaneous Squamous Cell Carcinoma and Targets Ras-Related Protein RRAS2

Cutaneous squamous cell carcinoma (cSCC) is one of the most common types of cancer with metastatic potential. MicroRNAs regulate gene expression at the post-transcriptional level. In this study, we report that miR-23b is downregulated in cSCCs and in actinic keratosis and that its expression is regulated by the MAPK signaling pathway. We show that miR-23b suppresses the expression of a gene network associated with key oncogenic pathways and that the miR-23b-gene signature is enriched in human cSCCs. miR-23b decreased the expression of FGF2 both at mRNA and protein levels and impaired the angiogenesis-inducing ability of cSCC cells. miR23b overexpression suppressed the capacity of cSCC cells to form colonies and spheroids, whereas the CRISPR/Cas9-mediated deletion of MIR23B resulted in increased colony and tumor sphere formation in vitro. In accordance with this, miR-23b-overexpressing cSCC cells formed significantly smaller tumors upon injection into immunocompromised mice with decreased cell proliferation and angiogenesis. Mechanistically, we verify RRAS2 as a direct target of miR-23b in cSCC. We show that RRAS2 is overexpressed in cSCC and that interference with its expression impairs angiogenesis and colony and tumorsphere formation. Taken together, our results suggest that miR-23b acts in a tumor-suppressive manner in cSCC, and its expression is decreased during squamous carcinogenesis.

Spatial Profiling of Circular RNAs in Cancer Reveals High Expression in Muscle and Stromal Cells

Circular RNAs (circRNA) are covalently closed molecules that can play important roles in cancer development and progression. Hundreds of differentially expressed circRNAs between tumors and adjacent normal tissues have been identified in studies using RNA sequencing or microarrays, emphasizing a strong translational potential. Most previous studies have been performed using RNA from bulk tissues and lack information on the spatial expression patterns of circRNAs. Here, we showed that the majority of differentially expressed circRNAs from bulk tissue analyses of colon tumors relative to adjacent normal tissues were surprisingly not differentially expressed when comparing cancer cells directly with normal epithelial cells. Manipulating the proliferation rates of cells grown in culture revealed that these discrepancies were explained by circRNAs accumulating to high levels in quiescent muscle cells due to their high stability; on the contrary, circRNAs were diluted to low levels in the fast-proliferating cancer cells due to their slow biogenesis rates. Thus, different subcompartments of colon tumors and adjacent normal tissues exhibited striking differences in circRNA expression patterns. Likewise, the high circRNA content in muscle cells was also a strong confounding factor in bulk analyses of circRNAs in bladder and prostate cancers. Together, these findings emphasize the limitations of using bulk tissues for studying differential circRNA expression in cancer and highlight a particular need for spatial analysis in this field of research.

SIGNIFICANCE

The abundance of circRNAs varies systematically between subcompartments of solid tumors and adjacent tissues, implying that differentially expressed circRNAs discovered in bulk tissue analyses may reflect differences in cell type composition between samples.

Circ_TNFRSF21 promotes cSCC metastasis and M2 macrophage polarization via miR-214-3p/CHI3L1

BACKGROUND

Cutaneous squamous cell carcinoma (cSCC) is a highly invasive disease with the potential to metastasize and cause fatality. Therefore, it is crucial to understand the mechanism behind cSCC in order to devise effective strategies to combat this disease.

OBJECTIVE

We investigated the function of circ_TNFRSF21/miR-214-3p/CHI3L1 axis in cSCC.

METHODS

The features of circ_TNFRSF21 was characterized using Sanger sequencing, and RNase R/actinomycin D treatment. Genes and M1/M2 markers levels were assessed by qRT-PCR and IHC. The proliferation, migration, and invasion of cells were evaluated by CCK-8, colony formation, EdU incorporation, and transwell assays. Tumor growth and metastasis in vivo were evaluated by nude mouse xenograft model. Interactions of circ_TNFRSF21/miR-214-3p and miR-214-3p/CHI3L1 were validated by RNA immunoprecipitation and dual luciferase assay.

RESULTS

Circ_TNFRSF21 and CHI3L1 expression were elevated in both human cSCC tissues and cells, whereas miR-214-3p was reduced. Circ_TNFRSF21 silencing or miR-214-3p overexpression suppressed cSCC cell proliferation, migration, invasion, and M2 macrophage polarization. Circ_TNFRSF21 functioned as a sponge for miR-214-3p while miR-214-3p directly targeted CHI3L1. Knockdown of miR-214-3p reversed the effects of circ_TNFRSF21 knockdown on cSCC development, while CHI3L1 upregulation reversed the effects of miR-214-3p overexpression. Furthermore, knockdown of circ_TNFRSF21 inhibited cSCC tumor growth and metastasis in vivo.

CONCLUSION

Circ_TNFRSF21 plays a significant role in cSCC progression by enhancing cell proliferation, migration, invasion, and M2 macrophage polarization through inhibiting miR-214-3p and subsequent disinhibition of CHI3L1. These findings deepen our understanding of the molecular mechanism of cSCC and propose the circ_TNFRSF21/miR-214-3p/CHI3L1 axis as promising diagnosis markers or therapeutic targets for cSCC.

Transcriptomic Study on Human Skin Samples: Identification of Two Subclasses of Actinic Keratoses

Actinic keratoses (AKs) are sun-damaged skin areas that affect 20% of the European adult population and more than 50% of people aged 70 years and over. There are currently no clinical or histological features allowing us to identify to which clinical class (i.e., regression or progression) an AK belongs. A transcriptomic approach seems to be a robust tool for AK characterization, but there is a need for additional studies, including more patients and elucidating the molecular signature of an AK. In this context, the present study, including the largest number of patients to date, is the first aiming at identifying biological features to objectively distinguish different AK signatures. We highlight two distinct molecular profiles: AKs featuring a molecular profile similar to squamous cell carcinomas (SCCs), which are called "lesional AKs" (AK_Ls), and AKs featuring a molecular profile similar to normal skin tissue, which are called "non-lesional AKs" (AK_NLs). The molecular profiles of both AK subclasses were studied, and 316 differentially expressed genes (DEGs) were identified between the two classes. The 103 upregulated genes in AK_L were related to the inflammatory response. Interestingly, downregulated genes were associated with keratinization. Finally, based on a connectivity map approach, our data highlight that the VEGF pathway could be a promising therapeutic target for high-risk lesions.

Molecular Alterations in Cutaneous Squamous Cell Carcinoma in Immunocompetent and Immunosuppressed Hosts-A Systematic Review

The characterization of cutaneous squamous cell carcinoma (cSCC) at the molecular level is lacking in the current literature due to the high mutational burden of this disease. Immunosuppressed patients afflicted with cSCC experience considerable morbidity and mortality. In this article, we review the molecular profile of cSCC among the immunosuppressed and immunocompetent populations at the genetic, epigenetic, transcriptomic, and proteometabolomic levels, as well as describing key differences in the tumor immune microenvironment between these two populations. We feature novel biomarkers from the recent literature which may serve as potential targets for therapy.

LINP1 represses unfolded protein response by directly inhibiting eIF2α phosphorylation to promote cutaneous squamous cell carcinoma

BACKGROUND

Endoplasmic reticulum stress (ER stress) may destroy endoplasmic reticulum homeostasis (ER homeostasis) and leads to programmable cell death. Unfolded protein response (UPR) originally stimulated by ER stress is critical for the survival of tumor cells through trying to re-establish ER homeostasis as an adaption to harsh microenvironment. However, mechanisms involving key regulators in modulating UPR remain underexplored.

METHODS

The expression of LINP1 in cutaneous squamous cell carcinoma (cSCC) tissues and cell lines was assessed. Subsequently, LINP1 was knocked out, knocked down or overexpressed in cSCC cells. CCK-8 assays, colony forming assays, transwell migration assays and invasiveness measurement by matrigel-coated transwell were performed to examine the role of LINP1 in cSCC development through gain-of-function and loss-of-function experiments. Transcriptomic sequencing (RNA-Seq) was conducted and indicated the key downstream signaling events regulated by LINP1 including UPR and apoptosis signaling. Furthermore, the direct interaction between LINP1 and eIF2α to modulate UPR and apoptosis was confirmed by RNA pulldown, RNA immunoprecipitation (RIP), ChIP-qPCR and in vitro phosphorylation assays.

RESULTS

In this study, LncRNA in non-homologous end joining pathway 1 (LINP1) was identified to be one of the top ten highest-expressed LncRNAs in cSCC, the second most common cancer in the world. Functional studies using in vitro and in vivo models revealed that LINP1 functions as an oncogene to promote cell proliferation, colony formation, migration and invasiveness while inhibiting cell apoptosis in cSCC. Transcriptomic sequencing after knockdown of LINP1 indicated LINP1 negatively regulates UPR-related pathways involving key effectors for activating UPR and the apoptosis following the prolonged UPR. Mechanistic study showed LINP1 physically interacts with eIF2α to inhibit its phosphorylation for avoiding unmitigated UPR. Loss of LINP1 followed by enhanced eIF2α phosphorylation led to overactivated UPR and induced DDIT3 expression, contributing to ER stress-induced apoptosis and suppression of cSCC development.

CONCLUSIONS

Our findings demonstrate a novel regulatory hierarchy of UPR by demonstrating LINP1 as a critical modulator for eIF2α phosphorylation and a suppressor of UPR-mediated apoptosis, which suggests a novel therapeutic target for cSCC treatment.

Tirbanibulin for Actinic Keratosis: Insights into the Mechanism of Action

Actinic keratosis (AK) is a common pre-neoplastic skin lesion constituted by uncontrolled proliferation of atypical keratinocytes that may evolve to squamous cell carcinoma. With global prevalence increasing, AK is expected to be the most common carcinoma of the skin. Tirbanibulin is a reversible tubulin polymerization inhibitor with potent anti-proliferative and anti-tumoral effects. In-vivo and in-vitro studies have shown that tirbanibulin significantly inhibits cell proliferation, tumor growth and downregulates Src signaling with no overt toxicity. Early phase and Phase III trials have shown high lesion clearance, compliance, and few side effects of once daily tirbanibulin treatment. This review discusses tirbanibulin anti-cancer activity, focusing on tubulin polymerization and Src signaling inhibitory effects, highlighting relevant literature and novel preclinical results from the ATNXUS-KX01-001 study. Furthermore, we address the relevant findings obtained in recent clinical trials to evaluate the safety, efficacy, pharmacokinetics, clearance efficacy, and side effects of the 1% tirbanibulin ointment applied once daily. In summary, we highlight preclinical and clinical evidence on the use of tirbanibulin as an effective and safe treatment option for AK.

pTINCR microprotein promotes epithelial differentiation and suppresses tumor growth through CDC42 SUMOylation and activation

The human transcriptome contains thousands of small open reading frames (sORFs) that encode microproteins whose functions remain largely unexplored. Here, we show that TINCR lncRNA encodes pTINCR, an evolutionary conserved ubiquitin-like protein (UBL) expressed in many epithelia and upregulated upon differentiation and under cellular stress. By gain- and loss-of-function studies, we demonstrate that pTINCR is a key inducer of epithelial differentiation in vitro and in vivo. Interestingly, low expression of TINCR associates with worse prognosis in several epithelial cancers, and pTINCR overexpression reduces malignancy in patient-derived xenografts. At the molecular level, pTINCR binds to SUMO through its SUMO interacting motif (SIM) and to CDC42, a Rho-GTPase critical for actin cytoskeleton remodeling and epithelial differentiation. Moreover, pTINCR increases CDC42 SUMOylation and promotes its activation, triggering a pro-differentiation cascade. Our findings suggest that the microproteome is a source of new regulators of cell identity relevant for cancer.

A cis-regulatory lexicon of DNA motif combinations mediating cell-type-specific gene regulation

Gene expression is controlled by transcription factors (TFs) that bind cognate DNA motif sequences in cis-regulatory elements (CREs). The combinations of DNA motifs acting within homeostasis and disease, however, are unclear. Gene expression, chromatin accessibility, TF footprinting, and H3K27ac-dependent DNA looping data were generated and a random-forest-based model was applied to identify 7,531 cell-type-specific cis-regulatory modules (CRMs) across 15 diploid human cell types. A co-enrichment framework within CRMs nominated 838 cell-type-specific, recurrent heterotypic DNA motif combinations (DMCs), which were functionally validated using massively parallel reporter assays. Cancer cells engaged DMCs linked to neoplasia-enabling processes operative in normal cells while also activating new DMCs only seen in the neoplastic state. This integrative approach identifies cell-type-specific cis-regulatory combinatorial DNA motifs in diverse normal and diseased human cells and represents a general framework for deciphering cis-regulatory sequence logic in gene regulation.

Whole-genome identification and construction of the lncRNA-mRNA co-expression network in patients with actinic keratosis

BACKGROUND

Actinic keratosis (AK) is a common premalignant lesion induced by chronic exposure to ultraviolet radiation and may develop into invasive cutaneous squamous carcinoma (cSCC). The identification of specific biomarkers in AK are still unclear. Long non-coding RNAs (lncRNAs), as transcripts of more than 200 nucleotides, significantly involving in multiple biologic processes, especially in the development of tumors.

METHODS

In our study, we obtained data from RNA-sequencing analysis using two AK lesion tissues and three normal cutaneous tissues to comparatively analyze the differentially expressed (DE) lncRNAs and messenger RNAs (mRNAs). Firstly, we used microarray analyses to identify DE lncRNAs and DE mRNAs. Secondly, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis to analyze the primary function and find out significant pathways of these DE mRNA and lncRNAs. Finally, we used the top ten DE lncRNAs to construct a lncRNA-mRNA co-expression network.

RESULTS

Our results showed that there were a total of 2,097 DE lncRNAs and 2,043 DE mRNAs identified. GO and KEGG analysis and the lncRNA-mRNA co-expression network (using the top 10 DE lncRNAs comprises 130 specific co-expressed mRNAs to construct) indicated that lncRNA uc011fnr.2 may negatively regulate SCIMP and Toll-like receptor 4 (TLR4) and play an important role in Janus kinase-signal transducer and activator of transcription 3 (JAK-STAT3) signaling pathway of AK.

CONCLUSIONS

lncRNA uc011fnr.2 may play an important role in JAK-STAT3 signaling pathway of AK by modulating SCIMP, TLR4 and IL-6. Further research is required to validate the value of lncRNA uc011fnr.2 in the progression of AK.

Unravelling the landscape of skin cancer through single-cell transcriptomics

The human skin is a complex organ that forms the first line of defense against pathogens and external injury. It is composed of a wide variety of cells that work together to maintain homeostasis and prevent disease, such as skin cancer. The exponentially rising incidence of skin malignancies poses a growing public health challenge, particularly when the disease course is complicated by metastasis and therapeutic resistance. Recent advances in single-cell transcriptomics have provided a high-resolution view of gene expression heterogeneity that can be applied to skin cancers to define cell types and states, understand disease evolution, and develop new therapeutic concepts. This approach has been particularly valuable in characterizing the contribution of immune cells in skin cancer, an area of great clinical importance given the increasing use of immunotherapy in this setting. In this review, we highlight recent skin cancer studies utilizing bulk RNA sequencing, introduce various single-cell transcriptomics approaches, and summarize key findings obtained by applying single-cell transcriptomics to skin cancer.

Knockdown of circRNA circEBF1 inhibits cutaneous squamous cell carcinoma progression by regulating the miR-1247-5p/CSF3 axis

BACKGROUND

Cutaneous squamous cell carcinoma (CSCC) is one of the causes of cancer-related death worldwide. Circular RNAs (circRNAs) play a vital role in the pathological process of many malignant tumors. This study aimed to explore the specific role and potential mechanism of circRNA EBF transcription factor 1 (circEBF1) in CSCC.

METHODS

The levels of circEBF1, microRNA-1247-5p (miR-1247-5p) and colony stimulating factor 3 (CSF3) were determined by quantitative real-time PCR and Western blot. Cell proliferation was assessed by colony formation assay. Cell migration, invasion and apoptosis were evaluated by Transwell, wound healing, flow cytometry and Western blot assays. Cellular glycolysis, including glucose consumption, lactate production, and ATP level, was detected by commercial kits. The binding relationship between miR-1247-5p and circEBF1 or CSF3 was verified by dual-luciferase reporter assay, RNA immunoprecipitation (RIP) assay and RNA pull-down assay. Xenograft experiment was conducted to analyze tumor growth in vivo.

RESULTS

circEBF1 and CSF3 were up-regulated, while miR-1247-5p was down-regulated in CSCC tissues and cells. Interference of circEBF1 hindered the proliferation, migration, invasion, and glycolysis of CSCC cells and promoted apoptosis. In addition, circEBF1 sponged miR-1247-5p to regulate CSCC cell development. Also, miR-1247-5p suppressed CSCC cell progression by inhibiting CSF3. Moreover, depletion of circEBF1 blocked CSCC tumor growth in vivo.

CONCLUSION

Down-regulation of circEBF1 impeded CSCC progression by regulating the miR-1247-5p/CSF3 pathway, suggesting that circEBF1 might be a promising therapeutic target for CSCC.

Multiplex Analysis of CircRNAs from Plasma Extracellular Vesicle-Enriched Samples for the Detection of Early-Stage Non-Small Cell Lung Cancer

BACKGROUND

The analysis of liquid biopsies brings new opportunities in the precision oncology field. Under this context, extracellular vesicle circular RNAs (EV-circRNAs) have gained interest as biomarkers for lung cancer (LC) detection. However, standardized and robust protocols need to be developed to boost their potential in the clinical setting. Although nCounter has been used for the analysis of other liquid biopsy substrates and biomarkers, it has never been employed for EV-circRNA analysis of LC patients.

METHODS

EVs were isolated from early-stage LC patients (n = 36) and controls (n = 30). Different volumes of plasma, together with different number of pre-amplification cycles, were tested to reach the best nCounter outcome. Differential expression analysis of circRNAs was performed, along with the testing of different machine learning (ML) methods for the development of a prognostic signature for LC.

RESULTS

A combination of 500 μL of plasma input with 10 cycles of pre-amplification was selected for the rest of the study. Eight circRNAs were found upregulated in LC. Further ML analysis selected a 10-circRNA signature able to discriminate LC from controls with AUC ROC of 0.86.

CONCLUSIONS

This study validates the use of the nCounter platform for multiplexed EV-circRNA expression studies in LC patient samples, allowing the development of prognostic signatures.

Role of Circular RNAs in Pulmonary Fibrosis

Pulmonary fibrosis is a chronic progressive form of interstitial lung disease, characterized by the histopathological pattern of usual interstitial pneumonia. Apart from aberrant alterations of protein-coding genes, dysregulation of non-coding RNAs, including microRNAs, long non-coding RNAs, and circular RNAs (circRNAs), is crucial to the initiation and progression of pulmonary fibrosis. CircRNAs are single-stranded RNAs that form covalently closed loops without 5′ caps and 3′ tails. Different from canonical splicing of mRNA, they are produced from the back-splicing of precursor mRNAs and have unique biological functions, as well as potential biomedical implications. They function as important gene regulators through multiple actions, including sponging microRNAs and proteins, regulating transcription, and splicing, as well as protein-coding and translation in a cap-independent manner. This review comprehensively summarizes the alteration and functional role of circRNAs in pulmonary fibrosis, with a focus on the involvement of the circRNA in the context of cell-specific pathophysiology. In addition, we discuss the diagnostic and therapeutic potential of targeting circRNA and their regulatory pathway mediators, which may facilitate the translation of recent advances from bench to bedside in the future.

MAB21L4 Deficiency Drives Squamous Cell Carcinoma via Activation of RET

Epithelial squamous cell carcinomas (SCC) most commonly originate in the skin, where they display disruptions in the normally tightly regulated homeostatic balance between keratinocyte proliferation and terminal differentiation. We performed a transcriptome-wide screen for genes of unknown function that possess inverse expression patterns in differentiating keratinocytes compared with cutaneous SCC (cSCC), leading to the identification of MAB21L4 (C2ORF54) as an enforcer of terminal differentiation that suppresses carcinogenesis. Loss of MAB21L4 in human cSCC organoids increased expression of RET to enable malignant progression. In addition to transcriptional upregulation of RET, deletion of MAB21L4 preempted recruitment of the CacyBP-Siah1 E3 ligase complex to RET and reduced its ubiquitylation. In SCC organoids and in vivo tumor models, genetic disruption of RET or selective inhibition of RET with BLU-667 (pralsetinib) suppressed SCC growth while inducing concomitant differentiation. Overall, loss of MAB21L4 early during SCC development blocks differentiation by increasing RET expression. These results suggest that targeting RET activation is a potential therapeutic strategy for treating SCC.

SIGNIFICANCE

Downregulation of RET mediated by MAB21L4-CacyBP interaction is required to induce epidermal differentiation and suppress carcinogenesis, suggesting RET inhibition as a potential therapeutic approach in squamous cell carcinoma.

M2 macrophage facilitated angiogenesis in cutaneous squamous cell carcinoma via circ_TNFRSF21/miR-3619-5p/ROCK axis

In recent years, the role of circular RNA in cancer cells has been studied broadly; however, the functional significance of circular RNA in the regulation of the tumor microenvironment (TME) is not fully understood. In this study, we aimed to reveal the role of circ_TNFRSF21 in M2 macrophage-induced cutaneous squamous cell carcinoma (cSCC) angiogenesis. Quantitative polymerase chain reaction and Western blotting were performed to determine the levels of the indicated genes. Direct binding between circ_TNFRSF21 and miR-3619-5p, miR-3619-5p, and ROCK2 was verified by dual-luciferase activity. The migration and invasion of human umbilical vein endothelial cells were evaluated by wound healing and transwell assays. Tube formation was performed to detect in vitro angiogenesis. Circ_TNFRSF21 and ROCK2 were upregulated in cSCC tissue, while miR-3619-5p was downregulated. Circ_TNFRSF21 negatively regulated the expression of miR-3619-5p, while miR-3619-5p negatively regulated the expression of ROCK2. miR-3619-5p suppressed tube formation by inhibiting ROCK signaling. M2 macrophages facilitated tube formation via the circ_TNFRSF21/miR-3619-5p/ROCK2 axis. Our present study revealed that circ_TNFRSF21 was elevated in M2 macrophages and mediated M2 macrophage-induced tube formation in vitro.

Non-Melanoma Skin Cancer: A Genetic Update and Future Perspectives

Simple Summary

Non-melanoma skin cancer (NMSC) is the main type of cancer in the Caucasian population, and the number of cases continues to rise. Research mostly focuses on clinical characteristics analysis, but genetic features are crucial to malignancies’ establishment and advance. We aim to explore the genetic basics of skin cancer, surrounding microenvironment interactions, and regulation mechanisms to provide a broader perspective for new therapies’ development.

Abstract

Skin cancer is one of the main types of cancer worldwide, and non-melanoma skin cancer (NMSC) is the most frequent within this group. Basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) are the most common types. Multifactorial features are well-known for cancer development, and new hallmarks are gaining relevance. Genetics and epigenetic regulation play an essential role in cancer susceptibility and progression, as well as the variety of cells and molecules that interact in the tumor microenvironment. In this review, we provide an update on the genetic features of NMSC, candidate genes, and new therapies, considering diverse perspectives of skin carcinogenesis. The global health situation and the pandemic have been challenging for health care systems, especially in the diagnosis and treatment of patients with cancer. We provide innovative approaches to overcome the difficulties in the current clinical dynamics.

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.

Comprehensive analysis of lncRNA-mRNAs co-expression network identifies potential lncRNA biomarkers in cutaneous squamous cell carcinoma

Background

Cutaneous squamous cell carcinoma (cSCC) is the second most common type of skin cancer, the prognosis for patients with metastatic cSCC remains relatively poor. Thus, there is an urgent need to identify new diagnostic, prognostic, and therapeutic targets and pathways in cSCC.

Results

It detected a total of 37,507 lncRNA probes and 32,825 mRNA probes and found 3593 differentially expressed lncRNAs and 3236 differentially expressed mRNAs. It has been found that mRNAs ACY3, NR1D1, MZB1 has co-expression relationship with six lncRNAs, GXYLT1P3, LINC00348, LOC101928131, A-33-p3340852, A-21-p0003442 and LOC644838.

Conclusions

The aim of this study is to identify cSCC-specific lncRNAs and indicated that six unstudied lncRNAs may serve an important role in endoplasmic reticulum stress apoptosis, autophagy and the progression of cSCC by modulating ACY3, NR1D1 and MZB1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08481-0.

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.

The emerging roles of circRNAs in cancer and oncology

Over the past decade, circular RNAs (circRNAs) have emerged as a large class of primarily non-coding RNA molecules, many of which have key roles in cancer development and progression through diverse mechanisms of action. CircRNAs often have tissue-restricted and cancer-specific expression patterns, and accumulating data suggest that these molecules are of potential clinical relevance and utility. In particular, circRNAs have strong potential as diagnostic, prognostic and predictive biomarkers, which is underscored by their detectability in liquid biopsy samples such as in plasma, saliva and urine. However, technical issues in the detection and assessment of circRNAs as well as biological knowledge gaps need to be addressed to move this relatively young field of research forward and bring circRNAs to the forefront of clinical practice. Herein, we review the current knowledge regarding circRNA biogenesis, regulation and functions in cancer as well as their clinical potential as biomarkers, therapeutic agents and drug targets.

The transcriptional landscape and biomarker potential of circular RNAs in prostate cancer

Background

Circular RNAs (circRNAs) constitute a largely unexplored source for biomarker discovery in prostate cancer (PC). Here, we characterize the biomarker potential of circRNAs in PC, where the need for novel diagnostic and prognostic tools to facilitate more personalized management is pressing.

Methods

We profiled the transcriptomic landscape of circRNAs in PC by total RNA sequencing of 31 adjacent-normal and 143 tumor samples from localized (radical prostatectomy (RP)) and metastatic PC patients (cohort 1, training). Diagnostic and prognostic potential was evaluated in cohort 1, and 39 top circRNA candidates were selected for validation in two additional PC cohorts (cohort 2, n = 111; RP cohort 3, n = 191) by NanoString-based expression analysis. Biochemical recurrence (BCR)-free survival was assessed using Kaplan-Meier, univariate, and multivariate Cox regression analyses. The circRNA candidates were further detected in extracellular vesicle (EV)-enriched plasma samples from PC patients and controls (cohort 4, n = 54).

Results

Expression of circABCC4, circFAT3, circATRNL1, and circITGA7 was highly cancer-specific (area under the curve 0.71–0.86), while low circITGA7 expression was significantly (P < 0.05) associated with BCR in univariate analysis in two RP cohorts. Moreover, we successfully trained and validated a novel 5-circRNA prognostic signature (circKMD1A/circTULP4/circZNF532/circSUMF1/circMKLN1) significantly associated with BCR beyond routine clinicopathological variables (RP cohort 1: P = 0.02, hazard ratio = 2.1; RP cohort 3: P < 0.001, hazard ratio = 2.1). Lastly, we provide proof-of-principle for detection of candidate circRNAs in EV-enriched plasma samples from PC patients.

Conclusions

circRNAs hold great biomarker potential in PC and display both high cancer specificity and association to disease progression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-01009-3.

CCT3-LINC00326 axis regulates hepatocarcinogenic lipid metabolism

OBJECTIVE

To better comprehend transcriptional phenotypes of cancer cells, we globally characterised RNA-binding proteins (RBPs) to identify altered RNAs, including long non-coding RNAs (lncRNAs).

DESIGN

To unravel RBP-lncRNA interactions in cancer, we curated a list of ~2300 highly expressed RBPs in human cells, tested effects of RBPs and lncRNAs on patient survival in multiple cohorts, altered expression levels, integrated various sequencing, molecular and cell-based data.

RESULTS

High expression of RBPs negatively affected patient survival in 21 cancer types, especially hepatocellular carcinoma (HCC). After knockdown of the top 10 upregulated RBPs and subsequent transcriptome analysis, we identified 88 differentially expressed lncRNAs, including 34 novel transcripts. CRISPRa-mediated overexpression of four lncRNAs had major effects on the HCC cell phenotype and transcriptome. Further investigation of four RBP-lncRNA pairs revealed involvement in distinct regulatory processes. The most noticeable RBP-lncRNA connection affected lipid metabolism, whereby the non-canonical RBP CCT3 regulated LINC00326 in a chaperonin-independent manner. Perturbation of the CCT3-LINC00326 regulatory network led to decreased lipid accumulation and increased lipid degradation in cellulo as well as diminished tumour growth in vivo.

CONCLUSIONS

We revealed that RBP gene expression is perturbed in HCC and identified that RBPs exerted additional functions beyond their tasks under normal physiological conditions, which can be stimulated or intensified via lncRNAs and affected tumour growth.

Long non-coding RNA PICSAR serves as a non-invasive biomarker for the diagnosis and prognosis of cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma (cSCC) is one of the most common malignant tumors in dermatology, and its incidence is increasing year by year. Long non-coding RNAs (lncRNAs) play vital roles in the processes of various malignant tumors. This study aimed to evaluate the expression of long non-coding RNA PICSAR and investigate whether serum PICSAR could serve as a biomarker for the diagnosis and prognosis of cSCC. The expression level of PICSAR was measured using quantitative Real-Time PCR. The diagnostic value of PICSAR was evaluated by receiver operating characteristic (ROC) analysis. Survival curves were established using the Kaplan-Meier method, and the log-rank test was used to compare differences between the two curves. Prognostic value of PICSAR was confirmed by Cox regression analysis. The expression of PICSAR was upregulated in serum of cSCC patients and tumor tissues of patients. Additionally, serum PICSAR expression had relatively high diagnostic accuracy for the screening of cSCC. Moreover, PICSAR expression was correlated with tumor size, grade of differentiation and TNM stages, and significantly increased in cSCC patients with poor tumor differentiation and cSCC patients with III-IV TNM stage. Furthermore, patients with high PICSAR expression had lower overall survival than the patients with low PICSAR expression, and PICSAR expression was an independent prognostic factors for cSCC patients. The results of this study indicated that PICSAR was upregulated in cSCC patients and tumor tissues and might serve as a non-invasive biomarker for the diagnosis and prognosis of cSCC patients.

Detection of viral gene expression in risk-stratified biopsies reveals no active HPV in cutaneous squamous cell carcinoma

BACKGROUND

Human papillomavirus (HPV) infection is known to promote the development of mucosal squamous cell carcinoma (mSCC), including pathologically high-grade lesions, but its role in cutaneous squamous cell carcinoma (cuSCC) remains unclear, particularly in lesions that are considered high risk.

OBJECTIVE

We aimed to determine whether enhanced HPV transcriptional activity can be detected in high-risk cuSCC samples compared with low-grade SCC samples or normal skin.

METHODS

We performed RNA sequencing of cuSCC across 23 risk-stratified skin lesions. A subset of samples was tested for the presence of HPV DNA. High-quality, non-human reads from each sample group were used for viral analysis using Microbiome Coverage Profiler.

RESULTS

None of the samples analysed had detectable expression of HPV RNA, while 64% of samples tested positive for HPV DNA. All samples were found to have expression of human endogenous retrovirus, and multiple samples showed expression of other viruses.

CONCLUSIONS

Viral and prophage gene expression can be monitored in cuSCC or normal skin biopsies, yet no sample in our study showed evidence of active HPV gene expression despite evidence of HPV genome presence. This suggests HPV transcription does not play a role in differentiating high-risk cuSCCs from low-risk cuSCCs or normal skin.

MicroRNA31 and MMP-1 contribute to the differentiated pathway of invasion -with enhanced epithelial-to-mesenchymal transition- in squamous cell carcinoma of the skin

Epithelial to mesenchymal transition (EMT) is an important mechanism of invasion in cutaneous squamous cell carcinomas (cSCCs) and has been found to be enhanced in tumors originated from actinic keratosis with transformation limited to the basal epithelial layer -differentiated pathway-, compared to cases with invasion subsequent to complete epidermal transformation -classical pathway-. Several microRNAs and proteins can contribute to EMT modulation in cSCCs. MicroRNA21 and microRNA31 are involved in posttranscriptional regulation of protein expression and could play a relevant role in EMT and cSCC progression. Throughout the EMT process upregulation of matrix metalloproteinases (MMPs) enhances invasiveness and MMP-1 and MMP-3 contribute to local invasion, angiogenesis and metastasis in cSCCs. Additionally, cSCC development is associated with PTEN loss and NF-κB, NOTCH-1 and p63 activation. The aim of this work is to identify differences in the expression of those molecules between both pathways of cSCCs development. Eight tissue microarrays from 80 consecutive cSCCs were analyzed using LNA-based miRNA in situ hybridization for miRNA21 and miRNA31 evaluation, and immunohistochemistry for MMP-1, MMP-3, PTEN, NOTCH-1, NF-κB, p63 and CD31. Significantly higher expression of miRNA31 (p < 0.0001) and MMP-1 (p = 0.0072) and angiogenesis (p = 0.0199) were found in the differentiated pathway, whereas PTEN loss (p = 0.0430) was more marked in the classical pathway. No significant differences were found for the other markers. Our findings support a contribution of miRNA31 and MMP-1 in the differentiated pathway, associated to EMT and increased microvascularization. The greater PTEN loss in the classical pathway indicate that its relevance in cSCC is not EMT-related.

Oh, the Mutations You'll Acquire! A Systematic Overview of Cutaneous Squamous Cell Carcinoma

Nearly two million cases of cutaneous squamous cell carcinoma (cSCC) are diagnosed every year in the United States alone. cSCC is notable for both its prevalence and its propensity for invasion and metastasis. For many patients, surgery is curative. However, patients experiencing immunosuppression or recurrent, advanced, and metastatic disease still face limited therapeutic options and significant mortality. cSCC forms after decades of sun exposure and possesses the highest known mutation rate of all cancers. This mutational burden complicates efforts to identify the primary factors driving cSCC initiation and progression, which in turn hinders the development of targeted therapeutics. In this review, we summarize the mutations and alterations that have been observed in patients’ cSCC tumors, affecting signaling pathways, transcriptional regulators, and the microenvironment. We also highlight novel therapeutic opportunities in development and clinical trials.

The Role of Circular RNAs in Keratinocyte Carcinomas

Keratinocyte carcinomas (KC) include basal cell carcinomas (BCC) and cutaneous squamous cell carcinomas (cSCC) and represents the most common cancer in Europe and North America. Both entities are characterized by a very high mutational burden, mainly UV signature mutations. Predominately mutated genes in BCC belong to the sonic hedgehog pathway, whereas, in cSCC, TP53, CDKN2A, NOTCH1/2 and others are most frequently mutated. In addition, the dysregulation of factors associated with epithelial to mesenchymal transition (EMT) was shown in invasive cSCC. The expression of factors associated with tumorigenesis can be controlled in several ways and include non-coding RNA molecules, such as micro RNAs (miRNA) long noncoding RNAs (lncRNA) and circular RNAs (circRNA). To update findings on circRNA in KC, we reviewed 13 papers published since 2016, identified in a PubMed search. In both BCC and cSCC, numerous circRNAs were identified that were differently expressed compared to healthy skin. Some of them were shown to target miRNAs that are also dysregulated in KC. Moreover, some studies confirmed the biological functions of individual circRNAs involved in cancer development. Thus, circRNAs may be used as biomarkers of disease and disease progression and represent potential targets of new therapeutic approaches for KC.

Characterization of circular RNA transcriptomes in psoriasis and atopic dermatitis reveals disease-specific expression profiles

Atopic dermatitis (AD) and psoriasis are two common chronic inflammatory skin diseases that are associated with various comorbidities. Circular RNA (circRNA) constitutes a major class of non-coding RNAs that have been implicated in many human diseases, although their potential involvement in inflammatory skin diseases remains elusive. Here, we compare and contrast the circRNA expression landscapes in paired lesional and non-lesional skin from psoriasis and AD patients relative to skin from unaffected individuals using high-depth RNA-seq data. CircRNAs and their cognate linear transcripts were quantified using the circRNA detection algorithm, CIRI2, and in situ hybridization and Sanger sequencing was used for validation purposes. We identified 39,286 circRNAs among all samples and found that psoriasis and AD lesional skin could be distinguished from non-lesional and healthy skin based on circRNA expression landscapes. In general, circRNAs were less abundant in lesional relative to non-lesional and healthy skin. Differential expression analyses revealed many significantly downregulated circRNAs, mainly in psoriasis lesional skin, and a strong correlation between psoriasis and AD-related circRNA expression changes was observed. Two individual circRNAs, ciRS-7 (also known as CDR1as) and circZRANB1, were specifically dysregulated in psoriasis and show promise as biomarkers for discriminating AD from psoriasis. In conclusion, the circRNA transcriptomes of psoriasis and AD share expression features, including a global downregulation relative to healthy skin, but this is most pronounced in psoriasis, and only psoriasis is characterized by several circRNAs being dysregulated independently of their cognate linear transcripts. Finally, specific circRNAs could potentially be used to distinguish AD from psoriasis.

Seven key hub genes identified by gene co-expression network in cutaneous squamous cell carcinoma

Background

Cutaneous squamous cell carcinoma (cSCC) often follows actinic keratosis (AK) and is the second most common skin cancer worldwide. To reduce metastasis risk, it is important to diagnose and treat cSCC early. This study aimed to identify hub genes associated with cSCC and AK.

Methods

This study used three datasets GSE45216, GSE98774, and GSE108008. We combined samples from the GSE45216 and GSE98774 datasets into the new dataset GSE45216–98774. We applied a weighted gene co-expression network analysis (WGCNA) to investigate key modules and hub genes associated with cSCC and AK. We considered the hub genes found in both the GSE45216–98774 and GSE108008 datasets as validated hub genes. We tested whether the expression of hub genes could predict patient survival outcomes in other cancers using TCGA pan-cancer data.

Results

We identified modules most relevant to cSCC and AK. Additionally, we identified and validated seven hub genes of cSCC: GATM, ARHGEF26, PTHLH, MMP1, POU2F3, MMP10 and GATA3. We did not find validated hub genes for AK. Each hub gene was significantly associated with the survival of various cancer types. Only GATA3 was significantly associated with melanoma survival.

Conclusions

We applied WGCNA to find seven hub genes that play important roles in cSCC tumorigenesis. These results provide new insights that help explain the pathogenesis of cSCC. These hub genes may become biomarkers or therapeutic targets for accurate diagnosis and treatment of cSCC in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08604-y.

Expression patterns and prognostic potential of circular RNAs in mantle cell lymphoma: a study of younger patients from the MCL2 and MCL3 clinical trials

Mantle cell lymphoma (MCL) is characterized by marked differences in outcome, emphasizing the need for strong prognostic biomarkers. Here, we explore expression patterns and prognostic relevance of circular RNAs (circRNAs), a group of endogenous non-coding RNA molecules, in MCL. We profiled the circRNA expression landscape using RNA-sequencing and explored the prognostic potential of 40 abundant circRNAs in samples from the Nordic MCL2 and MCL3 clinical trials, using NanoString nCounter Technology. We report a circRNA-based signature (circSCORE) developed in the training cohort MCL2 that is highly predictive of time to progression (TTP) and lymphoma-specific survival (LSS). The dismal outcome observed in the large proportion of patients assigned to the circSCORE high-risk group was confirmed in the independent validation cohort MCL3, both in terms of TTP (HR 3.0; P = 0.0004) and LSS (HR 3.6; P = 0.001). In Cox multiple regression analysis incorporating MIPI, Ki67 index, blastoid morphology and presence of TP53 mutations, circSCORE retained prognostic significance for TTP (HR 3.2; P = 0.01) and LSS (HR 4.6; P = 0.01). In conclusion, circRNAs are promising prognostic biomarkers in MCL and circSCORE improves identification of high-risk disease among younger patients treated with cytarabine-containing chemoimmunotherapy and autologous stem cell transplant.

Integrated tissue proteome and metabolome reveal key elements and regulatory pathways in cutaneous squamous cell carcinoma

Cutaneous squamous cell carcinoma (CSCC) is a widespread malignancy but has a very low long-term survival rate for patients at the metastatic stage. Therefore, it is urgent to identify prognostic biomarkers for CSCC and improve our understanding of disease progression. Here we took advantage of a data-independent acquisition (DIA)-based nano liquid chromatography equipped with an orbitrap mass spectrometry (nLC-MS/MS) and ultraperformance LC coupled to a time-of-flight tandem MS (UPLC-TOF-MS/MS) technique to analyze cancer and corresponding noncancerous tissues from 20 CSCC patients for integrated proteomic and metabolomic analyses. Overall, 6241 tissue proteins were detected, while 136 proteins were significantly expressed in CSCC tissues. Further functional analysis revealed that various biological processes were highly enriched and participated in the pathogenesis of CSCC, especially DNA damage responses. Moreover, 641 named metabolites in total were identified, among which 181 were significantly changed in CSCC tissues. A total of 101 pathways were significantly enriched including apoptosis, autophagy, PI3K-Akt and mTOR signaling pathways. Interestingly, two pathways, protein digestion & absorption and platelet activation were both enriched in proteomic and metabolomic studies involving 5 proteins and 11 metabolites. Accordingly, a four-metabolite panel consisting of arachidonate, glutamine, glutamic acid, and proline (all area under the curve (AUC) values more than 0.9) was developed with a high accuracy (0.971) to distinguish the 20 detected cancer tissues from their noncancerous tissues. Collectively, our work highlighted the key elements and regulatory pathways involved in the pathogenesis of CSCC. More importantly, the present study not only provided potential biomarkers for the early diagnosis of CSCC, but also expanded our knowledge of the physiopathology of the disease. SIGNIFICANCE: CSCC is the second most common human cancer but has few treatment options and few sensitive biomarkers for diagnosis. Here we comprehensively revealed its molecular characteristics by performing integrated tissue proteomic and metabolomic analyses. Significantly distinct profiles and certain enriched pathways including DNA damage responses were identified as associated with CSCC. Moreover, protein digestion & absorption and platelet activation were both enriched in the proteome and metabolome. These identified molecular changes probably play significant roles in CSCC development. Finally, we developed a four-metabolite panel to distinguish CSCC with high accuracy. Overall, our data not only provided potential diagnostic biomarkers, but also extended knowledge on the pathogenesis of CSCC.

Hsa_circ_0008234 facilitates proliferation of cutaneous squamous cell carcinoma through targeting miR-127-5p to regulate ADCY7

Cutaneous squamous cell carcinoma (cSCC) is the second most common skin malignant tumor with 25-50% of 5-year survival. There exist urgent needs for the identification of novel biomarkers for the diagnostic and therapeutic strategies of cSCC. The differentially expressed circRNAs in cSCC tissues and non-lesional skin tissues were obtained through analyzing the circular RNAs (circRNAs) microarray dataset GSE74758. The expression pattern of the indicated circRNAs in cSCC tissues was confirmed by qRT-PCR. FISH analysis was used to detect the location of hsa_circ_0008234 in cells. RIP experiment was used to detect the interaction between hsa_circ_0008234 and miR-127-5p. CCK-8 analysis and colony formation assay were used to detect the proliferation of cSCC cells. qRT-PCR and western blot were adopted to detect the expression of ACDY7. Three differential expressed circRNAs were obtained from the microarray data (GSE74758), and hsa_circ_0008234 was confirmed to be highly expressed in cSCC tissues by qRT-PCR. Hsa_circ_0008234 was mainly located in cytoplasm and stable in cSCC cells. RIP experiment revealed that hsa_circ_0008234 directly interacts with miR-127-5p in cSCC cells. Hsa_circ_0008234 increased the cell viability and colony formation of cSCC cells through acting as the sponge of miR-127-5p. MiR-127-5p inhibited the expression of ADCY7 in cSCC cells through binding the 3'UTR of ADCY7. Hsa_circ_0008234 was positively associated with ADCY7 expression in cSCC tissues. Hsa_circ_0008234 facilitates the proliferation of cSCC through targeting miR-127-5p to regulate ADCY7 expression and has the potential to be a novel therapeutic target for cSCC.

Cutaneous squamous cell carcinoma of head and neck region: A single center analysis of 1296 tumors with clinical characteristics, comorbidities, treatment, and sun-protection behavior

Cutaneous squamous cell carcinoma (SCC) is the second most common non-melanoma skin cancer among Caucasians. We analyzed our data on tumors been treated between 2006 and 2016 at the Skin Cancer Center of an academic teaching hospital. Thousand hundred and ninety-four patients with 1296 SCC have been identified including 804 males (67.3%) and 390 females (32.7%). The mean age among females was significantly higher with (83.9 ± 7.9) years compared with males ([79.3 ± 8.1] years; P < .001). The most common tumor localizations were outer ears (n = 227, 17.5%), scalp (216, 16.7%), and forehead (215, 16.6%). The majority of tumors was stage I. 31.9% of patients had previous SCC or Bowen's cancer, 29.1% had actinic keratoses. Major non-dermatological comorbidities were cardiovascular diseases and diabetes. SCC were treated by delayed Mohs surgery. First Mohs procedure resulted in R0 status in 83.9%. Most cases with R1 resection were located on nose and eyelids. For a subpopulation of patients (n = 105), a structured questionnaire study evaluated changes in sun protection behavior after skin cancer. After the confirmed diagnosis of a cutaneous SCC concerns about harm to the skin by intense sunlight was reported by 78.4% (76) of participants, with 92.0% among females versus 73.6% among males. Eighty percent of females and 76.8% of males wanted to care more about sun-protection in the future, spend less time 84.0% (females) and 72.5% (males) in the sun, and adapt their leisure activities for a better sun protection in 92.0% (females) and 81.9% (males). Secondary prevention in males needs improvement.

Identifying an lncRNA-Related ceRNA Network to Reveal Novel Targets for a Cutaneous Squamous Cell Carcinoma

Simple Summary

The exact functions and molecular mechanism of lncRNAs, acting as competitive endogenous RNAs in a cutaneous squamous cell carcinoma, remain unexplored. The present study was conducted to identify the differentially expressed lncRNAs and mRNAs and establish the lncRNA-related competing endogenous RNA networks associated with a cutaneous squamous cell carcinoma. A competing endogenous RNA network consisting of 137 miRNA-lncRNA and 221 miRNA-mRNA pairs was constructed. As for the functional analysis of the mRNAs in the network, a FoxO signaling pathway, an autophagy and cellular senescence were the top enrichment terms based on the Kyoto Encyclopedia of Genes and Genomes analysis. We identified five core mRNAs and built a core mRNA-associated competing endogenous RNA network. Finally, one lncRNA HLA-F-AS1 and three mRNAs named AGO4, E2F1 and CCND1 in the core mRNA-associated competing endogenous RNA network were validated with the same expression patterns. The core mRNAs and their associated lncRNAs may provide potential therapeutic targets for cutaneous squamous cell carcinomas.

Abstract

A cutaneous squamous cell carcinoma (cSCC) derived from keratinocytes is the second most common cause of non-melanoma skin cancer. The accumulation of the mutational burden of genes and cellular DNA damage caused by the risk factors (e.g., exposure to ultraviolet radiation) contribute to the aberrant proliferation of keratinocytes and the formation of a cSCC. A cSCC encompasses a spectrum of diseases that range from recursor actinic keratosis (AK) and squamous cell carcinoma (SCC) in situ (SCCIS) to invasive cSCCs and further metastatic SCCs. Emerging evidence has revealed that lncRNAs are involved in the biological process of a cSCC. According to the ceRNA regulatory theory, lncRNAs act as natural miRNA sponges and interact with miRNA response elements, thereby regulating the mRNA expression of their down-stream targets. This study was designed to search for the potential lncRNAs that may become potential therapeutic targets or biomarkers of a cSCC. Considering the spirit of the study to be adequately justified, we collected microarray-based datasets of 19 cSCC tissues and 12 normal skin samples from the GEO database (GSE42677 and GSE45164). After screening the differentially expressed genes via a limma package, we identified 24 differentially expressed lncRNAs (DElncRNAs) and 3221 differentially expressed mRNAs (DEmRNAs). The miRcode, miRTarBase, miRDB and TargetScan databases were used to predict miRNAs that could interact with DElncRNAs and DEmRNAs. A total of 137 miRNA-lncRNA and 221 miRNA-mRNA pairs were retained in the ceRNA network, consisting of 31 miRNAs, 11 DElncRNAs and 155 DEmRNAs. For the functional analysis, the top enriched biological process was enhancer sequence-specific DNA binding in Gene Ontology (GO) terms. The FoxO signaling pathway, autophagy and cellular senescence were the top enrichment terms based on a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The combination of a STRING tool and Cytoscape software (plug-in MCODE) identified five core mRNAs and built a core mRNA-associated ceRNA network. The expression for five identified core mRNAs and their related nine lncRNAs was validated using the external dataset GSE7553. Finally, one lncRNA HLA-F-AS1 and three mRNAs named AGO4, E2F1 and CCND1 were validated with the same expression patterns. We speculate that lncRNA HLA-F-AS1 may sponge miR-17-5p or miR-20b-5p to regulate the expression of CCND1 and E2F1 in the cSCC. The present study may provide potential diagnostic and therapeutic targets for cSCC patients.

Involvement of transcribed lncRNA uc.291 and SWI/SNF complex in cutaneous squamous cell carcinoma

While non-melanoma skin cancers (NMSCs) are the most common tumours in humans, only the sub-type cutaneous squamous cell carcinoma (cSCC), might become metastatic with high lethality. We have recently identified a regulatory pathway involving the lncRNA transcript uc.291 in controlling the expression of epidermal differentiation complex genes via the interaction with ACTL6A, a component of the chromatin remodelling complex SWI/SNF. Since transcribed ultra-conserved regions (T-UCRs) are expressed in normal tissues and are deregulated in tumorigenesis, here we hypothesize a potential role for dysregulation of this axis in cSCC, accounting for the de-differentiation process observed in aggressive poorly differentiated cutaneous carcinomas. We therefore analysed their expression patterns in human tumour biopsies at mRNA and protein levels. The results suggest that by altering chromatin accessibility of the epidermal differentiation complex genes, down-regulation of uc.291 and BRG1 expression contribute to the de-differentiation process seen in keratinocyte malignancy. This provides future direction for the identification of clinical biomarkers in cutaneous SCC. Analysis of publicly available data sets indicates that the above may also be a general feature for SCCs of different origins.

Tacrolimus inhibits oral carcinogenesis through cell cycle control

Tacrolimus (TAC, FK506) is a major calcineurin inhibitor and has been commonly used in treatments of patients with organ transplants and immune diseases. Moreover, tacrolimus is recommended by the treatment guidelines for oral potentially malignant disorders (OPMDs) such as oral lichen planus (OLP). However, whether tacrolimus increases the risk of cancer remains controversial. We observed that in a 4-Nitroquinoline N-oxide (4NQO)-induced oral carcinogenesis model, tacrolimus treatment was associated with a significantly lower ratio of cancer formation (52.94% vs. 90%) and a lower proportion of Ki67 and proliferation cell nuclear antigen (PCNA) -positive cells in lesion areas (P < 0.001). Liver, kidney, and lung functions of rats and the tumor immune microenvironment of the tongue were not affected. These observations suggest that tacrolimus blocked oral carcinogenesis through epithelial cell proliferation inhibition, independent of its immunosuppressive effects. As a processing factor, tacrolimus decreased tumor formation and cell proliferation in different stages of oral squamous cell carcinoma (OSCC) progression in vivo and in vitro. Furthermore, we investigated effects on the cell cycle and expression of related proteins. Tacrolimus induced G1/S phase arrest and significantly downregulated the expression of cyclinD1, cyclinE1, and c-Myc. These results suggest that tacrolimus induces G1/S phase arrest via inhibition of cyclinD1, cyclinE1, and c-Myc expression and retards oral cell carcinogenesis in vitro and in vivo. Thus, application of tacrolimus is a safe therapeutic strategy for treating OPMDs.

Epigenetic regulation of cellular functions in wound healing

Stringent spatiotemporal regulation of the wound healing process involving multiple cell types is associated with epigenetic mechanisms of gene regulation, such as DNA methylation, histone modification and chromatin remodelling, as well as non-coding RNAs. Here, we discuss the epigenetic changes that occur during wound healing and the rapidly expanding understanding of how these mechanisms affect healing resolution in both acute and chronic wound milieu. We provide a focussed overview of current research into epigenetic regulators that contribute to wound healing by specific cell type. We highlight the role of epigenetic regulators in the molecular pathophysiology of chronic wound conditions. The understanding of how epigenetic regulators can affect cellular functions during normal and impaired wound healing could lead to novel therapeutic approaches, and we outline questions that can provide guidance for future research on epigenetic-based interventions to promote healing. Dissecting the dynamic interplay between cellular subtypes involved in wound healing and epigenetic parameters during barrier repair will deepen our understanding of how to improve healing outcomes in patients affected by chronic non-healing wounds.

MiR-130a Acts as a Tumor Suppressor MicroRNA in Cutaneous Squamous Cell Carcinoma and Regulates the Activity of the BMP/SMAD Pathway by Suppressing ACVR1

Cutaneous squamous cell carcinoma (cSCC) is a malignant neoplasm of the skin resulting from the accumulation of somatic mutations due to solar radiation. cSCC is one of the fastest increasing malignancies, and it represents a particular problem among immunosuppressed individuals. MicroRNAs are short noncoding RNAs that regulate the expression of protein-coding genes at the post-transcriptional level. In this study, we identify miR-130a to be downregulated in cSCC compared to healthy skin and precancerous lesions (actinic keratosis). Moreoever, we show that its expression is regulated at the transcriptional level by HRAS and MAPK signaling pathway. We demonstrate that overexpession of miR-130a suppresses long-term capacity of growth, cell motility and invasion ability of human cSCC cell lines. We report that miR-130a suppresses the growth of cSCC xenografts in mice. Mechanistically, miR-130a directly targets ACVR1 (ALK2), and changes in miR-130a levels result in the decreased activity of the BMP/SMAD pathway through ACVR1. These data reveal a link between activated MAPK signaling and decreased expression of miR-130a, which acts as a tumor-suppressor microRNA in cSCC and contribute to a better understanding of the molecular processes during malignant transformation of epidermal keratinocytes.

CircRNA circ_0067772 aggravates the malignant progression of cutaneous squamous cell carcinoma by regulating miR-1238-3p/FOXG1 axis

BACKGROUND

Cutaneous squamous cell carcinoma (CSCC) is a severe malignancy derived from skin. Dysregulated circular RNAs (circRNAs) might play vital roles in tumor development.

OBJECTIVE

Here, we aimed to explore the function of a novel circRNA circ_0067772 in CSCC.

METHODS

Quantitative real-time PCR (qRT-PCR) or Western blot assay was performed to determine the expression of circ_0067772, microRNA (miR)-1238-3p and forkhead box protein G1 (FOXG1). Cell proliferation was assessed by Cell Counting Kit-8 (CCK-8) assay and colony formation assay. Transwell assay and wound healing assay were employed to examine cell metastasis. Flow cytometry was employed to monitor cell cycle and apoptosis. The target association between miR-1238-3p and circ_0067772 or FOXG1 was validated by dual-luciferase reporter assay. Moreover, role of circ_0067772 in vivo was investigated via xenograft model in nude mice.

RESULTS

Circ_0067772 and FOXG1 were upregulated, while miR-1238-3p was downregulated in CSCC tissues and cells. Circ_0067772 knockdown conferred inhibitory effects on cell proliferation, migration and invasion of CSCC cells. MiR-1238-3p served as a target of circ_0067772, whose silencing could reverse circ_0067772 knockdown-induced inhibitory impact on the malignant cellular behaviors. Circ_0067772 positively regulated FOXG1 expression by antagonizing miR-1238-3p. Additionally, miR-1238-3p could repress CSCC cell proliferation, migration and invasion by targeting FOXG1. Also, circ_0067772 knockdown hindered CSCC tumor growth in vivo.

CONCLUSION

Our study identified a novel oncogenic circRNA and the involvement of circ_0067772/miR-1238-3p/FOXG1 axis in CSCC development, providing a target for CSCC therapy.

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.

Profiling of circRNAs using an enzyme-free digital counting method

Circular RNAs (circRNAs) have recently gained substantial attention as potential biomarkers in several human diseases, most prominently in cancer. However, the detection and quantification of circRNA biomarkers pose specific challenges owing to their circular nature. In particular, the use of enzymes, such as reverse transcriptases, may lead to impaired quantification, poor interlaboratory reproducibility and even false-positive results. Therefore, the development of methods for accurate quantification of circRNA biomarkers is a high priority. An enzyme-free and digital quantification method, named NanoString nCounter, was recently adapted for circRNA detection. In this review, I describe the advantages of using this technology for circRNA quantification as well as methodological considerations, potential pitfalls and disadvantages.