Next Article in Journal
Management of Resectable Stage III-N2 Non-Small-Cell Lung Cancer (NSCLC) in the Age of Immunotherapy
Previous Article in Journal
Multi-Stage Harmonization for Robust AI across Breast MR Databases
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 19
10.3390/cancers13194810
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Targeted Sequencing Revealed Distinct Mutational Profiles of Ocular and Extraocular Sebaceous Carcinomas

by
Hee Young Na
1,2,†,
Jeong Hwan Park
1,3,†,
Sun Ah Shin
1,4,
Sejoon Lee
5,
Heonyi Lee
6,
Heejoon Chae
7,
HoKyung Choung
8,
Namju Kim
9,
Jin-Haeng Chung
1,2,* and
Ji Eun Kim
1,3,*
1
Department of Pathology, Seoul National University College of Medicine, Seoul 03080, Korea
2
Department of Pathology and Translational Medicine, Seoul National University Bundang Hospital, Seongnam 13620, Korea
3
Department of Pathology, Seoul Metropolitan Government-Seoul National University Boramae Medical Center, Seoul 07067, Korea
4
Department of Pathology, National Cancer Center, Goyang 10408, Korea
5
Precision Medicine Center, Seoul National University Bundang Hospital, Seongnam 13620, Korea
6
Bioinformatics Collaboration Unit, Department of Biomedical Systems Informatics, Yonsei University College of Medicine, Seoul 03722, Korea
7
Division of Computer Science, Sookmyung Women’s University, Seoul 04312, Korea
8
Department of Ophthalmology, Seoul Metropolitan Government-Seoul National University Boramae Medical Center, Seoul 07067, Korea
9
Department of Ophthalmology, Seoul National University Bundang Hospital, Seongnam 13620, Korea
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Cancers 2021, 13(19), 4810; https://doi.org/10.3390/cancers13194810
Submission received: 19 August 2021 / Revised: 10 September 2021 / Accepted: 22 September 2021 / Published: 26 September 2021
(This article belongs to the Section Molecular Cancer Biology)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

The molecular pathogenesis of sebaceous carcinoma is largely unknown so far. The main aim of our study was to investigate genetic alterations in ocular and extraocular sebaceous carcinoma in Korean patients. By applying targeted next-generation sequencing, we demonstrated a distinct mutational profile in ocular and extraocular sebaceous carcinoma. In addition, clinically actionable variants involving tyrosine kinase genes and homologous recombination deficiency-associated genes were detected, suggesting the possibility of targeted therapeutic trials in intractable sebaceous carcinoma cases.

Abstract

The biological behavior of sebaceous carcinoma (SeC) is relatively indolent; however, local invasion or distant metastasis is sometimes reported. Nevertheless, a lack of understanding of the genetic background of SeC makes it difficult to apply effective systemic therapy. This study was designed to investigate major genetic alterations in SeCs in Korean patients. A total of 29 samples, including 20 ocular SeCs (SeC-Os) and 9 extraocular SeCs (SeC-EOs), were examined. Targeted next-generation sequencing tests including 171 cancer-related genes were performed. TP53 and PIK3CA genes were frequently mutated in both SeC-Os and SeC-EOs with slight predominance in SeC-Os, whereas the NOTCH1 gene was more commonly mutated in SeC-EOs. In clinical correlation, mutations in RUNX1 and ATM were associated with development of distant metastases, and alterations in MSH6 and BRCA1 were associated with inferior progression-free survival (all p < 0.05). In conclusion, our study revealed distinct genetic alterations between SeC-Os and SeC-EOs and some important prognostic molecular markers. Mutations in potentially actionable genes, including EGFR, ERBB2, and mismatch repair genes, were noted, suggesting consideration of a clinical trial in intractable cases.

1. Introduction

Sebaceous carcinoma (SeC) is a rare but aggressive neoplasm that is relatively prevalent in the Asian population [1,2,3]. The cellular origins of SeC are secretory cells producing oily substance or ductal cells of the sebaceous gland [4]. Diagnosis of SeC is very challenging because of the anatomical proximity and histological similarity to hair follicle or sweat glands. It predominantly occurs in periocular areas, where Meibomian glands are enriched, but extraocular sites, including the head and neck region, and less frequently the trunk or breast, can also be involved [1]. Wide excision is the mainstay of treatment for both ocular sebaceous carcinomas (SeC-Os) and extraocular sebaceous carcinomas (SeC-EOs) [5]. However, locally aggressive SeC-O often requires extensive surgery, including orbital exenteration, leading to high morbidity [6,7]. Radiotherapy can be considered in unresectable, locally recurrent, or residual disease [5]. The cancer-specific mortality rate of SeC is 6–22% [8,9,10], and disease progression, including locoregional recurrence and distant metastasis, is reported in 4–28% [1,11,12]. Systemic therapy, including conventional chemotherapy, immunotherapy, and targeted therapies, can be considered in locally aggressive or metastatic SeCs, although the outcome data are scarce [5,13]. The use of conventional chemotherapy composed of Adriamycin and cyclophosphamide or platinum-based chemotherapeutic agents with either paclitaxel or 5-fluorouracil and folinic acid has been reported by some researchers [14,15,16]. Immunotherapy for microsatellite unstable SeCs or targeted therapies, such as retinoic acid receptor-β agonist [17], androgen receptor (AR) antagonist [18], and epidermal growth factor receptor inhibitor [19], are other potential options. However, no single effective targeting agent has been established, primarily due to the lack of understanding of the molecular pathogenesis of this tumor.
Next-generation sequencing (NGS) has opened a new era of high-throughput biology and has enabled researchers to identify large-scale genomic alterations in a relatively short time, contributing to the fine-tuning of cancer management. Accordingly, several outstanding studies of a comprehensive molecular analysis of SeC using NGS were recently published [20,21,22]. Some researchers have proposed that SeC shows genetic signatures similar to those of cutaneous squamous cell carcinoma, demonstrating alterations of TP53, NOTCH1 and NOTCH2 [20,21,22]. However, our previous reports revealed that the molecular pathogenesis of SeC could be compared with that of breast carcinoma because of its ontogenetic proximity and hormone-responsive growth, reflecting the similarity between sebaceous glands and mammary glands [23,24]. Moreover, most molecular studies were performed in Western cohorts, and molecular pathologic data of SeC in Asian populations require further investigation. In the present study, we investigated the genetic profile of SeCs in Korean patients, focusing on a comparison of cutaneous squamous cell carcinoma and mammary carcinoma, expecting to find a new therapeutic target.

2. Materials and Methods

2.1. Patients and Specimens

Surgically resected formalin-fixed paraffin-embedded (FFPE) tissue specimens of 29 SeC cases were obtained from 26 selected patients who were treated at Seoul Metropolitan Government Seoul National University (SMG-SNU) Boramae Medical Center, Seoul National University Bundang Hospital, and Seoul National University from 2000 to 2015. Among these 29 tumors, there were 26 primary lesions, 1 locally recurrent SeC-O, 1 locally recurrent SeC-EO, and 1 metastatic SeC-O. Among the 3 locally recurrent or metastatic lesions, locally recurrent SeC-EO and metastatic SeC-O had matched primary lesions. In addition, 2 primary lesions (1 SeC-O and 1 SeC-EO) from 1 suspected Muir–Torre syndrome (MTS) patient were included. Clinicopathologic information, including age, sex, and stage, was retrieved from the electronic medical records and pathology reports. This study was approved by the Institutional Review Board (IRB) at SMG-SNU Boramae Medical Center (IRB No. 16-2017-17-031).

2.2. DNA Extraction and NGS

Targeted DNA sequencing was performed using the Axen Cancer Panel 2 (Macrogen, Korea) composed of 171 cancer-associated genes (Table S1). After reviewing each slide, tumor areas were carefully chosen for manual macrodissection. The tumor purity ranged from 40% to 60%. To extract the genomic DNA, a QIAamp® DNA FFPE Tissue kit (Qiagen, Hilden, Germany) was used. A total of 200 ng of DNA was used for library generation. The hybrid capture method was used for DNA library preparation and target enrichment, according to Illumina’s standard protocol using an Agilent SureSelectXT Target Enrichment Kit (Agilent Technologies, Santa Clara, CA, USA). The target region was sequenced using the HiSeq 2500 system (Illumina, San Diego, CA, USA). The average mean depth was ×1100.

2.3. Alignment of Reads and Variant Calling

The cutadapt [25] was used to remove adapter sequences from the raw sequencing reads, followed by aligning the reads to the reference genome (GRCh37/hg19) using Burrows–Wheeler Aligner-MEM (BWA-MEM) [26]. Poorly mapped reads with a mapping quality (MAPQ) below 20 were removed using Samtools version 1.3.1 [27]. The MuTect2 algorithm [28] was used to detect somatic mutations, including single nucleotide variants (SNVs) and small insertions and deletions (INDELs). Finally, SnpEff and SnpSift v4.3i [29] with dbNSFP v2.9.3 [30] were used to annotate all of the detected variants.

2.4. Single Nucleotide Variant Detection

Since the present study did not include the matched germline samples, variants satisfying the following criteria were excluded to reduce the proportion of false-positive or possible germline variants: (i) variants with less than 5% allele frequency and <100× read depth; (ii) variants with an allele frequency >1% in the Exome Aggregation Consortium (ExAC) East Asian database [31]; (iii) all synonymous, intronic, 3′- and 5′-untranslated region (UTR) variants; and (iv) variants previously reported to be benign or likely benign with ascertained criteria in the ClinVar archive [32]. Finally, pathogenic mutations described in the COSMIC [33] database and truncation mutations of tumor suppressor gene were selected in analysis.

2.5. Immunohistochemistry (IHC) for Molecular Subtypes

A total of 22 out of 29 SeCs were classified into 5 molecular subtypes in a similar way to breast cancer, as described in our previous study (luminal 1, luminal 2, HER2, all-negative, and core basal) (Table S2) [24]. The remaining 7 cases had no available tissue for IHC study. An automated immunostainer (Ventana BenchMark XT, Tucson, AZ, USA) was used to perform the IHC for anti-HER2 (Dako, Carpinteria, CA, USA), antiandrogen receptor (anti-AR; Thermo Scientific, Carlsbad, CA, USA), anti-ER (Dako), anti-PR (Dako), anti-CK5/6 (Dako), and anti-Ki67 (Dako) antibodies. The interpretation of HER2 IHC results was based on the 2013 American Society of Clinical Oncology (ASCO)/College of American Pathologists (CAP) guidelines [34]. Finally, scores of 2+ and 3+ were considered positive for HER2. For hormone receptors (HRs), including AR, ER, and PR, it was considered positive when ≥1% of the tumor cell nuclei were stained. CK5/6 was scored when ≥50% of the tumor cells were strongly immunoreactive [24].

2.6. Bioinformatics Analysis

We used the maftools [35] package and R version 3.6.3 (http://www.r-progect.org/, accessed on 10 September 2021) to generate an oncoplot, which depicts the top 20 mutated genes in SeCs. We expanded the altered genes of each SeC by weighting frequent recurrent mutations on String with adjustment of settings (associated gene; except textmining, high confidence (0.700), no more than 5 interactors). Then, we evaluated protein–protein interaction (PPI) and conducted gene ontology (GO) pathway analysis using String v. 11.0 [36].

2.7. Statistical Analysis

Clinicopathologic variables were correlated with the presence or the frequency of genomic alterations in SeC-Os and SeC-EOs. The differences were analyzed by using Pearson’s chi-square test, Fisher’s exact test, or the Mann–Whitney U test. Progression-free survival (PFS) of SeC patients was analyzed by Kaplan–Meier survival curves. All the statistical analyses were performed using the SPSS (v22.0) software package (SPSS, Chicago, IL, USA). All p-values were two-sided, and a p-value of less than 0.05 was considered statistically significant.

3. Results

3.1. Patient Characteristics

The clinical profile of all patients is summarized in the Table 1. Among 26 patients, there were 18 patients with SeC-O and 8 patients with SeC-EO. The median age of the patients with SeC-O and SeC-EO was 64 years (range, 41–84 years) and 55 years (range, 42–70 years), respectively. A total of 15 (53.6%) and 4 (14.3%) SeC-Os were categorized as pathological T (pT) categories 1–2 and 3–4, respectively, while all 9 (100.0%) SeC-EOs were staged as pT1–2. Disease progression was noted in 8 (30.8%) patients with SeC-Os and 1 (12.5%) patient with SeC-EOs. None of the patients died during the follow-up period (median, 29 months; range, 1–132 months).

3.2. Genetic Alterations in SeC-Os and SeC-EOs

Targeted sequencing in 29 SeCs revealed 129 pathogenic mutations in 34 different genes (Table 2 and Figure 1). On average, 4.44 pathogenic mutations per tumor were observed (range, 0–15). The number of pathogenic mutations did not differ between the SeC-Os and SeC-EOs (p = 0.864).
In SeC-Os, the highest number of pathogenic mutations was detected in TP53, which harbored a total of 29 mutations in 13 cases (65.0%). Among 29 mutations, 21 were missense mutations, and 7 were nonsense mutations. The remaining 1 alteration occurred in the spliced acceptor site (c.920-1G > A). A total of 6 mutations in PIK3CA were observed in 4 cases (20.0%). SMARCA4, ATM, FBXW7, and RB1 were mutated in 3 cases (15.0%). In SeC-EOs, the most commonly mutated gene was NOTCH1, of which a total of 6 alterations were detected in 5 tumors (55.5%). There were 6 missense mutations, 3 nonsense mutations, and 1 frameshift mutation. In 3 (33.3%) tumors, a total of 8 TP53 alterations were observed, of which 7 were missense mutations, and the remaining 1 was a frameshift mutation. For BRCA1, MSH2, and SMARCA4, a total of 2 cases (22.2%) showed pathogenic variants, respectively (Table 2). Mutations in other potentially targetable genes, such as BRAF, ERBB2, and EGFR, were also detected (Table 2). Mutations in BRAF were observed in 2 (10.0%) SeC-O and 1 (11.1%) SeC-EO cases, all of which were p.N581S. All 3 cases harbored mutations in either ERBB2 (p.A775_G776insYVMA) or EGFR (p.L747_E749del). In addition, HRAS alterations were found in 1 (5.0%) SeC-O and 1 (11.1%) SeC-EO, and KMT2A alteration was observed in 2 (10.0%) SeC-O and 1 (11.1%) SeC-EO cases.
Compared with SeC-Os, mutations in NOTCH1 were significantly more enriched in SeC-EOs (p = 0.001). Although TP53 was more frequently mutated in SeC-Os than in SeC-EOs, and RB1 was mutated only in SeC-O cases, the difference was not statistically significant (p = 0.532). Likewise, the mutations in MSH2 were only observed in SeC-EOs, but the difference was not significant (p = 0.234).
Mutations in RUNX1 were observed in 2 SeC-Os, and mutations in ATM were identified in 3 SeC-Os and 1 SeC-EO. Although the frequency of mutations in these 2 genes were low, they were associated with the development of distant metastases during the follow-up period (p = 0.046, p = 0.028). However, cases harboring mutations in RUNX1 and ATM were not associated with shorter PFS in Kaplan–Meier analysis (both p > 0.05) (Figure S1). In SeC-Os, cases harboring MSH6 alterations showed significantly shorter PFS (p < 0.001), whereas in SeC-EOs, cases with BRCA1 mutations presented inferior PFS (p = 0.046) (Figure 2).

3.3. PPI in SeCs

Analysis of PPI based on genes with recurrent mutations and associated genes suggested that the TP3 pathway, SWItch/Sucrose Non-Fermentable (SWI/SNF) pathway, and MMR pathway were mainly responsible for the pathogenesis of both SeC-O and SeC-EO (Figure 3). Additionally, the PI3K/AKT/mTOR pathway (PIK3CA) and Notch signaling pathway (NOTCH1) were considered other tumorigenic pathways in SeC-O and Sec-EO, respectively (Figure 3). In GO analysis, mutated and associated genes were enriched in androgen receptor binding or signaling pathway, and p53 binding pathway in both SeC-O and SeC-EO (Table S3 and Figure S2B and Figure S3A,B). In SeC-O, genes associated with ERBB2 signaling pathway were enriched (Table S3 and Figure S2A, whereas genes associated with cellular response to ultraviolet (UV) were enriched in SeC-EO, suggesting a possible role of UV in SeC-EO tumorigenesis (Table S3 and Figure S3A).

3.4. Genomic Alterations according to the Molecular Subtypes

The molecular subtypes were determined in a total of 22 cases with available tissue for IHC and FISH analysis. In SeC-Os, 6 (42.9%), 2 (14.3%), 2 (14.3%), 4(28.6%), and 0 (0.0%) cases were classified as luminal 1, luminal 2, HER2, all-negative, and core basal type (Table 1). In SeC-EOs, 1 (14.3%), 1 (14.3%), 2 (28.6%), 1 (14.3%), and 3 (24.9%) cases were categorized as luminal 1, luminal 2, HER2, all-negative, and core basal type (Table 1). The genes that harbored pathogenic mutations in more than 2 cases are summarized according to the molecular subtypes in Table 3. TP53 was the most frequently mutated gene across the 5 subtypes. FBXW7 mutations were significantly enriched in the luminal 1 subtype compared with the other subtypes (p = 0.006). Likewise, mutations in NOTCH1 were observed only in the core basal subtype (p = 0.003). Although mutations in TP53 and SMARCA4 were more common in the luminal 1 subtype, the differences were not significant (p > 0.05). Mutations involving ERBB2 were detected only in the luminal 1 subtype.

3.5. Alterations in MMR Genes and Lynch Syndrome

One patient clinically diagnosed with Muir–Torre syndrome (MTS) was included in this study. This patient had a microsatellite instability (MSI)-high colon cancer, followed by early gastric carcinoma and SeC-O 2 years later, and then was diagnosed with metachronous SeC-EOs. The NGS study performed in 2 SeCs revealed a common mutation of MLH1 (p.M587Hfs*6) with a variant allele frequency (VAF) of 45.5% in SeC-O and 60.16% in SeC-EO (Table 4). Considering the VAF of the MLH1 mutation, this mutation was highly suggestive of a germline variant. Interestingly, his SeC-O and SeC-EO had no overlapping pathogenic variants other than the MLH1 mutation (Table 4).
In addition to the patient described above, 5 other patients had alterations in MMR genes (Table S4). The number of patients with mutations in MLH1, MSH2, and MSH6 was 1, 1, and 2, respectively. Mutations in both MSH2 and MSH6 were detected in the remaining 1 patient. Only 1 case revealed a possible germline variant in MSH6 (p.F1088Lfs*5) with a VAF of 47.66% but no other accompanying internal malignancies conforming to MTS. Except for that case, the VAFs of the other mutations were not suggestive of germline mutations. Alterations in MMR genes were associated with mutations involving KMT2A (p = 0.010) and NF1 (p = 0.052) with borderline significance. The number of pathogenic mutations was higher in tumors with MMR gene alterations than in others with borderline significance (p = 0.053).

3.6. Comparison among Primary and Locally Recurrent or Metastatic Lesions

There were 3 cases in which the mutation profile of the primary site and recurrent or metastatic tumors could be compared (Figure 4). Among them, 2 cases revealed common pathogenic mutations between the primary sites and others. Samples 4 (the original site) and 14 (the metastatic lesion) shared mutations involving PIK3CA and TP53. Likewise, in sample 7 and its recurrent lesion (sample 15), the same BRCA1 nonsense mutation was detected (Figure 4). Notably, the number of pathogenic mutations decreased in the recurrent lesions (samples 14 and 15) compared with the primary lesions (samples 4 and 7). However, the tumor population harboring TP53 p.Q192* and BRCA1 p.R1443* mutations increased in the recurrent lesions (Figure 4). On the other hand, in 1 patient out of 3, there was no common pathogenic mutation shared by both primary (sample 6) and recurrent tumors (sample 10), with the exception of the mutation in FLT4 (p.E1336K), which is a variant of undetermined significance.

4. Discussion

In the current study, we comparatively analyzed the mutational profile of SeC-Os and Sec-EOs using targeted NGS. Our results suggest that SeC-Os and SeC-EOs are related but distinct entities with different molecular pathogeneses. SeC showed frequent alterations involving TP53, NOTCH1, PIK3CA, and SMARCA4, whereas SeC-Os revealed mutations in TP53, PIK3CA, FBXW7, and RB1. Notably, NOTCH1 mutations were exclusively found in SeC-EOs. Given that most of these NOTCH1 mutations are related to a loss of function, NOTCH1 may play a role as tumor suppressor in SeC-EOs. On the other hand, FBXW7, which is known to negatively regulate NOTCH [37] and leads to activation of the Notch signal when mutated, was more commonly altered in SeC-Os. Therefore, direct or indirect activation of the Notch signaling pathway appears to be involved in the tumorigenesis of both SeC-Os and SeC-EOs. However, genes that are related to the PI3K/AKT/mTOR pathway, such as PIK3CA [38], FBXW7 [39], PTEN, and TSC1 [40], were more frequently mutated in SeC-Os. PPI analysis suggested that the TP53 pathway, SWI/SNF pathway, and MMR pathway were the main pathogenic pathways in SeCs.
Our results are in line with the previous studies, which have also demonstrated frequent activation mutations in the PI3K/AKT/mTOR pathway in SeCs [20,21,22]. The PI3K/AKT/mTOR pathway plays an important role in controlling cell growth and proliferation [41]. Activation of this pathway contributes to cell proliferation by stimulating downstream targets [41]. This phenomenon is also well known in other tumors, including breast [42], ovarian [43], and urothelial cancers [44]. North et al. also suggested a distinct pathogenesis for SeC-Os and SeC-EOs; mutations in TP53, NOTCH1, NOTCH2, ZNF750, and RREB1 are specific to SeCs, whereas truncation mutations in NOTCH1 are concentrated in SeC-EOs [20]. Notably, SeC-Os were enriched with mutations in ZNF750, while SeC-EOs showed a high mutation burden associated with MSI and UV damage [20]. Since ZNF750 is known as a transcription factor involved in keratinocyte differentiation [45,46], their findings provide support for the theory that the pathogenesis of SeC is closely related to epidermal carcinoma. However, ZNF750 was not included in our gene panel, and there are other studies that have published some controversial results [22]. Except for ZNF750, the mutational profile of SeCs in our cohort was similar to that of cutaneous squamous cell carcinoma (SqCC), in which TP53 and NOCTH1/2 are characteristically mutated [47,48]. Basal cell carcinoma (BCC) is a frequent mimicker of SeC because of its similar morphologic characteristics and cellular origin (folliculosebaceous germinative cells). However, BCC is known to show a different mutation pattern from SeC; sonic hedgehog pathway genes, such as PTCH1, SMO, or SUFU, were the mainly altered genes in BCC in addition to TP53, unlike SeCs [49]. These results suggest that the molecular characteristics of SeCs are more akin to cutaneous SqCCs than to BCCs.
Our group previously proposed five molecular subtypes of SeCs according to the HRs, HER2, and CK5/6, in a similar way to breast cancer: luminal 1, luminal 2, HER2, all-negative, and core basal subtypes [24]. While SeC-EOs were more frequently classified as core basal subtypes, SeC-Os prevailed among the other four subtypes. When comparing mutational profiles, FBXW7 mutations were significantly enriched in the luminal 1 subtype, and mutations in NOTCH1 were observed only in the core basal subtype (p = 0.003). Although SeCs and breast cancer are similar in that PIK3CA and TP53 are the most common mutations, there was little overlap between the two tumors in terms of the types of alterations according to each molecular subtype [50,51]. For example, in general, basal-like breast carcinoma has been known to harbor mutations involving BRCA1, CCNE1, AKT3, or MYC [50,51]; however, these genes were largely wild type, and NOTCH1 mutations were relatively common in our SeC-EOs. Wang et al. demonstrated that approximately 13% of triple-negative breast cancer cases carry NOTCH1 mutations, which are centered in the PEST domain, and ultimately lead to the activation of the Notch signaling pathway [52]. In contrast, none of the six mutations detected in SeC-EOs occurred in the PEST domain, and most of the NOTCH1 mutations were predicted to damage protein function. Although SeCs can be classified into five immunohistochemical subtypes, their mutational profile did not completely match with these subtypes adopted from breast cancer.
Mutations in MMR genes were found in a total of six patients, including one with classical MTS harboring pathogenic MLH1 mutations in both SeC-O and SeC-EO. Interestingly, SeC-EO did not share other pathogenic mutations with SeC-O, suggesting multiple individual primary tumors rather than metastasis. In the other five patients with MMR defects, MLH1, MSH2, and MSH6 mutations were observed mostly in low VAF, suggesting a somatic origin. Since various types of mutations were found in these cases in addition to MMR defects, it is difficult to determine which of them is the driving mutation. Although we could not confirm the MSI status of these tumors due to tissue unavailability, it is possible that DNA instability induced by MMR defects resulted in numerous mutations. Therefore, we suggest that SeC should also be screened for MSI as for colorectal carcinoma and endometrial carcinoma. Previous studies have also demonstrated that somatic alterations in MMR genes occur especially in SeC-EOs, which show an MSI-high phenotype [20,21].
We had three pairs of primary and matched recurrent or metastatic lesions. In two patients, some of the pathogenic variants of the primary lesions remained in the recurrent lesions. The VAFs of certain variants in the primary lesions were significantly increased (TP53 p.Q192* in sample 4 and BRCA1 p.R1443* in sample 7) or decreased (TP53 p.R342P in sample 4) in recurrent lesions, suggesting an expansion of a specific tumoral subclone that conveys a survival advantage. Notably, samples 6 and 10 did not share any pathogenic variants. Instead, a variant of undetermined significance (FLT4 p.E1336K) was detected in both lesions. FLT4, also known as VEGFR3, is an oncogene that mediates cell proliferation, survival, and invasion [53]. Although somatic variants of FLT4 in solid malignancies are not commonly reported, mutations in colorectal cancer [54] and non-small cell lung cancer [55] have been documented. FLT4 p.E1336K has not previously been reported; however, the possibility of an oncogenic role can be considered under the above circumstances.
This study identified some potentially targetable alterations in SeCs. First, PIK3CA mutations account for 20% (4/20) of SeC-Os, and it can be targeted by PI3K inhibitors similar to hormone receptor-positive breast cancer [56]. Notably, alterations in homologous recombination deficiency (HRD)-related genes, such as BRCA1, BRCA2, ATM, and KMT2A, were detected in both SeC-Os and SeC-EOs, suggesting a possible role of poly (ADP-ribose) polymerase inhibitors [57,58]. Second, three patients showed exon 20 insertion mutations in ERBB2 (p.A775_G776insYVMA), for which a response to pan-HER kinase inhibitors has been reported in various types of cancer [59]. In addition, activating mutations in other kinases, such as EGFR (exon 19 deletion) and BRAF (p.N581S), were detected, which can also provide grounds for potential targeted therapy [60,61]. Finally, alterations in MMR genes can be an important target for immunotherapy, although further validation of MSI status is required in these tumors [62,63]. Although NOTCH1 mutations were prominent in SeC-EOs, most alterations were associated with a loss of protein function; therefore, NOTCH1 inhibitors do not seem promising based on our findings.
One of the novel findings of the present study was that mutations in RUNX1 (2/29, 6.9%) and ATM (4/29, 13.8%) were associated with the development of distant metastasis during the follow-up period. However, cases with mutations in RUNX1 or ATM were not associated with shorter PFS. This might be due to the fact that local recurrences also contribute to disease progression, especially in SeC-Os. RUNX1 is frequently altered in hematopoietic malignancies [64,65]. Albeit with a lower frequency, mutations in RUNX1 are also found in solid malignancies, such as breast cancer [66,67,68,69]. Notably, researchers demonstrated that RUNX1 knockdown induces the activation of transforming growth factor β and WNT signaling pathways, thereby promoting the epithelial–mesenchymal transition (EMT) process [70,71]. Mutations in ATM have been reported in various solid tumors, including endometrial cancer, bladder cancer, and colorectal cancer [66,67]. Loss of ATM protein expression was associated with poor overall survival in breast cancer [72], cervical cancer [73], and colorectal cancer [74], although the prognostic implications of ATM mutation are generally unknown to date. ATM signaling is known to induce tumor progression in a nuclear factor kappa B-dependent manner that promotes EMT [75]. Moreover, ATM signaling upregulated the alphavbeta3 integrin pathway, which can trigger invasive and metastatic potential [76]. As the EMT process is critical in the invasion and development of distant metastasis [77], mutations in RUNX1 and ATM might contribute to developing metastasis. Our findings suggest that mutations in RUNX1 and ATM can serve as prognostic markers in SeCs. In addition, MSH6 and BRCA1 mutations were associated with inferior PFS in SeC-Os and SeC-EOs, respectively. For MMR gene defects and BRCA mutations, contradictory results have been reported [78,79,80,81]. Further validation is required to determine the prognostic implications of MSH6 and BRCA1 mutations in SeC.

5. Conclusions

This is the first study performing NGS in Korean SeC-O and SeC-EO patients. Sequencing revealed a distinct pattern of genetic alterations in SeC-Os and SeC-EOs. Additionally, various potentially targetable alterations were detected, including PIK3CA, HRD-related genes, and MMR genes, which could lead to clinical benefit in patients with advanced SeCs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13194810/s1, Figure S1: Kaplan–Meier analysis according to RUNX1 and ATM mutation status, Figure S2: Gene ontology analysis in ocular sebaceous carcinoma, Figure S3: Gene ontology analysis in extraocular sebaceous carcinoma, Table S1: Information of targeted next-generation sequencing test (Axen Cancer Panel 2), Table S2: Determination of molecular subtypes of sebaceous carcinoma, Table S3: Gene ontology analysis of each ocular and extraocular sebaceous carcinoma, Table S4: Cases showing alterations in mismatch repair genes.

Author Contributions

Conceptualization, H.Y.N., J.-H.C., and J.E.K.; methodology, H.Y.N., J.H.P., S.A.S., S.L., H.L., and H.C. (Heejoon Chae); validation, H.Y.N., J.H.P., S.L., H.L., and H.C. (Heejoon Chae); formal analysis, H.Y.N., J.H.P., S.L., H.L., and H.C. (Heejoon Chae); investigation, H.Y.N., J.H.P., and S.A.S.; resources, S.A.S., H.C. (HoKyung Choung), N.K., J.-H.C., and J.E.K.; data curation, H.Y.N., J.H.P., S.A.S., S.L., H.L., and H.C. (Heejoon Chae); writing—original draft preparation, H.Y.N. and J.H.P.; writing—review and editing, H.Y.N., J.H.P., S.A.S., H.C. (HoKyung Choung), N.K., J.-H.C., and J.E.K.; visualization, H.Y.N., J.H.P., S.L., H.L., and H.C. (Heejoon Chae); supervision, H.C. (HoKyung Choung), N.K., J.-H.C., and J.E.K.; project administration, H.C. (HoKyung Choung), N.K., J.-H.C., and J.E.K.; funding acquisition, J.-H.C. and J.E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Seoul National University Bundang Hospital Research Fund (grant number: 13-2017-012) and National Research Foundation of Korea (NRF) Grant funded by the Korean government (MSIT) (grant number: 2017R1A5A1015626).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (IRB) at SMG-SNU Boramae Medical Center (IRB No. 16-2017-17-031).

Informed Consent Statement

Patient consent was waived because this study was a retrospective study using FFPE blocks only.

Data Availability Statement

The main research data supporting the results of this study are included in Table 1 and Table 2 and Figure 1 and Figure 3. Other data can be made available upon reasonable request from the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dasgupta, T.; Wilson, L.D.; Yu, J. A retrospective review of 1349 cases of sebaceous carcinoma. Cancer 2009, 115, 158–165. [Google Scholar] [CrossRef]
  2. Kumar, R. Clinicopathologic study of malignant eyelid tumours. Clin. Exp. Optom. 2010, 93, 224–227. [Google Scholar] [CrossRef] [PubMed]
  3. Ni, C.; Searl, S.S.; Kuo, P.K.; Chu, F.R.; Chong, C.S.; Albert, D.M. Sebaceous cell carcinomas of the ocular adnexa. Int. Ophthalmol. Clin. 1982, 22, 23–61. [Google Scholar] [CrossRef] [PubMed]
  4. Requena, L.; Sangueza, O. Cutaneous Adenexal Neoplasms, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2017; pp. 755–764. [Google Scholar]
  5. Owen, J.L.; Kibbi, N.; Worley, B.; Kelm, R.C.; Wang, J.V.; Barker, C.A.; Behshad, R.; Bichakjian, C.K.; Bolotin, D.; Bordeaux, J.S.; et al. Sebaceous carcinoma: Evidence-based clinical practice guidelines. Lancet Oncol. 2019, 20, e699–e714. [Google Scholar] [CrossRef]
  6. Spencer, J.M.; Nossa, R.; Tse, D.T.; Sequeira, M. Sebaceous carcinoma of the eyelid treated with Mohs micrographic surgery. J. Am. Acad. Dermatol. 2001, 44, 1004–1009. [Google Scholar] [CrossRef] [PubMed]
  7. Snow, S.N.; Larson, P.O.; Lucarelli, M.J.; Lemke, B.N.; Madjar, D.D. Sebaceous carcinoma of the eyelids treated by mohs micrographic surgery: Report of nine cases with review of the literature. Dermatol. Surg. 2002, 28, 623–631. [Google Scholar] [CrossRef] [PubMed]
  8. Esmaeli, B.; Nasser, Q.J.; Cruz, H.; Fellman, M.; Warneke, C.L.; Ivan, D. American Joint Committee on Cancer T category for eyelid sebaceous carcinoma correlates with nodal metastasis and survival. Ophthalmology 2012, 119, 1078–1082. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Rao, N.A.; Hidayat, A.A.; McLean, L.C.I.W.; Zimmerman, L.E. Sebaceous carcinomas of the ocular adnexa: A clinicopathologic study of 104 cases, with five-year follow-up data. Hum. Pathol. 1982, 13, 113–122. [Google Scholar] [CrossRef]
  10. Shields, J.A.; Demirci, H.; Marr, B.P.; Eagle, R.C., Jr.; Shields, C.L. Sebaceous carcinoma of the eyelids: Personal experience with 60 cases. Ophthalmology 2004, 111, 2151–2157. [Google Scholar] [CrossRef] [PubMed]
  11. Kyllo, R.L.; Brady, K.L.; Hurst, E.A. Sebaceous carcinoma: Review of the literature. Dermatol. Surg. 2015, 41, 1–15. [Google Scholar] [CrossRef]
  12. Panda, B.B.; Parija, S.; Pujahari, S.; Mallick, J. Sebaceous gland carcinoma of eyelid—A tarnished masquerade. J. Clin. Diagn. Res. 2016, 10, ND03–ND05. [Google Scholar] [CrossRef] [PubMed]
  13. Kibbi, N.; Worley, B.; Owen, J.L.; Kelm, R.C.; Bichakjian, C.K.; Chandra, S.; Demirci, H.; Kim, J.; Nehal, K.S.; Thomas, J.R.; et al. Sebaceous carcinoma: Controversies and their evidence for clinical practice. Arch. Dermatol. Res. 2019, 312, 25–31. [Google Scholar] [CrossRef] [PubMed]
  14. Kumar, V.; Xu, Y. Unusual presentation of metastatic sebaceous carcinoma and its response to chemotherapy: Is genotyping a right answer for guiding chemotherapy in rare tumours? Curr. Oncol. 2015, 22, 316–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kaliki, S.; Ayyar, A.; Nair, A.G.; Mishra, D.K.; Reddy, V.A.P.; Naik, M.N. Neoadjuvant systemic chemotherapy in the management of extensive eyelid sebaceous gland carcinoma: A study of 10 cases. Ophthalmic Plast. Reconstr. Surg. 2016, 32, 35–39. [Google Scholar] [CrossRef]
  16. Joshi, P.; Joshi, A.; Norohna, V.; Prabhash, K.; Kane, S.; D’Cruz, A.K. Sebaceous carcinoma and systemic chemotherapy. Indian J. Med. Paediatr. Oncol. 2012, 33, 239–241. [Google Scholar] [CrossRef] [Green Version]
  17. Yu, C.; Shahsavari, M.; Stevens, G.; Liskanich, R.; Horowitz, D. Isotretinoin as monotherapy for sebaceous hyperplasia. J. Drugs Dermatol. 2010, 9, 699–701. [Google Scholar]
  18. Mulay, K.; Shah, S.J.; Aggarwal, E.; White, V.A.; Honavar, S.G. Periocular sebaceous gland carcinoma: Do androgen receptor (NR3C4) and nuclear survivin (BIRC5) have a prognostic significance? Acta Ophthalmol. 2014, 92, e681–e687. [Google Scholar] [CrossRef] [PubMed]
  19. Ivan, D.; Prieto, V.G.; Esmaeli, B.; Wistuba, I.I.; Tang, X.; Lazar, A.J.F. Epidermal growth factor receptor (EGFR) expression in periocular and extraocular sebaceous carcinoma. J. Cutan. Pathol. 2010, 37, 231–238. [Google Scholar] [CrossRef]
  20. North, J.P.; Golovato, J.; Vaske, C.J.; Sanborn, J.Z.; Nguyen, A.; Wu, W.; Goode, B.; Stevers, M.; McMullen, K.; White, B.E.P.; et al. Cell of origin and mutation pattern define three clinically distinct classes of sebaceous carcinoma. Nat. Commun. 2018, 9, 1894. [Google Scholar] [CrossRef]
  21. Tetzlaff, M.T.; Singh, R.R.; Seviour, E.G.; Curry, J.L.; Hudgens, C.W.; Bell, D.; Wimmer, D.A.; Ning, J.; Czerniak, B.A.; Zhang, L.; et al. Next-generation sequencing identifies high frequency of mutations in potentially clinically actionable genes in sebaceous carcinoma. J. Pathol. 2016, 240, 84–95. [Google Scholar] [CrossRef]
  22. Xu, S.; Moss, T.J.; Rubin, M.L.; Ning, J.; Eterovic, K.; Yu, H.; Jia, R.; Fan, X.; Tetzlaff, M.T.; Esmaeli, B. Whole-exome sequencing for ocular adnexal sebaceous carcinoma suggests PCDH15 as a novel mutation associated with metastasis. Mod. Pathol. 2020, 33, 1256–1263. [Google Scholar] [CrossRef]
  23. Oftedal, O.T. The mammary gland and its origin during synapsid evolution. J. Mammary Gland. Biol. Neoplasia 2002, 7, 225–252. [Google Scholar] [CrossRef]
  24. Na, H.Y.; Choe, J.-Y.; Shin, S.A.; Choung, H.-K.; Oh, S.; Chung, J.-H.; Park, M.; Kim, J.E. Proposal of a provisional classification of sebaceous carcinoma based on hormone receptor expression and HER2 status. Am. J. Surg. Pathol. 2016, 40, 1622–1630. [Google Scholar] [CrossRef] [PubMed]
  25. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  26. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 2013, arXiv:1303.3997v2. [Google Scholar]
  27. Li, H.; Handsaker, R.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The sequence alignment/map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [Green Version]
  28. Cibulskis, K.; Lawrence, M.S.; Carter, S.L.; Sivachenko, A.; Jaffe, D.B.; Sougnez, C.; Gabriel, S.B.; Meyerson, M.L.; Lander, E.S.; Getz, G. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat. Biotechnol. 2013, 31, 213–219. [Google Scholar] [CrossRef]
  29. Cingolani, P.; Patel, V.M.; Coon, M.; Nguyen, T.; Land, S.J.; Ruden, D.M.; Lu, X. Using Drosophila melanogaster as a model for genotoxic chemical mutational studies with a new program, SnpSift. Front. Genet. 2012, 3, 35. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, X.; Wu, C.; Li, C.; Boerwinkle, E. dbNSFP v3.0: A one-stop database of functional predictions and annotations for human nonsynonymous and splice-site SNVs. Hum. Mutat. 2015, 37, 235–241. [Google Scholar] [CrossRef] [Green Version]
  31. Lek, M.; Karczewski, K.J.; Minikel, E.V.; Samocha, K.; Banks, E.; Fennell, T.; O’Donnell-Luria, A.; Ware, J.; Hill, A.J.; Cummings, B.B.; et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016, 536, 285–291. [Google Scholar] [CrossRef] [Green Version]
  32. Landrum, M.J.; Lee, J.M.; Benson, M.; Brown, G.R.; Chao, C.; Chitipiralla, S.; Gu, B.; Hart, J.; Hoffman, D.; Jang, W.; et al. ClinVar: Improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2017, 46, D1062–D1067. [Google Scholar] [CrossRef] [Green Version]
  33. Tate, J.G.; Bamford, S.; Jubb, H.C.; Sondka, Z.; Beare, D.M.; Bindal, N.; Boutselakis, H.; Cole, C.G.; Creatore, C.; Dawson, E.; et al. COSMIC: The catalogue of somatic mutations in cancer. Nucleic Acids Res. 2018, 47, D941–D947. [Google Scholar] [CrossRef] [Green Version]
  34. Wolff, A.; Hammond, M.E.H.; Hicks, D.G.; Dowsett, M.; McShane, L.M.; Allison, K.H.; Allred, D.C.; Bartlett, J.M.; Bilous, M.; Fitzgibbons, P.; et al. Recommendations for human epidermal growth factor receptor 2 testing in breast cancer: American Society of Clinical Oncology/College of American Pathologists clinical practice guideline update. J. Clin. Oncol. 2013, 31, 3997–4013. [Google Scholar] [CrossRef] [PubMed]
  35. Mayakonda, A.; Koeffler, H.P. Maftools: Efficient analysis, visualization and summarization of MAF files from large-scale cohort based cancer studies. bioRxiv 2016. [Google Scholar] [CrossRef] [Green Version]
  36. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Mo, J.-S.; Kim, M.-Y.; Han, S.-O.; Kim, I.-S.; Ann, E.-J.; Lee, K.S.; Seo, M.-S.; Kim, J.-Y.; Lee, S.-C.; Park, J.-W.; et al. Integrin-linked kinase controls Notch1 signaling by down-regulation of protein stability through Fbw7 ubiquitin ligase. Mol. Cell. Biol. 2007, 27, 5565–5574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Engelman, J.A.; Luo, J.; Cantley, L. The evolution of phosphatidylinositol 3-kinases as regulators of growth and metabolism. Nat. Rev. Genet. 2006, 7, 606–619. [Google Scholar] [CrossRef] [PubMed]
  39. Mao, J.-H.; Kim, I.-J.; Wu, D.; Climent, J.; Kang, H.C.; DelRosario, R.; Balmain, A. FBXW7 targets mTOR for degradation and cooperates with PTEN in tumor suppression. Science 2008, 321, 1499–1502. [Google Scholar] [CrossRef] [Green Version]
  40. Nathan, N.; Keppler-Noreuil, K.M.; Biesecker, L.G.; Moss, J.; Darling, T.N. Mosaic disorders of the PI3K/PTEN/AKT/TSC/mTORC1 signaling pathway. Dermatol. Clin. 2016, 35, 51–60. [Google Scholar] [CrossRef] [Green Version]
  41. Laplante, M.; Sabatini, D.M. mTOR signaling in growth control and disease. Cell 2012, 149, 274–293. [Google Scholar] [CrossRef] [Green Version]
  42. Raphael, J.; Desautels, D.; Pritchard, K.I.; Petkova, E.; Shah, P.S. Phosphoinositide 3-kinase inhibitors in advanced breast cancer: A systematic review and meta-analysis. Eur. J. Cancer 2018, 91, 38–46. [Google Scholar] [CrossRef]
  43. Aziz, A.U.R.; Farid, S.; Qin, K.; Wang, H.; Liu, B. PIM kinases and their relevance to the PI3K/AKT/mTOR pathway in the regulation of ovarian cancer. Biomolecules 2018, 8, 7. [Google Scholar] [CrossRef] [Green Version]
  44. Serra, V.; Markman, B.; Scaltriti, M.; Eichhorn, P.; Valero, V.; Guzman, M.; Botero, M.L.; Llonch, E.; Atzori, F.; Di Cosimo, S.; et al. NVP-BEZ235, a dual PI3K/mTOR inhibitor, prevents PI3K signaling and inhibits the growth of cancer cells with activating PI3K mutations. Cancer Res. 2008, 68, 8022–8030. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Boxer, L.; Barajas, B.; Tao, S.; Zhang, J.; Khavari, P.A. ZNF750 interacts with KLF4 and RCOR1, KDM1A, and CTBP1/2 chromatin regulators to repress epidermal progenitor genes and induce differentiation genes. Genes Dev. 2014, 28, 2013–2026. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Sen, G.L.; Boxer, L.; Webster, D.; Bussat, R.T.; Qu, K.; Zarnegar, B.J.; Johnston, D.; Siprashvili, Z.; Khavari, P.A. ZNF750 is a p63 target gene that induces KLF4 to drive terminal epidermal differentiation. Dev. Cell 2012, 22, 669–677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Pickering, C.; Zhou, J.H.; Lee, J.J.; Drummond, J.A.; Peng, S.A.; Saade, R.E.; Tsai, K.; Curry, J.L.; Tetzlaff, M.T.; Lai, S.; et al. Mutational landscape of aggressive cutaneous squamous cell carcinoma. Clin. Cancer Res. 2014, 20, 6582–6592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. South, A.P.; Purdie, K.J.; Watt, S.A.; Haldenby, S.; Breems, N.Y.D.; Dimon, M.; Arron, S.; Kluk, M.J.; Aster, J.C.; McHugh, A.; et al. NOTCH1 mutations occur early during cutaneous squamous cell carcinogenesis. J. Investig. Dermatol. 2014, 134, 2630–2638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Bonilla, X.; Parmentier, L.; King, B.; Bezrukov, F.; Kaya, G.; Zoete, V.; Seplyarskiy, V.B.; Sharpe, H.; 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]
  50. The Cancer Genome Atlas Network. Comprehensive molecular portraits of human breast tumours. Nature 2012, 490, 61–70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Pereira, B.; Chin, S.-F.; Rueda, O.M.; Vollan, H.-K.M.; Provenzano, E.; Bardwell, H.A.; Pugh, M.; Jones, L.; Russell, R.; Sammut, S.-J.; et al. The somatic mutation profiles of 2433 breast cancers refines their genomic and transcriptomic landscapes. Nat. Commun. 2016, 7, 11479. [Google Scholar] [CrossRef] [Green Version]
  52. Wang, K.; Zhang, Q.; Li, D.; Ching, K.; Zhang, C.; Zheng, X.; Ozeck, M.; Shi, S.; Li, X.; Wang, H.; et al. PEST domain mutations in Notch receptors comprise an oncogenic driver segment in triple-negative breast cancer sensitive to a gamma-secretase inhibitor. Clin. Cancer Res. 2015, 21, 1487–1496. [Google Scholar] [CrossRef] [Green Version]
  53. Su, J.-L.; Yang, P.-C.; Shih, J.-Y.; Yang, C.-Y.; Wei, L.-H.; Hsieh, C.-Y.; Chou, C.-H.; Jeng, Y.-M.; Wang, M.-Y.; Chang, K.-J.; et al. The VEGF-C/Flt-4 axis promotes invasion and metastasis of cancer cells. Cancer Cell 2006, 9, 209–223. [Google Scholar] [CrossRef] [Green Version]
  54. Oliveira, D.M.; Grillone, K.; Mignogna, C.; De Falco, V.; Laudanna, C.; Biamonte, F.; Locane, R.; Corcione, F.; Fabozzi, M.; Sacco, R.; et al. Next-generation sequencing analysis of receptor-type tyrosine kinase genes in surgically resected colon cancer: Identification of gain-of-function mutations in the RET proto-oncogene. J. Exp. Clin. Cancer Res. 2018, 37, 84. [Google Scholar] [CrossRef] [Green Version]
  55. Ding, L.; Getz, G.; Wheeler, D.A.; Mardis, E.R.; McLellan, M.D.; Cibulskis, K.; Sougnez, C.; Greulich, H.; Muzny, D.M.; Morgan, M.B.; et al. Somatic mutations affect key pathways in lung adenocarcinoma. Nature 2008, 455, 1069–1075. [Google Scholar] [CrossRef]
  56. Andre, F.; Ciruelos, E.; Rubovszky, G.; Campone, M.; Loibl, S.; Rugo, H.S.; Iwata, H.; Conte, P.; Mayer, I.A.; Kaufman, B.; et al. Alpelisib for PIK3CA-mutated, hormone receptor-positive advanced breast cancer. N. Engl. J. Med. 2019, 380, 1929–1940. [Google Scholar] [CrossRef]
  57. Audeh, M.W.; Carmichael, J.; Penson, R.T.; Friedlander, M.; Powell, B.; Bell-McGuinn, K.M.; Scott, C.; Weitzel, J.N.; Oaknin, A.; Loman, N.; et al. Oral poly(ADP-ribose) polymerase inhibitor olaparib in patients with BRCA1 or BRCA2 mutations and recurrent ovarian cancer: A proof-of-concept trial. Lancet 2010, 376, 245–251. [Google Scholar] [CrossRef]
  58. Mateo, J.; Carreira, S.; Sandhu, S.; Miranda, S.; Mossop, H.; Perez-Lopez, R.; Rodrigues, D.N.; Robinson, D.; Omlin, A.; Tunariu, N.; et al. DNA-repair defects and olaparib in metastatic prostate cancer. N. Engl. J. Med. 2015, 373, 1697–1708. [Google Scholar] [CrossRef] [PubMed]
  59. Hyman, D.M.; Piha-Paul, S.; Won, H.; Rodon, J.; Saura, C.; Shapiro, G.I.; Juric, D.; Quinn, D.I.; Moreno, V.; Doger, B.; et al. HER kinase inhibition in patients with HER2- and HER3-mutant cancers. Nature 2018, 554, 189–194. [Google Scholar] [CrossRef] [PubMed]
  60. Bracht, J.W.P.; Karachaliou, N.; Bivona, T.; Lanman, R.B.; Faull, I.; Nagy, R.J.; Drozdowskyj, A.; Berenguer, J.; Fernandez-Bruno, M.; Molina-Vila, M.A.; et al. BRAF mutations classes I, II, and III in NSCLC patients included in the SLLIP trial: The need for a new pre-clinical treatment rationale. Cancers 2019, 11, 1381. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Yang, J.C.-H.; Wu, Y.-L.; Schuler, M.; Sebastian, M.; Popat, S.; Yamamoto, N.; Zhou, C.; Hu, C.-P.; O’Byrne, K.; Feng, J.; et al. Afatinib versus cisplatin-based chemotherapy for EGFR mutation-positive lung adenocarcinoma (LUX-Lung 3 and LUX-Lung 6): Analysis of overall survival data from two randomised, phase 3 trials. Lancet Oncol. 2015, 16, 141–151. [Google Scholar] [CrossRef] [Green Version]
  62. Marcus, L.; Lemery, S.J.; Keegan, P.; Pazdur, R. FDA approval summary: Pembrolizumab for the treatment of microsatellite instability-high solid tumors. Clin. Cancer Res. 2019, 25, 3753–3758. [Google Scholar] [CrossRef] [Green Version]
  63. Banerjee, N.; Hossain, F.; Wirtschafter, E.; Fathizadeh, P. Pembrolizumab in the treatment of microsatellite instability—High sebaceous carcinoma: A case report with review of the literature. JCO Precis. Oncol. 2020, 4, 61–65. [Google Scholar] [CrossRef]
  64. Mendler, J.; Maharry, K.; Radmacher, M.D.; Mrózek, K.; Becker, H.; Metzeler, K.; Schwind, S.; Whitman, S.P.; Khalife, J.; Kohlschmidt, J.; et al. RUNX1 mutations are associated with poor outcome in younger and older patients with cytogenetically normal acute myeloid leukemia and with distinct gene and MicroRNA expression signatures. J. Clin. Oncol. 2012, 30, 3109–3118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Haferlach, T.; Stengel, A.; Eckstein, S.; Perglerová, K.; Alpermann, T.; Kern, W.; Meggendorfer, M. The new provisional WHO entity ‘RUNX1 mutated AML’ shows specific genetics but no prognostic influence of dysplasia. Leukemia 2016, 30, 2109–2112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio cancer genomics portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Ellis, M.J.; Ding, L.; Shen, D.; Luo, J.; Suman, V.J.; Wallis, J.W.; Van Tine, B.A.; Hoog, J.; Goiffon, R.; Goldstein, T.C.; et al. Whole-genome analysis informs breast cancer response to aromatase inhibition. Nature 2012, 486, 353–360. [Google Scholar] [CrossRef]
  69. Banerji, S.; Cibulskis, K.; Rangel-Escareno, C.; Brown, K.; Carter, S.L.; Frederick, A.M.; Lawrence, M.S.; Sivachenko, A.Y.; Sougnez, C.; Zou, L.; et al. Sequence analysis of mutations and translocations across breast cancer subtypes. Nature 2012, 486, 405–409. [Google Scholar] [CrossRef]
  70. Wang, L.; Brugge, J.S.; Janes, K.A. Intersection of FOXO- and RUNX1-mediated gene expression programs in single breast epithelial cells during morphogenesis and tumor progression. Proc. Natl. Acad. Sci. USA 2011, 108, E803–E812. [Google Scholar] [CrossRef] [Green Version]
  71. Chimge, N.-O.; Little, G.H.; Baniwal, S.K.; Adisetiyo, H.; Xie, Y.; Zhang, T.; O’Laughlin, A.; Liu, Z.Y.; Ulrich, P.; Martin, A.; et al. RUNX1 prevents oestrogen-mediated AXIN1 suppression and beta-catenin activation in ER-positive breast cancer. Nat. Commun. 2016, 7, 10751. [Google Scholar] [CrossRef] [Green Version]
  72. Feng, X.; Li, H.; Dean, M.; E Wilson, H.; Kornaga, E.; Enwere, E.K.; Tang, P.; Paterson, A.; Lees-Miller, S.P.; Magliocco, A.M.; et al. Low ATM protein expression in malignant tumor as well as cancer-associated stroma are independent prognostic factors in a retrospective study of early-stage hormone-negative breast cancer. Breast Cancer Res. 2015, 17, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Ho, C.; Kornaga, E.; Klimowicz, A.; Enwere, E.; Dean, M.; Bebb, G.; Phan, T.; Ghatage, P.; Magliocco, A.; Lees-Miller, S.; et al. Expression of DNA damage response proteins in cervical cancer patients treated with radical chemoradiotherapy. Gynecol. Oncol. 2017, 145, 176–184. [Google Scholar] [CrossRef]
  74. Lu, Y.; Lu, Y.; Gao, J. Downregulated Ku70 and ATM associated to poor prognosis in colorectal cancer among Chinese patients. OncoTargets Ther. 2014, 7, 1955–1961. [Google Scholar] [CrossRef] [Green Version]
  75. Pazolli, E.; Alspach, E.; Milczarek, A.; Prior, J.; Piwnica-Worms, D.; Stewart, S.A. Chromatin remodeling underlies the senescence-associated secretory phenotype of tumor stromal fibroblasts that supports cancer progression. Cancer Res. 2012, 72, 2251–2261. [Google Scholar] [CrossRef] [Green Version]
  76. Stucci, S.; Tucci, M.; Passarelli, A.; Silvestris, F. Avβ3 integrin: Pathogenetic role in osteotropic tumors. Crit. Rev. Oncol. 2015, 96, 183–193. [Google Scholar] [CrossRef]
  77. Kang, Y.; Massagué, J. Epithelial-mesenchymal transitions: Twist in development and metastasis. Cell 2004, 118, 277–279. [Google Scholar] [CrossRef] [Green Version]
  78. Navarro, E.B.; López, E.V.; Quijano, Y.; Caruso, R.; Ferri, V.; Durand, H.; Cabrera, I.F.; Reques, E.D.; Ielpo, B.; Glagolieva, A.Y.; et al. Impact of BRCA1/2 gene mutations on survival of patients with pancreatic cancer: A case-series analysis. Ann. Hepatbiliary Pancreat. Surg. 2019, 23, 200–205. [Google Scholar] [CrossRef] [Green Version]
  79. Nagle, C.M.; O’Mara, T.A.; Tan, Y.; Buchanan, D.D.; Obermair, A.; Blomfield, P.; Quinn, M.A.; Webb, P.M.; Spurdle, A.B. Endometrial cancer risk and survival by tumor MMR status. J. Gynecol. Oncol. 2018, 29, e39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Kim, S.I.; Lee, M.; Kim, H.S.; Chung, H.H.; Kim, J.-W.; Park, N.H.; Song, Y.-S. Effect of BRCA mutational status on survival outcome in advanced-stage high-grade serous ovarian cancer. J. Ovarian Res. 2019, 12, 1–10. [Google Scholar] [CrossRef] [PubMed]
  81. Hou, J.-T.; Zhao, L.-N.; Zhang, D.-J.; Lv, D.-Y.; He, W.-L.; Chen, B.; Li, H.-B.; Li, P.-R.; Chen, L.-Z.; Chen, X.-L. Prognostic value of mismatch repair genes for patients with colorectal cancer: Meta-analysis. Technol. Cancer Res. Treat. 2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Cancers 13 04810 g001 550
Figure 1. Oncoplot of ocular and extraocular sebaceous carcinomas. TP53 and PIK3CA genes were frequently mutated in ocular sebaceous carcinomas, while the NOTCH1 gene was more commonly mutated in extraocular sebaceous carcinomas.
Figure 1. Oncoplot of ocular and extraocular sebaceous carcinomas. TP53 and PIK3CA genes were frequently mutated in ocular sebaceous carcinomas, while the NOTCH1 gene was more commonly mutated in extraocular sebaceous carcinomas.
Cancers 13 04810 g001
Cancers 13 04810 g002 550
Figure 2. Kaplan–Meier analysis revealed (A) MSH6 mutation, and (B) BRCA1 mutation was associated with inferior progression-free survival in ocular and extraocular sebaceous carcinomas, respectively.
Figure 2. Kaplan–Meier analysis revealed (A) MSH6 mutation, and (B) BRCA1 mutation was associated with inferior progression-free survival in ocular and extraocular sebaceous carcinomas, respectively.
Cancers 13 04810 g002
Cancers 13 04810 g003 550
Figure 3. Analysis of PPI based on genes with recurrent mutation and associated genes. (A) Ocular sebaceous carcinoma (SeC-O). (B) Extraocular sebaceous carcinoma (SeC-EO). The TP3 pathway, SWItch/Sucrose Non-Fermentable (SWI/SNF) pathway, and MMR pathway were shared by both SeC-O and SeC-EO. The PI3K/AKT/mTOR pathway (PIK3CA) and Notch signaling pathway (NOTCH1) were unique to SeC-O and Sec-EO, respectively.
Figure 3. Analysis of PPI based on genes with recurrent mutation and associated genes. (A) Ocular sebaceous carcinoma (SeC-O). (B) Extraocular sebaceous carcinoma (SeC-EO). The TP3 pathway, SWItch/Sucrose Non-Fermentable (SWI/SNF) pathway, and MMR pathway were shared by both SeC-O and SeC-EO. The PI3K/AKT/mTOR pathway (PIK3CA) and Notch signaling pathway (NOTCH1) were unique to SeC-O and Sec-EO, respectively.
Cancers 13 04810 g003
Cancers 13 04810 g004 550
Figure 4. Comparison between primary and recurrent lesions. (A) Two patients showed common pathogenic mutations in TP53, PIK3CA, and BRCA1. Sample 6 and 10 from other patient harbored a common FLT4 mutation, which is of unknown clinical significance. (B) The frequency of mutated TP53 and BRCA1 alleles increased in recurrent lesion in sample 14 and 15.
Figure 4. Comparison between primary and recurrent lesions. (A) Two patients showed common pathogenic mutations in TP53, PIK3CA, and BRCA1. Sample 6 and 10 from other patient harbored a common FLT4 mutation, which is of unknown clinical significance. (B) The frequency of mutated TP53 and BRCA1 alleles increased in recurrent lesion in sample 14 and 15.
Cancers 13 04810 g004
Table
Table 1. Clinicopathological parameters of ocular and extraocular sebaceous carcinomas.
Table 1. Clinicopathological parameters of ocular and extraocular sebaceous carcinomas.
ParametersOcular SeCExtraocular SeCp-Value
Age (mean, range, years) (n = 26)64 (41–84)55 (42–70)0.663
Gender (n = 26) 0.197
Male4 (22.2%)4 (50.0%)
Female14 (77.8%)4 (50.0%)
Lymphatic invasion (n = 29)3 (15.0%)1 (11.1%)>0.999
Neural invasion (n = 29)3 (15.0%)0 (0.0%)0.532
pT (n = 28) 0.273
pT1 and pT215 (53.6%)9 (100.0%)
pT3 and pT44 (14.3%)0 (0.0%)
pN (n = 28) >0.999
pN017 (60.7%)9 (100.0%)
pN12 (7.1%)0 (0.0%)
pM (n = 28) >0.999
pM018 (64.3%)9 (100.0%)
pM11 (3.6%)0 (0.0%)
Molecular classification (n = 22) 0.067
Luminal 16 (42.9%)1 (14.3%)
Luminal 22 (14.3%)1 (14.3%)
HER22 (14.3%)2 (28.6%)
All-negative4 (28.6%)1 (14.3%)
Core basal0 (0.0%)3 (42.9%)
Local recur (n = 26)4 (15.4%)1 (12.5%)>0.999
Distant metastasis during F/U (n = 26)6 (23.1%)0 (0.0%)0.132
Disease progression (n = 26)8 (30.8%)1 (12.5%)0.190
Progression-free survival (months)24 (2–114)24 (1–36)0.322
F/U duration (months)31 (2–132)26 (1–69)0.083
SeC, sebaceous carcinoma; F/U, follow-up.
Table
Table 2. List of genes that harbor mutations in ocular and extraocular sebaceous carcinomas.
Table 2. List of genes that harbor mutations in ocular and extraocular sebaceous carcinomas.
GeneOcular SeC
(n = 20)		GeneExtraocular SeC
(n = 9)
TP5313 (65.0%)NOTCH15 (55.5%)
PIK3CA4 (20.0%)TP533 (33.3%)
SMARCA43 (15.0%)MSH22 (22.2%)
ATM3 (15.0%)SMARCA42 (22.2%)
FBXW73 (15.0%)ATM1 (11.1%)
RB13 (15.0%)BAP11 (11.1%)
BRAF2 (10.0%)BRAF1 (11.1%)
ERBB22 (10.0%)BRCA11 (11.1%)
MLH12 (10.0%)CREBBP1 (11.1%)
RUNX12 (10.0%)DNMT3A1 (11.1%)
VHL2 (10.0%)ERBB21 (11.1%)
MSH62 (10.0%)FBXW71 (11.1%)
RNF432 (10.0%)FGFR21 (11.1%)
KMT2A2 (10.0%)HRAS1 (11.1%)
NF12 (10.0%)KMT2A1 (11.1%)
APC1 (5.0%)MLH11 (11.1%)
BRCA11 (5.0%)MSH61 (11.1%)
BRCA21 (5.0%)NF11 (11.1%)
CREBBP1 (5.0%)NOTCH21 (11.1%)
EGFR1 (5.0%)PIK3CA1 (11.1%)
GNAS1 (5.0%)PIK3R11 (11.1%)
HRAS1 (5.0%)SMAD41 (11.1%)
IDH11 (5.0%)VHL1 (11.1%)
MSH21 (5.0%)
NOTCH21 (5.0%)
PTEN1 (5.0%)
SMAD41 (5.0%)
TSC11 (5.0%)
WT11 (5.0%)
SeC, sebaceous carcinoma.
Table
Table 3. List of genes with frequent pathogenic mutations according to the molecular subtypes.
Table 3. List of genes with frequent pathogenic mutations according to the molecular subtypes.
GeneLuminal 1
(n = 7)		Luminal 2
(n = 3)		HER2
(n = 4)		All-Negative
(n = 5)		Core Basal
(n = 3)
TP535 (38.5%)1 (7.7%)3 (23.1%)3 (23.1%)1 (7.7%)
SMARCA43 (60%)0 (0%)1 (20%)0 (0%)1 (20%)
FBXW74 (100%) *0 (0%)0 (0%)0 (0%)0 (0%)
NOTCH10 (0%)0 (0%)1 (25%)0 (0%)3 (75%) †
ATM1 (33.3%)0 (0%)1 (33.3%)1 (33.3%)0 (0%)
BRAF2 (66.7%)0 (0%)1 (33.3%)0 (0%)0 (0%)
PIK3CA2 (66.7%)0 (0%)1 (33.3%)0 (0%)0 (0%)
BRCA11 (50%)0 (0%)0 (0%)1 (50%)0 (0%)
ERBB22 (100%)0 (0%)0 (0%)0 (0%)0 (0%)
HRAS1 (50%)1 (50%)0 (0%)0 (0%)0 (0%)
MLH11 (50%)0 (0%)0 (0%)0 (0%)1 (50%)
MSH20 (0%)0 (0%)1 (50%)1 (50%)0 (0%)
MSH60 (0%)0 (0%)1 (50%)1 (50%)0 (0%)
NOTCH21 (50%)0 (0%)1 (50%)0 (0%)0 (0%)
VHL2 (100%)0 (0%)0 (0%)0 (0%)0 (0%)
* Significantly enriched in luminal A than in other subtypes. † Significantly enriched in core basal than in other subtypes.
Table
Table 4. Mutations in suspected Muir–Torre syndrome patient.
Table 4. Mutations in suspected Muir–Torre syndrome patient.
Sample 29 SeC-OSample 12 SeC-EO
MLH1 mutationVAF (%)MLH1 mutationVAF (%)
p.M587Hfs * 645.5p.M587Hfs * 660.16
Other pathogenic mutationsOther pathogenic mutations
FBXW7p.Q303 *NOTCH1p.W745 *
TP53p.R342P
TP53p.R175H
SeC-O, ocular sebaceous carcinoma; SeC-EO, extraocular sebaceous carcinoma; VAF, variant allele frequency.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Na, H.Y.; Park, J.H.; Shin, S.A.; Lee, S.; Lee, H.; Chae, H.; Choung, H.; Kim, N.; Chung, J.-H.; Kim, J.E. Targeted Sequencing Revealed Distinct Mutational Profiles of Ocular and Extraocular Sebaceous Carcinomas. Cancers 2021, 13, 4810. https://doi.org/10.3390/cancers13194810

AMA Style

Na HY, Park JH, Shin SA, Lee S, Lee H, Chae H, Choung H, Kim N, Chung J-H, Kim JE. Targeted Sequencing Revealed Distinct Mutational Profiles of Ocular and Extraocular Sebaceous Carcinomas. Cancers. 2021; 13(19):4810. https://doi.org/10.3390/cancers13194810

Chicago/Turabian Style

Na, Hee Young, Jeong Hwan Park, Sun Ah Shin, Sejoon Lee, Heonyi Lee, Heejoon Chae, HoKyung Choung, Namju Kim, Jin-Haeng Chung, and Ji Eun Kim. 2021. "Targeted Sequencing Revealed Distinct Mutational Profiles of Ocular and Extraocular Sebaceous Carcinomas" Cancers 13, no. 19: 4810. https://doi.org/10.3390/cancers13194810

APA Style

Na, H. Y., Park, J. H., Shin, S. A., Lee, S., Lee, H., Chae, H., Choung, H., Kim, N., Chung, J. -H., & Kim, J. E. (2021). Targeted Sequencing Revealed Distinct Mutational Profiles of Ocular and Extraocular Sebaceous Carcinomas. Cancers, 13(19), 4810. https://doi.org/10.3390/cancers13194810

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