Skin basal cell carcinomas assemble a pro-tumorigenic spatially organized and self-propagating Trem2+ myeloid niche.

Cancer immunotherapies have revolutionized treatment but have shown limited success as single-agent therapies highlighting the need to understand the origin, assembly, and dynamics of heterogeneous tumor immune niches. Here, we use single-cell and imaging-based spatial analysis to elucidate three microenvironmental neighborhoods surrounding the heterogeneous basal cell carcinoma tumor epithelia. Within the highly proliferative neighborhood, we find that TREM2+ skin cancer-associated macrophages (SCAMs) support the proliferation of a distinct tumor epithelial population through an immunosuppression-independent manner via oncostatin-M/JAK-STAT3 signaling. SCAMs represent a unique tumor-specific TREM2+ population defined by VCAM1 surface expression that is not found in normal homeostatic skin or during wound healing. Furthermore, SCAMs actively proliferate and self-propagate through multiple serial tumor passages, indicating long-term potential. The tumor rapidly drives SCAM differentiation, with intratumoral injections sufficient to instruct naive bone marrow-derived monocytes to polarize within days. This work provides mechanistic insights into direct tumor-immune niche dynamics independent of immunosuppression, providing the basis for potential combination tumor therapies.

LY6D marks pre-existing resistant basosquamous tumor subpopulations.

Improved response to canonical therapies requires a mechanistic understanding of dynamic tumor heterogeneity by identifying discrete cellular populations with enhanced cellular plasticity. We have previously demonstrated distinct resistance mechanisms in skin basal cell carcinomas, but a comprehensive understanding of the cellular states and markers associated with these populations remains poorly understood. Here we identify a pre-existing resistant cellular population in naive basal cell carcinoma tumors marked by the surface marker LY6D. LY6D+ tumor cells are spatially localized and possess basal cell carcinoma and squamous cell carcinoma-like features. Using computational tools, organoids, and spatial tools, we show that LY6D+ basosquamous cells represent a persister population lying on a central node along the skin lineage-associated spectrum of epithelial states with local environmental and applied therapies determining the kinetics of accumulation. Surprisingly, LY6D+ basosquamous populations exist in many epithelial tumors, such as pancreatic adenocarcinomas, which have poor outcomes. Overall, our results identify the resistant LY6D+ basosquamous population as an important clinical target and suggest strategies for future therapeutic approaches to target them.

Gibbin mesodermal regulation patterns epithelial development.

Proper ectodermal patterning during human development requires previously identified transcription factors such as GATA3 and p63, as well as positional signalling from regional mesoderm1-6. However, the mechanism by which ectoderm and mesoderm factors act to stably pattern gene expression and lineage commitment remains unclear. Here we identify the protein Gibbin, encoded by the Xia-Gibbs AT-hook DNA-binding-motif-containing 1 (AHDC1) disease gene7-9, as a key regulator of early epithelial morphogenesis. We find that enhancer- or promoter-bound Gibbin interacts with dozens of sequence-specific zinc-finger transcription factors and methyl-CpG-binding proteins to regulate the expression of mesoderm genes. The loss of Gibbin causes an increase in DNA methylation at GATA3-dependent mesodermal genes, resulting in a loss of signalling between developing dermal and epidermal cell types. Notably, Gibbin-mutant human embryonic stem-cell-derived skin organoids lack dermal maturation, resulting in p63-expressing basal cells that possess defective keratinocyte stratification. In vivo chimeric CRISPR mouse mutants reveal a spectrum of Gibbin-dependent developmental patterning defects affecting craniofacial structure, abdominal wall closure and epidermal stratification that mirror patient phenotypes. Our results indicate that the patterning phenotypes seen in Xia-Gibbs and related syndromes derive from abnormal mesoderm maturation as a result of gene-specific DNA methylation decisions.

c-FOS drives reversible basal to squamous cell carcinoma transition.

While squamous transdifferentiation within subpopulations of adenocarcinomas represents an important drug resistance problem, its underlying mechanism remains poorly understood. Here, using surface markers of resistant basal cell carcinomas (BCCs) and patient single-cell and bulk transcriptomic data, we uncover the dynamic roadmap of basal to squamous cell carcinoma transition (BST). Experimentally induced BST identifies activator protein 1 (AP-1) family members in regulating tumor plasticity, and we show that c-FOS plays a central role in BST by regulating the accessibility of distinct AP-1 regulatory elements. Remarkably, despite prominent changes in cell morphology and BST marker expression, we show using inducible model systems that c-FOS-mediated BST demonstrates reversibility. Blocking EGFR pathway activation after c-FOS induction partially reverts BST in vitro and prevents BST features in both mouse models and human tumors. Thus, by identifying the molecular basis of BST, our work reveals a therapeutic opportunity targeting plasticity as a mechanism of tumor resistance.

AP-1 and TGFß cooperativity drives non-canonical Hedgehog signaling in resistant basal cell carcinoma.

Tumor heterogeneity and lack of knowledge about resistant cell states remain a barrier to targeted cancer therapies. Basal cell carcinomas (BCCs) depend on Hedgehog (Hh)/Gli signaling, but can develop mechanisms of Smoothened (SMO) inhibitor resistance. We previously identified a nuclear myocardin-related transcription factor (nMRTF) resistance pathway that amplifies noncanonical Gli1 activity, but characteristics and drivers of the nMRTF cell state remain unknown. Here, we use single cell RNA-sequencing of patient tumors to identify three prognostic surface markers (LYPD3, TACSTD2, and LY6D) which correlate with nMRTF and resistance to SMO inhibitors. The nMRTF cell state resembles transit-amplifying cells of the hair follicle matrix, with AP-1 and TGFß cooperativity driving nMRTF activation. JNK/AP-1 signaling commissions chromatin accessibility and Smad3 DNA binding leading to a transcriptional program of RhoGEFs that facilitate nMRTF activity. Importantly, small molecule AP-1 inhibitors selectively target LYPD3+/TACSTD2+/LY6D+ nMRTF human BCCs ex vivo, opening an avenue for improving combinatorial therapies.