Integrated multi-omics reveals cellular and molecular interactions governing the invasive niche of basal cell carcinoma

Contribution of Keratinocytes in Skin Cancer Initiation and Progression

Keratinocytes are major cellular components of the skin and are strongly involved in its homeostasis. Oncogenic events, starting mainly from excessive sun exposure, lead to the dysregulation of their proliferation and differentiation programs and promote the initiation and progression of non-melanoma skin cancers (NMSCs). Primary melanomas, which originate from melanocytes, initiate and develop in close interaction with keratinocytes, whose role in melanoma initiation, progression, and immune escape is currently being explored. Recent studies highlighted, in particular, unexpected modes of communication between melanocytic cells and keratinocytes, which may be of interest as sources of new biomarkers in melanomagenesis or potential therapeutic targets. This review aims at reporting the various contributions of keratinocytes in skin basal cell carcinoma (BCC), cutaneous squamous cell carcinoma (cSCC), and melanoma, with a greater focus on the latter in order to highlight some recent breakthrough findings. The readers are referred to recent reviews when contextual information is needed.

Distinctive Phenotypic and Microenvironmental Characteristics of Neuroendocrine Carcinoma and Adenocarcinoma Components in Gastric Mixed Adenoneuroendocrine Carcinoma

This study aimed to conduct an in-depth examination of gene expression and microenvironmental profiles of gastric neuroendocrine carcinoma (NEC) and mixed adeno-NEC (MANEC). Tissue microarrays from 55 patients with gastric MANEC (N = 32) or NEC (N = 23) were analyzed using digital spatial profiling (GeoMx DSP, NanoString Technologies). Representative regions of interest were selected from the adenocarcinoma (ADC) portion (ADC-MANEC) and the NEC portion (NEC-MANEC) of the MANEC cores, and pure NEC (pNEC) cores. All regions of interest were separated into epithelial components and stromal components using the masking procedure in the GeoMx platform, followed by transcriptome analysis. Comparison of gene expression between ADC-MANEC and NEC-MANEC/pNEC identified several differentially expressed genes in the epithelial (including PEG10, MAP1B, STMN3, and AKT3) and stromal (FN1, COL1A1, SPARC, and BGN) components. Gene set enrichment analysis revealed that pathways related to the E2F target and G2M checkpoint were more enriched in NEC-MANEC and pNEC than in ADC-MANEC. Deconvolution analysis showed that the microenvironmental profile varied according to histologic differentiation. In ADC-MANEC, intraepithelial infiltrating immune cells were relatively more numerous, whereas fibroblasts in the stroma were more abundant in NEC-MANEC and pNEC. This study confirmed the distinct expression profile of each histologic component of MANEC according to its tumor vs stromal compartment using the DSP platform. Although each component of MANEC shares the same genetic origin, distinctive phenotypes should not be overlooked when managing patients with MANEC. This study provides a useful validation data set for future studies.

Markers of Dermal Fibroblast Subpopulations for Viable Cell Isolation via Cell Sorting: A Comprehensive Review

Fibroblasts are among the most abundant cell types in the human body, playing crucial roles in numerous physiological processes, including the structural maintenance of the dermis, production of extracellular matrix components, and mediation of inflammatory responses. Despite their importance, fibroblasts remain one of the least characterized cell populations. The advent of single-cell analysis techniques, particularly single-cell RNA sequencing (scRNA-seq) and fluorescence-activated cell sorting (FACS), has enabled detailed investigations into fibroblast biology. In this study, we present an extensive analysis of fibroblast surface markers suitable for cell sorting and subsequent functional studies. We reviewed over three thousand research articles describing fibroblast populations and their markers, characterizing and comparing subtypes based on their surface markers, as well as their intra- and extracellular proteins. Our detailed analysis identified a variety of distinct fibroblast subpopulations, each with unique markers, characteristics dependent on their location, and the physiological or pathophysiological environment. These findings underscore the diversity of fibroblasts as a cellular population and could lead to the development of novel diagnostic and therapeutic tools.

High-throughput transcriptome profiling indicates ribosomal RNAs to be associated with resistance to immunotherapy in non-small cell lung cancer (NSCLC)

BACKGROUND

Despite the impressive outcomes with immune checkpoint inhibitor (ICI) in non-small cell lung cancer (NSCLC), only a minority of the patients show long-term benefits from ICI. In this study, we used retrospective cohorts of ICI treated patients with NSCLC to discover and validate spatially resolved protein markers associated with resistance to programmed cell death protein-1 (PD-1) axis inhibition.

METHODS

Pretreatment samples from 56 patients with NSCLC treated with ICI were collected and analyzed in a tissue microarray (TMA) format in including four different tumor regions per patient using the GeoMx platform for spatially informed transcriptomics. 34 patients had assessable tissue with tumor compartment in all 4 TMA spots, 22 with leukocyte compartment and 12 with CD68 compartment. The patients' tissue that was not assessable in fourfold redundancy in each compartment was designated as the validation cohort; cytokeratin (CK) (N=22), leukocytes CD45 (N=31), macrophages, CD68 (N=43). The human whole transcriptome, represented by~18,000 individual genes assessed by oligonucleotide-tagged in situ hybridization, was sequenced on the NovaSeq platform to quantify the RNAs present in each region of interest.

RESULTS

54,000 gene variables were generated per case, from them 25,740 were analyzed after removing targets with expression lower than a prespecified frequency. Cox proportional-hazards model analysis was performed for overall and progression-free survival (OS, PFS, respectively). After identifying genes significantly associated with limited survival benefit (HR>1)/progression per spot per patient, we used the intersection of them across the four TMA spots per patient. This resulted in a list of 12 genes in the tumor-cell compartment (RPL13A, GNL3, FAM83A, CYBA, ACSL4, SLC25A6, EPAS1, RPL5, APOL1, HSPD1, RPS4Y1, ADI1). RPL13A, GNL3 in tumor-cell compartment were also significantly associated with OS and PFS, respectively, in the validation cohort (CK: HR, 2.48; p=0.02 and HR, 5.33; p=0.04). In CD45 compartment, secreted frizzled-related protein 2, was associated with OS in the discovery cohort but not in the validation cohort. Similarly, in the CD68 compartment ARHGAP and PNN interacting serine and arginine rich protein were significantly associated with PFS and OS, respectively, in the majority but not all four spots per patient.

CONCLUSION

This work highlights RPL13A and GNL3 as potential indicative biomarkers of resistance to PD-1 axis blockade that might help to improve precision immunotherapy strategies for lung cancer.

Targeting PDGF signaling of cancer-associated fibroblasts blocks feedback activation of HIF-1α and tumor progression of clear cell ovarian cancer

Ovarian clear cell carcinoma (OCCC) is a gynecological cancer with a dismal prognosis; however, the mechanism underlying OCCC chemoresistance is not well understood. To explore the intracellular networks associated with the chemoresistance, we analyze surgical specimens by performing integrative analyses that combine single-cell analyses and spatial transcriptomics. We find that a chemoresistant OCCC subpopulation with elevated HIF activity localizes mainly in areas populated by cancer-associated fibroblasts (CAFs) with a myofibroblastic phenotype, which is corroborated by quantitative immunostaining. CAF-enhanced chemoresistance and HIF-1α induction are recapitulated in co-culture assays, which show that cancer-derived platelet-derived growth factor (PDGF) contributes to the chemoresistance and HIF-1α induction via PDGF receptor signaling in CAFs. Ripretinib is identified as an effective receptor tyrosine kinase inhibitor against CAF survival. In the co-culture system and xenograft tumors, ripretinib prevents CAF survival and suppresses OCCC proliferation in the presence of carboplatin, indicating that combination of conventional chemotherapy and CAF-targeted agents is effective against OCCC.

The Reign of Follistatin in Tumors and Their Microenvironment: Implications for Drug Resistance

Follistatin (FST) is a potent neutralizer of the transforming growth factor-β superfamily and is associated with normal cellular programs and various hallmarks of cancer, such as proliferation, migration, angiogenesis, and immune evasion. The aberrant expression of FST by solid tumors is a well-documented observation, yet how FST influences tumor progression and therapy response remains unclear. The recent surge in omics data has revealed new insights into the molecular foundation underpinning tumor heterogeneity and its microenvironment, offering novel precision medicine-based opportunities to combat cancer. In this review, we discuss these recent FST-centric studies, thereby offering an updated perspective on the protean role of FST isoforms in shaping the complex cellular ecosystem of tumors and in mediating drug resistance.

MHC-I upregulation safeguards neoplastic T cells in the skin against NK cell-mediated eradication in mycosis fungoides

Cancer-associated immune dysfunction is a major challenge for effective therapies. The emergence of antibodies targeting tumor cell-surface antigens led to advancements in the treatment of hematopoietic malignancies, particularly blood cancers. Yet their impact is constrained against tumors of hematopoietic origin manifesting in the skin. In this study, we employ a clonality-supervised deep learning methodology to dissect key pathological features implicated in mycosis fungoides, the most common cutaneous T-cell lymphoma. Our investigations unveil the prominence of the IL-32β-major histocompatibility complex (MHC)-I axis as a critical determinant in tumor T-cell immune evasion within the skin microenvironment. In patients' skin, we find MHC-I to detrimentally impact the functionality of natural killer (NK) cells, diminishing antibody-dependent cellular cytotoxicity and promoting resistance of tumor skin T-cells to cell-surface targeting therapies. Through murine experiments in female mice, we demonstrate that disruption of the MHC-I interaction with NK cell inhibitory Ly49 receptors restores NK cell anti-tumor activity and targeted T-cell lymphoma elimination in vivo. These findings underscore the significance of attenuating the MHC-I-dependent immunosuppressive networks within skin tumors. Overall, our study introduces a strategy to reinvigorate NK cell-mediated anti-tumor responses to overcome treatment resistance to existing cell-surface targeted therapies for skin lymphoma.

Inhibition of anti-tumor immunity by melanoma cell-derived Activin-A depends on STING

The transforming growth factor-β (TGF-β) family member activin A (hereafter Activin-A) is overexpressed in many cancer types, often correlating with cancer-associated cachexia and poor prognosis. Activin-A secretion by melanoma cells indirectly impedes CD8+ T cell-mediated anti-tumor immunity and promotes resistance to immunotherapies, even though Activin-A can be proinflammatory in other contexts. To identify underlying mechanisms, we here analyzed the effect of Activin-A on syngeneic grafts of Braf mutant YUMM3.3 mouse melanoma cells and on their microenvironment using single-cell RNA sequencing. We found that the Activin-A-induced immune evasion was accompanied by a proinflammatory interferon signature across multiple cell types, and that the associated increase in tumor growth depended at least in part on pernicious STING activity within the melanoma cells. Besides corroborating a role for proinflammatory signals in facilitating immune evasion, our results suggest that STING holds considerable potential as a therapeutic target to mitigate tumor-promoting Activin-A signaling at least in melanoma.

Multiscale spatial mapping of cell populations across anatomical sites in healthy human skin and basal cell carcinoma

Our understanding of how human skin cells differ according to anatomical site and tumour formation is limited. To address this, we have created a multiscale spatial atlas of healthy skin and basal cell carcinoma (BCC), incorporating in vivo optical coherence tomography, single-cell RNA sequencing, spatial global transcriptional profiling, and in situ sequencing. Computational spatial deconvolution and projection revealed the localisation of distinct cell populations to specific tissue contexts. Although cell populations were conserved between healthy anatomical sites and in BCC, mesenchymal cell populations including fibroblasts and pericytes retained signatures of developmental origin. Spatial profiling and in silico lineage tracing support a hair follicle origin for BCC and demonstrate that cancer-associated fibroblasts are an expansion of a POSTN+ subpopulation associated with hair follicles in healthy skin. RGS5+ pericytes are also expanded in BCC suggesting a role in vascular remodelling. We propose that the identity of mesenchymal cell populations is regulated by signals emanating from adjacent structures and that these signals are repurposed to promote the expansion of skin cancer stroma. The resource we have created is publicly available in an interactive format for the research community.

Spatially Resolved Transcriptomics Technology Facilitates Cancer Research

Single cell RNA sequencing (scRNA-seq) provides a great convenience for studying tumor occurrence and development for its ability to study gene expression at the individual cell level. However, patient-derived tumor tissues are composed of multiple types of cells including tumor cells and adjacent non-malignant cells such as stromal cells and immune cells. The spatial locations of various cells in situ tissues plays a pivotal role in the occurrence and development of tumors, which cannot be elucidated by scRNA-seq alone. Spatially resolved transcriptomics (SRT) technology emerges timely to explore the unrecognized relationship between the spatial background of a particular cell and its functions, and is increasingly used in cancer research. This review provides a systematic overview of the SRT technologies that are developed, in particular the more widely used cutting-edge SRT technologies based on next-generation sequencing (NGS). In addition, the main achievements by SRT technologies in precisely unveiling the underappreciated spatial locations on gene expression and cell function with unprecedented high-resolution in cancer research are emphasized, with the aim of developing more effective clinical therapeutics oriented to a deeper understanding of the interaction between tumor cells and surrounding non-malignant cells.

epiAneufinder identifies copy number alterations from single-cell ATAC-seq data

Single-cell open chromatin profiling via scATAC-seq has become a mainstream measurement of open chromatin in single-cells. Here we present epiAneufinder, an algorithm that exploits the read count information from scATAC-seq data to extract genome-wide copy number alterations (CNAs) for individual cells, allowing the study of CNA heterogeneity present in a sample at the single-cell level. Using different cancer scATAC-seq datasets, we show that epiAneufinder can identify intratumor clonal heterogeneity in populations of single cells based on their CNA profiles. We demonstrate that these profiles are concordant with the ones inferred from single-cell whole genome sequencing data for the same samples. EpiAneufinder allows the inference of single-cell CNA information from scATAC-seq data, without the need of additional experiments, unlocking a layer of genomic variation which is otherwise unexplored.

Molecular Classifiers in Skin Cancers: Challenges and Promises

Skin cancers are common and heterogenous malignancies affecting up to two in three Australians before age 70. Despite recent developments in diagnosis and therapeutic strategies, the mortality rate and costs associated with managing patients with skin cancers remain high. The lack of well-defined clinical and histopathological features makes their diagnosis and classification difficult in some cases and the prognostication difficult in most skin cancers. Recent advancements in large-scale "omics" studies, including genomics, transcriptomics, proteomics, metabolomics and imaging-omics, have provided invaluable information about the molecular and visual landscape of skin cancers. On many occasions, it has refined tumor classification and has improved prognostication and therapeutic stratification, leading to improved patient outcomes. Therefore, this paper reviews the recent advancements in omics approaches and appraises their limitations and potential for better classification and stratification of skin cancers.

Functionally distinct cancer-associated fibroblast subpopulations establish a tumor promoting environment in squamous cell carcinoma

Cutaneous squamous cell carcinoma (cSCC) is a serious public health problem due to its high incidence and metastatic potential. It may progress from actinic keratosis (AK), a precancerous lesion, or the in situ carcinoma, Bowen's disease (BD). During this progression, malignant keratinocytes activate dermal fibroblasts into tumor promoting cancer-associated fibroblasts (CAFs), whose origin and emergence remain largely unknown. Here, we generate and analyze >115,000 single-cell transcriptomes from healthy skin, BD and cSCC of male donors. Our results reveal immunoregulatory and matrix-remodeling CAF subtypes that may derive from pro-inflammatory and mesenchymal fibroblasts, respectively. These CAF subtypes are largely absent in AK and interact with different cell types to establish a pro-tumorigenic microenvironment. These findings are cSCC-specific and could not be recapitulated in basal cell carcinomas. Our study provides important insights into the potential origin and functionalities of dermal CAFs that will be highly beneficial for the specific targeting of the cSCC microenvironment.

Reprogramming of cancer-associated fibroblasts combined with immune checkpoint inhibitors: A potential therapeutic strategy for cancers

Activated fibroblasts, namely cancer-associated fibroblasts (CAFs), are highly heterogeneous in phenotypes, functions, and origins. CAFs originated from varieties of cell types, including local resident fibroblasts, epithelial cells, mesenchymal stromal cells, or others. These cells participate in tumor angiogenesis, mechanics, drug access, and immune suppression, with the latter being particularly important. It was difficult to distinguish CAFs by subsets due to their complex origins until the use of scRNA-seq. Reprogramming CAFs with TGFβ-RI inhibitor, a CXCR4 blocker, or other methods increases T cells activation and infiltration, together with a decrease in CAFs recruitment, thus improving the prognosis. As depletion of CAFs can't bring clinical benefit, the combination of reprogramming CAFs and immune checkpoint inhibitors (ICIs) come into consideration. It has shown better outcomes compared with monotherapy respectively in basic/preclinical researches, and needs more data on clinical trials. Combination therapy may be a promising and expecting method for treatment of cancer.

Activin A-Mediated Polarization of Cancer-Associated Fibroblasts and Macrophages Confers Resistance to Checkpoint Immunotherapy in Skin Cancer

PURPOSE

Cemiplimab is approved for the treatment of locally advanced basal cell carcinomas (BCC), although with mitigated results. We sought to interrogate the cellular and molecular transcriptional reprogramming underlying BCC resistance to immunotherapy.

EXPERIMENTAL DESIGN

Here, we combined spatial and single-cell transcriptomics to deconvolute the spatial heterogeneity of the tumor microenvironment in regard with response to immunotherapy, in a cohort of both naïve and resistant BCCs.

RESULTS

We identified subsets of intermingled cancer-associated fibroblasts (CAF) and macrophages contributing the most to CD8 T-cell exclusion and immunosuppression. Within this spatially resolved peritumoral immunosuppressive niche, CAFs and adjacent macrophages were found to display Activin A-mediated transcriptional reprogramming towards extracellular matrix remodeling, suggesting active participation to CD8 T-cell exclusion. In independent datasets of human skin cancers, Activin A-conditioned CAFs and macrophages were associated with resistance to immune checkpoint inhibitors (ICI).

CONCLUSIONS

Altogether, our data identify the cellular and molecular plasticity of tumor microenvironment (TME) and the pivotal role of Activin A in polarizing the TME towards immune suppression and ICI resistance.

The Use of Single-Cell RNA-Sequencing and Spatial Transcriptomics in Understanding the Pathogenesis and Treatment of Skin Diseases

The development of multiomic profiling tools has rapidly expanded in recent years, along with their use in profiling skin tissues in various contexts, including dermatologic diseases. Among these tools, single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST) have emerged as widely adopted and powerful assays for elucidating key cellular components and their spatial arrangement within skin disease. In this paper, we review the recent biological insights gained from the use of scRNA-seq and ST and the advantages of combining both for profiling skin diseases, including aberrant wound healing, inflammatory skin diseases, and cancer. We discuss the role of scRNA-seq and ST in improving skin disease treatments and moving toward the goal of achieving precision medicine in dermatology, whereby patients can be optimally matched to treatments that maximize therapeutic response.

Non-Canonical Activin A Signaling Stimulates Context-Dependent and Cellular-Specific Outcomes in CRC to Promote Tumor Cell Migration and Immune Tolerance

We have shown that activin A (activin), a TGF-β superfamily member, has pro-metastatic effects in colorectal cancer (CRC). In lung cancer, activin activates pro-metastatic pathways to enhance tumor cell survival and migration while augmenting CD4+ to CD8+ communications to promote cytotoxicity. Here, we hypothesized that activin exerts cell-specific effects in the tumor microenvironment (TME) of CRC to promote anti-tumoral activity of immune cells and the pro-metastatic behavior of tumor cells in a cell-specific and context-dependent manner. We generated an Smad4 epithelial cell specific knockout (Smad4-/-) which was crossed with TS4-Cre mice to identify SMAD-specific changes in CRC. We also performed IHC and digital spatial profiling (DSP) of tissue microarrays (TMAs) obtained from 1055 stage II and III CRC patients in the QUASAR 2 clinical trial. We transfected the CRC cells to reduce their activin production and injected them into mice with intermittent tumor measurements to determine how cancer-derived activin alters tumor growth in vivo. In vivo, Smad4-/- mice displayed elevated colonic activin and pAKT expression and increased mortality. IHC analysis of the TMA samples revealed increased activin was required for TGF-β-associated improved outcomes in CRC. DSP analysis identified that activin co-localization in the stroma was coupled with increases in T-cell exhaustion markers, activation markers of antigen presenting cells (APCs), and effectors of the PI3K/AKT pathway. Activin-stimulated PI3K-dependent CRC transwell migration, and the in vivo loss of activin lead to smaller CRC tumors. Taken together, activin is a targetable, highly context-dependent molecule with effects on CRC growth, migration, and TME immune plasticity.

Spatial transcriptomics in human skin research

Spatial transcriptomics is a revolutionary technique that enables researchers to characterise tissue architecture and localisation of gene expression. A plethora of technologies that map gene expression are currently being developed, aiming to facilitate spatially resolved, high-dimensional assessment of gene transcription in the context of human skin research. Knowing which gene is expressed by which cell and in which location within skin, facilitates understanding of skin function and dysfunction in both health and disease. In this review, we summarise the available spatial transcriptomic methods and we describe their application to a broad spectrum of dermatological diseases.

Recurrence or de novo? Intradermal Basal Cell Carcinoma of the Scrotum: A Report of Two Cases

Basal cell carcinoma (BCC) is the most common cutaneous malignancy, usually occurring in sun-exposed areas. Although BCC in the scrotal region is uncommon, it carries a higher risk of metastasis than BCC at other sites. Here, we report two cases of BCC that developed in the scrotal region: Case 1 presented as a superficial nodule and Case 2 as a subcutaneous nodule. Histopathologically, both tumors lacked continuity with the surface epidermis and formed an intradermal nodule. In Case 1, BCC occurred for the first time and presumably developed de novo. Case 2 underwent excision of a scrotal BCC 5 years previously, and the histopathological diagnosis at that time was nodular BCC. However, when the original specimen was re-examined, it was determined that, although the tumor had been completely resected, part of the lesion had moved away from the nodular area to represent a micronodular phenotype, an aggressive BCC subtype. We hypothesized that partial evolution from a nodular to a micronodular phenotype may have contributed to the recurrence of BCC in Case 2.

SRTsim: spatial pattern preserving simulations for spatially resolved transcriptomics

Spatially resolved transcriptomics (SRT)-specific computational methods are often developed, tested, validated, and evaluated in silico using simulated data. Unfortunately, existing simulated SRT data are often poorly documented, hard to reproduce, or unrealistic. Single-cell simulators are not directly applicable for SRT simulation as they cannot incorporate spatial information. We present SRTsim, an SRT-specific simulator for scalable, reproducible, and realistic SRT simulations. SRTsim not only maintains various expression characteristics of SRT data but also preserves spatial patterns. We illustrate the benefits of SRTsim in benchmarking methods for spatial clustering, spatial expression pattern detection, and cell-cell communication identification.

Bodywide ecological interventions on cancer

Historically, cancer research and therapy have focused on malignant cells and their tumor microenvironment. However, the vascular, lymphatic and nervous systems establish long-range communication between the tumor and the host. This communication is mediated by metabolites generated by the host or the gut microbiota, as well by systemic neuroendocrine, pro-inflammatory and immune circuitries-all of which dictate the trajectory of malignant disease through molecularly defined biological mechanisms. Moreover, aging, co-morbidities and co-medications have a major impact on the development, progression and therapeutic response of patients with cancer. In this Perspective, we advocate for a whole-body 'ecological' exploration of malignant disease. We surmise that accumulating knowledge on the intricate relationship between the host and the tumor will shape rational strategies for systemic, bodywide interventions that will eventually improve tumor control, as well as quality of life, in patients with cancer.