William John Hatchett1
Marta Brunetti1
Kristin Andersen1
Maren Randi Tandsæther1
Ingvild Lobmaier2Marius Lund-Iversen2Thomas Lien-Dahl2Francesca Micci1Ioannis Panagopoulos1*- 1Section for Cancer Cytogenetics, Institute for Cancer Genetics and Informatics, The Norwegian Radium Hospital, Oslo University Hospital, Oslo, Norway
- 2Department of Pathology, The Norwegian Radium Hospital, Oslo University Hospital, Oslo, Norway
Introduction: Intramuscular myxomas are benign tumors that are challenging to diagnose, especially on core needle biopsies. Acquired chromosomal aberrations and pathogenic variants in codon 201 or codon 227 in GNAS complex locus gene (GNAS) have been reported in these tumors. Here we present our genetic findings in a series of 22 intramuscular myxomas.
Materials and methods: The tumors were investigated for the presence of acquired chromosomal aberrations using G-banding and karyotyping. Pathogenic variants in codon 201 or codon 227 of GNAS were assessed using direct cycle Sanger sequencing and Ion AmpliSeq Cancer Hotspot Panel v2 methodologies.
Results: Eleven tumors carried chromosomal abnormalities. Six tumors had numerical, four had structural, and one had both numerical and structural chromosomal aberrations. Gains of chromosomes 7 and 8 were the most common abnormalities being found in five and four tumors respectively. Pathogenic variants in GNAS were detected in 19 myxomas (86%) with both methodologies. The detected pathogenic variants were p.R201H in nine cases (seven with abnormal and two with normal karyotypes), p.R201C in five cases, all with normal karyotypes, p.R201S in three cases (two with abnormal and one with normal karyotype), p.R201G in one case with a normal karyotype, and p.Q227E in one case with a normal karyotype.
Conclusion: Firstly, our data indicate a possible association between chromosomal abnormalities and GNAS pathogenic variants in intramuscular myxomas. Secondly, the presence of the rare pathogenic variants R201S, p.R201G and p.Q227E in 26% (5 out of 19) of myxomas with GNAS pathogenic variants shows that methodologies designed to detect only the common “hotspot” of p.R201C and p.R201H will give false negative results. Finally, a comparison between Ion AmpliSeq Cancer Hotspot Panel v2 and direct cycle Sanger sequencing showed that direct cycle Sanger sequencing provides a quick, reliable, and relatively cheap method to detect GNAS pathogenic variants, matching even the most cutting-edge sequencing methods.
Introduction
Intramuscular myxomas are rare benign soft tissue tumors, characterized by stellate and/or spindle-shaped cells embedded in a hypo-vascular, abundant myxoid stroma. They were first pathologically described by Enzinginer in 1965 [1] and the description has since been validated in later studies [2–6]. The reported annual incidence rate of intramuscular myxoma is 0.10–0.13/100,000 and mostly occurs in adults, with women having a predisposition for the condition (∼60–70%) [7–9]. Intramuscular myxoma is a somatic disorder with the majority of patients presenting asymptomatically with painless, palpable, well-defined tumors that slowly enlarge to around 2–15 cm in diameter [9]. They are mostly located in the thigh, and do not recur after a simple excision [8]. Multiple intramuscular myxomas are rare, generally occurring as an isolated lesion not associated with other clinical abnormalities [8]. However, multiple intramuscular myxomas in combination with polyostotic fibrous dysplasia can be associated with Mazabraud’s syndrome [7, 9]. Intramuscular myxoma, Mazabraud’s syndrome and the closely related McCune-Albright syndrome (also associated with fibrous dysplasia, café au lait macules and endocrine disorders) are mostly caused by the activating missense mutations in codon 201 or codon 227 of the GNAS complex locus gene (GNAS) [10–14].
Pathogenic variants in codon 201 of the GNAS gene on chromosome 20 are one of the primary cancer-causing variants in heterotrimeric G proteins and is one of the leading causes of oncogenesis in a variety of low-grade malignant and benign tumors [15]. The GNAS gene encodes the stimulatory G-alpha subunit of the heterotrimeric G-protein complex, which regulates the activation of adenylyl cyclase that converts adenosine triphosphate into cyclic adenosine monophosphate. When a cell has a GNAS pathogenic variant, there is an overproduction of cyclic adenosine monophosphate and activation of downstream signaling pathways [16–22]. The first description of a somatic GNAS pathogenic variant in intramuscular myxomas was published by Okamoto et al. in 2000 [10]. Thereafter, many studies have shown GNAS pathogenic variants occurring frequently in sporadic intramuscular myxomas [10, 18–22]. Additionally, GNAS pathogenic variants are absent in low-grade myxofibrosarcoma, which can be used as a malignant differential diagnosis to intramuscular myxomas [23, 24]. A precise diagnosis between these two is essential for the correct treatment; low-grade myxofibrosarcomas commonly recur and require wide excisions, whereas a more precise or marginal excision is sufficient to treat intramuscular myxomas [24, 25].
Intramuscular myxomas are challenging to diagnose, especially on core needle biopsies [26]. Also, they have a broad differential diagnosis at a morphological level and no clear karyotypic diagnosis [22, 27]. Additionally, there are many different types of GNAS pathogenic variants found in intramuscular myxomas, located on either codon 201, in exon 8, or codon 227, in exon 9. GNAS also has a low range of pathogenic allele frequencies in intramuscular myxomas, ranging between ∼ 5–30%. The most common pathogenic variants are p.R201H (c.602G>A) and p.R201C (c.601C>T), making up approximately 80%–90% of all cases [10, 19, 20, 22]. This has influenced the design of some diagnostic methods to only detect common pathogenic variants, and not the rarer ones [28, 29]. Depending on the method, GNAS pathogenic variants were found in ∼30–80% of intramuscular myxomas [20–22].
The gold standard for sensitive screening for genes that can predispose a patient to particular cancer types, such as pathogenic variants in GNAS, was PCR amplification, followed by sequencing, e.g., Sanger sequencing (also known as dideoxy chain termination method) of the amplified product [30–34]. However, in the last two decades, next-generation sequencing (NGS, also known as high throughput sequencing and massive parallel sequencing) has become common practice in cancer genetic screening [35]. Methods such as whole genome sequencing, exome sequencing, and transcriptome sequencing can provide more sensitive and accurate results. The main difference between Sanger sequencing and NGS is the sequencing volume. The Sanger methodology can only sequence a single DNA fragment at a time, whereas NGS simultaneously sequences millions of fragments per run. Thus, NGS results in hundreds to thousands of genes being sequenced simultaneously. Additionally, NGS can start with relatively low input DNA compared to that of Sanger methodology and provides enhanced detection capabilities to detect novel or rare variants through deep sequencing. Comparing the results obtained from NGS and Sanger sequencing there is close to 100% concordance between these methodologies [36–38]. NGS-based gene panels have also been developed to target specific genes or pathogenic variants (hotspots) that are associated to a particular cancer type [39–41]. These panels are nowadays used in clinical cancer genetics, providing accurate information on oncogenic drivers and actionable genetic alterations for a variety of genes [39–41]. The commercially available Ion AmpliSeq Cancer Hotspot Panel v2 (ThermoFisher Scientific, Waltham, MA, United States) covers approximately 2,800 pathogenic variants which are found in the Catalogue Of Somatic Mutations In Cancer (COSMIC) [42]. The pathogenic variants detected by Ion AmpliSeq Cancer Hotspot Panel v2 are from 50 oncogenes and tumor suppressor genes, including those at codons 201 and 227 of GNAS.
In the present study, we cytogenetically analyzed intramuscular myxomas and compared a PCR based/Sanger sequencing methodology with the above-mentioned panel for the detection of GNAS pathogenic variants in codons 201 and 227.
Materials and methods
Patients
Our materials consisted of 22 intramuscular myxomas (Table 1) and all samples were collected at the Norwegian Radium Hospital between 2015 and 2023. The patients were 17 females and five males between the ages of 39 and 80. The study was approved by the Regional Committees for Medical Research Ethics South East Norway.
TABLE 1. Clinicopathological data, karyotypes, and status of pathogenic variants of the GNAS gene on the intramuscular myxomas.
G-banding and karyotyping
A representative tumor area was investigated cytogenetically as previously described [27, 43]. Material for cytogenetic examination was available from all 22 samples. After mechanical and enzymatic disaggregation of the tissue sample, the resulting cells were short-term cultured, harvested, and processed for cytogenetic examination. To obtain G-banding of chromosomes, Wright’s stain was used (Sigma Aldrich; St Louis, MO, United States). The subsequent cytogenetic analysis and karyotype description followed the recommendations of the International System for Human Cytogenomic Nomenclature (ISCN) 2020 guidelines [44].
DNA extraction from tumor and FFPE tissue samples
Frozen (−80°C) tumor tissue adjacent to that used for cytogenetic analysis and histologic examination was available for 21 myxomas (Table 1; Supplementary Table S1) and formalin-fixed, paraffin-embedded (FFPE) sample was used for case 14 (Supplementary Table S1). DNA was extracted from the frozen tumor samples, using the Maxwell RSC Tissue DNA Kit and the Maxwell RSC Instrument (Promega, Madison, WI, United States). For FFPE sample, DNA was extracted using the QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany). The DNA concentration was then estimated for each sample, using the QuantiFluor ONE dsDNA System using Quantus™ Fluorometer (Promega).
Sequencing of GNAS pathogenic variants using Sanger sequencing
The primers used for PCR amplification and Sanger sequencing are listed in Table 2. The BigDye Direct Cycle Sequencing Kit was used for PCR/cycle (Sanger) sequencing according to the company’s recommendations (ThermoFisher Scientific, Waltham, MA, United States). The primer combination M13For-GNASint7-F1 and M13Rev-GNASint8-R1 was used for amplification and sequencing of exon 8 of GNAS to detect possible pathogenic variants at the codon 201 (p.R201) (Table 2). The primer combination M13For-GNASint8-F1 and M13Rev-GNASint9-R1 was used for amplification and sequencing of exon 9 of GNAS to detect possible pathogenic variants at codon 227 (p.Q227) (Table 2). Sequencing was run on the Applied Biosystems SeqStudio Genetic Analyzer system (ThermoFisher Scientific). For analysis of sequence data, the basic local alignment search tool (BLAST) software1 was used [45]. The obtained sequences were aligned against the GNAS reference sequences NM_000516.7 and NG_016194.2 corresponding to GNAS complex locus (GNAS) transcript variant 1 and RefSeqGene on chromosome 20, respectively. Codons 201 (p.R201) and 227 (p.Q227) are based on the NCBI reference sequence NP_000507.1 which corresponds to protein GNAS isoform GNASL.
TABLE 2. Designation, sequence (5′->3′), and position in reference sequence of the forward (F) and reverse (R) primers used for BigDyeTM Direct Cycle Sequencing of exons 8 and 9 of the GNAS gene.
Sequencing of GNAS pathogenic variants using Ion AmpliSeq Cancer Hotspot Panel v2
The Ion Ampliseq™ Cancer Hotspot Panel v2 Chef-Ready Kit was used to detect GNAS pathogenic variants following the company’s recommendations (ThermoFisher Scientific). DNA samples were diluted to between 0.4 and 1.6 ng/μL or used as is when the sample concentration was below 0.6 ng/μL. All runs on the Ion Chief and Ion torrent were created in TorrenSuite™ (ThermoFisher Scientific). Libraries were prepared using Ion Ampliseq™ Cancer Hotspot Research Panel v2 primers and by using the Ion AmpliSeq™ Kit for Chef DL8 and were run using recommended setting from the company (ThermoFisher Scientific). All samples were templated via the automated template preparation protocol Ion 540™ Kit—Chef (ThermoFisher Scientific), which was then sequenced on the Ion GeneStudio S5 following the company’s recommendations (ThermoFisher Scientific).
Results
Cytogenetics
The cytogenetic data are presented in Table 1. Normal karyotypes were found in 11 intramuscular myxomas and were abnormal in 11 tumors. Six abnormal karyotypes had numerical chromosomal aberrations, involving one or more of chromosomes X, 7, 8, 9, 20, and four tumors had structural aberration, and one had both numerical and structural aberrations (Table 1). Gains of chromosome 7 were the most frequent abnormalities found in five myxomas, followed by gains of chromosome 8, which were found in four tumors. No recurrent structural aberration was found (Table 1). The structural aberrations were t (6; 13) (p23; q14), inv (12) (q13q21), t (10; 22) (q24; q13), add (15) (p11∼13), and add (8) (p?11) (Figure 1; Table 1).
FIGURE 1. Cytogenetic analysis of myxomas 2, 3, and 4. Partial Karyotypes showing the abnormal chromosomes together with the corresponding normal chromosome homologs. For myxoma 2, der (6) and der (13), together with chromosomes 6 and 13 are shown. For myxoma 3, inv (12) together with normal 12 are shown. For myxoma 4, der (10) and der (22), together with chromosomes 10 and 22 are shown.
Sanger sequencing
DNA was successfully extracted from all tumor samples. Pathogenic variants in the GNAS gene were detected in 19 out of 22 myxomas (86%) (Supplementary Table S1). The detected pathogenic variants were p.R201H in nine cases which were found in seven myxomas with abnormal karyotypes and two with a normal karyotype, p.R201C in five cases all of them with normal karyotypes, p.R201S in three cases which two of them had abnormal karyotype and one had a normal karyotype, p.R201G in one case with a normal karyotype, and p.Q227E in one case with an abnormal karyotype (Table 1; Figure 2).
FIGURE 2. The detected GNAS pathogenic variants. Partial chromatograms showing the GNAS pathogenic variants in cases 10, 11, 13, 16, 18, and 19. Chromatograms from both forward and reverse sequencing are shown. The frequencies of mutated alleles, found with the Ion Ampliseq™ Cancer Hotspot Panel v2, are also shown.
Ion AmpliSeq Cancer Hotspot Panel v2
The Ion torrent results were the same as that of Sanger sequencing (Supplementary Table S1) and are therefore both represented in Table 1 for each case. Pathogenic variants in the GNAS gene were detected in the same 19 samples in which were detected by Sanger sequencing. Ion torrent also provided the allele frequencies of the mutated alleles ranging from 4.3%–33%, with an average of 13% (Table 1).
Discussion
Over 40,000 new GNAS p.R201 pathogenic variant cases were reported in 2021 in the United States alone, illustrating the high prevalence of this genetic abnormality and importance of detecting this genetic abnormality [15]. Diagnosing myxoma is difficult as there are no useful immunohistochemical markers to distinguish intramuscular myxomas from low-grade myxofibrosarcomas, both have an abundant myxoid matrix, and both may express the CD34 molecule [22, 46, 47]. Detection of acquired genetic abnormalities in combination with histopathology, can be a powerful diagnostic tool to help distinguish between morphologically similar neoplasms, such as intramuscular myxomas and low-grade myxofibrosarcomas [18, 23, 24]. Morphologically, myxofibrosacomas have more prominent blood vessels and nuclear atypia, but cellular myxomas share similar features, so morphology alone is not always sufficient to make a definite diagnosis, especially in small biopsies. Additionally, small biopsies may also create issues with pathogenic variants detection which can be hampered by low frequency of pathogenic variants alleles, i.e., sequences of mutated alleles “diluted” by wild-type sequences. For pathogenic variant detection it is recommended to analyses at least 50% lesioned tissue. This may require painstaking and time-consuming micro-dissection of the relevant tissue and can be difficult to achieve in samples from needle biopsies [18, 21]. It has been suggested that a minimum of 10% tumor cell fraction is required to detect GNAS pathogenic variants [22], but the tumor cell fraction may have not been estimated accurately [48, 49]. This is especially difficult in intramuscular myxoma tumors, even when using fully excised tumors, as there is low amounts of DNA available for extraction and a generally low level of pathogenic variants which could easily become diluted by wild-type sequences [18].
Another complication is that there are no pathognomonic chromosomal aberrations for intermuscular myxomas [27]. The abnormal karyotypes, which were found in 11 out of 22 myxomas (Table 1), confirmed the result and conclusion of our previous study [27]. Abnormalities of chromosomes 7 and 8, both mainly trisomy, were the most common aberrations. These abnormalities lack specificity and are found in a variety of different neoplasms [50–55] and trisomy 7 is also seen in non-neoplastic lesions [27, 56] Furthermore, the observed structural chromosomal aberrations were not recurrent (Table 1; [27]). Even though “acquisition of clonal chromosome aberrations is an integral part of the disease process” in intermuscular myxomas, the cytogenetic information cannot be used for the diagnosis of myxomas (Table 1, [27]).
A useful tool for the diagnosis of intramuscular myxomas is that intramuscular myxomas often carry a somatic GNAS pathogenic variants. Since GNAS pathogenic variants are absent in low-grade myxofibrosarcoma, detection of these pathogenic variants can be a useful differential diagnosis [10, 18–24]. Detecting pathogenic variants within the GNAS gene in intramuscular myxoma can be difficult due to the poor cellularity of the tumor and the low frequency of mutated alleles, which was found to be between 5% and 30%, [3, 18, 20–22]. A further complication for GNAS pathogenic variant detection is that PCR-based methods can produce varying results. Okamato et al. (2000) detected GNAS pathogenic variants in five out of six (83%) intramuscular myxomas (60% p.R201H and 40% p.R201C) using PCR methodology together with single strand conformation polymorphism [10, 57]. Delaney et al. (2009) detected pathogenic variants in eight out of 28 (29%) intramuscular myxomas using PCR followed by pathogenic variant-specific restriction enzyme digestion [18]. However, they found GNAS pathogenic variants in 17 out of 28 (61%) intramuscular myxomas with COLD-PCR followed by pathogenic variant-specific restriction enzyme digestion [18]. Walther et al. (2014) detected pathogenic variants in 23 out of 63 (36%) intramuscular myxomas (52% p.R201C and 48% p.R201H) using PCR and direct sequencing [19]. In the present work, we detected GNAS pathogenic variants in 19 out of 22 samples (86%) using direct cycle Sanger sequencing with primers tailed on their 5′-end with M13 universal primer sequences (Table 2). We detected a total of five different pathogenic variants: p.R201H (47.3%), p.R201C (26.3%), p.R201S (16%), p.Q227E (5.2%), and the rare pathogenic variant p.R201G (5.2%) (Table 1; Figure 2). The p.R201G pathogenic variants has so far only been reported in one patient with McCune-Albright Syndrome over 20 years ago [58]. One complication when using Sanger sequencing is that it can be hard to detect allele frequencies below 15%. Intramuscular myxomas tumors have low allele frequencies of GNAS pathogenic variants and in our investigation, we found an average allele frequency of 13% (Table 1). To combat this a positive result from both bidirectional Sanger sequencing is mandatory, and reruns are used to confirm if only one direction is positive. Additionally, Sanger sequencing cannot accurately determine the allele frequency and is only used as a qualitative method (Yes/No) whereas NGS can be used as both qualitative and quantitative (Yes/No and frequency of various alleles) (Table 1; Figure 2).
NGS has also been used to detect GNAS pathogenic variants in intramuscular myxomas [20–22]. Sunitsch et al. (2018) detected GNAS pathogenic variants in 12 out of 13 (92%) specimens (four cases with p.R201C, six cases with p.R201H and two cases with p.Q227E) using Ion Ampliseq primer panel and Ion torrent sequencing [21]. Bekers et al. (2019) detected GNAS pathogenic variants in 16 out of 28 (57%) samples with frequency of the mutant allele between 5% and 27% using single molecule tagged molecular inversion probe assay and Illumina NextSeq sequencing [20]. In addition to p.R201H and p.R201C, which were detected in 44% and 31% of the samples respectively, four other GNAS pathogenic variants were found in intramuscular myxomas; p.R201S, p.R201L, p.R201P and p.Q227R [20]. Libbrecht et al. (2019) detected eight out of nine (89%) intramuscular myxomas to have a GNAS pathogenic variants (62.5% p.R201H and 37.5% p. R201C) with frequency of the mutant allele ranging between 5% and 28% using TruSight Tumor 26 panel [22]. In the present study, we detected five different GNAS pathogenic variants in 19 out of 22 samples 83%, with frequency of the mutated allele to ranging between 4.3% and 28% (Table 1) using the Ion Ampliseq™ Cancer Hotspot Panel v2. Our data obtained with the NGS panel were similar to the above-mentioned studies. We also found that results obtained from the NGS panel were identical to those obtained by direct cycle Sanger sequencing, indicating that the latter methodology is a highly sensitive technique to detect GNAS pathogenic variants. The present study (Table 1) together with the studies published by Bekers et al. (2019) and Sunitsch et al. (2019) showed that approximately 23% of myxomas with GNAS pathogenic variants carried variants in other locations other than the “hotspot” variants p.R201C and p.R201H. Thus, assays targeting only these pathogenic variants, p.R201C and p.R201H, such as pyrosequencing, quantitative real-time PCR, and digital PCR should be used with caution for the detection of GNAS pathogenic variants in intramuscular myxomas [20, 21, 59, 60].
Even though cytogenetic information cannot be used for the diagnosis of myxomas, we found an association between karyotypes and specific GNAS pathogenic variant. The pathogenic variant p.R201H was detected in seven tumors carrying abnormal karyotype and in only two tumors with normal karyotype. The detected p.R201C pathogenic variants were found in five tumors with only normal karyotype (Table 1). To the best of our knowledge, this association between cytogenetic and molecular aberrations in intramuscular myxomas is reported for the first time in the present study. Although the number of cases in this report make up to few to draw any conclusions, our data may imply functional variations between the different pathogenic variants. Very little is known about possible functional differences between GNAS pathogenic variants. All constitutive pathogenic variants activate Gsα adenylate cyclase with overproduction of cyclic adenosine monophosphate (cAMP) and are therefore often treated as one and the same [61–63]. GNAS pathogenic variants are associated with upregulation of Wnt/beta-catening signaling [64–66] and were found in many tumors [15, 67]. In murine models, either p.R201C or p. R201H resulted in fibrous dysplasia [64, 68]. In patients with fibrous dysplasia/McCune-Albright syndrome carrying p.R201H or p.R201C, no clear genotype-phenotype correlation was found [69, 70]. Additionally, in patients with fibrous dysplasia no association was found between age, site, size, specimen type and GNAS mutational status [71]. A pan-cancer cohort of patients from 1050 GNAS mutant tumors showed a heavy disposition to the p.R201 codon, with a classic gain-of-function mutation. Both p.R201H and p.R201C pathogenic variants were found to drive tumor cell growth both in vitro and in vivo [15]. GNAS was overexpressed in the LS174T cell line using a doxycycline-inducible promoter, and the resulting effect on clonogenic capacity was assessed [15]. In another study, the basal and maximal adenylyl cyclase activity (cAMP accumulation) dose response under isoproterenol stimulation was found to be higher in p.R201H than that of p.R201C [72]. However, the authors have suggested that the GNAS expression was not stringently controlled, and this result could be a technical issue. This warranted further experimentation with both pathogenic variants and wildtype alleles to further elucidate functional variation like our karyotype results suggest.
It has been suggested that there can be notable variability found in phenotypes, pathogenicity and oncogenic effects, based on the specific substitution at a single location [73–75]. For example, the pathogenic missense variants which alter codon 336 in the GARS1 gene that results in different phenotypic expression in a range of genetic neuropathies [73]. Two pathogenic variants resulting in a missense modifying amino acid at codon 336 in the catalytic domain of GARS1, were found in two unrelated patients - one a female with infantile spinal muscular atrophy; and the second, a male with Charcot–Marie–Tooth disease type 2D [73]. Exchanges in amino acids that change the pH of the residue have also been linked to cancers [75]. White et al. (2017) found that the pathogenic variants resulting in a change from arginine residue to a histidine residue produced a rise in intracellular pH, which in turn, conferred these mutants with oncogenic effects [75]. This is the same change in amino acid residue we see in the amino acid pathogenic variant p.R201H which could suggest functional variations between the different pathogenic variants in GNAS.
In conclusion, acquired genomic abnormalities, both at chromosomal and molecular levels, are found in intramuscular myxomas. However, the detected chromosomal aberrations are not pathognomonic for myxomas since they lack specificity and are found in various neoplasms or even in non-neoplastic cells. In our experiment, we reliably detected somatic GNAS pathogenic variants at the molecular level. Together with H&E stain morphology and immunohistochemistry, we can reliably diagnose and detect the majority of intramuscular myxomas, which in turn leads to the appropriate treatment. The two methods used, i.e., direct cycle Sanger sequencing (BigDye Direct Cycle Sequencing Kit) and NGS panel (Ion AmpliSeq Cancer Hotspot Panel v2) gave identical results detected GNAS pathogenic variant in 83% of the examined myxomas. In 26% of the myxomas which carried GNAS pathogenic variants, both methods detected pathogenic variants other than the common “hotspot” of p.R201C and p.R201H. This shows that methods designed to detect only “hotspot” variants might give false negative results. Direct cycle Sanger sequencing is a quick, reliable, and relatively cheap method to detect the GNAS pathogenic variants, matching even the most cutting-edge of sequencing methods while only lacking the detection of frequency of mutant allele, not required for diagnostic purposes. Ultimately, this will allow laboratories that do not have the funding for expensive NGS methods to accurately diagnose GNAS pathogenic variants. This will find application not only in intramuscular myxomas but other tumor types as well.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The studies involving humans were approved by the Regional Committees for Medical Research Ethics South East Norway. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
IP planned and designed the experiment, WH, MB, KA, MT, and IP conducted genetic/cytogenetic diagnostics and IL, ML-I, and TL-D conducted pathological examinations. WH, MB, KA, and IP did the experiment, WH and IP did data analysis, and wrote the manuscript. Finally, WH, MB, KA, MT, IL, ML-I, TL-D, FM, and IP edited the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.por-journal.com/articles/10.3389/pore.2024.1611553/full#supplementary-material
Footnotes
1https://blast.ncbi.nlm.nih.gov/Blast.cgi
References
1. Enzinger, FM. Intramuscular myxoma: a review and follow-up study of 34 cases. Am J Clin Pathol (1965) 43:104–13. doi:10.1093/ajcp/43.2.104
2. Miettinen, M, and Markku, M. Modern soft tissue pathology: tumors and non-neoplastic conditions. Cambridge University Press (2016). p. 1062. Available from: https://www.cambridge.org/core/product/716C3345EDB53691CEE30DB640C84E7E (Accessed June 13, 2023).
3. Nielsen, GP, O’Connell, JX, and Rosenberg, AE. Intramuscular myxoma: a clinicopathologic study of 51 cases with emphasis on hypercellular and hypervascular variants. Am J Surg Pathol (1998) 22:1222–7. doi:10.1097/00000478-199810000-00007
4. Kindblom, L-G, Stener, B, and Angervall, L. Intramuscular myxoma. Cancer (1974) 34:1737–44. doi:10.1002/1097-0142(197411)34:5<1737::AID-CNCR2820340523>3.0.CO;2-U
5. Miettinen, M, Höckerstedt, K, Reitamo, J, and Tötterman, S. Intramuscular myxoma—a clinicopathological study of twenty-three cases. Am J Clin Pathol (1985) 84:265–72. doi:10.1093/ajcp/84.3.265
6. Hashimoto, H, Tsuneyoshi, M, Daimaru, Y, Enjoji, M, and Shinohara, N. Intramuscular myxoma: a clinicopathologic, immunohistochemical, and electron microscopic study. Cancer (1986) 58:740–7. doi:10.1002/1097-0142(19860801)58:3<740::AID-CNCR2820580322>3.0.CO;2-K
7. Heymans, O, Gebhart, M, Alexiou, J, de Saint Aubain, N, and Larsimont, D. Intramuscular myxoma. Acta Chir Belg (1998) 98:120–2. doi:10.1080/00015458.1998.12098394
8. Allen, PW. Myxoma is not a single entity: a review of the concept of myxoma. Ann Diagn Pathol (2000) 4:99–123. doi:10.1016/S1092-9134(00)90019-4
9. Murphey, MD, McRae, GA, Fanburg-Smith, JC, Temple, HT, Levine, AM, and Aboulafia, AJ. Imaging of soft-tissue myxoma with emphasis on CT and MR and comparison of radiologic and pathologic findings. Radiology (2002) 225:215–24. doi:10.1148/radiol.2251011627
10. Okamoto, S, Hisaoka, M, Ushijima, M, Nakahara, S, Toyoshima, S, and Hashimoto, H. Activating Gs(alpha) mutation in intramuscular myxomas with and without fibrous dysplasia of bone. Virchows Archiv (2000) 437:133–7. doi:10.1007/s004280000217
11. Schwindinger, WF, Francomano, CA, and Levine, MA. Identification of a mutation in the gene encoding the alpha subunit of the stimulatory G protein of adenylyl cyclase in McCune-Albright syndrome. Proc Natl Acad Sci (1992) 89:5152–6. doi:10.1073/pnas.89.11.5152
12. Shenker, A, Weinstein, LS, Moran, A, Pescovitz, OH, Charest, NJ, Boney, CM, et al. Severe endocrine and nonendocrine manifestations of the McCune-Albright syndrome associated with activating mutations of stimulatory G protein Gs. J Pediatr (1993) 123:509–18. doi:10.1016/S0022-3476(05)80943-6
13. Shenker, A, Weinstein, LS, Sweet, DE, and Spiegel, AM. An activating Gs alpha mutation is present in fibrous dysplasia of bone in the McCune-Albright syndrome. J Clin Endocrinol Metab (1994) 79:750–5. doi:10.1210/jcem.79.3.8077356
14. Weinstein, LS, Shenker, A, Gejman, PV, Merino, MJ, Friedman, E, and Spiegel, AM. Activating mutations of the stimulatory G protein in the McCune–albright syndrome. New Engl J Med (1991) 325:1688–95. doi:10.1056/NEJM199112123252403
15. More, A, Ito, I, Haridas, V, Chowdhury, S, Gu, Y, Dickson, P, et al. Oncogene addiction to GNAS in GNASR201 mutant tumors. Oncogene (2022) 41:4159–68. doi:10.1038/s41388-022-02388-6
16. O’Hayre, M, Vázquez-Prado, J, Kufareva, I, Stawiski, EW, Handel, TM, Seshagiri, S, et al. The emerging mutational landscape of G proteins and G-protein-coupled receptors in cancer. Nat Rev Cancer (2013) 13:412–24. doi:10.1038/nrc3521
17. Weinstein, LS, Chen, M, Xie, T, and Liu, J. Genetic diseases associated with heterotrimeric G proteins. Trends Pharmacol Sci (2006) 27:260–6. doi:10.1016/j.tips.2006.03.005
18. Delaney, D, Diss, TC, Presneau, N, Hing, S, Berisha, F, Idowu, BD, et al. GNAS1 mutations occur more commonly than previously thought in intramuscular myxoma. Mod Pathol (2009) 22:718–24. doi:10.1038/modpathol.2009.32
19. Walther, I, Walther, BM, Chen, Y, and Petersen, I. Analysis of GNAS1 mutations in myxoid soft tissue and bone tumors. Pathol Res Pract (2014) 210:1–4. doi:10.1016/j.prp.2013.09.003
20. Bekers, EM, Eijkelenboom, A, Rombout, P, van Zwam, P, Mol, S, Ruijter, E, et al. Identification of novel GNAS mutations in intramuscular myxoma using next-generation sequencing with single-molecule tagged molecular inversion probes. Diagn Pathol (2019) 14:15. doi:10.1186/s13000-019-0787-3
21. Sunitsch, S, Gilg, MM, Kashofer, K, Gollowitsch, F, Leithner, A, and Liegl-Atzwanger, B. Detection of GNAS mutations in intramuscular/cellular myxomas as diagnostic tool in the classification of myxoid soft tissue tumors. Diagn Pathol (2018) 13:52. doi:10.1186/s13000-018-0734-8
22. Libbrecht, L, Bempt, IV, Schubert, T, Sciot, R, and Galant, C. Next generation sequencing for GNAS uncovers CD34 as a sensitive marker for intramuscular myxoma. Ann Diagn Pathol (2019) 43:151409. doi:10.1016/j.anndiagpath.2019.151409
23. Bartuma, H, Möller, E, Collin, A, Domanski, HA, Von Steyern, FV, Mandahl, N, et al. Fusion of the FUS and CREB3L2 genes in a supernumerary ring chromosome in low-grade fibromyxoid sarcoma. Cancer Genet Cytogenet (2010) 199:143–6. doi:10.1016/j.cancergencyto.2010.02.011
24. Mertens, F, Fletcher, CDM, Antonescu, CR, Coindre, J-M, Colecchia, M, Domanski, HA, et al. Clinicopathologic and molecular genetic characterization of low-grade fibromyxoid sarcoma, and cloning of a novel FUS/CREB3L1 fusion gene. Lab Invest (2005) 85:408–15. doi:10.1038/labinvest.3700230
25. Mentzel, T, Calonje, E, Wadden, C, Camplejohn, RS, Beham, A, Smith, MA, et al. Myxofibrosarcoma. Clinicopathologic analysis of 75 cases with emphasis on the low-grade variant. Am J Surg Pathol (1996) 20:391–405. doi:10.1097/00000478-199604000-00001
26. Ogilvie, CM, Torbert, JT, Finstein, JL, Fox, EJ, and Lackman, RD. Clinical utility of percutaneous biopsies of musculoskeletal tumors. Clin Orthop Relat Res (2006) 450:95–100. doi:10.1097/01.blo.0000229302.52147.c7
27. Panagopoulos, I, Gorunova, L, Lobmaier, I, Bjerkehagen, B, and Heim, S. Karyotyping and analysis of GNAS locus in intramuscular myxomas. Oncotarget (2017) 8:22086–94. doi:10.18632/oncotarget.14986
28. Fecteau, RE, Lutterbaugh, J, Markowitz, SD, Willis, J, and Guda, K. GNAS mutations identify a set of right-sided, RAS mutant, villous colon cancers. PLoS One (2014) 9:e87966. doi:10.1371/journal.pone.0087966
29. Dal Molin, M, Matthaei, H, Wu, J, Blackford, A, Debeljak, M, Rezaee, N, et al. Clinicopathological correlates of activating GNAS mutations in intraductal papillary mucinous neoplasm (IPMN) of the pancreas. Ann Surg Oncol (2013) 20:3802–8. doi:10.1245/s10434-013-3096-1
30. Gal, S. Sequencing of double-stranded PCR products. In: CA Graham, and AJM Hill, editors. DNA sequencing protocols. New Jersey: Humana Press (1993). p. 183–90. doi:10.1385/0-89603-248-5:183
31. Dorit, RL, Ohara, O, Hwang, CB-C, Kim, JB, and Blackshaw, S. Direct DNA sequencing of PCR products. In: FM Ausubel, R Brent, RE Kingston, DD Moore, JG Seidman, and K Struhl, editors. Current protocols in molecular biology. Hoboken, NJ, USA: John Wiley and Sons, Inc. (2001). doi:10.1002/0471142727.mb1502s56
32. Sanger, F, Nicklen, S, and Coulson, AR. DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci U S A (1977) 74:5463–7. doi:10.1073/pnas.74.12.5463
33. Bevan, IS, Rapley, R, and Walker, MR. Sequencing of PCR-amplified DNA. PCR Methods Appl (1992) 1:222–8. doi:10.1101/gr.1.4.222
34. Kretz, K, Callen, W, and Hedden, V. Cycle sequencing. PCR Methods Appl (1994) 3:S107–12. doi:10.1101/gr.3.5.s107
35. Mardis, ER. The impact of next-generation sequencing on cancer genomics: from discovery to clinic. Cold Spring Harb Perspect Med (2019) 9:a036269. doi:10.1101/cshperspect.a036269
36. Baudhuin, LM, Lagerstedt, SA, Klee, EW, Fadra, N, Oglesbee, D, and Ferber, MJ. Confirming variants in next-generation sequencing panel testing by sanger sequencing. J Mol Diagn (2015) 17:456–61. doi:10.1016/j.jmoldx.2015.03.004
37. Mu, W, Lu, H-M, Chen, J, Li, S, and Elliott, AM. Sanger confirmation is required to achieve optimal sensitivity and specificity in next-generation sequencing panel testing. J Mol Diagn (2016) 18:923–32. doi:10.1016/j.jmoldx.2016.07.006
38. De Cario, R, Kura, A, Suraci, S, Magi, A, Volta, A, Marcucci, R, et al. Sanger validation of high-throughput sequencing in genetic diagnosis: still the best practice? Front Genet (2020) 11:592588. doi:10.3389/fgene.2020.592588
39. Nagahashi, M, Shimada, Y, Ichikawa, H, Kameyama, H, Takabe, K, Okuda, S, et al. Next generation sequencing-based gene panel tests for the management of solid tumors. Cancer Sci (2019) 110:6–15. doi:10.1111/cas.13837
40. Colomer, R, Miranda, J, Romero-Laorden, N, Hornedo, J, González-Cortijo, L, Mouron, S, et al. Usefulness and real-world outcomes of next generation sequencing testing in patients with cancer: an observational study on the impact of selection based on clinical judgement. EClinicalMedicine (2023) 60:102029. doi:10.1016/j.eclinm.2023.102029
41. Mauer, CB, Pirzadeh-Miller, SM, Robinson, LD, and Euhus, DM. The integration of next-generation sequencing panels in the clinical cancer genetics practice: an institutional experience. Genet Med (2014) 16:407–12. doi:10.1038/gim.2013.160
42. Tate, JG, Bamford, S, Jubb, HC, Sondka, Z, Beare, DM, Bindal, N, et al. COSMIC: the Catalogue of somatic mutations in cancer. Nucleic Acids Res (2019) 47:D941–D947. doi:10.1093/nar/gky1015
43. Panagopoulos, I, Gorunova, L, Lund-Iversen, M, Andersen, K, Andersen, HK, Lobmaier, I, et al. Cytogenetics of spindle cell/pleomorphic lipomas: karyotyping and FISH analysis of 31 tumors. Cancer Genomics Proteomics (2018) 15:193–200. doi:10.21873/cgp.20077
44. J McGowan-Jordan, R Hastings, and S Moore, editors. ISCN 2020: an international System for human cytogenomic nomenclature. Karger (2020). p. 164. doi:10.1159/isbn.978-3-318-06867-2
45. Altschul, SF, Gish, W, Miller, W, Myers, EW, and Lipman, DJ. Basic local alignment search tool. J Mol Biol (1990) 215:403–10. doi:10.1016/S0022-2836(05)80360-2
46. Smith, SC, Poznanski, AA, Fullen, DR, Ma, L, McHugh, JB, Lucas, DR, et al. CD34-positive superficial myxofibrosarcoma: a potential diagnostic pitfall. J Cutan Pathol (2013) 40:639–45. doi:10.1111/cup.12158
47. Sugiura, Y, Machinami, R, Matsumoto, S, Kanda, H, Ae, K, Takazawa, Y, et al. Prognostic value of CD34 expression status in patients with myxofibrosarcomas and undifferentiated pleomorphic sarcomas. Sci Rep (2021) 11:15494. doi:10.1038/s41598-021-94834-w
48. Smits, AJJ, Kummer, JA, de Bruin, PC, Bol, M, van den Tweel, JG, Seldenrijk, KA, et al. The estimation of tumor cell percentage for molecular testing by pathologists is not accurate. Mod Pathol (2014) 27:168–74. doi:10.1038/modpathol.2013.134
49. Lhermitte, B, Egele, C, Weingertner, N, Ambrosetti, D, Dadone, B, Kubiniek, V, et al. Adequately defining tumor cell proportion in tissue samples for molecular testing improves interobserver reproducibility of its assessment. Virchows Archiv (2017) 470:21–7. doi:10.1007/s00428-016-2042-6
50. Sreekantaiah, C, Karakousis, CP, Leong, SP, and Sandberg, AA. Trisomy 8 as a nonrandom secondary change in myxoid liposarcoma. Cancer Genet Cytogenet (1991) 51:195–205. doi:10.1016/0165-4608(91)90132-e
51. Bridge, JA, Swarts, SJ, Buresh, C, Nelson, M, Degenhardt, JM, Spanier, S, et al. Trisomies 8 and 20 characterize a subgroup of benign fibrous lesions arising in both soft tissue and bone. Am J Pathol (1999) 154:729–33. doi:10.1016/S0002-9440(10)65319-9
52. Paulsson, K, and Johansson, B. Trisomy 8 as the sole chromosomal aberration in acute myeloid leukemia and myelodysplastic syndromes. Pathol Biol (Paris) (2007) 55:37–48. doi:10.1016/j.patbio.2006.04.007
53. Hemsing, AL, Hovland, R, Tsykunova, G, and Reikvam, H. Trisomy 8 in acute myeloid leukemia. Expert Rev Hematol (2019) 12:947–58. doi:10.1080/17474086.2019.1657400
54. Broberg, K, Toksvig-Larsen, S, Lindstrand, A, and Mertens, F. Trisomy 7 accumulates with age in solid tumors and non-neoplastic synovia. Genes Chromosomes Cancer (2001) 30:310–5. doi:10.1002/1098-2264(2000)9999:9999<::aid-gcc1096>3.0.co;2-c
55. Swansbury, J. Cytogenetic and genetic studies in solid tumors: background. Methods Mol Biol (2003) 220:125–33. doi:10.1385/1-59259-363-1:125
56. Johansson, B, Heim, S, Mandahl, N, Mertens, F, and Mitelman, F. Trisomy 7 in nonneoplastic cells. Genes Chromosomes Cancer (1993) 6:199–205. doi:10.1002/gcc.2870060402
57. Hayashi, K. PCR-SSCP: a simple and sensitive method for detection of mutations in the genomic DNA. PCR Methods Appl (1991) 1:34–8. doi:10.1101/gr.1.1.34
58. Riminucci, M, Fisher, LW, Majolagbe, A, Corsi, A, Lala, R, De Sanctis, C, et al. A novel GNAS1 mutation, R201G, in McCune-albright syndrome. J Bone Mineral Res (1999) 14:1987–9. doi:10.1359/jbmr.1999.14.11.1987
59. Liang, Q, Wei, M, Hodge, L, Fanburg-Smith, JC, Nelson, A, Miettinen, M, et al. Quantitative analysis of activating alpha subunit of the G protein (Gsα) mutation by pyrosequencing in fibrous dysplasia and other bone lesions. J Mol Diagn (2011) 13:137–42. doi:10.1016/j.jmoldx.2010.10.003
60. Romanet, P, Philibert, P, Fina, F, Cuny, T, Roche, C, Ouafik, L, et al. Using digital droplet polymerase chain reaction to detect the mosaic GNAS mutations in whole blood DNA or circulating cell-free DNA in fibrous dysplasia and McCune-albright syndrome. J Pediatr (2019) 205:281–5. doi:10.1016/j.jpeds.2018.09.070
61. Lin, Y-L, Ma, R, and Li, Y. The biological basis and function of GNAS mutation in pseudomyxoma peritonei: a review. J Cancer Res Clin Oncol (2020) 146:2179–88. doi:10.1007/s00432-020-03321-8
62. Weinstein, LS, Liu, J, Sakamoto, A, Xie, T, and Chen, M. Minireview: GNAS: normal and abnormal functions. Endocrinology (2004) 145:5459–64. doi:10.1210/en.2004-0865
63. Weinstein, LS, Yu, S, Warner, DR, and Liu, J. Endocrine manifestations of stimulatory G protein alpha-subunit mutations and the role of genomic imprinting. Endocr Rev (2001) 22:675–705. doi:10.1210/edrv.22.5.0439
64. Khan, SK, Yadav, PS, Elliott, G, Hu, DZ, Xu, R, and Yang, Y. Induced GnasR201H expression from the endogenous Gnas locus causes fibrous dysplasia by up-regulating Wnt/β-catenin signaling. Proc Natl Acad Sci U S A (2018) 115:E418–E427. doi:10.1073/pnas.1714313114
65. Nomura, R, Saito, T, Mitomi, H, Hidaka, Y, Lee, S, Watanabe, S, et al. GNAS mutation as an alternative mechanism of activation of the Wnt/β-catenin signaling pathway in gastric adenocarcinoma of the fundic gland type. Hum Pathol (2014) 45:2488–96. doi:10.1016/j.humpath.2014.08.016
66. Wilson, CH, McIntyre, RE, Arends, MJ, and Adams, DJ. The activating mutation R201C in GNAS promotes intestinal tumourigenesis in Apc(Min/+) mice through activation of Wnt and ERK1/2 MAPK pathways. Oncogene (2010) 29:4567–75. doi:10.1038/onc.2010.202
67. Tirosh, A, Jin, DX, De Marco, L, Laitman, Y, and Friedman, E. Activating genomic alterations in the Gs alpha gene (GNAS) in 274 694 tumors. Genes Chromosomes Cancer (2020) 59:503–16. doi:10.1002/gcc.22854
68. Zhao, X, Deng, P, Iglesias-Bartolome, R, Amornphimoltham, P, Steffen, DJ, Jin, Y, et al. Expression of an active Gαs mutant in skeletal stem cells is sufficient and necessary for fibrous dysplasia initiation and maintenance. Proc Natl Acad Sci U S A (2018) 115:E428–E437. doi:10.1073/pnas.1713710115
69. Zhadina, M, Roszko, KL, Geels, RES, de Castro, LF, Collins, MT, and Boyce, AM. Genotype-phenotype correlation in fibrous dysplasia/McCune-albright syndrome. J Clin Endocrinol Metab (2021) 106:1482–90. doi:10.1210/clinem/dgab053
70. Lumbroso, S, Paris, F, and Sultan, CEuropean Collaborative Study. Activating Gsalpha mutations: analysis of 113 patients with signs of McCune-Albright syndrome--a European Collaborative Study. J Clin Endocrinol Metab (2004) 89:2107–13. doi:10.1210/jc.2003-031225
71. Jour, G, Oultache, A, Sadowska, J, Mitchell, T, Healey, J, Nafa, K, et al. GNAS mutations in fibrous dysplasia: a comparative study of standard sequencing and locked nucleic acid PCR sequencing on decalcified and nondecalcified formalin-fixed paraffin-embedded tissues. Appl Immunohistochem Mol Morphol (2016) 24:660–7. doi:10.1097/PAI.0000000000000242
72. Wheeler, EC, Vora, S, Mayer, D, Kotini, AG, Olszewska, M, Park, SS, et al. Integrative RNA-omics discovers GNAS alternative splicing as a phenotypic driver of splicing factor–mutant neoplasms. Cancer Discov (2022) 12:836–55. doi:10.1158/2159-8290.CD-21-0508
73. Meyer, AP, Forrest, ME, Nicolau, S, Wiszniewski, W, Bland, MP, Tsao, C, et al. Pathogenic missense variants altering codon 336 of GARS1 lead to divergent dominant phenotypes. Hum Mutat (2022) 43:869–76. doi:10.1002/humu.24372
74. Roeber, S, Grasbon-Frodl, E-M, Windl, O, Krebs, B, Xiang, W, Vollmert, C, et al. Evidence for a pathogenic role of different mutations at codon 188 of PRNP. PLoS One (2008) 3:e2147. doi:10.1371/journal.pone.0002147
Keywords: intramuscular myxoma, chromosomal aberrations, GNAS, codon 201, codon 227, Sanger sequencing, Ion torrent
Citation: Hatchett WJ, Brunetti M, Andersen K, Tandsæther MR, Lobmaier I, Lund-Iversen M, Lien-Dahl T, Micci F and Panagopoulos I (2024) Genetic characterization of intramuscular myxomas. Pathol. Oncol. Res. 30:1611553. doi: 10.3389/pore.2024.1611553
Received: 16 October 2023; Accepted: 10 January 2024;
Published: 22 January 2024.
Edited by:
Andrea Ladányi, National Institute of Oncology (NIO), HungaryCopyright © 2024 Hatchett, Brunetti, Andersen, Tandsæther, Lobmaier, Lund-Iversen, Lien-Dahl, Micci and Panagopoulos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ioannis Panagopoulos, aW9hbm5pcy5wYW5hZ29wb3Vsb3NAcnItcmVzZWFyY2gubm8=