Skip to main content

Genetic profiling of patients with adenoid cystic carcinoma of the Bartholin’s glands reveals potential new routes for targeted therapies: a case report



Bartholin gland carcinomas (BGCs) are rare tumor types, for which no molecular analyses including genomic sequencing have been reported to date. Adenoid cystic carcinomas (ACCs) of the Bartholin’s glands are an atypical histological type of BGC, and currently nothing is known regarding their genetic profiles or similarity to ACC carcinogenesis in other organs including the salivary glands, thereby limiting possible therapeutic options using precision medicine.

Case presentation

We used targeted gene sequencing to analyze the occurrence of 160 cancer-related genes in two patients with BG-ACC. KRAS and KDM6A mutations were detected in tumor samples collected from each patient. No KRAS mutations have been previously reported in salivary gland ACCs, indicating that the carcinogenesis of BG-ACC differs from that of the salivary gland ACCs. KDM6A mutations are often reported in salivary gland ACCs and facilitate novel gene-targeted therapy, including the use of BET and HDAC inhibitors.


A better understanding of the underlying genetic mechanisms will help to clarify the carcinogenesis of BG-ACC. In turn, this will enable treatment with novel targeting agents, as well as the initial exploration of gene-based precision oncological therapies, which aim to improve treatment outcomes for patients with this disease.


Bartholin gland carcinomas (BGCs) are extremely rare tumors, comprising less than 1% of total malignancies in female genitals [1]. They have been categorized histologically as adenocarcinomas and squamous, adenosquamous, and adenoid cystic carcinomas (ACCs) [1]. Among these, adenocarcinomas and squamous carcinomas are the most common, accounting for approximately 40% of cases, whereas adenosquamous and ACCs comprise 5 and 15%, respectively [1]. BGCs have no distinguishing clinical manifestations; therefore, they are difficult to suspect or diagnose during clinical practice. Consequently, most primary BGCs are misdiagnosed and erroneously treated as benign Bartholin gland cysts or abscesses. Even with early correct diagnosis and successful local management, most patients die because of metastases within 10 years [1]. BGCs often metastasize to neighboring nerves, a mechanism known as perineural invasion, which is indicative of poor prognosis and increases the chance of local recurrence [2, 3]. Patients with BGC have poor prognosis owing to a lack of standard treatments, including chemotherapy.

The pathogenesis and carcinogenesis of BGCs remains unclear. The histological classifications, including adenocarcinoma, squamous, and adenosquamous, are similar to those of cervical cancer. In previous studies, it has been shown that a proportion of the squamous type might be related to human papillomavirus (HPV) infection, particularly HPV type 16, as well as to cervical cancer [2, 3]. However, due to a paucity of clinical trials, there is currently no consensus regarding the efficacy of radiation therapy, one of the approaches used to treat cervical cancer in patients with squamous-type BGC. Most ACCs occur in the salivary glands and it is an atypical histological type among genital cancers, although extremely rare cases have been reported in the breast, skin, respiratory system, and Bartholin glands [4]. Currently, given the rarity of this disease, the lack of prospective trials, and unknown carcinogenesis, valid chemotherapies do not exist for the long-term management for patients with BGCs.

Precision medicine, such as that using comprehensive genome sequencing, is a possible treatment approach for a number of cancers. In particular, recent molecular and whole-exome sequencing is beginning to reveal the genetic basis of this uncommon disease [5,6,7,8].

To date, genetic profiling for BG-ACCs has not been reported. In the present study, we analyzed 160 cancer-related genes from two patients with BG-ACCs as an attempt to bridge the current knowledge gap regarding the genetic characteristics of this disease. This information will help to clarify the carcinogenesis of BG-ACCs and will be useful to identify potential therapeutic targets for future precision medicine.

Case presentation

Patient 1 was a 63-year-old woman undergoing excision of a vulvar lesion based on the diagnosis of a Bartholin gland cyst. The resected specimens were reviewed by a pathologist, who diagnosed it as a stage I adenoid cystic carcinoma in the Bartholin glands (BG-ACC). The patient was not provided with adjuvant therapy because there were no residual tumor tissues. However, metastasis to the lung was observed 4 years later. Although the patient had undergone video-assisted thoracoscopic surgery, recurrence at a site in the vulva was detected 5 years after this surgery. She subsequently underwent excision of the recurrent vulvar lesion, which was used for molecular analysis.

Patient 2 was a 61-year-old woman undergoing biopsy for a vulvar lesion and was diagnosed as having stage IV BG-ACC. She subsequently underwent total pelvic exenteration with resection of the inferior ramus. The patient had no adjuvant therapy because there were no residual tumors. However, metastases to the lung and liver were observed 2 years later. A biopsy obtained from the recurrent lesion in the liver was used for molecular analysis to establish an appropriate treatment approach.

Similar pathological findings were obtained for the tumor in case patient 1 and one of the tumors in the case patient 2. Low-power views showed tumor cells organized in nests, sheets, and trabeculae, with a cribriform pattern (Fig. 1; case 1: upper left, case 2: upper right). Higher power views showed nests and sheets of relatively uniform cells with a cribriform pattern. The cribriform regions comprised unstructured eosinophilic basement membrane-like components (Fig. 1; case 1: lower left, case 2: lower right). Expression of p16 was not detected in either tumor. Based on a pathologist’s assessment of the resected samples, the masses were determined to be BG-ACCs (case 1: pT1NxM0, case 2: pT1NxM0).

Fig. 1
figure 1

Histology of Bartholin gland adenoid cystic carcinomas. Cells were found to be organized in nests, sheets, and trabeculae, with a cribriform pattern (low-power view; case 1: upper left, case 2: upper right). We also observed comparatively homogeneous cells in nests and sheets with a cribriform pattern containing unstructured eosinophilic basement membrane-like components (high-power view; case 1: lower left, case 2: lower right)

Genomic DNA sequences were obtained for the two samples classified as BG-ACCs, with average sequencing depths of 706.5× and 695.3× for patient 1 and patient 2, respectively. The average tumor cellularity was 30 and 80%, respectively, as determined histologically. Tumor cellularities estimated using a variant allele frequencies were 30 (case and 80% as well. Profiles of the gene alterations detected in each sample are shown in Table 1.

Table 1 Gene alterations in tumor samples

Actionable gene alterations were detected in each sample. A KRAS point mutation was detected in the tumor from patient 1, and a KDM6A splicing alteration was detected in the tumor from patient 2 (Table 1). Details of these variants are described in Supplemental Figs. S1 and S2. No gene amplification or loss was detected in either sample. Of these two variants, the KDM6A mutation is potentially druggable. For both samples, the tumor mutation burden calculated from our pipeline was 1.3 single-nucleotide variants per megabase. Copy-number variation box and variant allele frequency plots are provided in Supplementary Figs. S3 and S4. With respect to the secondary germline, no ACMG-recommended genes for testing were observed in either patient.

Discussion and conclusions

We sequenced 160 cancer-related genes in tumor samples obtained from two patients diagnosed with BG-ACC, and accordingly, detected KRAS and KDM6A mutations. Given that p16 expression was detected in neither of the samples, we hypothesize that HPV infection, which has been associated with squamous cell carcinomas [1, 9], was probably not involved in carcinogenesis in either of the two patients.

In these cases, we focused on the association between histological types and genetic events. The histological type of both cases was ACC, an uncommon malignancy that can arise in more than one organ site, despite being observed most often in the salivary glands [4]. For salivary gland ACCs, a few alterations have been identified in known cancer-related genes implicated in chromatin regulation, Notch signaling, and a number of other pathways, including PIK3CA, ATM, and TSC1 [7, 10, 11]. Furthermore, recent studies have demonstrated a recurrent t(6;9)(q22–23;p23–24) translocation arising from the fusion of the v-myb myeloblastosis viral oncogene homolog (MYB) gene located on chromosome 6 with the nuclear factor I/B (NFIB) gene located on chromosome 9 in 44% of ACC cases [5, 11, 12]. However, given that analyses of the MYB/NFIB fusion gene have been performed for ACC mainly in the salivary gland, it remains unclear whether this gene is associated with the carcinogenesis of BG-ACC. The degree of contribution to the disease by other genes, and thereby the utility of the genes as possible therapeutic targets, is uncertain. Further, the extent to which other genes contribute to this disease and might constitute additional targets for potential therapeutic exploitation has not been well established, as previous genetic investigations have focused on salivary gland ACCs [10], and no similar sequencing has been performed with respect to BG-ACC.

In the report of Stephens et al., multiple mutations were identified in half of the examined cases, which collectively implicated chromatin deregulation [7]. In addition, somatic gene mutations were discerned in previously identified cancer-associated genes, including PIK3CA, ATM, CDKN2A, SF3B1, SUFU, TSC1, and CYLD. Togashi et al. detected MYB or MYBL1 locus rearrangements in nearly all ACCs examined, suggesting that these might constitute good diagnostic markers for ACCs [12]. However, these authors found that the fusion transcript-specific RT-PCR for MYB-NFIB and MYBL1-NFIB and ordinary split FISH assays for MYB and MYBL1 were less sensitive [12]. In our BG-ACC patient 1, we identified a KRAS mutation, which, to the best of our knowledge, has not been previously detected in salivary gland ACCs; this suggests differences between the carcinogenesis of BG-ACC and salivary gland ACCs [7, 11].

In our BG-ACC patient 2, we identified a mutation in KDM6A, which is related to chromatin remodeling. Previous studies have found alterations in genes implicated in histone modification and chromatin remodeling in cancer development [13]. Studies have reported somatic mutations in ARID1A in clear-cell, renal, transitional-cell bladder, and gastric carcinomas [13,14,15,16], as well as mutations in CREBBP and EP300 in non-Hodgkin lymphoma [17] and in KDM6A in renal and transitional cell bladder carcinomas [13, 14]. However, to date, only truncating mutations of these genes have been reported. These findings, indicating that chromatin regulators are consistently altered in ACC, suggest an important role for the epigenetic control of gene expression in ACC and that the genetic profile of BG-ACC could be similar to ACCs of the salivary gland [6, 7]. Furthermore, it has been previously demonstrated that the loss of KDM6A triggers super-enhancers that cause sex-specific squamous-like pancreatic cancer with sensitivity to bromodomain and extra-terminal motif (BET) protein inhibitors [18]. Because histone pathology has a marked role in this pathway, novel compounds that target histone-related biomarkers might merit further study. For example, romidepsin, a histone deacetylase (HDAC) inhibitor, is an approved therapy for cutaneous and peripheral T cell lymphomas [19] with potential applications for ACC therapy. Given that a number of other pathways have also been linked to ACC development, grouping patients by molecular and histologic subtype will be key to increase our understanding of the clinical, pathological, and molecular correlations among them. However, it is not known whether these shared pathways are indicative of the potential effectiveness of combining current therapies. If this turns out to be the case, fundamental downstream pathways linked to proliferation, cell cycle regulation, angiogenesis, and cell adhesion also suggest possible therapeutic routes. Alternatively, the use of mutations with a lower prevalence in ACC than in other cancers (e.g. SMARCA2) as drug targets for BG-ACC might require further investigation and confirmation. Moreover, there is increasing evidence from human cancer genome sequencing that KDM6A participates in oncogenesis, suggesting it as a therapeutic target for BG-ACCs.

The recent and growing interest in the unique genetic pathways involved in neoplasia is partially a consequence of the promising developments in targeted therapies, which contribute to the increasing potential of personalized-medicine strategies. Previous studies have discovered a number of genetic abnormalities associated with oncogenes and tumor suppressor genes in the tumorigenesis of multiple cancer types. Given that no alterations in the RAS/ERK signaling pathway, which includes KRAS, have been detected in salivary gland ACCs to date, the detection of a KRAS mutation in the tumor of one of our patients suggests that we might have discovered a novel pathway related to BG-ACC carcinogenesis, which could in turn provide a novel therapeutic target.

Our study is limited by both its retrospective nature, which lacked information on fusion genes including MYB/NFIB, as well as by the limited number of patients with BG-ACC. Nonetheless, this is the first study reporting the genomic profiling of BG-ACC. Further genetic analyses of BGCs will undoubtedly elucidate the carcinogenesis of BGC and BG-ACC and will thereby provide a basis for novel gene-targeted therapy as a step toward developing precision medicine strategies.

Availability of data and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.



Bartholin gland carcinoma


Adenoid cystic carcinoma


Human papilloma virus


Bromodomain and extra-terminal motif


Histone deacetylase


  1. Hill DS, Butterfield A. Bartholin’s gland squamous cell carcinoma, a rare vulvar neoplasm. J Diag Med Sonogr. 2010;26:296–8.

    Article  Google Scholar 

  2. Fordice J, Kershaw C, El-Naggar A, Goepfert H. Adenoid cystic carcinoma of the head and neck: predictors of morbidity and mortality. Arch Otolaryngol Head Neck Surg. 1999;125:149–52.

    Article  CAS  Google Scholar 

  3. Marcinow A, Ozer E, Teknos T, Wei L, Hurtuk A, Old M, et al. Clinicopathologic predictors of recurrence and overall survival in adenoid cystic carcinoma of the head and neck: a single institutional experience at a tertiary care center. Head Neck. 2014;36:1705–11.

    Article  Google Scholar 

  4. Brahmer JR, Tykodi SS, Chow LQ, Hwu WJ, Topalian SL, Hwu P, et al. Safety and activity of anti-PD-L1 antibody in patients with advanced cancer. N Engl J Med. 2012;366:2455–65.

    Article  CAS  Google Scholar 

  5. Persson M, Andrén Y, Mark J, Horlings HM, Persson F, Stenman G. Recurrent fusion of MYB and NFIB transcription factor genes in carcinomas of the breast and head and neck. Proc Natl Acad Sci U S A. 2009;106:18740–4.

    Article  CAS  Google Scholar 

  6. Ho AS, Kannan K, Roy DM, Morris G, Ganly I, Katabi N, et al. The mutational landscape of adenoid cystic carcinoma. Nat Genet. 2013;45:791–8.

    Article  CAS  Google Scholar 

  7. Stephens PJ, Davies HR, Mitani Y, Van Loo P, Shlien A, Tarpey PS, et al. Whole exome sequencing of adenoid cystic carcinoma. J Clin Invest. 2013;123:2965–8.

    Article  CAS  Google Scholar 

  8. Ross JS, Wang K, Rand JV, Sheehan CE, Jennings TA, Al-Rohil RN, et al. Comprehensive genomic profiling of relapsed and metastatic adenoid cystic carcinomas by next-generation sequencing reveals potential new routes to targeted therapies. Am J Surg Pathol. 2014;38:235–8.

    Article  Google Scholar 

  9. Felix JC, Cote RJ, Kramer EE, Saigo P, Goldman GH. Carcinomas of Bartholin’s gland. Histogenesis and the etiological role of human papillomavirus. Am J Pathol. 1993;142:925–33.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Wellcome Trust Sanger Institute. Catalogue of somatic mutations in cancer. COSMIC web site Accessed 10 May 2013.

  11. Rettig EM, Talbot CC, Sausen M, Jones S, Bishop JA, Wood LD, et al. Whole-genome sequencing of salivary gland adenoid cystic carcinoma. Cancer Prev Res (Phila). 2016;9:265–74.

    CAS  Google Scholar 

  12. Togashi Y, Dobashi A, Sakata S, Sato Y, Baba S, Seto A, et al. MYB and MYBL1 in adenoid cystic carcinoma: diversity in the mode of genomic rearrangement and transcripts. Mod Pathol. 2018;31:934–46.

    Article  CAS  Google Scholar 

  13. Varela I, Tarpey P, Raine K, Huang D, Ong CK, Stephens P, et al. Exome sequencing identifies frequent mutation of the SWI/SNF complex gene PBRM1 in renal carcinoma. Nature. 2011;469:539–42.

    Article  CAS  Google Scholar 

  14. Gui Y, Guo G, Huang Y, Hu X, Tang A, Gao S, et al. Frequent mutations of chromatin remodeling genes in transitional cell carcinoma of the bladder. Nat Genet. 2011;43:875–8.

    Article  CAS  Google Scholar 

  15. Wang K, Kan J, Yuen ST, Shi ST, Chu KM, Law S, et al. Exome sequencing identifies frequent mutation of ARID1A in molecular subtypes of gastric cancer. Nat Genet. 2011;43:1219–23.

    Article  CAS  Google Scholar 

  16. Wiegand KC, Shah SP, Al-Agha OM, Zhao Y, Tse K, Zeng T, et al. ARID1A mutations in endometriosis-associated ovarian carcinomas. N Engl J Med. 2010;363:1532–43.

    Article  CAS  Google Scholar 

  17. Pasqualucci L, Dominguez-Sola D, Chiarenza A, Fabbri G, Grunn A, Trifonov V, et al. Inactivating mutations of acetyltransferase genes in B-cell lymphoma. Nature. 2011;471:189–95.

    Article  CAS  Google Scholar 

  18. Andricovich J, Perkail S, Kai Y, Casasanta N, Peng W, Tzatsos A. Loss of KDM6A activates super-enhancers to induce gender-specific squamous-like pancreatic cancer and confers sensitivity to BET inhibitors. Cancer Cell. 2018;33:512–526.e8.

    Article  CAS  Google Scholar 

  19. Foss F, Coiffier B, Horwitz S, Pro B, Prince HM, Sokol L, et al. Tolerability to romidepsin in patients with relapsed/refractory T-cell lymphoma. Biomark Res. 2014;2:16.

    Article  Google Scholar 

Download references


We thank Emmy Yanagida and Hiroshi Yamada for providing technical assistance.


This study was supported by JSPS KAKENHI (Grant-in-Aid for Young Scientists B) grant 20 K18232, and Japan Agency for Medical Research and Development (AMED) grant JP 16kk0205006. AMED had no control over the interpretation, writing, or publication of this work.

Author information

Authors and Affiliations



Drafting of manuscript, KN.; molecular genetic studies, EA, MI and HH; clinical information, HN1; sequence alignment and data curation, ST; conceptualization, and participation in design and coordination of study and draft preparation, HN2. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Kohei Nakamura.

Ethics declarations

Ethics approval and consent to participate

Acquisition of tissue specimens and genomic analyses were approved by the Ethics Committee of Keio University, Tokyo, Japan.

Consent for publication

Written informed consent was obtained from the patients for publication of this report and accompanying images. A copy of the written consent is available for review by the editor of this journal.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Details of the KRAS point mutation. Figure S2. Details of the KDM6A alteration. Figure S3. (Copy-number alteration and variant allele frequency (VAF) in case 1. The horizontal axis corresponds to each examined gene and the vertical axis corresponds to (A) copy number or (B) VAF. Figure S4. Copy-number alteration and variant allele frequency (VAF) in case 2. The horizontal axis corresponds to each examined gene and the vertical axis corresponds to the (A) copy number or (B) VAF.

Additional file 2: Table S1.

Genes (160) examined in the PleSSision test

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nakamura, K., Aimono, E., Tanishima, S. et al. Genetic profiling of patients with adenoid cystic carcinoma of the Bartholin’s glands reveals potential new routes for targeted therapies: a case report. Diagn Pathol 15, 64 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: