Basal cell carcinomas (BCCs) rely on Hedgehog (HH) pathway growth signal amplification by the microtubule-based organelle, the primary cilium. Despite naive tumor responsiveness to Smoothened inhibitors (Smoⁱ), resistance in advanced tumors remains common. Although the resistant BCCs usually maintain HH pathway activation, squamous cell carcinomas with Ras/MAPK pathway activation also arise, and the molecular basis of tumor type and pathway selection are still obscure. Here, we identify the primary cilium as a critical determinant controlling tumor pathway switching. Strikingly, Smoothened inhibitor-resistant BCCs have an increased mutational load in ciliome genes, resulting in reduced primary cilia and HH pathway activation compared with naive or Gorlin syndrome patient BCCs. Gene set enrichment analysis of resistant BCCs with a low HH pathway signature showed increased Ras/MAPK pathway activation. Tissue analysis confirmed an inverse relationship between primary cilia presence and Ras/MAPK activation, and primary cilia removal in BCCs potentiated Ras/MAPK pathway activation. Moreover, activating Ras in HH-responsive cell lines conferred resistance to both canonical (vismodegib) and noncanonical (atypical protein kinase C and MRTF inhibitors) HH pathway inhibitors and conferred sensitivity to MAPK inhibitors. Our results provide insights into BCC treatment and identify the primary cilium as an important lineage gatekeeper, preventing HH-to-Ras/MAPK pathway switching.

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.

Rationale

Previous studies identified an interaction between HLA and oral peanut exposure. HLA-DQA1*01:02 had a protective role with the induction of Ara h 2 epitope-specific IgG4 associated with peanut consumption during the LEAP clinical trial for prevention of peanut allergy, while it was a risk allele for peanut allergy in the peanut avoidance group. We have now evaluated this gene-environment interaction in two subsequent peanut oral immunotherapy (OIT) trials - IMPACT and POISED - to better understand the potential for the HLA-DQA1*01:02 allele as an indicator of higher likelihood of desensitization, sustained unresponsiveness, and peanut allergy remission.

Methods

We determined HLA-DQA1*01:02 carrier status using genome sequencing from POISED (N=118, age: 7-55yr) and IMPACT (N=126, age: 12-<48mo). We tested for association with remission, sustained unresponsiveness (SU), and desensitization in the OIT groups, as well as peanut component specific IgG4 (psIgG4) using generalized linear models and adjusting for relevant covariates and ancestry.

Results

While not quite statistically significant, a higher proportion of HLA-DQA1*01:02 carriers receiving OIT in IMPACT were desensitized (93%) compared to non-carriers (78%); odds ratio (OR)=5.74 (p=0.06). In this sample we also observed that a higher proportion of carriers achieved remission (35%) compared to non-carriers (22%); OR=1.26 (p=0.80). In POISED, carriers more frequently attained continued desensitization (80% versus 61% among non-carriers; OR=1.28, p=0.86) and achieved SU (52% versus 31%; OR=2.32, p=0.19). psIgG4 associations with HLA-DQA1*01:02 in the OIT arm of IMPACT which included younger study subjects recapitulated patterns noted in LEAP, but no associations of note were observed in the older POISED study subjects.

Conclusions

Findings across three clinical trials show a pattern of a gene environment interaction between HLA and oral peanut exposure. Age, and prior sensitization contribute additional determinants of outcomes, consistent with a mechanism of restricted antigen recognition fundamental to driving protective immune responses to OIT.

Hedgehog pathway-dependent cancers can escape Smoothened (SMO) inhibition through mutations in genes encoding canonical hedgehog pathway components; however, around 50% of drug-resistant basal cell carcinomas (BCCs) lack additional variants of these genes. Here we use multidimensional genomics analysis of human and mouse drug-resistant BCCs to identify a noncanonical hedgehog activation pathway driven by the transcription factor serum response factor (SRF). Active SRF along with its coactivator megakaryoblastic leukemia 1 (MKL1) binds DNA near hedgehog target genes and forms a previously unknown protein complex with the hedgehog transcription factor glioma-associated oncogene family zinc finger-1 (GLI1), causing amplification of GLI1 transcriptional activity. We show that cytoskeletal activation through Rho and the formin family member Diaphanous (mDia) is required for SRF-MKL-driven GLI1 activation and for tumor cell viability. Remarkably, nuclear MKL1 staining served as a biomarker in tumors from mice and human subjects to predict tumor responsiveness to MKL inhibitors, highlighting the therapeutic potential of targeting this pathway. Thus, our study illuminates, for the first time, cytoskeletal-activation-driven transcription as a personalized therapeutic target for combatting drug-resistant malignancies.

Background

Eczema herpeticum (EH) is a rare complication of atopic dermatitis (AD) caused by disseminated herpes simplex virus (HSV) infection. The role of rare and/or deleterious genetic variants in disease etiology is largely unknown. This study aimed to identify genes that harbor damaging genetic variants associated with HSV infection in AD with a history of recurrent eczema herpeticum (ADEH+).

Methods

Whole genome sequencing (WGS) was performed on 49 recurrent ADEH+ (≥3 EH episodes), 491 AD without a history of eczema herpeticum (ADEH-) and 237 non-atopic control (NA) subjects. Variants were annotated, and a gene-based approach (SKAT-O) was used to identify genes harboring damaging genetic variants associated with ADEH+. Genes identified through WGS were studied for effects on HSV responses and keratinocyte differentiation.

Results

Eight genes were identified in the comparison of recurrent ADEH+to ADEH-and NA subjects: SIDT2, CLEC7A, GSTZ1, TPSG1, SP110, RBBP8NL, TRIM15, and FRMD3. Silencing SIDT2 and RBBP8NL in normal human primary keratinocytes (NHPKs) led to significantly increased HSV-1 replication. SIDT2-silenced NHPKs had decreased gene expression of IFNk and IL1b in response to HSV-1 infection. RBBP8NL-silenced NHPKs had decreased gene expression of IFNk, but increased IL1b. Additionally, silencing SIDT2 and RBBP8NL also inhibited gene expression of keratinocyte differentiation markers keratin 10 (KRT10) and loricrin (LOR).

Conclusion

SIDT2 and RBBP8NL participate in keratinocyte's response to HSV-1 infection. SIDT2 and RBBP8NL also regulate expression of keratinocyte differentiation genes of KRT10 and LOR.