Extramammary Paget’s disease (EMPD) is a rare adenocarcinoma typically presenting as an intraepithelial carcinoma [1]. Areas with high apocrine gland concentration such as the vulva, scrotum, penis, perineum, perianal region and axillae are predominantly affected [1]. Although there is favourable prognosis for the early disease stage, in certain cases, it can advance to invasive EMPD. Invasive EMPD is characterized by infiltration into deeper tissues and is associated to nodal and distant metastasis [1]. EMPD usually occurs in elderly individuals aged from 60 to 80 years, with an annual incidence ranging from 0.1 to 2.4 new cases per million people [1]. Only about 20% of patients with EMPD present with distant metastatic disease [1]. As a result, there is little guidance on treatment strategies, resulting in dismal outcomes [2].
Managing EMPD is complex due to the lack of standardized treatment protocols [3]. Surgical intervention, such as wide local excision or radical vulvectomy, is common for localized disease, yet recurrence rates range from 20 to 70% [3]. Invasive cases may require inguinofemoral lymphadenectomy or sentinel lymph node biopsy, while non-surgical options like radiation therapy (RT) are considered for unresectable, recurrent, or metastatic disease [3]. Despite these treatment modalities, improving patient outcomes through personalized targeted therapy continues to pose a challenge.
EMPD share similar histological and clinical features as Mammary Paget’s Disease (MPD) [4]. MPD is a skin condition marked by eczematous lesions usually presenting on the nipple or areolar regions and is frequently linked to underlying in situ or invasive breast cancer [4]. Human epidermal growth factor receptor 2 (HER2) is a transmembrane receptor tyrosine kinase encoded by the ERBB2 gene. Within the epidermal growth factor receptor family, HER2 is widely recognized for its oncogenic role in a subset of invasive breast cancers [5]. In MPD, ERBB2 amplification and consequent HER2 overexpression has been demonstrated in up to 90% of the cases [6, 7]. In EMPD, however, the frequency of HER2 overexpression has been generally inconsistent, with small-scale studies reporting lower frequencies of HER2 overexpression ranging from 15 to 65% [7,8,9,10], while ERBB2 gene amplification ranged from 13 to 43% [8, 11,12,13]. Despite the anecdotal occurrences of metastatic EMPD cases, many studies have posited a correlation between HER2 overexpression and the invasive nature of EMPD [7, 9, 12, 13]. Recent advancements made towards developing HER2-targeted therapies have greatly expanded the treatment options of HER2-expressing breast cancer [11, 13]. Such HER2-directed approaches have also been described in a few EMPD cases [2, 14,15,16]. However, the treatment efficacy has yet to be ascertained given differing patient responses within the limited number of cases studied. Furthermore, none of the studies have described the clinical observations in the context of the underlying tumor genomic landscape.
Herein we present a case report of a patient with de novo metastatic EMPD who responded rapidly to HER2-directed therapy combined with other agents in a comprehensive treatment regimen. Whole genome sequencing (WGS) revealed a distinctive copy number landscape with potential clinical implications.
Materials and methods
Patient data and biospecimen collection
All clinical information was retrieved from electronic medical records. Verification of demographic data including sex, age and ethnicity of the affected patient was corroborated by National Registry Identification Card. All immunohistochemistry was performed at the pathology services laboratory of the Singapore General Hospital and histological parameters were reviewed by an expert dermatopathologist (B.I.).
Whole genome sequencing (WGS) and variant calling
DNA isolated from snap-frozen tumor tissue and matched whole blood was selected for whole-genome sequencing. Whole-genome sequencing was performed on the Illumina NovaSeq 6000 platform as paired-end 150 bp reads, using DNA inserts averaging 350 bp (NovogeneAIT Genomics). The raw sequencing data obtained from the WGS was subjected to a series of data pre-processing steps using the nextflow sarek workflow [17]. This workflow incorporates various tools and pipelines for quality control, alignment, and variant calling. The following steps were performed: Sequencing quality control (FastQC) [18]; Map Reads to Reference GRCh38 (GATK) (BWA mem) [19]; Mark Duplicates (GATK MarkDuplicates) [20]; Base (Quality Score) [21]; Recalibration (GATK BaseRecalibrator, GATK ApplyBQSR) [22]; Preprocessing quality control (samtools stats) [23]; Preprocessing quality control (mosdepth) [24]; Overall pipeline run summaries (MultiQC) [25]; variant calling (Mutect2) [26]; variant annotation (VEP) [27]. After annotating the variants, a rigorous selection process was applied to retain solely non-synonymous variants with deleterious effects, resulting in a refined subset of 43 variants.
Mutational signatures
The mSigAct tool [28] was employed to investigate the presence and activity of mutational signatures within our WGS data. mSigAct employs a conservative maximum likelihood approach to determine the presence of specific mutational signatures within a spectrum. Additionally, it identifies the minimum subset of signatures necessary to reconstruct the observed spectrum. The sparse assign signatures functionality favors using the fewest signatures possible. It also enables maximum a posteriori estimation of signature activity by considering the proportion of tumors with a particular signature in a given type and the likelihood of a specific signature combination generating the observed spectrum.
Copy number variations (CNV) analysis
Whole-genome sequencing data for the matched tumor-normal pairs was analyzed to detect copy number variations (CNVs) using CNVkit [29]. Based on the assumption that the normal sample has a diploid genome with no major copy number alterations, CNVkit compares the read counts in the cancer sample with those in the normal sample. The baseline for copy number calling on the tumor cohort was applied using the mean value of overall coverage. Copy counts were computed based on the depth of coverage in each location after GC bias and regional bias were corrected. The threshold of a notable CNV gain was set as log2FC > 0.9; likewise, a notable CNV loss was reported if log2FC < -0.9. Only genes that met the threshold of CNV gain and loss were considered for subsequent pathway enrichment analysis. The WEB-based GEne SeT AnaLysis Toolkit (Web Gestalt) tool was used for Over Representation Analysis (ORA) [30]. Enriched gene sets with a false discovery rate (FDR) of < 0.05 were considered significant. Pathway Enrichment via Functional Annotation Clustering was conducted with the Database for Annotation, Visualization, and Integrated Discovery (DAVID) [31]. Annotation clusters with Benjamini values < 0.05 were considered significant.
The median VAF was inferred from the SNV analysis. To obtain the purity-normalized absolute copy number, the sample was assumed to be copy-neutral (ploidy = 2) and a cutoff log2(copy_ratio) of 0.9 was applied to identify regions of significant amplification. AmpliconArchitect (AA) [32] was used to detect and characterize the amplicon structures. AA used an absolute copy number cutoff of 5 and a minimum span of 100 kilobase pairs (kbp) for the detection of any circular DNA structures that can carry amplified genetic material. This criterion was used to identify the seeding intervals and potentially seek out amplicon structures. The annotations provided by AA revealed the amplified genes associated with cancer within the identified amplicon. To further analyze the structure of amplicon, AA_classifier was used. Through the analysis, it was determined that the identified amplicon lacked cyclic edges, ruling out the possibility of extrachromosomal DNA (ecDNA) or associations to a break-fusion-bridge event and concluding that a simple-linear amplification event is observed in the sample.
TSGene, an evidence-based online repository dedicated to tumor suppressor genes (TSGs) [33], catalogs data pertaining to 1217 human TSGs (including 1018 genes and 199 long non-coding RNAs). Venn analysis was conducted on genes that met the specified threshold for CNV loss (log2FC < -0.9), revealing the subset of genes that shared overlap with known TSGs.
Results
Clinical case presentation
The patient is a 68-year-old man of Chinese ethnicity, who initially presented with progressive low back pain associated with unintended weight loss of 5 kg over two months. Physical examination revealed an erythematous hyperkeratotic rash with central ulceration over the scrotum (Fig. 1A). He also had a fever measuring 38.6℃. Fluorodeoxyglucose-Positron-Emission Tomography/Computed Tomography (FDG-PET/CT) imaging showed an FDG-avid scrotal wall mass measuring 2.8 cm x 2.5 cm (Maximum Standardized Uptake Value, SUVmax 6.9), extensive FDG-avid mixed lytic-sclerotic bony lesions involving the spine and ribs (SUVmax 9.3), as well as the right iliac and inguinal lymph nodes (SUVmax 6.7) (Fig. 1B). Histologic examination of the scrotal skin showed thickened epidermis replaced by atypical epithelioid cells with enlarged hyperchromatic nuclei and abundant clear to eosinophilic or vacuolated cytoplasm, as well as regions of invasive carcinoma. Bone marrow biopsy showed moderate haematopoietic effacement. Primary tumor cells were diffusely positive for CK7, CK20, Cam5.2 and negative for SOX10, TTF1, PSA, and CDX2. Metastatic tumor cells exhibited a similar staining pattern. Additional immunostaining on the metastatic tumor was positive for GATA3, androgen receptor and estrogen receptor, but negative for progesterone receptor. HER2 was positive in over 90% of the tumors cells with at least 40% showing 3 + intensity (Fig. 2). Given the known association of Paget’s disease with the genitourinary and gastrointestinal tracts [34], oesophagogastroduodenoscopy and colonoscopy, as well as urine cytology were performed, all of which gave negative results for malignancy.
Over the course of investigation, the patient developed an emergent syndrome indicative of rapid disease progression. Serial blood work over 14 days revealed worsening microangiopathic haemolytic anaemia (MAHA) with rapidly falling haemoglobin levels (nadir, 5.2 g/dL), platelet counts (nadir, 33 × 109/L), raised serum lactate dehydrogenase (LDH) levels (from 1,908 to > 15,000 U/L), serum hyperferritinaemia (from 5,300 to 31,354 µg/L), decreased serum haptoglobin levels (< 0.10 g/L) and negative Direct Coomb’s test. The patient subsequently responded promptly to a combination of intravenous paclitaxel (80 mg/m2) and intravenous trastuzumab (8 mg/kg). Subcutaneous denosumab (120 mg) was administered on days 1, 8 and 15. Further doses of paclitaxel on days 8 and 15 were omitted due to the development of pneumonia requiring antibiotic therapy. By the time he was due for cycle 2 of treatment, ferritin (3,465 µg/L), LDH (983 U/L), haemoglobin (8.8 g/dL) and platelet counts (318 × 109/L) had improved significantly. The patient was continued on HER2-directed therapy along with an anti-androgen (oral bicalutamide 150 mg daily, subcutaneous trastuzumab 600 mg and denosumab 120 mg once every 4 weeks) for another 7 cycles until disease progression (partial response at best). Subsequently, bicalutamide was replaced by paclitaxel (175 mg/m2, given once every 4 weeks) for a total of 6 cycles until disease progression, achieving partial response (PR) as best response. Thereafter, he was treated with capecitabine (1,000 mg twice daily for 2 weeks) and lapatinib (1,250 mg daily) given every 3-weeks for nine cycles (stable disease as best response), followed by a single dose of TDM-1 (no response). Then, the patient’s condition worsened with rapid elevation of ferritin and LDH levels, leading to his eventual demise about 22 months after diagnosis. The treatment strategy and serial PET/CT images are shown in Fig. 3.
Copy number variants (CNV) analysis
The WGS data revealed several chromosomal regions exhibiting copy number gains and losses in comparison to the reference normal sample. Among these regions, 81 unique genes on genomic loci chromosomes 7 and 8 demonstrated notable copy number gains (log2FC > 0.9), of which 92.6% of them (75 genes) are located on chromosome 8 (Fig. 4). Web-based databases were used to functionally annotate the 75 genes on chromosome 8 that harbored copy number gains (Table 1). This gene list showed an enrichment for the Response to transforming growth factor-beta (TGFβ) Pathway (GO:0071559, FDR = 0.0376, Enrichment Ratio = 8.12). Genes involved in the TGFβ pathway included ADAM9, FNTA, HTRA4, SFRP1, STAR and ZNF703. Additionally, these genes were enriched in a FGFR1 fusion signaling pathway (REACTOME Pathway ID: R-HSA-8,853,336, FDR = 0.0082, Enrichment Ratio = 2.3). The genes identified in this pathway were BAG4, ERLIN2 and FGFR1 (Fig. 5). Other significant pathways that were enriched include the Benzene-containing Compound Metabolic process (GO:0042537, FDR = 0.0376, Enrichment Ratio = 40.2). The genes involved are IDO1, IDO2 and STAR (Fig. 5). ERBB2 gene on chromosome 17 presented a copy number log2FC value of 0.4 despite strong HER2 protein overexpression observed from immunohistochemistry.
A tumor purity of about 44% was estimated using the median Variant Allele Frequency (VAF) from SNV analysis, resulting in a purity-normalized absolute copy number of 5.937. Among five seeding intervals identified on chromosomes 7 and 8 using the AmpliconArchitect (AA) tool, only one amplicon showed gene amplification, including ZNF703, HOOK3, DDHD2, LSM1, WHSC1L1, ADAM9, BRF2, MYST3, and FGFR1. These genes were also detected in the CNV analysis. The lack of cyclic edges in this amplicon suggests that the amplification event is a simple-linear one.
Genes with notable copy number loss (log2FC < -0.9) were primarily located on the X and Y sex chromosomes, as well as chromosomes 7, 8 and 9. Among these genes, CDKN2A, CDKN2B, CUX1, LRP1B and PTPRD were retrieved as tumor suppressors from the TSGene database.
Somatic variants and mutational signatures
Variant calling analysis identified 43 genes harbouring non-synonymous variants with potentially deleterious effects (Table 2), including PRAMEF2 and PTPRD. From WGS data, COSMIC mutational signatures were extracted, with 49.90%, 10.91%, and 10.13% of the data resembling COSMIC signatures ID5, ID1, and ID8, respectively. Additionally, 34.25%, 15.76% and 11.36% resembled signature SBS5, SBS2, and SBS40 respectively. For DBS signatures, 34.85%, 20.71% and 13.13% resembled DBS11, DBS7, and DBS10 respectively (Fig. 6). These signatures primarily relate to aging (ID5, ID1, ID8, SBS5, SBS40) and defective DNA repair (DBS7, DBS10), with SBS2 and DBS11 linked to APOBEC. SBS2’s co-occurrence with SBS13 [35] aligns with our findings; 10% of the data resembled SBS13.
Discussion
Most patients with intraepithelial EMPD achieve remission with an anticipated 5-year survival of over 60% [2]. However, when the disease invades the dermis, this percentage drops to 10%, with about one-third of patients developing lymph node involvement and distant metastases [2]. Spinal metastases in EMPD are rare but have been reported in a few cases. In two of such described cases, compression fractures in the vertebral bodies were described, and both cases resulted in patient demise within 5 to 6.2 months of spinal metastasis detection [36, 37]. Several studies have investigated the potential correlation between EMPD invasiveness and HER2 protein overexpression, although conclusive findings are hindered by small sample sizes. Invasive EMPD typically have a larger percentage of 2 + or 3 + HER2 IHC scores compared to in situ EMPD. However, IHC HER2 positivity ranges widely across reports between 18 − 71% and between 10 − 62% in invasive EMPD and in situ EMPD, respectively [7, 10, 12, 38, 39]. Therefore, despite serving as a useful biomarker, it cannot be assumed that HER2 positivity signifies EMPD invasiveness. In our case, the patient’s tumor was metastatic and exhibited strong IHC HER2 positivity. However, ERBB2 presented minimal CNV gain from our CNV analysis. Previous studies have also reported IHC HER2 positive invasive EMPD cases with an absence of ERBB2 gene amplification via FISH [8] and CISH [10]. While gene amplification detected by CNV gains using WGS cannot be directly compared against ISH methods, there seems to be a weak correlation between ERBB2 gene amplification and HER2 protein overexpression. Exploring other signaling pathways associated with HER2 protein overexpression instead could be worthwhile.
The patient was administered with a unique EMPD treatment regimen combining HER2-targeting trastuzumab with paclitaxel together with denosumab, followed by the subsequent replacement of paclitaxel with anti-androgen bicalutamide. The patient exhibited a PR to this treatment for 7 cycles. A study by Sekiguchi et al. administered trastuzumab and paclitaxel combination chemotherapy as well. Of the four cases reported, one case exhibited PR with progression free survival (PFS) of 12 months, while the other 3 cases met with early disease progression [16]. Contrastingly, a study by Takahagi et al. administering the same combination chemotherapy eradicated HER2-expressing tumor cells in the patient’s dermal nests and lymph vessels 6 months into treatment. However, metastases in the CNS could not be reduced and eventually led to the demise of the patient [15]. In these studies, because the patients were previously treated with docetaxel, cross-resistance between docetaxel and paclitaxel was a concern [16]. The response of our patient was not influenced by any prior use of docetaxel.
Estrogen receptor, androgen receptor and progesterone receptor expression has been reported to be positive in 13%, 40% and 8% of patients with EMPD, respectively [40]. Previous EMPD treatments combining bicalutamide and leuprorelin with trastuzumab [2] and other anti-androgen drugs [41] have failed to prevent rapid disease progression. Only one case presented reduced bone and lymph node metastases upon bicalutamide and leuprorelin acetate treatment [42]. It has also been hypothesized that the androgen receptor overexpression might just be an inherent characteristic of glandular cells with apocrine differentiation [2]. With the limited anti-androgen treatments administered to EMPD, it is difficult to ascertain whether EMPD tumorigenesis depends on the androgen receptor signaling pathway. Finding an effective personalized treatment approach for each EMPD case remains a significant challenge even with the presence of certain biomarkers, and this is further complicated by the limited number of previous studies available as references.
An enrichment of genes on chromosome 8 associated with TGFβ was also reported. Although the role of TGFβ is not well studied in EMPD, a study by Hirakawa et al. [43] has investigated the onset of epithelial-mesenchymal transition (EMT) in EMPD, possibly enhanced by TGFβ. EMT is a process whereby polarized epithelial cells undergo biochemical remodelling to develop into mesenchymal cells, thus contributing to metastasis [44]. Our patient experienced metastases in the right iliac and inguinal lymph nodes, spine, and ribs, suggesting that lymphoangiogenesis and subsequent lymphatic invasions contributed to the metastatic spread. Histological examination of EMPD lesions have also revealed notable difference in the development and enlargement of lymphatic and blood vessels in the dermis compared to healthy skin or even other skin cancers including melanoma [43]. Previous studies have observed that invasive Paget cells demonstrated lower E-cadherin expression and higher expression of mesenchymal markers including N-cadherin and/or vimentin than normal epidermal keratinocytes and in situ carcinoma Paget cells, and the expression of EMT-related markers correlated with the incidence of lymphatic invasion in primary skin tumor [43, 45]. These findings substantiate the involvement of EMT-like mechanisms in the invasion of tumor cells in EMPD, particularly lymphatic invasion, and this is further enhanced by TGFβ. Additionally, some cell-based models have investigated the interaction between TGFβ and HER2 signalling pathways as well. For instance, TGFβ modulates the actin cytoskeleton, inducing cell migration and EMT in HER2-overexpressing breast cancer cells [44]. Exogenous TGFβ could also induces trastuzumab resistance in HER2-positive breast cancer cell lines [5]. Taken together, these findings highlight the crucial roles of HER2 and TGFβ in breast cancer metastasis. The functional crosstalk between TGFβ and HER2 could be of potential relevance to EMPD.
There was also an enrichment of the pathways associated with FGFR1 signalling. Fibroblast growth factor receptor (FGFR) is a target that has been explored as part of metastatic breast cancer therapies [46]. Various studies have shown that gene-translocations can lead to the expression of fusion proteins with FGFR, resulting to aberrant FGFR signalling which can act as oncogenic drivers of cancer [46, 47]. Despite limited research on the role of FGFR in EMPD, one study by Ishida et al. described that out of 87 cases, only 3 cases harbored FGFR1 amplification [48]. Interestingly, the study postulated that amplifications or mutations in ERBB2, mutations in ERBB3 and FGFR1 amplification were mutually exclusive, suggesting that alterations in these 3 genes represents distinct disease subsets [48]. In our case, we observed HER2 overexpression via IHC and FGFR1 amplification via notable CNV gains. Further research is necessary to establish the significance of our preliminary WGS data.
Among the genes that were showing notable CNV gains, IDO1 and IDO2 (log2 FC = 1.22256) were particularly interesting. IDO is an immunoregulatory enzyme; upon inflammatory stimulus, IDO converts tryptophan to kynurenine, inducing immunosuppression and inhibiting immune activation [49]. Similar to checkpoint molecules such as PDL1, cancer cells exploit these pathways to evade the immune system and suppress T-cell-mediated antitumor responses [49]. Promising IDO-inhibitors including indoximod (1-methyl-D-tryptophan) and epacadostat have shown promising results in clinical trials, although often in combination with chemotherapy and checkpoint inhibitors. Notably, there is a metastatic EMPD case reported with durable response to checkpoint inhibitor immunotherapy; 4 cycles of ipilimumab plus nivolumab treatment resulted in 7 months PR [50]. Collectively, these findings highlight the potential of future EMPD treatment regimens involving IDO inhibitors.
Among the genes with notable CNV loss, 5 tumor suppressors (CDKN2A, CDKN2B, CUX1, LRP1B and PTPRD) were identified. A prior EMPD study also revealed deletions of areas harbouring CDKN2A via CNV analysis [48]. CUX1 was identified as a potential EMPD driver mutation, exhibiting deletions and truncating mutations in 19.5% of the EMPD cases studied, with 9.2% showing nonsense or splice-site mutations as well [48].. The paralleled findings suggest that the CNV loss of these tumor suppressor genes may have contributing roles to EMPD development.
Conclusion
In conclusion, the presented case exemplifies the extrapolation of information from treating common cancers– in this case, breast cancer– to rare cancers by identifying their similar histological and molecular profile, thereby improving the treatment approaches. Despite being confronted with poor prognostic factors including spinal metastases, associated MAHA and HER2 positivity, the patient responded remarkably and survived for another 22 months post diagnosis under the administration of HER2-based treatment in combination with other therapeutic agents that were introduced sequentially over multiple treatment cycles. Through this case, both clinical information and genomic analysis were integrated to attain a more comprehensive understanding of the disease landscape.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request. The datasets supporting the conclusions of this article are available in the European Genome-Phenome Archive repository (accession number EGAS50000000243).
References
Ohara K, Fujisawa Y, Yoshino K, Kiyohara Y, Kadono T, Murata Y, et al. A proposal for a TNM staging system for extramammary Paget disease: retrospective analysis of 301 patients with invasive primary tumors. J Dermatol Sci. 2016;83:234–9.
Zattarin E, Nichetti F, Ligorio F, Mazzeo L, Lobefaro R, Fucà G et al. Case Report: prolonged clinical benefit with sequential trastuzumab-containing treatments in a patient with advanced extramammary paget disease of the groin. Front Oncol 2022;12.
Caruso G, Barcellini A, Mazzeo R, Gallo R, Vitale MG, Passarelli A, et al. Vulvar Paget’s Disease: a systematic review of the MITO Rare Cancer Group. Cancers. 2023;15:1803.
St. Claire K, Hoover A, Ashack K, Khachemoune A. Extramammary Paget disease. Dermatol Online J 2019;25.
Chow A, Arteaga CL, Wang SE. When Tumor suppressor TGFβ meets the HER2 (ERBB2) Oncogene. J Mammary Gland Biol Neoplasia. 2011;16:81–8.
Bianco MK, Vasef MA. HER-2 gene amplification in Paget Disease of the Nipple and Extramammary Site: a chromogenic in situ hybridization study. Diagn Mol Pathol. 2006;15:131.
Gatalica Z, Vranic S, Krušlin B, Poorman K, Stafford P, Kacerovska D, et al. Comparison of the biomarkers for targeted therapies in primary extra-mammary and mammary Paget’s disease. Cancer Med. 2020;9:1441–50.
Lu X, Zhang P, Zhu Y, Ye D. Human epidermal growth factor receptor 2 amplification as a biomarker for treatment in patients with lymph node-metastatic penoscrotal extramammary Paget’s disease. Oncol Lett. 2019;17:2677–86.
Takata M, Fujimoto A, Aoki H, Hatta N, Ooi A, Takehara K. erbB-2 overexpression but no activation of beta-catenin gene in extramammary Paget’s disease. J Invest Dermatol. 1999;113:258–62.
Tanskanen M, Jahkola T, Asko-Seljavaara S, Jalkanen J, Isola J. HER2 oncogene amplification in extramammary Paget’s disease. Histopathology. 2003;42:575–9.
Klos KS, Zhou X, Lee S, Zhang L, Yang W, Nagata Y, et al. Combined trastuzumab and paclitaxel treatment better inhibits ErbB-2-mediated angiogenesis in breast carcinoma through a more effective inhibition of akt than either treatment alone. Cancer. 2003;98:1377–85.
Maeda T, Kitamura S, Nishihara H, Yanagi T. Extramammary Paget’s disease patient-derived xenografts harboring ERBB2 S310F mutation show sensitivity to HER2-targeted therapies. Oncogene. 2020;39:5867–75.
Rassy E, Rached L, Pistilli B. Antibody drug conjugates targeting HER2: clinical development in metastatic breast cancer. Breast off J Eur Soc Mastology. 2022;66:217–26.
Hanawa F, Inozume T, Harada K, Kawamura T, Shibagaki N, Shimada S. A case of metastatic Extramammary Paget’s Disease responding to Trastuzumab plus Paclitaxel Combination Therapy. Case Rep Dermatol. 2011;3:223–7. h.
Takahagi S, Noda H, Kamegashira A, Madokoro N, Hori I, Shindo H, et al. Metastatic extramammary Paget’s disease treated with paclitaxel and trastuzumab combination chemotherapy. J Dermatol. 2009;36:457–61.
Sekiguchi N, Kubota S, Noguchi T, Fukushima T, Kobayashi T, Kanda S, et al. Experiences of trastuzumab plus paclitaxel combination therapy in metastatic human epidermal growth factor receptor 2-positive extramammary Paget’s disease: four cases and a review. J Dermatol. 2020;47:1276–9.
Garcia M, Juhos S, Larsson M, Olason PI, Martin M, Eisfeldt J, et al. Sarek: a portable workflow for whole-genome sequencing analysis of germline and somatic variants. F1000Research. 2020;9:63.
Simon A. FastQC: a quality control tool for high throughput sequence data. 2010.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.
van der Auwera G, O’Connor BD. Genomics in the cloud: using Docker, GATK, and WDL in Terra. First edition. Sebastopol, CA: O’Reilly Media; 2020.
Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using Phred I. Accuracy Assessment. Genome Res. 1998;8:175–85.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Pedersen BS, Quinlan AR. Mosdepth: quick coverage calculation for genomes and exomes. Bioinformatics. 2018;34:867–8.
Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinforma Oxf Engl. 2016;32:3047–8.
Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31:213–9.
McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GRS, Thormann A, et al. Ensembl Variant Effect Predictor Genome Biol. 2016;17:122.
Rozen SG. mSigAct 2023.
Talevich E, Shain AH, Botton T, Bastian BC. CNVkit: genome-wide Copy Number Detection and visualization from targeted DNA sequencing. PLOS Comput Biol. 2016;12:e1004873.
Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019;47:W199–205.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Deshpande V, Luebeck J, Nguyen N-PD, Bakhtiari M, Turner KM, Schwab R, et al. Exploring the landscape of focal amplifications in cancer using AmpliconArchitect. Nat Commun. 2019;10:392.
Zhao M, Sun J, Zhao Z. TSGene: a web resource for tumor suppressor genes. Nucleic Acids Res. 2013;41:D970–6.
Jones RE, Austin C, Ackerman AB. Extramammary Paget’s disease. A critical reexamination. Am J Dermatopathol. 1979;1:101–32. https://doi.org/10.1097/00000372-197900120-00002.
Sundaramoorthy R, Murugan A, Nakoaka H, Inoue I. APOBEC: a molecular driver in cervical cancer pathogenesis. Cancer Lett. 2021;496:104–16.
Beleznay KM, Levesque MA, Gill S. Response to 5-fluorouracil in metastatic extramammary paget disease of the scrotum presenting as pancytopenia and back pain. Curr Oncol. 2009;16:81–3.
Little JT, Lehman VT, Morris JM, Lehman JS, Diehn FE. Spinal metastases of Extramammary Paget Disease with radiologic-pathologic correlation. J Radiol Case Rep. 2016;10:1–8.
Richter CE, Hui P, Buza N, Silasi D-A, Azodi M, Santin AD, et al. HER-2/NEU overexpression in vulvar paget disease: the Yale experience. J Clin Pathol. 2010;63:544–7.
Tanaka R, Sasajima Y, Tsuda H, Namikawa K, Tsutsumida A, Otsuka F, et al. Human epidermal growth factor receptor 2 protein overexpression and gene amplification in extramammary Paget disease. Br J Dermatol. 2013;168:1259–66.
Angelico G, Santoro A, Inzani F, Straccia P, Arciuolo D, Mulè A, et al. Hormonal environment and HER2 status in Extra-mammary Paget’s Disease (eMPD): a systematic literature review and Meta-analysis with clinical considerations. Diagnostics. 2020;10:1040.
Yoneyama K, Kamada N, Kinoshita K, Kawashima T, Otani M, Endo H, et al. Androgen-deprivation regimen for multiple bone metastases of extramammary Paget disease. Br J Dermatol. 2005;153:853–5.
Liegl B, Horn L-C, Moinfar F. Androgen receptors are frequently expressed in mammary and extramammary Paget’s disease. Mod Pathol. 2005;18:1283–8.
Hirakawa S, Detmar M, Kerjaschki D, Nagamatsu S, Matsuo K, Tanemura A, et al. Nodal lymphangiogenesis and metastasis: role of tumor-induced lymphatic vessel activation in extramammary Paget’s disease. Am J Pathol. 2009;175:2235–48.
Derynck R, Muthusamy BP, Saeteurn KY. Signaling pathway cooperation in TGF-β-induced epithelial–mesenchymal transition. Curr Opin Cell Biol. 2014;31:56–66.
Cho Z, Konishi E, Kanemaru M, Isohisa T, Arita T, Kawai M, et al. Podoplanin expression in peritumoral keratinocytes predicts aggressive behavior in extramammary Paget’s disease. J Dermatol Sci. 2017;87:29–35.
De Luca A, Esposito Abate R, Rachiglio AM, Maiello MR, Esposito C, Schettino C, et al. FGFR fusions in Cancer: from Diagnostic approaches to therapeutic intervention. Int J Mol Sci. 2020;21:6856.
Krook MA, Reeser JW, Ernst G, Barker H, Wilberding M, Li G, et al. Fibroblast growth factor receptors in cancer: genetic alterations, diagnostics, therapeutic targets and mechanisms of resistance. Br J Cancer. 2021;124:880–92.
Ishida Y, Kakiuchi N, Yoshida K, Inoue Y, Irie H, Kataoka TR, et al. Unbiased detection of driver mutations in Extramammary Paget Disease. Clin Cancer Res off J Am Assoc Cancer Res. 2021;27:1756–65.
Peyraud F, Guegan J-P, Bodet D, Cousin S, Bessede A, Italiano A. Targeting Tryptophan Catabolism in Cancer Immunotherapy Era: challenges and perspectives. Front Immunol. 2022;13:807271.
Guercio BJ, Iyer G, Kidwai WZ, Lacouture ME, Ghafoor S, Rossi AM, et al. Treatment of metastatic Extramammary Paget Disease with Combination Ipilimumab and Nivolumab: a Case Report. Case Rep Oncol. 2021;14:430–8.
Acknowledgments
The authors would like to thank the patient for participating in this study.
Funding
This work was supported by the Singapore Ministry of Health’s National Medical Research Council under its Transition Award (TA21jun-0005), RTF Seed Fund (SEEDFD21jun-0002) and TETRAD II Collaborative Centre Grant (CG21APR2002), as well as the SingHealth Duke-NUS AM/ACP-Designated Philanthropic Fund Grant Award (08/FY2023/EX/27-A65).
Author information
Authors and Affiliations
Cancer Discovery Hub, National Cancer Centre Singapore, Singapore, Singapore
Boon Yee Lim, Zexi Guo, Tun Kiat Ko, Elizabeth Chun Yong Lee, Bavani Kannan, Jing Yi Lee, Abner Herbert Lim, Zhimei Li, Cedric Chuan-Young Ng & Jason Yongsheng Chan
Division of Cellular and Molecular Research, National Cancer Centre Singapore, Singapore, Singapore
Jing Quan Lim
Duke-NUS Medical School, Singapore, Singapore
Jing Quan Lim & Jason Yongsheng Chan
Department of Anatomical Pathology, Singapore General Hospital, Singapore, Singapore
Inny Busmanis
Divison of Medical Oncology, National Cancer Centre Singapore, Singapore, Singapore
Jason Yongsheng Chan
Contributions
B.Y.L. analyzed the data and drafted the manuscript; Z.G. A.H.L., J.Q.L. and T.K.K. performed the bioinformatics analyses; I.B. provided pathological assessment of tissues; J.Y.C provided patient samples and clinical data; B.K., E.C.Y.L, J.Y.L, Z.L. and C.C.N. provided technical expertise and performed various experiments. J.Y.C. conceived the study, interpreted the results, and revised the manuscript; and all authors read and approved the final version of the manuscript.
Corresponding author
Correspondence to Jason Yongsheng Chan.
Ethics declarations
Ethics approval and consent to participate
Written informed consent from the patient for use of biospecimens and clinical data was obtained in accordance with the Declaration of Helsinki. Tissue collection and consent protocols were under ethics approval from the SingHealth Centralized Institution Review Board (2010/426/B).
Consent for publication
Written informed consent has been obtained from the patient to publish this paper.
Competing interests
The authors declare no conflict of interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.