Next Article in Journal
Prevalence of Germline BRCA1/2 Variants in Ashkenazi and Non-Ashkenazi Prostate Cancer Populations: A Systematic Review and Meta-Analysis
Previous Article in Journal
Protoporphyrin IX (PpIX) Fluorescence during Meningioma Surgery: Correlations with Histological Findings and Expression of Heme Pathway Molecules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Basal Cell Carcinoma (BCC) Tumour Microenvironment to Other Solid Malignancies

by
Eliana-Ruobing Zhang
1,†,
Sarah Ghezelbash
1,2,3,†,
Pingxing Xie
1,4,
Misha Fotovati
1,
Ivan V. Litvinov
1,2,4 and
Philippe Lefrançois
1,2,3,4,5,*
1
Faculty of Medicine, McGill University, Montreal, QC H3T 1E2, Canada
2
Division of Experimental Medicine, Department of Medicine, McGill University, Montreal, QC H3T 1E2, Canada
3
Lady Davis Institute for Medical Research, Montreal, QC H3T 1E2, Canada
4
Division of Dermatology, Department of Medicine, McGill University, Montreal, QC H3T 1E2, Canada
5
Division of Dermatology, Department of Medicine, Jewish General Hospital, Montreal, QC H3T 1E2, Canada
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2023, 15(1), 305; https://doi.org/10.3390/cancers15010305
Submission received: 30 November 2022 / Revised: 28 December 2022 / Accepted: 30 December 2022 / Published: 2 January 2023
(This article belongs to the Section Tumor Microenvironment)

Abstract

:

Simple Summary

Basal cell carcinoma is the most common human cancer. Most BCCs are low-risk and are easily treated; however, 1–2% are aggressive and highly destructive to the surrounding skin, called advanced BCC. It was discovered that this subtype has a different immune profile than routine BCC, and they contain a special type of stem cell population that helps them grow and spread. There is currently no reliable laboratory model for advanced BCC, making it hard to further study it and find new treatments. For these reasons, this project was conducted using genomic data from 11,000 tumours coming from 33 non-BCC cancer types. Using computational biology, we have compared the immune cell makeup of the tumour microenvironments to determine the top three most similar cancers, which we will call BCC “relatives”. We will examine how these “relatives” develop and grow, as well as current existing treatments and their response to such treatments.

Abstract

Basal cell carcinoma (BCC) is the most common form of skin cancer, contributing to nearly a third of new cancer cases in Western countries. Most BCCs are considered low risk “routine” lesions that can either be excised through surgery or treated with chemotherapeutic agents. However, around 1–2% of BCC cases are locally aggressive, present a high risk of metastasis, and often develop chemoresistance, termed advanced BCC. There currently exists no animal model or cell line that can recapitulate advanced BCC, let alone intermediate-risk and high-risk early BCC. We previously found that aggressive BCC tumours presented a Th2 cytokine inflammation profile, mesenchymal stem cell properties, and macrophage-induced tumoral inflammation. In this study, we aimed to identify potential BCC “relatives” among solid-organ malignancies who present similar immune cell proportions in their microenvironment compositions. Using immune cell type deconvolution by CIBERSORTx, and cell type enrichment by xCell, we determined three cancers with the most similar tumour microenvironments as compared to BCC. Specifically, chromophobe renal cell carcinoma, sarcoma, and skin cutaneous melanoma presented significance in multiple cell types, namely in CD4+ T lymphocytes, gammadelta T lymphocytes, and NK cell populations. Consequently, further literature analysis was conducted to understand similarities between BCC and its “relatives”, as well as investigating novel treatment targets. By identifying cancers most like BCC, we hope to propose prospective druggable pathways, as well as to gain insight on developing a reliable animal or cell line model to represent advanced BCC.

1. Introduction

Basal cell carcinoma (BCC) is the most common type of skin cancer and accounts for roughly 80% of all skin cancers, with increased risk in older, fair-skinned individuals [1,2]. In particular, the risk of Caucasians being diagnosed with BCC reaches up to 40% [3]. Major risk factors include age, UV exposure (particularly UVB), and immunosuppression. Current understanding of BCC progression and tumorigenesis revolves around the Sonic hedgehog (Shh) pathway, an important regulator of cell differentiation and tumorigenesis [4,5]. Its constitutive activation, specifically through PTCH1 inactivation [6], has been shown to be crucial for both sporadic and familial forms of BCC [7]. Small molecule inhibitors such as vismodegib and sonidegib, which target SMO (Smoothened) in the Shh pathway, have been FDA-approved for treatment of advanced BCC [8]; however, they are limited by their severe side effects and failure of complete patient response [9].
Although most BCCs are considered “routine”, or low-risk, and are simple to treat, around 1–2% of BCC cases are locally aggressive and present a high risk of metastasis [10]. These advanced cases are often resistant to the targeted therapies mentioned above. Currently, there exists no reliable lab models in either animal or cell lines to study advanced BCC, owing to their lack of important features in experimental models [11], or their difficulty of manipulation [12], respectively. We recently determined that, in comparison to non-advanced BCC, aggressive BCC tumours presented an inflammatory cytokine profile that shifted to a Th2 cytokine inflammation profile, which is thought to be more tumour permissive [13] and may include pathways such as Toll-like receptor, PDGFR, and extracellular matrix remodeling [14]. Furthermore, advanced BCC had a higher enrichment for mesenchymal stem cells [13].
Given the lack of reliable experimental models for BCC, especially more aggressive, higher risk, or advanced disease, we aimed to use our findings related to tumour microenvironment to identify other cancers/cancer subtypes most resembling BCC. Hence, this manuscript aims to identify non-BCC cancers that present comparable tumour microenvironment features to that of BCC (“relatives”), in hopes to contribute to both the development of a reliable model for advanced BCC and novel potential therapeutic avenues.

2. Materials and Methods

2.1. Data Acquisition and RNA-Seq Processing

Whole-genome RNA-sequencing data for 75 BCC samples was obtained and processed as previously described [13,14]. Non-BCC whole-genome RNA-sequencing data was extracted from all 33 cancer types in The Cancer Genome Atlas consortium, and processed as similarly described [15].

2.2. Cell-Type Enumeration Using RNA Deconvolution

To identify the closest BCC “relatives”, RNA deconvolution in CIBERSORTx was performed, using the standard LM22 leukocyte signature matrix obtained from 22 pure immune cell lines [16] and 100 permutations, to estimate the relative number of each cell subtype. The following analyses were performed: all B cells, all CD4+ T cells, CD8+ T cells, Treg cells, all NK T cells, T γδ (gamma delta) cells, total lymphocytes, and total macrophages, as previously performed [17]. Primary analyses were conducted using CIBERSORTx’s standard relative score. xCell analysis was also performed to determine cell type abundance scores using xCell’s standard 64 cell type signatures [18]. To consider relatedness of cancers, we tested, for a specific feature, whether the distributions of scores for that feature were similar between BCC and the non-BCC malignancy. We considered non-rejection of the null hypothesis using asymptotic two-sample Kolmogorov–Smirnov test to compare and validate significant distributions. The top three (if available) cancers with the highest p-values over 0.05 were considered as a potential BCC “relative”. CIBERSORTx and xCell scores were separately Bonferroni-corrected for multiple hypothesis testing. Figures were generated using Graphpad Prism (violin plots) and Rtsne R package (tSNE plot).

3. Results

We compared the microenvironment of 75 BCC samples to over 11,000 tumours across 33 cancer subtypes from the TGCA using xCell (Figure 1) and CIBERSORTx (Figure 2), with immune cell types selected according to previous findings [13]. Using these guiding cell type populations, we identified the top three most recurring cancers. Chromophobe renal cell carcinoma (KICH) presented similar distributions in its scores for CD4+ cells (p = 0.17), macrophages (p = 0.4904), Th2 cells (p = 0.151), and Tgammadelta cells (p = 0.5525). Sarcoma (SARC) showed significant p-values in its total lymphocyte (p = 0.2290), NK cells (p = 0.1979), Tgammadelta cells (p = 0.07291), and CD4+ (p = 0.6344) cell populations. Bladder urothelial carcinoma (BLCA) demonstrated significant distributions in its NK (p = 0.1718), CD4+ (p = 0.2715), and macrophage (p = 0.3333) cell scores. Cholangiocarcinoma (CHOL) also showed high significance in its Tgammadelta (p = 0.6099), NK (p = 0.05944), and CD4+ (p = 0.6959) cell populations compared to BCC. Lymphoid neoplasm diffuse large B-cell lymphoma (DLBC) presented moderately significant p-values in scores for macrophages (p = 0.08976), Tgammadelta (p = 0.08334), and NK cells (p = 0.1221). Glioblastoma multiforme (GBM) also demonstrated high significance in its CD8+ (p = 0.2741), B cells (p = 0.2374), and Tgammadelta (p = 0.244) cell populations to BCC. Kidney renal clear cell carcinoma (KIRC) showed similar distributions for its NK (p = 0.2311), Th2 (p = 0.1368), and Tgammadelta (p = 0.2002) cell populations. Skin cutaneous melanoma (SKCM) had significantly similar cell populations of MSCs (p = 0.3872), NK cells (p = 0.05057), and CD4+ cells (p = 0.5715) to BCC. Finally, thyroid carcinoma (THCA) showed significance in its Tgammadelta (p = 0.7752), CD4+ (p = 0.4649), and NK (p = 0.0512) cell populations compared to BCC.
Given that seven cancers had three similarly distributed features compared to BCC (BLCA, CHOL, DLBC, GBM, KIRC, SKCM, THCA), we favoured one of the three potential close relatives over the others. SKCM and BCC both originate in the skin and are UV-driven. A summary of the top three most similar cancers by cell type population is presented in Table 1, with an expanded version in Supplementary Table S1. Nonlinear dimensional reduction (tSNE) was used to visualize clustering results, using a perplexity of 50, depicted in Figure 3. Uveal melanoma (UVM) and thyroid carcinomas (THCAs) cluster in close spatial proximity.

4. Discussion

Similar to advanced BCC, kidney tumours such as KICH carry a heterogeneous microenvironment comprising both malignant and solid stromal cells, and a particularly elevated number of macrophages, and a Th2 cytokine shift. In particular, tumour-associated macrophages (TAMs), such as M2 macrophages, promote tumour growth through angiogenesis, which is enhanced when Th2-associated cytokines such as IL-4 and IL-10 are upregulated. Angiogenesis consists of the formation of new capillaries following ECM degradation, aided by matrix metalloproteinases (MMPs) and urokinase-like plasminogen activators (uPas) and their regulators. This process provides the tumour with oxygen, nutrients, and a pathway into circulation to metastasize [19]. Another similarity with BCC is that only a very low percentage of kidney chromophobe metastasizes: some studies estimate the risk of distant spread in KICH as low as 1.3%, compared to 0.0028–0.55% for BCC (average 0.04%) [20,21,22]. Thirdly, KICH differentiates itself from other kidney tumour cancers through its high mutation rate in the TERT gene promoter, in PTEN, and in TP53 [23]. TP53 is the second most commonly mutated gene in BCC after Sonic hedgehog alterations such as PTCH1 loss of function [24]. Moreover, studies focusing on mutations in the TERT promoter have been gaining much traction, with over 50% of BCC samples harboring the genomic aberration, specifically with C > T or CC > TT changes, distinctive for UV exposure [25,26,27]. BCC has the highest mutational burden of any human cancer, in contrast to KICH which has a very low mutational burden [24,28].
Currently, there are no existing precision oncology therapies for metastatic kidney chromophobe cancers due to its rarity and lack of genomic data [29]. However, it is particularly interesting to observe how both advanced BCC and KICH show contributions from tumour-associated macrophages (TAMs) and elevated Th1/Th2 immune cytokine scores [13,19]. Potential targeted therapies to reprogram or eliminate TAMs in KICH may also apply to advanced BCC [19]. In the future, the relationships between large amounts of TAMs, low metastatic rate, and mutations in TP53 and TERT promoter should be further examined.
Next on the list, sarcoma (SARC) was selected as the second BCC “relative”, specifically regarding its total lymphocytes, Tgammadelta, NK, and CD4 cells scores. Sarcomas have over 80 subtypes and comprise malignant tumours concerning non-epithelial connective tissue [30]. Molecular findings across various types of sarcomas are quite heterogeneous [31]. For most subtypes, first-line management includes a combination of surgery and radiation therapy, the latter either used adjuvantly and/or neoadjuvantly. Some sarcomas, such as dermatofibrosarcoma protuberans (DFSP), have actionable targets. In the case of DFSP, imatinib and other tyrosine kinase inhibitors can alter the overexpressed COL1A1-PDGFBB fusion proteins resulting in constitutive activation of PDGFR signaling [32]. Unfortunately, 25–50% of sarcoma patients still develop recurrent and/or metastatic disease after the surgery, which has prompted research to be focused on immunotherapy as a future treatment, in the form of immune checkpoint inhibitors, adoptive T cell transfer, or cancer vaccine [30].
Another common subtype of soft tissue sarcoma, myxofibrosarcoma, could be of great interest to BCC given the peculiar fibromyxoid stroma characterizing BCC [33]. For BCC, this stroma is believed to protect tumour cells and play a role in cancer promotion [34]. Myxofibrosarcoma tumour cells display complex karyotypes without a clear recurrent structural variant and molecular events are still poorly understood. However, the most frequent structural variation in this sarcoma is loss of chromosome arm 13q, which leads to the loss of well-characterized tumour suppressor RB in sarcomatous cells [35]. Circumventing dysfunctional Rb is a potentially actionable target [36].
The third “relative” from this analysis is skin cutaneous melanoma (SKCM), not surprising given that cutaneous melanoma is also a skin cancer and shares similar risk factors to BCC, such as ultraviolet radiation, fair skin complexion, familial history of skin cancer, genetic susceptibility, and numbers of nevi (acquired common and atypical nevi). Driver mutations of SKCM include mutations in BRAF, NRAS, and KIT oncogenes, and in NF1, PTEN, and TP53 tumour suppressors, among others [37]. Some driver mutations are sex-specific, such as DDX3X in males [15]. The current standard of treatment is surgical resection, followed by adjuvant treatment with immune checkpoint inhibitors (ICIs) or targeted therapy drugs, depending on the tumour staging [38]. It has been shown that BRAF mutations are often absent in non-melanoma skin cancers [39]; however, other less-frequently occurring mutations may provide insight on BCC. Mutations in KIT, a known tumour marker and proto-oncogene, accounts for about 3% of all melanomas, but are more prevalent in melanomas on chronically sun-damaged skin [40]. There have been conflicting results on whether c-KIT is implicated in BCC tumorigenesis, some in support [41,42,43], and others proving otherwise [44,45]. Despite this, the use of tyrosine inhibitors on BCC remains unfavourable, especially since its use in AML patients mainly caused BCC as a secondary cancer [46]. A more promising, yet under-researched, target may be JNK. JNK signaling is highly responsive to UV radiation and inflammatory cytokines, among others. Its downstream effects include upregulating AP1, a JNK effector responsible for extracellular matrix (ECM) remodeling [47]. In melanoma, JNK has been found to be implicated in melanoma growth and progression, initiated through activation by IL-1β, whose receptors are highly expressed on Th2 cells [48,49]. Concerning BCC, it has been reported that JNK is a master mediator for numerous key pathways, such as Wnt, Shh, and YAP [50]. In addition, the aforementioned involvement in the TERT promoter gene, which is mutated in over 70% of melanoma cases, and over 50% in BCC, may be of interest as a therapeutic target [51]. TERT mutations have been shown to be positively correlated with a higher tumour mutational burden (TMB) value, neoantigen load, and tumour purity. Its expression as a tumour-associated antigen has primed it to be an ideal target of immunotherapy, with numerous studies and trials currently in development and underway [52]. As such, along with its crucial role in melanoma resistance, more studies on JNK and TERT promoter in BCC may uncover more information on advanced and treatment-resistant BCC.
Omics profiling has yielded important clues into tumour initiation, promotion, and progression which resulted in development of targeted therapies [53]. For BCC, Sonic hedgehog inhibitors that selectively inhibit SMO include vismodegib (Erivedge) [8] and sonidegib (Odomzo) [54]. Vismodegib has been shown to lead to the recruitment of cytotoxic T cells into the tumour microenvironment, with concurrent upregulation of type I MHC in BCC tumour cells [55]. Treatment of SMO inhibitors leads to at least a partial response in ~2/3 of advanced BCC patients, but >25% of patients discontinue therapy due to frequent, hard-to-tolerate side effects such as alopecia, dysgeusia, and muscle spasms [9]. BCC tumour cells often acquire resistance to SMO inhibitors, which may result from a cell identity switch towards a mesenchymal-stem-cell-like phenotype [56]. Despite its very high tumour mutational burden, anti-PD-1 therapy such as cemiplimab fails in ~70% of advanced BCC patients [57].

5. Conclusions

We have identified three potential BCC “relatives” sharing similar tumour microenvironment findings. We also reviewed underlying literature for current management therapies and potential novel approaches in these three relatives which may apply to advanced BCC patients. Kidney chromophobe particularly stood out as a promising BCC “relative”, demonstrating strong similarities in TAM and Th2 cytokine profiles. Targeting the RB pathway in myxofibrosarcoma may be of interest regarding modulating the fibromyxoid stroma of BCC, rendering it less suitable to promote BCC growth. Exploring lesser-known contributors to melanoma, such as c-KIT and JNK, may also be of interest. Although there still does not exist a definitive treatment for advanced BCC, new targeted therapeutics used in other malignancies with closely related microenvironment changes might provide novel management strategies to be studied and potentially employed in the future.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15010305/s1, Table S1: Summary of all significant and top three (if available, in bold) closest BCC “relatives” by cell type.

Author Contributions

Conceptualization, E.-R.Z., S.G., P.X. and P.L.; methodology, S.G., P.X. and P.L.; software, P.L., E.-R.Z. and S.G.; validation, S.G.; formal analysis, E.-R.Z., S.G., P.X. and P.L.; investigation, E.-R.Z., S.G. and M.F.; resources, P.L.; data curation, P.X. and P.L.; writing—original draft preparation, E.-R.Z., S.G., P.X., M.F., I.V.L. and P.L.; writing—review and editing, E.-R.Z., S.G. and P.L.; visualization, S.G.; supervision, S.G., P.X., I.V.L. and P.L.; project administration, P.L.; funding acquisition, P.L. and E.-R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

P.L. received grants from the Lady Davis Institute for Medical Research, from the Jewish General Hospital Foundation, from the Jewish General Hospital Department of Medicine, from the Fonds de recherche du Québec—Santé (#312768 and #324151), and from the Marathon of Hope Cancer Centre Network—Terry Fox Research Institute (#3253) for this work. This research was enabled in part by support provided by Calcul Québec (www.calculquebec.ca), Compute Ontario (www.computeontario.ca), WestGrid (www.westgrid.ca), and Compute Canada (www.computecanada.ca) (RRG #4318). E.R.Z. received support from the following scholarship: Dr. Raymond and Dr. Leora Kahn Research Bursary in Memory of Mrs. Marion Kahn.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Previously published whole-genome RNA-Seq data originated from: Atwood et al. (GEO accession number GSE58377) [58], Bonilla et al. (EGA accession number EGAS00001001540) [59], and Sharpe et al. (EGA accession number EGAS00001000845) [60]. The pan-cancer datasets analyzed are available from The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/ accessed on 19 June 2019).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Muzic, J.G.; Schmitt, A.R.; Wright, A.C.; Alniemi, D.T.; Zubair, A.S.; Lourido, J.M.O.; Seda, I.M.S.; Weaver, A.L.; Baum, C.L. Incidence and Trends of Basal Cell Carcinoma and Cutaneous Squamous Cell Carcinoma. Mayo Clin. Proc. 2017, 92, 890–898. [Google Scholar] [CrossRef] [PubMed]
  2. Chlebicka, I.; Stefaniak, A.; Matusiak, Ł.; Szepietowski, J. Basal cell carcinoma: What new can be learned about the most common human cancer? A cross-sectional prospective study of 180 cases in a single centre. Adv. Dermatol. Allergol. 2021, 38, 1086–1091. [Google Scholar] [CrossRef] [PubMed]
  3. Miller, D.L.; Weinstock, M.A. Nonmelanoma skin cancer in the United States: Incidence. J. Am. Acad. Dermatol. 1994, 30 Pt 1, 774–778. [Google Scholar] [CrossRef]
  4. Carballo, G.B.; Honorato, J.R.; De Lopes, G.P.F.; de Sampaio e Spohr, T.C.L. A highlight on Sonic hedgehog pathway. Cell Commun. Signal. 2018, 16, 11. [Google Scholar] [CrossRef]
  5. Marzuka, A.G.; Book, S.E. Basal Cell Carcinoma: Pathogenesis, Epidemiology, Clinical Features, Diagnosis, Histopathology, and Management. Yale J. Boil. Med. 2015, 88, 167–179. [Google Scholar]
  6. Quinn, A.G.; Epstein, J.E. Patched, Hedgehog, and Skin Cancer; Humana Press: Totowa, NJ, USA, 2003; pp. 85–95. [Google Scholar]
  7. Rubin, A.I.; Chen, E.H.; Ratner, D. Basal-Cell Carcinoma. N. Engl. J. Med. 2005, 353, 2262–2269. [Google Scholar] [CrossRef] [PubMed]
  8. Sekulic, A.; Migden, M.R.; Oro, A.E.; Dirix, L.; Lewis, K.D.; Hainsworth, J.D.; Solomon, J.A.; Yoo, S.; Arron, S.T.; Friedlander, P.A.; et al. Efficacy and Safety of Vismodegib in Advanced Basal-Cell Carcinoma. N. Engl. J. Med. 2012, 366, 2171–2179. [Google Scholar] [CrossRef] [Green Version]
  9. Xie, P.; Lefrançois, P. Efficacy, safety, and comparison of sonic hedgehog inhibitors in basal cell carcinomas: A systematic review and meta-analysis. J. Am. Acad. Dermatol. 2018, 79, 1089–1100.e17. [Google Scholar] [CrossRef]
  10. Mohan, S.V.; Chang, A.L.S. Advanced Basal Cell Carcinoma: Epidemiology and Therapeutic Innovations. Curr. Dermatol. Rep. 2014, 3, 40–45. [Google Scholar] [CrossRef] [Green Version]
  11. Jiang, L.-P.; Shen, Q.-S.; Yang, C.-P.; Chen, Y.-B. Establishment of basal cell carcinoma animal model in Chinese tree shrew (Tupaia belangeri chinensis). Zool. Res. 2017, 38, 180–190. [Google Scholar] [CrossRef] [Green Version]
  12. Noubissi, F.K.; Kim, T.; Kawahara, T.N.; Aughenbaugh, W.D.; Berg, E.; Longley, B.; Athar, M.; Spiegelman, V.S. Role of CRD-BP in the Growth of Human Basal Cell Carcinoma Cells. J. Investig. Dermatol. 2014, 134, 1718–1724. [Google Scholar] [CrossRef] [PubMed]
  13. Lefrançois, P.; Xie, P.; Gunn, S.; Gantchev, J.; Villarreal, A.M.; Sasseville, D.; Litvinov, I.V. In silico analyses of the tumor microenvironment highlight tumoral inflammation, a Th2 cytokine shift and a mesenchymal stem cell-like phenotype in advanced in basal cell carcinomas. J. Cell Commun. Signal. 2020, 14, 245–254. [Google Scholar] [CrossRef] [PubMed]
  14. Litvinov, I.V.; Xie, P.; Gunn, S.; Sasseville, D.; Lefrançois, P. The transcriptional landscape analysis of basal cell carcinomas reveals novel signalling pathways and actionable targets. Life Sci. Alliance 2021, 4, e202000651. [Google Scholar] [CrossRef] [PubMed]
  15. Alkallas, R.; Lajoie, M.; Moldoveanu, D.; Hoang, K.V.; Lefrançois, P.; Lingrand, M.; Ahanfeshar-Adams, M.; Watters, K.; Spatz, A.; Zippin, J.H.; et al. Multi-omic analysis reveals significantly mutated genes and DDX3X as a sex-specific tumor suppressor in cutaneous melanoma. Nat. Cancer 2020, 1, 635–652. [Google Scholar] [CrossRef]
  16. Newman, A.M.; Liu, C.L.; Green, M.R.; Gentles, A.J.; Feng, W.; Xu, Y.; Hoang, C.D.; Diehn, M.; Alizadeh, A.A. Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 2015, 12, 453–457. [Google Scholar] [CrossRef] [Green Version]
  17. Lefrançois, P.; Xie, P.; Wang, L.; Tetzlaff, M.T.; Moreau, L.; Watters, A.K.; Netchiporouk, E.; Provost, N.; Gilbert, M.; Ni, X.; et al. Gene expression profiling and immune cell-type deconvolution highlight robust disease progression and survival markers in multiple cohorts of CTCL patients. OncoImmunology 2018, 7, e1467856. [Google Scholar] [CrossRef] [Green Version]
  18. Aran, D.; Hu, Z.; Butte, A.J. xCell: Digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017, 18, 220. [Google Scholar] [CrossRef] [Green Version]
  19. Kovaleva, O.V.; Samoilova, D.V.; Shitova, M.S.; Gratchev, A. Tumor Associated Macrophages in Kidney Cancer. Anal. Cell. Pathol. 2016, 2016, 9307549. [Google Scholar] [CrossRef] [Green Version]
  20. Goldenberg, G.; Karagiannis, T.; Palmer, J.B.; Lotya, J.; O’Neill, C.; Kisa, R.; Herrera, V.; Siegel, D.M. Incidence and prevalence of basal cell carcinoma (BCC) and locally advanced BCC (LABCC) in a large commercially insured population in the United States: A retrospective cohort study. J. Am. Acad. Dermatol. 2016, 75, 957–966.e2. [Google Scholar] [CrossRef]
  21. Volpe, A.; Novara, G.; Antonelli, A.; Bertini, R.; Billia, M.; Carmignani, G.; Cunico, S.C.; Longo, N.; Martignoni, G.; Minervini, A.; et al. Chromophobe renal cell carcinoma (RCC): Oncological outcomes and prognostic factors in a large multicentre series. BJU Int. 2011, 110, 76–83. [Google Scholar] [CrossRef]
  22. Piva de Freitas, P.; Senna, C.G.; Tabai, M.; Chone, C.T.; Altemani, A. Metastatic Basal Cell Carcinoma: A Rare Manifestation of a Common Disease. Case Rep. Med. 2017, 2017, 8929745. [Google Scholar] [CrossRef] [PubMed]
  23. Shea, C.R.; McNutt, N.S.; Volkenandt, M.; Lugo, J.; Prioleau, P.G.; Albino, A.P. Overexpression of p53 protein in basal cell carcinomas of human skin. Am. J. Pathol. 1992, 141, 25–29. [Google Scholar] [PubMed]
  24. Jayaraman, S.S.; Rayhan, D.J.; Hazany, S.; Kolodney, M.S. Mutational Landscape of Basal Cell Carcinomas by Whole-Exome Sequencing. J. Investig. Dermatol. 2014, 134, 213–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Griewank, K.G.; Murali, R.; Schilling, B.; Schimming, T.; Möller, I.; Moll, I.; Schwamborn, M.; Sucker, A.; Zimmer, L.; Schadendorf, D.; et al. TERT Promoter Mutations Are Frequent in Cutaneous Basal Cell Carcinoma and Squamous Cell Carcinoma. PLoS ONE 2013, 8, e80354. [Google Scholar] [CrossRef]
  26. Pópulo, H.; Boaventura, P.; Vinagre, J.; Batista, R.; Mendes, A.; Caldas, R.; Pardal, J.; Azevedo, F.; Honavar, M.; Guimarães, I.; et al. TERT Promoter Mutations in Skin Cancer: The Effects of Sun Exposure and X-Irradiation. J. Investig. Dermatol. 2014, 134, 2251–2257. [Google Scholar] [CrossRef] [Green Version]
  27. Scott, G.A.; Laughlin, T.S.; Rothberg, P.G. Mutations of the TERT promoter are common in basal cell carcinoma and squamous cell carcinoma. Mod. Pathol. 2014, 27, 516–523. [Google Scholar] [CrossRef] [Green Version]
  28. Garje, R.; Elhag, D.; A Yasin, H.; Acharya, L.; Vaena, D.; Dahmoush, L. Comprehensive review of chromophobe renal cell carcinoma. Crit. Rev. Oncol./Hematol. 2021, 160, 103287. [Google Scholar] [CrossRef]
  29. Papanikolaou, D.; Ioannidou, P.; Koukourikis, P.; Moysidis, K.; Meditskou, S.; Koutsoumparis, D.; Hatzimouratidis, K.; Hatzivassiliou, E. Systemic therapy for chromophobe renal cell carcinoma: A systematic review. Urol. Oncol. Semin. Orig. Investig. 2020, 38, 137–149. [Google Scholar] [CrossRef]
  30. Clemente, O.; Ottaiano, A.; Di Lorenzo, G.; Bracigliano, A.; Lamia, S.; Cannella, L.; Pizzolorusso, A.; Di Marzo, M.; Santorsola, M.; De Chiara, A.; et al. Is immunotherapy in the future of therapeutic management of sarcomas? J. Transl. Med. 2021, 19, 173. [Google Scholar] [CrossRef]
  31. Dancsok, A.R.; Asleh-Aburaya, K.; Nielsen, T.O. Advances in sarcoma diagnostics and treatment. Oncotarget 2016, 8, 7068–7093. [Google Scholar] [CrossRef] [Green Version]
  32. Navarrete-Dechent, C.; Mori, S.; Barker, C.A.; Dickson, M.A.; Nehal, K.S. Imatinib Treatment for Locally Advanced or Metastatic Dermatofibrosarcoma Protuberans. JAMA Dermatol. 2019, 155, 361–369. [Google Scholar] [CrossRef] [PubMed]
  33. Kaporis, H.G.; Guttman-Yassky, E.; A Lowes, M.; Haider, A.S.; Fuentes-Duculan, J.; Darabi, K.; Whynot-Ertelt, J.; Khatcherian, A.; Cardinale, I.; Novitskaya, I.; et al. Human Basal Cell Carcinoma Is Associated with Foxp3+ T cells in a Th2 Dominant Microenvironment. J. Investig. Dermatol. 2007, 127, 2391–2398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Mackiewicz-Wysocka, M.; Bowszyc-Dmochowska, M.; Strzelecka-Węklar, D.; Dańczak-Pazdrowska, A.; Adamski, Z. Reviews Basal cell carcinoma—Diagnosis. Contemp. Oncol./Współczesna Onkol. 2013, 4, 337–342. [Google Scholar] [CrossRef] [PubMed]
  35. Nishio, J.; Iwasaki, H.; Nabeshima, K.; Naito, M. Cytogenetics and Molecular Genetics of Myxoid Soft-Tissue Sarcomas. Genet. Res. Int. 2011, 2011, 497148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Witkiewicz, A.K.; Chung, S.; Brough, R.; Vail, P.; Franco, J.; Lord, C.; Knudsen, E.S. Targeting the Vulnerability of RB Tumor Suppressor Loss in Triple-Negative Breast Cancer. Cell Rep. 2018, 22, 1185–1199. [Google Scholar] [CrossRef] [Green Version]
  37. Van Allen, E.M.; Miao, D.; Schilling, B.; Shukla, S.A.; Blank, C.; Zimmer, L.; Sucker, A.; Hillen, U.; Foppen, M.H.G.; Goldinger, S.M.; et al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science 2015, 350, 207–211. [Google Scholar] [CrossRef] [Green Version]
  38. I Daud, A.; Daud; Soon, C.W. Treatment of cutaneous melanoma: Current approaches and future prospects. Cancer Manag. Res. 2010, 2, 197–211. [Google Scholar] [CrossRef]
  39. Libra, M.; Malaponte, G.; Bevelacqua, V.; Siciliano, R.; Castrogiovanni, P.; Fulvi, A.; Micali, G.; Ligresti, G.; Mazzarino, M.C.; Stivala, F.; et al. Absence of BRAF Gene Mutation in Non-Melanoma Skin Tumors. Cell Cycle 2006, 5, 968–970. [Google Scholar] [CrossRef] [Green Version]
  40. Curtin, J.A.; Busam, K.; Pinkel, D.; Bastian, B.C. Somatic Activation of KIT in Distinct Subtypes of Melanoma. J. Clin. Oncol. 2006, 24, 4340–4346. [Google Scholar] [CrossRef]
  41. Terada, T. Expression of NCAM (CD56), chromogranin A, synaptophysin, c-KIT (CD117) and PDGFRA in normal non-neoplastic skin and basal cell carcinoma: An immunohistochemical study of 66 consecutive cases. Med. Oncol. 2013, 30, 444. [Google Scholar] [CrossRef]
  42. Leon, A.; Ceauşu, Z.; Ceauşu, M.; Ardeleanu, C.; Mehedinţi, R. Mast cells and dendritic cells in basal cell carcinoma. Rom. J. Morphol. Embryol. 2009, 50, 85–90. [Google Scholar] [PubMed]
  43. Castillo, J.-M.; Knol, A.-C.; Nguyen, J.-M.; Khammari, A.; Saint-Jean, M.; Dreno, B. Immunohistochemical markers of advanced basal cell carcinoma: CD56 is associated with a lack of response to vismodegib. Eur. J. Dermatol. 2016, 26, 452–459. [Google Scholar] [CrossRef] [PubMed]
  44. Dessauvagie, B.F.; Wood, B.A. CD117 and CD43 are useful adjuncts in the distinction of adenoid cystic carcinoma from adenoid basal cell carcinoma. Pathology 2015, 47, 130–133. [Google Scholar] [CrossRef] [PubMed]
  45. Yang, D.T.; Holden, J.A.; Florell, S.R. CD117, CK20, TTF-1, and DNA topoisomerase II-α antigen expression in small cell tumors. J. Cutan. Pathol. 2004, 31, 254–261. [Google Scholar] [CrossRef]
  46. Verma, D.; Kantarjian, H.; Strom, S.S.; Rios, M.B.; Jabbour, E.; Quintas-Cardama, A.; Verstovsek, S.; Ravandi, F.; O’Brien, S.; Cortes, J. Malignancies occurring during therapy with tyrosine kinase inhibitors (TKIs) for chronic myeloid leukemia (CML) and other hematologic malignancies. Blood 2011, 118, 4353–4358. [Google Scholar] [CrossRef]
  47. Jin, Y.-J.; Park, I.; Hong, I.-K.; Byun, H.-J.; Choi, J.; Kim, Y.-M.; Lee, H. Fibronectin and vitronectin induce AP-1-mediated matrix metalloproteinase-9 expression through integrin α5β1/αvβ3-dependent Akt, ERK and JNK signaling pathways in human umbilical vein endothelial cells. Cell. Signal. 2011, 23, 125–134. [Google Scholar] [CrossRef]
  48. Qin, Y.; Milton, D.R.; Oba, J.; Ding, Z.; Lizée, G.; Ekmekcioglu, S.; Grimm, E.A. Inflammatory IL-1β-driven JNK activation in stage III melanoma. Pigment Cell Melanoma Res. 2015, 28, 236–239. [Google Scholar] [CrossRef]
  49. Santarlasci, V.; Cosmi, L.; Maggi, L.; Liotta, F.; Annunziato, F. IL-1 and T Helper Immune Responses. Front. Immunol. 2013, 4, 182. [Google Scholar] [CrossRef] [Green Version]
  50. Hammouda, M.B.; Ford, A.E.; Liu, Y.; Zhang, J.Y. The JNK Signaling Pathway in Inflammatory Skin Disorders and Cancer. Cells 2020, 9, 857. [Google Scholar] [CrossRef] [Green Version]
  51. Huang, F.W.; Hodis, E.; Xu, M.J.; Kryukov, G.V.; Chin, L.; Garraway, L.A. Highly Recurrent TERT Promoter Mutations in Human Melanoma. Science 2013, 339, 957–959. [Google Scholar] [CrossRef] [Green Version]
  52. Negrini, S.; De Palma, R.; Filaci, G. Anti-cancer Immunotherapies Targeting Telomerase. Cancers 2020, 12, 2260. [Google Scholar] [CrossRef] [PubMed]
  53. Scatena, C.; Murtas, D.; Tomei, S. Cutaneous Melanoma Classification: The Importance of High-Throughput Genomic Technologies. Front. Oncol. 2021, 11, 635488. [Google Scholar] [CrossRef] [PubMed]
  54. Migden, M.R.; Guminski, A.; Gutzmer, R.; Dirix, L.; Lewis, K.D.; Combemale, P.; Herd, R.M.; Kudchadkar, R.; Trefzer, U.; Gogov, S.; et al. Treatment with two different doses of sonidegib in patients with locally advanced or metastatic basal cell carcinoma (BOLT): A multicentre, randomised, double-blind phase 2 trial. Lancet Oncol. 2015, 16, 716–728. [Google Scholar] [CrossRef] [PubMed]
  55. Passarelli, A.; Galdo, G.; Aieta, M.; Fabrizio, T.; Villonio, A.; Conca, R. A Vismodegib Experience in Elderly Patients with Basal Cell Carcinoma: Case Reports and Review of the Literature. Int. J. Mol. Sci. 2020, 21, 8596. [Google Scholar] [CrossRef] [PubMed]
  56. Biehs, B.; Dijkgraaf, G.J.P.; Piskol, R.; Alicke, B.; Boumahdi, S.; Peale, F.; Gould, S.E.; de Sauvage, F.J. A cell identity switch allows residual BCC to survive Hedgehog pathway inhibition. Nature 2018, 562, 429–433. [Google Scholar] [CrossRef] [PubMed]
  57. Stratigos, A.J.; Sekulic, A.; Peris, K.; Bechter, O.; Prey, S.; Kaatz, M.; Lewis, K.D.; Basset-Seguin, N.; Chang, A.L.S.; Dalle, S.; et al. Cemiplimab in locally advanced basal cell carcinoma after hedgehog inhibitor therapy: An open-label, multi-centre, single-arm, phase 2 trial. Lancet Oncol. 2021, 22, 848–857. [Google Scholar] [CrossRef]
  58. Atwood, S.X.; Sarin, K.Y.; Whitson, R.J.; Li, J.R.; Kim, G.; Rezaee, M.; Ally, M.S.; Kim, J.; Yao, C.; Chang, A.L.S.; et al. Smoothened Variants Explain the Majority of Drug Resistance in Basal Cell Carcinoma. Cancer Cell 2015, 27, 342–353. [Google Scholar] [CrossRef] [Green Version]
  59. Bonilla, X.; Parmentier, L.; King, B.; Bezrukov, F.; Kaya, G.; Zoete, V.; Seplyarskiy, V.B.; Sharpe, H.J.; McKee, T.; Letourneau, A.; et al. Genomic analysis identifies new drivers and progression pathways in skin basal cell carcinoma. Nat. Genet. 2016, 48, 398–406. [Google Scholar] [CrossRef]
  60. Sharpe, H.J.; Pau, G.; Dijkgraaf, G.J.; Basset-Seguin, N.; Modrusan, Z.; Januario, T.; Tsui, V.; Durham, A.B.; Dlugosz, A.A.; Haverty, P.M.; et al. Genomic Analysis of Smoothened Inhibitor Resistance in Basal Cell Carcinoma. Cancer Cell 2015, 27, 327–341. [Google Scholar] [CrossRef]
Figure 1. BCC vs. BCC “relatives”: xCell. In silico cell type enrichment scores obtained using xCell, comparing BCC (dark red) to other BCC “relatives” (white), with statistically similar cancers denoted (grey), with the top three in colour (dark blue, purple, violet) with their corresponding p-values above. Truncated violin plots of Th2 immune cytokine (a), and mesenchymal stem cells (MSCs) (b) are shown in order of ascending median. BCC median, 1st quartile, and 3rd quartile are delineated in red (where applicable).
Figure 1. BCC vs. BCC “relatives”: xCell. In silico cell type enrichment scores obtained using xCell, comparing BCC (dark red) to other BCC “relatives” (white), with statistically similar cancers denoted (grey), with the top three in colour (dark blue, purple, violet) with their corresponding p-values above. Truncated violin plots of Th2 immune cytokine (a), and mesenchymal stem cells (MSCs) (b) are shown in order of ascending median. BCC median, 1st quartile, and 3rd quartile are delineated in red (where applicable).
Cancers 15 00305 g001
Figure 2. BCC vs. BCC “relatives”: CIBERSORTx. In silico immune cell fractions obtained using CIBERSORTx, comparing BCC (dark red) to other BCC “relatives” (white), with statistically similar cancers denoted (grey), with the top three in colour (dark blue, purple, violet) with their corresponding p-values above. Truncated violin plots of total CD8 T cells (a), T gamma delta cells (b), total B-cells (c), total lymphocytes (d), total macrophages (e), total NK cells (f), regulatory T cells (Treg) (g), and total CD4+ cells (h) are shown in order of ascending median. BCC median, 1st quartile, and 3rd quartile are delineated in red (where applicable).
Figure 2. BCC vs. BCC “relatives”: CIBERSORTx. In silico immune cell fractions obtained using CIBERSORTx, comparing BCC (dark red) to other BCC “relatives” (white), with statistically similar cancers denoted (grey), with the top three in colour (dark blue, purple, violet) with their corresponding p-values above. Truncated violin plots of total CD8 T cells (a), T gamma delta cells (b), total B-cells (c), total lymphocytes (d), total macrophages (e), total NK cells (f), regulatory T cells (Treg) (g), and total CD4+ cells (h) are shown in order of ascending median. BCC median, 1st quartile, and 3rd quartile are delineated in red (where applicable).
Cancers 15 00305 g002
Figure 3. t-distributed stochastic neighbour embedding plot of the tumour microenvironment scores (from CIBERSORTx and xCell) of 33 TGCA cancers with BCC, coloured by cancer type.
Figure 3. t-distributed stochastic neighbour embedding plot of the tumour microenvironment scores (from CIBERSORTx and xCell) of 33 TGCA cancers with BCC, coloured by cancer type.
Cancers 15 00305 g003
Table 1. Summary of the top three significant (if available) closest BCC “relatives” by cell type.
Table 1. Summary of the top three significant (if available) closest BCC “relatives” by cell type.
Cell Type BCC Score (Median), Q1, Q3BCC “Relatives” Score (Median)Q1Q3p-Value
Th24.188 × 10−18, 0, 0.01980129PRAD = 1.07 × 10−1700.017065650.6102
KICH = 0.00105670500.0073810450.151
KIRC = 2.23 × 10−1700.030514090.1368
MSC0.188433346, 0.08336783, 0.28205374SKCM = 0.160968670.072024250.254464260.3872
Total CD8+0.00445073, 0, 0.043739228GBM = 0.0104815500.030218790.2741
LGG = 000.019965480.05719
Tgd0, 0, 0.008231755THCA = 000.0080639120.7752
CHOL = 000.012820030.6099
LIHC = 000.0067202860.5979
B-cells 0.04857795, 0.020757612, 0.086950544THYM = 0.0551831040.031081080.084553590.3715
GBM = 0.048157240.0283188760.0753449910.2374
Total Lymphocytes0.39140159, 0.3132104, 0.4688854SARC = 0.4169986540.31243020.51537580.2290
Macrophages0.25033106, 0.1763854, 0.3784326ESCA = 0.23312470.16757340.33521330.5397
KICH = 0.2767290.18433250.39498210.4904
BLCA = 0.25906530.166937270.337628170.3333
Total NK0.0325587, 0.005046464, 0.053041073THYM = 0.026824970.0094086740.0491247710.6696
ESCA = 0.025239730.011523860.046047260.4214
KIRP = 0.028787490.013828300.048743960.3381
Tregs 0, 0, 0ACC = 0000.9991
LAML = 000.003150340.3533
Total CD4+0.15296126, 0.1031139, 0.2082608CHOL = 0.1561890.13211290.19526960.6959
SARC = 0.15625860.11542910.20837010.6334
SKCM = 14126040.087283840.191687630.5715
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, E.-R.; Ghezelbash, S.; Xie, P.; Fotovati, M.; Litvinov, I.V.; Lefrançois, P. Comparison of the Basal Cell Carcinoma (BCC) Tumour Microenvironment to Other Solid Malignancies. Cancers 2023, 15, 305. https://doi.org/10.3390/cancers15010305

AMA Style

Zhang E-R, Ghezelbash S, Xie P, Fotovati M, Litvinov IV, Lefrançois P. Comparison of the Basal Cell Carcinoma (BCC) Tumour Microenvironment to Other Solid Malignancies. Cancers. 2023; 15(1):305. https://doi.org/10.3390/cancers15010305

Chicago/Turabian Style

Zhang, Eliana-Ruobing, Sarah Ghezelbash, Pingxing Xie, Misha Fotovati, Ivan V. Litvinov, and Philippe Lefrançois. 2023. "Comparison of the Basal Cell Carcinoma (BCC) Tumour Microenvironment to Other Solid Malignancies" Cancers 15, no. 1: 305. https://doi.org/10.3390/cancers15010305

APA Style

Zhang, E. -R., Ghezelbash, S., Xie, P., Fotovati, M., Litvinov, I. V., & Lefrançois, P. (2023). Comparison of the Basal Cell Carcinoma (BCC) Tumour Microenvironment to Other Solid Malignancies. Cancers, 15(1), 305. https://doi.org/10.3390/cancers15010305

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop