Skip to main content
  • Letter to the Editor
  • Open access
  • Published:

Targeted next-generation sequencing of circulating free DNA enables non-invasive tumor detection in myxoid liposarcomas

Main text

Myxoid liposarcoma (MLS), a malignant soft-tissue tumor derived from lipocytes, is characterized by specific genetic translocations t(12;16) and t(12;22) on a background of few additional chromosomal changes [1]. About 30% of patients with localized high-grade MLS will eventually develop distant metastases [2]. Unlike other soft tissue sarcomas (STS), MLS exhibit a distinct pulmonary and extra-pulmonary metastatic pattern. Imaging for follow-up is thus extensive and requires whole-body Magnetic Resonance Imaging (MRI) or a combination of various imaging modalities [3].

Previously, we investigated the potential of circulating tumor DNA (ctDNA) to detect tumor recurrence and monitor treatment response [4]. However, due to the unique nature of each translocation and few hotspot mutations, quantification of ctDNA was technically demanding and assays were not suitable for routine diagnostics.

To overcome these limitations, we developed targeted next generation sequencing (NGS) - based approaches, which allow ultrasensitive detection of MLS DNA in a routine diagnostic setting and prospective clinical trials. As standard NGS panels don’t cover common genetic alterations in MLS, we designed a lockdown panel which encompasses genes with a reported mutation frequency of at least 5%. The 36,541 base pair (bp) standard panel covers the introns of DDIT3, FUS and EWS where the t (12;16) and t (12;22) translocations occur, the TERT promoter region and mutation hotspots within exons from seven genes (Supplementary Fig. 1) [5,6,7,8]. Applying molecular barcodes for digital error correction [9] allowed the detection of mutations with a variant allele frequency (VAF) of 0.05% (Supplementary Fig. 2). 51 MLS tumors and two MLS cell lines (402-91 and 1765-92) were sequenced. Matched normal DNA was available for 23 tumors. Breakpoints could be identified in 49 tumors, and both cell lines. Translocations occurred in 87.7% between DDIT3 and FUS and in 8.3% between DDIT3 and EWSR1 (Fig. 1A). No translocations were observed in matched leukocyte DNA (sensitivity of 96% and specificity of 100%).

Fig. 1
figure 1

Mutational profiling of myxoid liposarcomas. A 51 MLS tumors and two MLS cell lines (402-91 and 1765-92) were sequenced with an MLS specific lockdown panel. Chromosomal translocations could be detected in 49/51 tumors (1 - 49) and both cell lines (C1, C2) and occurred between DDIT3 and FUS in 87.7% and between DDIT3 and EWSR1 in 8.3%. Breakpoints clustered to several distinct regions within DDIT3 and FUS but were not restricted to a single site. For 36 tumors, both breakpoints of the reciprocal translocation could be determined. Arrows depict the sites where the chromosomal breaks occurred. Areas with an increased likelihood of chromosomal breaks are colored in blue. Breakpoints were annotated to homo sapiens (human) genome assembly GRCh38 (hg38). B For 36 tumors, where both breakpoints could be sequenced, loss or gain of DNA during the translocation event could be determined and specifically mapped to either DDIT3 or FUS. A mean loss of 7 bp (SD 64 bp) occurred on chromosome 12 (DDIT3) and of 11 bp (SD 73 bp) on chromosome 16 (FUS). There was considerable intertumor variability. Each symbol of the graph represents one tumor, with the largest gain of 282 bp on chromosome 16 and the biggest deletion of 138 bp on chromosome 12. C Point mutations occurred most commonly in the TERT promoter region (73%) and PIK3CA (33%). Only one tumor showed an additional point mutation in TET2 and one tumor analyzed for intratumor heterogeneity (tumor 2 in Fig. 1 D) displayed two mutations in PTEN. TERT promoter mutations were all detected at the well-known hotspot locations C228T and C250T, with a prevalence of 61% and 12% of all analyzed tumors respectively. PIK3CA mutations occurred at well-known hotspot mutations in exon 9 (c.1624G > A, c.1633G > A, c.1633G > C, c1634A > G) and exon 20 (c.3140A > G) but also at less commonly annotated positions in exon 5 (c.1035 T > A) and exon 8 (c. 1345C > A). D To determine intratumor heterogeneity of MLS, 20 individual samples of tumor 1 were taken at uniform distances and 10 samples were taken from tumor 2. Each sample was analyzed separately with the standard panel. The individual tumor-specific breakpoints identified before (Fig. 1 A) were detectable in all samples. There were no TERT promoter mutations in tumor 1, however PIK3CA mutations were present in 3/20 samples. Two showed a hotspot mutation in exon 9 (chr3:179,218,294; c.1624G > A) and one samples had a hotspot mutation in exon 20 (chr3:179,234,297; c.3140A > G). In tumor 2, TERT promoter mutations were present in all samples, however nine samples contained the C250T and one sample the C228T mutations. PIK3CA mutations were even more diverse. 6/10 samples contained a PIK3CA mutation. Of these samples three showed the well-known hotspot mutation in exon 9 (chr3:179,218,304; c.1634A > G) and the remaining three the hotspot mutation in exon 20 (chr3:179,234,297; c.3140A > G). A deletion (c.388del) and a point mutation in (c.377C > A) were identified in PTEN. The intratumor heterogeneity of PIK3CA, PTEN and the TERT promoter point towards multiple subclones which emerge from a tumor ancestor that initially acquired the characteristic t (12;16) driver translocation. Depicted is the tumor (red), which was reconstructed from the patients’ MRI scans. The black lines depict how the tumor was sectioned for histopathologic assessment. Although the samples were taken at uniform distances throughout the tumor, the exact location of each sample within the tumor cannot be determined due to the retrospective nature of the study

We could determine both breakpoints of the balanced translocations in 36 of the 51 tumors and both cell lines. This allowed us to determine if deletions or insertions occurred during the translocation event. We observed a mean loss of 7 bp (SD 64 bp) on chromosome 12 (DDIT3) and of 11 bp (SD 73 bp) on chromosome 16 (FUS) with a high intertumor variability (Fig. 1B).

Point mutations were detectable in 74.5% of tumors. Mutations in the TERT promoter region were prevalent in 73% of tumors. Thereof the C228T mutation occurred in 61% and the C250T mutation in 12% of analyzed tumors. PIK3CA mutations were found in 33% of MLS samples. Besides the well-known hotspot mutations in exon 9 (c.1624G > A, c.1633G > A, c.1633G > C, c1634A > G) and exon 20 (c.3140A > G) we identified less commonly annotated mutations in exon 5 (c.1035 T > A) and exon 8 (c. 1345C > A). Only two additional genes were mutated at low frequency in our cohort, TET2 (2%) and PTEN (4%) (Fig. 1C). Taken together, the lockdown panel could detect at least one mutation in 96% of all the tumor samples. Combining the breakpoints and point mutations, the panel detected an average 2.8 somatic mutations (1.7 breakpoints and 1.1 point mutations) per tumor, which can be targeted in circulating free DNA (cfDNA).

To determine the impact of tumor heterogeneity on ctDNA detection, spatially separate samples from two tumors (tumor 1 and 2 from Fig. 1A) were analyzed. Results were compared to matched white blood cells. Mutations with a VAF of at least 5 % were recorded. The two t (12;16) breakpoints were detected in all specimens. The presence of the major MLS driver translocation confirms its importance from tumor initiation to promotion. The consistency of patient individual breakpoints was previously reported for multifocal MLS [10]. In contrast, there was marked intratumor heterogeneity for TERT promoter, PIK3CA and PTEN mutations (Fig. 1D). This supports the hypothesis that additional mutations seem to occur secondarily and evolve within different tumor subclones [11]. Thus, tracking of breakpoint fragments in cfDNA promises detection of the primary tumor and all its potential metastases. In contrast, point mutations identified in the primary tumor might or might not be present in its metastases, depending of their clones of origin.

The assay was subsequently employed to quantify ctDNA in plasma samples of MLS patients. Nine plasma samples, collected during the 2-year treatment of patient 1 (see Fig. 1D) with a localized MLS who later developed metastatic disease were analyzed with the standard panel. The two PIK3CA mutations (c.1624G > A and c.3140A > G) were additionally quantified by droplet digital PCR. ctDNA levels decreased after tumor resection and increased when metastatic disease was detected. We observed a decline in ctDNA when radio/chemotherapy was initiated. However, with increasing tumor burden, concentrations rose again after several months (Fig. 2A and Supplementary Fig. 3 A). Serial ctDNA testing promises monitoring of treatment response in metastatic MLS. It might be especially beneficial for patients treated with immunotherapeutic agents, that challenge established imaging-based response assessment criteria [12, 13].

Fig. 2
figure 2

Quantification of ctDNA in patients’ plasma samples. A Quantification of ctDNA in 9 plasma samples of patient 1 collected during 2 years of treatment. ctDNA was determined by NGS (standard panel) and PIK3CA mutations (c.1624G > A and c.3140A > G) were additionally quantified by ddPCR. He initially presented with a localized MLS of the thigh which was completely resected. Soon after, he developed metastatic disease with predominantly osseous lesions. He then received radiotherapy of bone metastasis and several courses of chemotherapy. Repeated imaging during follow-up showed numerous new skeletal lesions and the patient again received radiotherapy to selected metastasis. He succumbed to his disease 1.5 years after removal of the primary tumor. ctDNA increased to 172 copies/ml (sample 6) when metastatic disease was detected and decreased during radio/chemotherapy to 14 copies/ml (sample 7) and 7 copies/ml (sample 8). There was a rapid incline in ctDNA when multiple new metastases were detected (sample 9: 552 copies/ml). Standard imaging, which reflects the total mass of viable and necrotic tumor cells, showed a steady increase (blue area). The irradiated tumor volume is depicted as a surrogate marker for the necrotic tumor mass (green area). t (12;16) ctDNA levels were higher than PIK3CA ctDNA concentrations. This most likely reflect intertumor heterogeneity with only a fraction of metastases carrying PIK3CA mutations (Fig. 1 D and Supplementary Fig. 4). B Additional target mutations from exome sequencing increase sensitivity of ctDNA detection. Tumor 2 was subjected to exome sequencing to identify additional target mutations. Together with breakpoints and mutations from the standard panel, a 7320 bp hybrid exome panel targeting 15 genomic regions was designed. ctDNA in plasma obtained during treatment was determined by the standard and exome panel. He initially received neo-adjuvant radiotherapy to an MLS of his right thigh and subsequently the tumor was completely resected. Two plasma samples were collected prior to commencement of radiotherapy, a third sample before surgery and a fourth sample after tumor resection. ctDNA quantified by the exome panel (red line) was present in similar amounts at the two time points before treatment, declined after radiotherapy and was not detectable after tumor resection. The standard panel (dashed black line) could detect ctDNA only in the first sample, showing reduced sensitivity compared to the exome panel. The blue area represents the tumor volume as calculated from the MRI scans. C Comparison of different assays in detecting MLS tumor-DNA. Dilution series of MLS tumor-DNA from two tumors (patient 2 and 3) in matched normal DNA were analyzed by ddPCR (PIK3CA mutations p.N345K, c.1035 T > A and p.E545G, c.1634A > G), the NGS standard panel and respective exome panels. Depicted are mean values and linear regression of n = 2 tumors for ddPCR, n = 2 for the standard panel and n = 2 for the exome panels. We observed a similar performance for ddPCR and the standard panel, whereas detection of tumor-DNA with exome panels was clearly superior. D Patient 3 presented with two small localized tumors of his legs (red) after numerous prior resections at another hospital. The tumors were completely resected, but he repeatedly developed local recurrences at both locations in the following years. E These recurrences were subsequently resected at four consecutive operations before a small lung metastasis (0.3 cm3) was detected and subsequently removed. Exome panels were obtained from sequencing one of the primary lesions and the lung metastasis. ctDNA was subsequently quantified with both panels in 15 plasma samples obtained during the course of his treatment. During multifocal localized disease, ctDNA values undulated at low concentrations depending on the presence of viable tumor tissue. The exome panel from the primary tumor best reflected the clinical course (enlarged image section). There was one outlier (circle). Despite complete tumor resections, ctDNA values never reached the baseline indicating MRD. In contrast the plasma sample obtained shortly before resection of the lung metastasis showed markedly increased ctDNA with a decline after its resection. The red line represents ctDNA values measured by the exome panel from the primary lesion and the dashed line ctDNA measured by the exome panel obtained from the lung metastasis. The blue area depicts the tumor volume

From the two PIK3CA mutations identified in the primary tumor, only c.3140A > G ctDNA correlated with t (12;16) ctDNA. This indicates that most metastases originated from a clone in the primary tumor, which harbored this mutation (Supplementary Fig. 4). As new treatment opportunities which target PIK3CA are only effective in PIK3CA mutated cells [14], liquid biopsy may enable us to identify patients with druggable metastases without the need of repeated biopsies.

Further mutations obtained by exome sequencing of individual tumors were added to the mutations already identified by the standard panel to lower the limit of detection (LoD) (Supplementary Fig. 5). These hybrid panels (exome panels) allowed us to monitor ctDNA of localized tumors as exemplified in the following scenario. Patient 2 received neoadjuvant radiotherapy and subsequently complete resection (tumor necrosis rate > 90%) of a localized MLS. Two plasma samples obtained at an interval of 9 days before initiation of radiotherapy showed similar concentrations of ctDNA (Fig. 2B). The third specimen obtained after radiotherapy demonstrated markedly decreased ctDNA. No ctDNA was detectable in the fourth sample collected after surgery. Analysis with the standard panel detected ctDNA in the initial plasma sample only. Calculating relative amounts, the fraction of ctDNA was between 0% and 0.05% (Supplementary Fig. 3 B). Detection of these minute amounts in limited plasma samples requires extremely sensitive assays. This was accomplished by exome panels, which target multiple mutations simultaneously, thus detecting more mutant copies in the same amounts of tumor DNA than ddPCR and the standard panel (Fig. 2C and Supplementary Fig. 5).

The impact of intertumor heterogeneity on ctDNA quantification was assessed in a patient who was initially presented with multifocal disease of his legs after previous resections at another hospital (Fig. 2D). In the course of his treatment, he suffered from four local recurrences before a small lung metastasis (0.3 cm3) was identified and resected. Exome sequencing was conducted from the initial leg tumor and the lung metastasis. The exome panels targeted 22 and 17 genomic regions respectively, and 6 mutations were identical in both panels. ctDNA in 15 plasma samples was quantified with both panels. During localized disease ctDNA fluctuated depending on the amount of viable tumor mass. The panel from the primary lesion performed superior during localized disease. The small lung metastasis however, led to markedly increased ctDNA concentrations which again decreased after resection. Both panels performed similarly during metastatic disease (Fig. 2E). Comparatively high ctDNA concentrations of lung metastases in MLS might enable detection of recurrence earlier by liquid biopsy than with imaging-based approaches [4, 15].


In this study, we present an approach for ctDNA monitoring in MLS patients in a routine diagnostic setting using a disease and patient-specific hybrid capture NGS technique. Quantification of ctDNA on the basis of cancer genomic profiling could help to predict tumor recurrence, and monitor tumor heterogeneity and treatment response in metastatic disease with minimal invasiveness and at affordable cost. The assay can easily be adapted to other translocation driven tumors, e.g. synovial sarcomas. Given our promising results, the methods we have described warrant investigations in prospective trials with larger cohorts, so they can timely be translated into clinical practice.

Availability of data and materials

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



Myxoid liposarcoma


Soft tissue sarcoma


Circulating tumor DNA


Circulating free DNA


Next generation sequencing


Base pair


Formalin-Fixed Paraffin-Embedded tissue


Homo sapiens (human) genome assembly GRCh38


Droplet digital PCR


Variant allele frequency


Magnetic Resonance Imaging


Standard Deviation


Limit of Detection


  1. Hofvander J, Viklund B, Isaksson A, Brosjö O, Vult von Steyern F, Rissler P, et al. Different patterns of clonal evolution among different sarcoma subtypes followed for up to 25 years. Nat Commun. 2018;9:3662.

    Article  Google Scholar 

  2. Muratori F, Bettini L, Frenos F, Mondanelli N, Greto D, Livi L, et al. Myxoid Liposarcoma: prognostic factors and metastatic pattern in a series of 148 patients treated at a single institution. Int J Surg Oncol. 2018;2018(8928706):9.

  3. Dürr HR, Rauh J, Baur-Melnyk A, Knösel T, Lindner L, Roeder F, et al. Myxoid liposarcoma: local relapse and metastatic pattern in 43 patients. BMC Cancer. 2018;18:304.

    Article  Google Scholar 

  4. Braig D, Becherer C, Bickert C, Braig M, Claus R, Eisenhardt AE, et al. Genotyping of circulating cell-free DNA enables noninvasive tumor detection in myxoid liposarcomas. Int J Cancer. 2019;145:1148–61.

    Article  CAS  Google Scholar 

  5. Xiang H, Wang J, Hisaoka M, Zhu X. Characteristic sequence motifs located at the genomic breakpoints of the translocation t (12;16) and t (12;22) in myxoid liposarcoma. Pathology. 2008;40:547–52.

    Article  CAS  Google Scholar 

  6. Koelsche C, Renner M, Hartmann W, Brandt R, Lehner B, Waldburger N, et al. TERT promoter hotspot mutations are recurrent in myxoid liposarcomas but rare in other soft tissue sarcoma entities. J Exp Clin Cancer Res. 2014;33:33.

    Article  Google Scholar 

  7. Joseph CG, Hwang H, Jiao Y, Wood LD, Kinde I, Wu J, et al. Exomic analysis of myxoid liposarcomas, synovial sarcomas, and osteosarcomas. Genes Chromosomes Cancer. 2014;53:15–24.

    Article  CAS  Google Scholar 

  8. Movva S, Wen W, Chen W, Millis SZ, Gatalica Z, Reddy S, et al. Multi-platform profiling of over 2000 sarcomas: identification of biomarkers and novel therapeutic targets. Oncotarget. 2015;6:12234–47.

    Article  Google Scholar 

  9. MacConaill LE, Burns RT, Nag A, Coleman HA, Slevin MK, Giorda K, et al. Unique, dual-indexed sequencing adapters with UMIs effectively eliminate index cross-talk and significantly improve sensitivity of massively parallel sequencing. BMC Genomics. 2018;19:30.

    Article  Google Scholar 

  10. de Vreeze R, de Jong D, Nederlof P, Ruijter HJ, Boerrigter L, Haas R, et al. Multifocal myxoid liposarcoma--metastasis or second primary tumor?: a molecular biological analysis. J Mol Diagn. 2010;12:238–43.

    Article  Google Scholar 

  11. Ferreira MSV, Crysandt M, Braunschweig T, Jost E, Voss B, Bouillon A-S, et al. Presence of TERT promoter mutations is a secondary event and associates with elongated telomere length in Myxoid Liposarcomas. Int J Mol Sci. 2018;19(2):608.

    Article  Google Scholar 

  12. Unterrainer M, Ruzicka M, Fabritius MP, Mittlmeier LM, Winkelmann M, Rübenthaler J, et al. PET/CT imaging for tumour response assessment to immunotherapy: current status and future directions. Eur Radiol Exp. 2020;4:63.

    Article  Google Scholar 

  13. Bratman SV, Yang SYC, Iafolla MAJ, Liu Z, Hansen AR, Bedard PL, et al. Personalized circulating tumor DNA analysis as a predictive biomarker in solid tumor patients treated with pembrolizumab. Nat Cancer. 2020;1:873–81.

    Article  Google Scholar 

  14. André F, Ciruelos E, Rubovszky G, Campone M, Loibl S, Rugo HS, et al. Alpelisib for PIK3CA-mutated, hormone receptor-positive advanced breast Cancer. N Engl J Med. 2019;380:1929–40.

    Article  Google Scholar 

  15. Avanzini S, Kurtz DM, Chabon JJ, Moding EJ, Hori SS, Gambhir SS, et al. A mathematical model of ctDNA shedding predicts tumor detection size. Sci Adv. 2020;6:eabc4308.

    Article  Google Scholar 

Download references


The authors thank Dr. Pierre Åman, University of Gothenburg and Dr. Marcus Renner, University of Heidelberg for the kind donation of MLS cell lines (402-91 and 1765-92) and FFPE tissue samples. Dr. Dietmar Pfeifer and his team for the use of their NGS facility, technical assistance and advice. Marie Follo and the team of the Lighthouse Core Facility for their assistance with ddPCR. Prof. S. Laßmann and her team for use of the Fragment Analyzer and NGS facility. Prof. Börries and her team for bioinformatic support. Dr. A. Flörcken, principal investigator of the soft tissue sarcoma biobank (ZeBanC) and head of the interdisciplinary sarcoma board of the Comprehensive Cancer Center Charité (CCCC) for her support and contribution. This research has been conducted using biological samples and data obtained from the Central Biobank Charité (ZeBanC).


This study was supported by the German Research Foundation (DFG) (Grant number 396168587; BR 5712/1-1, CL 427/4-1, EI 866/7-1) and the Fördergesellschaft Forschung Tumorbiologie (Grant “Liquid-Biopsy-Initiative”). S.U.E. is a Heisenberg Professor of the DFG (EI 866/9-1). SUE is also supported by a personal DFG project grant (EI 866/10-1) that is not related to the study.

Author information

Authors and Affiliations



Anja E. Eisenhardt: study design; sample acquisition; data acquisition; data analysis and interpretation; drafting of manuscript. Adrian Schmid: data acquisition, analysis and interpretation of data; revision of manuscript. Julia Esser: data acquisition and analysis. Zacharias Brugger: data acquisition and analysis. Ute Lausch: data acquisition and analysis. Jurij Kiefer: sample acquisition; data interpretation; revision of manuscript. Moritz Braig: analysis of MRI data; revision of manuscript. Alexander Runkel: data acquisition and analysis. Julius Wehrle: data acquisition and analysis. Rainer Claus: data interpretation; revision of manuscript. Peter Bronsert: sample acquisition; data interpretation; revision of manuscript. Andreas Leithner: sample acquisition; revision of manuscript. Bernadette Liegl-Atzwanger: sample acquisition; revision of manuscript. Johannes Zeller: sample acquisition; revision of manuscript. Remo Papini: sample acquisition, revision of manuscript. Maximilian von Laffert: sample acquisition, revision of manuscript. Berit Pfitzner: sample acquisition, revision of manuscript. Georgios Koulaxouzidis: sample acquisition, revision of manuscript. Riccardo E. Giunta: data interpretation; revision of manuscript. Steffen U. Eisenhardt: study design; sample acquisition; data analysis and interpretation; drafting of manuscript. David Braig: study design; sample acquisition; data acquisition; data analysis and interpretation; drafting of manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to David Braig.

Ethics declarations

Ethics approval and consent to participate

The Ethics Committee of the Albert-Ludwigs-University of Freiburg, Germany, approved the study (study number 243/13 and 236/16). The design and performance of the study are in accordance with the Declaration of Helsinki. Signed informed consent was obtained from all participants before inclusion, allowing analysis of tumor tissue, blood samples and clinical data.

Consent for publication

Not applicable.

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: Supplementary Figure 1.

Specifications of standard lockdown panel for MLS. Supplementary Figure 2. Evaluation of limit of detection (LoD). Supplementary Figure 3. Comparison of absolute and relative ctDNA quantification. Supplementary Figure 4. Impact of tumor heterogeneity on ctDNA detection. Supplementary Figure 5. Additional target mutations from exome sequencing increase sensitivity of tumor DNA detection. Supplementary Figure 6. Determination of sensitivity and specificity of standard panel.

Additional file 2:

Materials and Methods.

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

Eisenhardt, A.E., Schmid, A., Esser, J. et al. Targeted next-generation sequencing of circulating free DNA enables non-invasive tumor detection in myxoid liposarcomas. Mol Cancer 21, 50 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: