- Research article
- Open Access
- Open Peer Review
Distinct molecular subtypes of uterine leiomyosarcoma respond differently to chemotherapy treatment
BMC Cancer volume 17, Article number: 639 (2017)
Uterine leiomyosarcoma (ULMS) is an aggressive form of soft tissue tumors. The molecular heterogeneity and pathogenesis of ULMS are not well understood.
Expression profiling data were used to determine the possibility and optimal number of ULMS molecular subtypes. Next, clinicopathological characters and molecular pathways were analyzed in each subtype to prospect the clinical applications and progression mechanisms of ULMS.
Two distinct molecular subtypes of ULMS were defined based on different gene expression signatures. Subtype I ULMS recapitulated low-grade ULMS, the gene expression pattern of which resembled normal smooth muscle cells, characterized by overexpression of smooth muscle function genes such as LMOD1, SLMAP, MYLK, MYH11. In contrast, subtype II ULMS recapitulated high-grade ULMS with higher tumor weight and invasion rate, and was characterized by overexpression of genes involved in the pathway of epithelial to mesenchymal transition and tumorigenesis, such as CDK6, MAPK13 and HOXA1.
We identified two distinct molecular subtypes of ULMS responding differently to chemotherapy treatment. Our findings provide a better understanding of ULMS intrinsic molecular subtypes, and will potentially facilitate the development of subtype-specific diagnosis biomarkers and therapy strategies for these tumors.
Uterine leiomyosarcoma (ULMS) is a type of malignant soft tissue tumors showing distinctive morphologic features and molecular signatures . ULMS has poor prognosis and high recurrence rate [2,3,4]. Currently, the treatment of ULMS is mostly by surgery with some adjuvant therapies, such as cytotoxic chemotherapy and radiotherapy [5,6,7]. Due to the complex molecular heterogeneity of ULMS and unavailability of targeted therapeutic methods, the five-year survival rate of ULMS is still low . By gene expression profiling methods, a number of malignant tumors, including breast cancer, gastric cancer, and uterine carcinosarcoma [9,10,11,12,13], have been categorized into different molecular subtypes. Based on the subtype information, patients may be able to receive better diagnosis and more effective therapeutic options . Thus, it is important to classify molecular subgroups of ULMS, which may provide a better understanding of disease mechanism and guide future precision treatment. Previously, we analyzed leiomyosarcoma (LMS) cases from uterine and extra-uterine sites to classify molecular subtypes . Nevertheless, the treatments for LMS were different regarding LMS locations, which were between uterine and extra-uterine LMS patients in clinical practice. Italiano et al. demonstrated the molecular heterogeneity of LMS from extra-uterus by genetic profiling . But the molecular heterogeneity of ULMS was less investigated until now.
In this study, by analyzing gene expression data sets, we identified and defined two molecular subtypes of ULMS, each of which presents subtype-specific gene expression patterns. Genes and pathways enriched in subtype I ULMS were associated with smooth muscle function, while genes and pathways involved in epithelial to mesenchymal transition (EMT) and tumorigenesis were enriched in subtype II ULMS. Our findings will provide a better understanding of ULMS pathogenesis and facilitate the development of more effective and individualized therapies.
Bioinformatic analyses and immunohistochemistry staining
To identify the molecular subtypes of ULMS, we collected and analyzed expression profile data sets from TCGA database and determined the optimal number of molecular subtypes of ULMS by Consensus Clustering (R package, ConsensusClusteringPlus ). The accuracy of subtype assignments was evaluated by Silhouette analysis (R package cluster ). The subtype-specific gene expression patterns were investigated by Gene Set Enrichment Analysis (GSEA) and Significance Analysis of Microarrays (SAM-seq). To identify the pathways that enriched in each subtype, the KEGG pathway analysis was performed online (https://david.ncifcrf.gov/). Cluster 3.0 and TreeView were used to perform hierarchical clustering to view the top 500 significantly over-expressed genes from each subtype. To develop subtype specific diagnostic biomarkers, genes overexpressed in subtype I (LMOD1, 1:20, Sigma, CAT#HPA028325) and subtype II (ARL4C, 1:120, Sigma, CAT#HPA028927) ULMS were selected for immunohistochemistry staining (IHC) based on SAM-seq result and the antibody availability. The procedure and scoring of IHC were performed as described previously .
Statistical significance was assessed by the chi-square and Fisher exact tests. For all statistical analyses, p value less than 0.05 was considered statistically significant.
Consensus clustering of gene expression profiles revealed two molecular subtypes of uterine leiomyosarcoma
Level 3 RNAseq expression data of 29 ULMS cases were collected from The Cancer Genome Atlas (TCGA) and used to determine the molecular heterogeneity of ULMS by consensus clustering (Fig. 1a), a method that estimates cluster stability by iterative resampling of genes and samples . The consensus clustering demonstrated that two subtypes were the optimal number for ULMS, as indicated by the empirical cumulative distribution plots, showing the greatest increase in the area under CDF curve (Additional file 1: Figure S1A and B). Next, the confidence of subtype assignment from Consensus Clustering was evaluated by silhouette analysis (Fig. 1b), which showed that all cases from both subtypes have a positive silhouette value, confirming the two molecular ULMS subtypes.
Clinicopathologic features of ULMS molecular subtypes
Next, we compared the clinicopathologic features between subtype I and subtype II ULMS patients. As shown in Table 1, the ULMS subtype is significantly associated with clinical treatment response. Specifically, subtype I patients were significantly more responded to chemotherapy treatment than subtype II. However, there is no significant association between molecular subtypes with other clinicopathologic characteristics, including tumor weight, metastasis status, invasion and necrosis (Table 1).
Distinct molecular subtypes of ULMS have different gene expression patterns
We next investigated the subtype-specific gene expression patterns of ULMS by Gene Set Enrichment Analysis (GSEA) . GSEA analysis showed that 2669 gene sets were enriched. Among the 2669 gene sets, 1568 gene sets were over-expressed in subtype II while the other 1101 gene sets were over-expressed in subtype I (Fig. 2a). Gene sets associated with leiomyosarcoma and myogenic targets were enriched in subtype I (Fig. 2b), whereas gene sets involved in EMT and tumorigenesis were enriched in subtype II (Fig. 2c). Interestingly, genes over-expressed in subtype II were also associated with cell cycle, proliferation, organ development and tumorigenesis. These genes include CDK6, BMP1, MAPK13, PDGFRL and HOXA1 (Fig. 3). Subtype I ULMS was enriched with genes involved in smooth muscle function (Fig. 3), including LMOD1, SLMAP, MYLK, and MYH11, all of which are the smooth muscle-specific markers [20,21,22].
To further explore the different gene expression patterns of the two subtypes of ULMS, we performed Significance Analysis by SAM-seq. Our results showed that 1947 genes were significantly differently expressed between two subtypes. Among these genes, 1050 were over-expressed in subtype II ULMS while 897 genes were over-expressed in subtype I. Next, we analyzed the top 500 over-expressed genes in each subtype by hierarchical clustering. As the heatmap shows, those genes were significantly over-expressed in subtype I and subtype II, respectively (Additional file 2: Figure S2). Consistent with the GSEA results, the pathways enriched in subtype I were associated with smooth muscle function, such as vascular smooth muscle contraction, calcium signaling pathway, and regulation of actin cytoskeleton (Table 2). Pathways involved in tumorigenesis associated with subtype II, such as pathways in cancer, TGF-β and Hedgehog signaling pathway (Table 2).
Immunohistochemistry staining of ULMS subtype specific biomarkers
To further validate ULMS subtypes and develop subtype specific diagnostic biomarkers, we selected ULMS subtype specific genes for IHC staining by SAM-seq result and availability of commercial antibodies. The IHC results showed that subtype I biomarker LMOD1 was positively stained in 51 of 68 ULMS cases (75%), while 13 subtype II ULMS cases were positive for ARL4C (13/68, 19%) (Fig. 4a), and the correlation coefficient of staining results of these two genes across all the ULMS cases is −0.43 (Fig. 4b).
Uterine sarcomas are composed of leiomyosarcoma, endometrial stromal sarcoma and carcinosarcoma. Among these, leiomyosarcoma is the most common subclass, mainly found in postmenopausal women [1, 23]. Although early diagnosis could improve the survival rate of ULMS patients, there are still challenges for treating late stage ULMS patients due to its high invasiveness and relatively high resistance to radiotherapy and chemotherapy . Molecular subtyping of tumors based on their gene expression profiling have guided subtype-specific diagnosis, prognosis, and aided to develop subtype targeted therapies . In our study, we identified two molecular subtypes of ULMS and found that these two subtypes exhibited significantly different gene expression patterns and distinct sensitivities to chemotherapy treatment. Nevertheless, please kindly noted that the number of patients with clinical treatment response information in this study is low, and additional large scale validation of treatment responses on distinct molecular subtypes will be needed.
Among the genes and pathways enriched, subtype I ULMS showed overexpression of smooth muscle-specific markers, including LMOD1, SLMAP, MYLK, MYH11. LMOD1, also known as Leiomodin-1, could be activated by serum response factor (SRF) or myocardin (MYOCD) and functions in smooth muscle cell differentiation . SLMAP, or sarcolemmal membrane-associated protein, is involved in microtubule organization , excitation-contraction coupling  and myoblast fusion . Myosin light chain kinase and Myosin-11 protein, encoded by MYLK and MYH11, respectively, are components of smooth muscle cells SMC contractile apparatus . MYLK is a Ca2+/CaM-dependent kinase and involved in smooth muscle contraction by promoting the interaction between myosin and actin filaments . Myosin-11 belongs to the myosin heavy chain family, and has contractile activity by hydrolyzing ATP .
Subtype II ULMS showed overexpression of genes enriched in pathways including cancer, TGF-β and Hedgehog signaling. Particularly, over-expression was found in the following genes, including CDK6, MAPK13 and HOXA1. CDK6 (Cyclin-dependent kinase 6) is a cell cycle regulator and forms a complex with cyclin D to initiate G1 to S phase transition by phosphorylating and inactivating Rb [31,32,33]. MAPK13 (Mitogen-activated protein kinase 13) belongs to the MAP kinase family and functions in cell proliferation, differentiation and development. MAPK13 is also involved in cell motility and invasion, serving as a diagnostic marker for cholangiocarcinoma . MAPK13 is highly expressed in uterine and ovary tumor tissues, especially in gynecological cancer stem cells, and has tumor-initiating activity that was involved in tumorigenesis . Highly expressed in many types of tumor cells, HOXA1 (Homeobox A1) is a DNA-binding protein and involved in facilitating cell proliferation, invasion, metastasis, and tumor progression [36, 37]. These proto-oncogenes were all over-expressed in subtype II, suggesting that the extent of malignancy of subtype II ULMS may be higher than that of subtype I, and subtype II may represent high-grade ULMS.
In conclusion, we characterized distinct intrinsic molecular subtypes of ULMS with different gene signatures. Our findings provide new insights into the understanding of the pathogenesis of ULMS, facilitate the development of subtype-specific diagnostic biomarkers and targeted treatment for ULMS. Furthermore, our finding may provide valuable information to develop individualized medicine for ULMS patients.
Cyclin-dependent kinase 6
Epithelial to mesenchymal transition
Gene Set Enrichment Analysis
- MAP K13:
Mitogen-activated protein kinase 13
Myosin light chain kinase
Significance Analysis of Microarrays
Sarcolemmal membrane-associated protein
Serum response factor
The Cancer Genome Atlas
Miyata T, Sonoda K, Tomikawa J, Tayama C, Okamura K, Maehara K, et al. Genomic, Epigenomic, and Transcriptomic profiling towards identifying Omics features and specific biomarkers that distinguish uterine Leiomyosarcoma and Leiomyoma at molecular levels. Sarcoma. 2015;2015:412068.
Gauthe M, Testart Dardel N, Nascimento C, Trassard M, Banal A, Alberini JL. Uterine leiomyosarcoma metastatic to thyroid shown by 18F-FDG PET/CT imaging. Revista espanola de medicina nuclear e imagen molecular. 2017;36(2):113–5.
Dandapani M, Seagle BL, Abdullah A, Hatfield B, Samuelson R, Shahabi S. Metastatic uterine Leiomyosarcoma involving bilateral ovarian Stroma without capsular involvement implies a local route of Hematogenous dissemination. Case Rep Obstet Gynecol. 2015;2015:950373.
Artioli G, Borgato L, Calamelli S, Azzarello G. Unusual cardiac metastasis of uterine leiomyosarcoma: case report and literature review. Tumori. 2016;102(Suppl. 2):90–2.
Bogani G, Fuca G, Maltese G, Ditto A, Martinelli F, Signorelli M, et al. Efficacy of adjuvant chemotherapy in early stage uterine leiomyosarcoma: a systematic review and meta-analysis. Gynecol Oncol. 2016;143(2):443–7.
Ducie JA, Leitao MM Jr. The role of adjuvant therapy in uterine leiomyosarcoma. Expert Rev Anticancer Ther. 2016;16(1):45–55.
Monti STP, Mesirov J, Golub T. Consensus clustering: a Resampling-based method for class discovery and visualization of gene expression microarray data. Mach Learn. 2003;52:91–118.
Zivanovic O, Jacks LM, Iasonos A, Leitao MM Jr, Soslow RA, Veras E, et al. A nomogram to predict postresection 5-year overall survival for patients with uterine leiomyosarcoma. Cancer. 2012;118(3):660–9.
Bertucci F, Finetti P, Rougemont J, Charafe-Jauffret E, Cervera N, Tarpin C, et al. Gene expression profiling identifies molecular subtypes of inflammatory breast cancer. Cancer Res. 2005;65(6):2170–8.
Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, et al. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001;98(19):10869–74.
Lei Z, Tan IB, Das K, Deng N, Zouridis H, Pattison S, et al. Identification of molecular subtypes of gastric cancer with different responses to PI3-kinase inhibitors and 5-fluorouracil. Gastroenterology. 2013;145(3):554–65.
Cristescu R, Lee J, Nebozhyn M, Kim KM, Ting JC, Wong SS, et al. Molecular analysis of gastric cancer identifies subtypes associated with distinct clinical outcomes. Nat Med. 2015;21(5):449–56.
An Y, Wang H, Jie J, Tang Y, Zhang W, Ji S, et al. Identification of distinct molecular subtypes of uterine carcinosarcoma. Oncotarget. 2017;8:15878.
Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thurlimann B, Senn HJ, et al. Strategies for subtypes--dealing with the diversity of breast cancer: highlights of the St. Gallen international expert consensus on the primary therapy of early breast cancer 2011. Ann Oncol. 2011;22(8):1736–47.
Guo X, Jo VY, Mills AM, Zhu SX, Lee CH, Espinosa I, et al. Clinically relevant molecular subtypes in Leiomyosarcoma. Clin Cancer Res. 2015;21(15):3501–11.
Italiano A, Lagarde P, Brulard C, Terrier P, Lae M, Marques B, et al. Genetic profiling identifies two classes of soft-tissue leiomyosarcomas with distinct clinical characteristics. Clin Cancer Res. 2013;19(5):1190–6.
Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010;26(12):1572–3.
Peter RJ. Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. J Comput Appl Math. 1987;20:53–65.
Hung JH, Yang TH, Hu Z, Weng Z, DeLisi C. Gene set enrichment analysis: performance evaluation and usage guidelines. Brief Bioinform. 2012;13(3):281–91.
Nanda V, Miano JM. Leiomodin 1, a new serum response factor-dependent target gene expressed preferentially in differentiated smooth muscle cells. J Biol Chem. 2012;287(4):2459–67.
Yeung KK, Bogunovic N, Keekstra N, Beunders AA, Pals J, van der Kuij K, et al. Transdifferentiation of human dermal fibroblasts to smooth muscle-like cells to study the effect of MYH11 and ACTA2 mutations in aortic aneurysms. Hum Mutat. 2017;38:439.
Okumu LA, Bruinton S, Braden TD, Simon L, Goyal HO. Estrogen-induced maldevelopment of the penis involves down-regulation of myosin heavy chain 11 (MYH11) expression, a biomarker for smooth muscle cell differentiation. Biol Reprod. 2012;87(5):109.
Skorstad M, Kent A, Lieng M. Uterine leiomyosarcoma - incidence, treatment, and the impact of morcellation. A nationwide cohort study. Acta Obstet Gynecol Scand. 2016;95(9):984–90.
Samarnthai N, Hall K, Yeh IT. Molecular profiling of endometrial malignancies. Obstet Gynecol Int. 2010;2010:162363.
Guzzo RM, Sevinc S, Salih M, Tuana BS. A novel isoform of sarcolemmal membrane-associated protein (SLMAP) is a component of the microtubule organizing centre. J Cell Sci. 2004;117(Pt 11):2271–81.
Nader M, Westendorp B, Hawari O, Salih M, Stewart AF, Leenen FH, et al. Tail-anchored membrane protein SLMAP is a novel regulator of cardiac function at the sarcoplasmic reticulum. Am J Physiol Heart Circ Physiol. 2012;302(5):H1138–45.
Guzzo RM, Wigle J, Salih M, Moore ED, Tuana BS. Regulated expression and temporal induction of the tail-anchored sarcolemmal-membrane-associated protein is critical for myoblast fusion. Biochem J. 2004;381(Pt 3):599–608.
Renard M, Callewaert B, Baetens M, Campens L, MacDermot K, Fryns JP, et al. Novel MYH11 and ACTA2 mutations reveal a role for enhanced TGFbeta signaling in FTAAD. Int J Cardiol. 2013;165(2):314–21.
Yin F, Hoggatt AM, Zhou J, Herring BP. 130-kDa smooth muscle myosin light chain kinase is transcribed from a CArG-dependent, internal promoter within the mouse mylk gene. Am J Physiol Cell Physiol. 2006;290(6):C1599–609.
Kuang SQ, Kwartler CS, Byanova KL, Pham J, Gong L, Prakash SK, et al. Rare, nonsynonymous variant in the smooth muscle-specific isoform of myosin heavy chain, MYH11, R247C, alters force generation in the aorta and phenotype of smooth muscle cells. Circ Res. 2012;110(11):1411–22.
Meyerson M, Harlow E. Identification of G1 kinase activity for cdk6, a novel cyclin D partner. Mol Cell Biol. 1994;14(3):2077–86.
Bertoli C, Skotheim JM, de Bruin RA. Control of cell cycle transcription during G1 and S phases. Nat Rev Mol Cell Biol. 2013;14(8):518–28.
Ezhevsky SA, Ho A, Becker-Hapak M, Davis PK, Dowdy SF. Differential regulation of retinoblastoma tumor suppressor protein by G(1) cyclin-dependent kinase complexes in vivo. Mol Cell Biol. 2001;21(14):4773–84.
Tan FL, Ooi A, Huang D, Wong JC, Qian CN, Chao C, et al. p38delta/MAPK13 as a diagnostic marker for cholangiocarcinoma and its involvement in cell motility and invasion. Int J Cancer. 2010;126(10):2353–61.
Yasuda K, Hirohashi Y, Kuroda T, Takaya A, Kubo T, Kanaseki T, et al. MAPK13 is preferentially expressed in gynecological cancer stem cells and has a role in the tumor-initiation. Biochem Biophys Res Commun. 2016;472(4):643–7.
Yuan C, Zhu X, Han Y, Song C, Liu C, Lu S, et al. Elevated HOXA1 expression correlates with accelerated tumor cell proliferation and poor prognosis in gastric cancer partly via cyclin D1. J Exp Clin Cancer Res. 2016;35:15.
Wang H, Liu G, Shen D, Ye H, Huang J, Jiao L, et al. HOXA1 enhances the cell proliferation, invasion and metastasis of prostate cancer cells. Oncol Rep. 2015;34(3):1203–10.
We thank Feng Selina Ji, who read and improved the manuscript.
This work was supported by the National Natural Science Foundation of China (No. 81602362, to XG), the program for Science and Technology Development in Henan Province (No.162102310391, to XG), the program for Young Key Teacher of Henan Province (No.2016GGJS-214, to XG), the key program for Science and Technology Research of Henan Department of Education (No.15A310001, to YA), the supporting grants of Henan University (No.2015YBZR048, to XG; No.B2015151, to XG; No.2015YBZR046, to YA), the program for Innovative Talents of Science and Technology in Henan Province (No. 18HASTIT048, to XG), and Yellow River Scholar Program (No.H2016012, to XG). The funding bodies were not involved in the study design, data collection, analysis and interpretation of data, or in writing of this manuscript.
Availability of data and materials
The Level 3 RNAseq expression data of uterine leimyosarcoma cases were obtained from The Cancer Genome Atlas (TCGA) database.
Ethics approval and consent to participate
Paraffin blocks of uterine leiomyosarcomas were obtained with Health Insurance Portability and Accountability Act (HIPAA)-compliant Stanford University Medical Center Institutional Review Board approval and a waiver of consent due to the archival nature of the specimens.
Consent for publication
The authors confirm that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Delineation of two distinct molecular subtypes of ULMS. (A) Empirical cumulative distribution plots. (B) The increased area under the CDF curve along with increased number of molecular subtypes. (TIFF 215 kb)
Heatmap of the top 500 genes over-expressed in distinct molecular subtypes. Based on SAM-seq result, the top 500 genes over-expressed in each subtype were selected to form the TOP500 genes. Hierarchical clustering of the TOP500 genes was performed by Cluster 3.0 using centroid linkage method. Each row denotes a gene and each column corresponds to a case of ULMS. Red, over-expressed genes; Green, down-expressed genes. (TIFF 2173 kb)