A genomic case study of desmoplastic small round cell tumor: comprehensive analysis reveals insights into potential therapeutic targets and development of a monitoring tool for a rare and aggressive disease

Background Genome-wide profiling of rare tumors is crucial for improvement of diagnosis, treatment, and, consequently, achieving better outcomes. Desmoplastic small round cell tumor (DSRCT) is a rare type of sarcoma arising from mesenchymal cells of abdominal peritoneum that usually develops in male adolescents and young adults. A specific translocation, t(11;22)(p13;q12), resulting in EWS and WT1 gene fusion is the only recurrent molecular hallmark and no other genetic factor has been associated to this aggressive tumor. Here, we present a comprehensive genomic profiling of one DSRCT affecting a 26-year-old male, who achieved an excellent outcome. Methods We investigated somatic and germline variants through whole-exome sequencing using a family based approach and, by array CGH, we explored the occurrence of genomic imbalances. Additionally, we performed mate-paired whole-genome sequencing for defining the specific breakpoint of the EWS-WT1 translocation, allowing us to develop a personalized tumor marker for monitoring the patient by liquid biopsy. Results We identified genetic variants leading to protein alterations including 12 somatic and 14 germline events (11 germline compound heterozygous mutations and 3 rare homozygous polymorphisms) affecting genes predominantly involved in mesenchymal cell differentiation pathways. Regarding copy number alterations (CNA) few events were detected, mainly restricted to gains in chromosomes 5 and 18 and losses at 11p, 13q, and 22q. The deletions at 11p and 22q indicated the presence of the classic translocation, t(11;22)(p13;q12). In addition, the mapping of the specific genomic breakpoint of the EWS-WT1 gene fusion allowed the design of a personalized biomarker for assessing circulating tumor DNA (ctDNA) in plasma during patient follow-up. This biomarker has been used in four post-treatment blood samples, 3 years after surgery, and no trace of EWS-WT1 gene fusion was detected, in accordance with imaging tests showing no evidence of disease and with the good general health status of the patient. Conclusions Overall, our findings revealed genes with potential to be associated with risk assessment and tumorigenesis of this rare type of sarcoma. Additionally, we established a liquid biopsy approach for monitoring patient follow-up based on genomic information that can be similarly adopted for patients diagnosed with a rare tumor. Electronic supplementary material The online version of this article (doi:10.1186/s40246-016-0092-0) contains supplementary material, which is available to authorized users.


Methods:
We investigated somatic and germline variants through whole-exome sequencing using a family based approach and, by array CGH, we explored the occurrence of genomic imbalances. Additionally, we performed mate-paired whole-genome sequencing for defining the specific breakpoint of the EWS-WT1 translocation, allowing us to develop a personalized tumor marker for monitoring the patient by liquid biopsy. Results: We identified genetic variants leading to protein alterations including 12 somatic and 14 germline events (11 germline compound heterozygous mutations and 3 rare homozygous polymorphisms) affecting genes predominantly involved in mesenchymal cell differentiation pathways. Regarding copy number alterations (CNA) few events were detected, mainly restricted to gains in chromosomes 5 and 18 and losses at 11p, 13q, and 22q. The deletions at 11p and 22q indicated the presence of the classic translocation, t(11;22)(p13;q12). In addition, the mapping of the specific genomic breakpoint of the EWS-WT1 gene fusion allowed the design of a personalized biomarker for assessing circulating tumor DNA (ctDNA) in plasma during patient follow-up. This biomarker has been used in four post-treatment blood samples, 3 years after surgery, and no trace of EWS-WT1 gene fusion was detected, in accordance with imaging tests showing no evidence of disease and with the good general health status of the patient.
(Continued on next page) (Continued from previous page) Conclusions: Overall, our findings revealed genes with potential to be associated with risk assessment and tumorigenesis of this rare type of sarcoma. Additionally, we established a liquid biopsy approach for monitoring patient follow-up based on genomic information that can be similarly adopted for patients diagnosed with a rare tumor.
Keywords: Desmoplastic small round cell tumor, Genomic profiling, Whole-exome sequencing, EWS-WT1 gene fusion, Personalized biomarker, Liquid biopsy Background Comprehensive molecular profiling is an especially important tool to gain insights on the biological pathways involved in tumor onset and to improve the management and treatment of rare tumors. Desmoplastic small round cell tumor (DSRCT) is a very rare type of sarcoma, with an age-adjusted incidence rate of 0.3 cases/million [1], which typically arises from the abdominal or pelvic peritoneum and occurs mainly in male adolescents and young adults (peak incidence at 20-24 years of age) [1]. Current therapeutic approaches involve the use of multimodal therapeutic regimen, including aggressive polychemotherapy, debulking surgery, and whole abdominal radiation [2].
Apart from this translocation, no other recurrent genomic alteration has been reported in DSRCT cases. Silva et al. [6] detected a somatic amplification involving AURKB and MCL1 genes in one patient, and La Starza et al. [7] found specific genomic imbalances, including gain at chromosome 3 reported in two cases and chromosome 5 polysomy in one case. In terms of point mutations, the data is even scarcer. Variants of unknown clinical significance were reported in ARID1A and RUNX1 genes in one patient [6], whereas in another study, no mutations were detected in a panel of 29 genes evaluated in a cohort of 24 DSRCT cases [8]. This limited genomic information about DSRCT impairs new and more efficient therapeutic opportunities for the young patients affected with this rare tumor.
Here, aiming to contribute with the knowledge of the genomic abnormalities that underlies DSRCT, we performed a comprehensive genomic profiling using a family based approach in one case of DSRCT diagnosed in a young male patient with pelvic tumor, who presented excellent outcome sustained for above 3 years. We identified somatic mutations in a genomic background of rare germline variants either homozygous or as compound heterozygous inheritance, which can improve the understanding of the genetic basis of this rare tumor. We have also generated the profile of genomic imbalances, which was confirmed by whole-exome sequencing. In addition, as we were able to define the precise genomic breakpoints of the EWS-WT1 translocation by whole-genome and Sanger sequencing, we managed to establish a personalized strategy for tracking DNA tumor traces in plasma, allowing an accurate monitoring of tumor recurrence.
Overall, our analysis revealed potential genes and pathways associated with this rare sarcoma and demonstrated the feasibility of using genomic profiling for the benefit of patients affected by rare tumors by developing a personalized monitoring strategy.

DNA extraction
Tumor tissues and blood samples were collected following the technical and ethical procedures of A.C. Camargo Tumor Bank, registered at National Council for Ethics in Research by the number B001 [9]. Genomic DNA and plasma DNA were extracted in DNA and RNA Bank [10] using QIASymphony DNA Mini kit (QIAGEN, Hilden, Germany) for tumor and leukocyte DNA and QIAamp DNA blood Midi kit (QIAGEN, Hilden, Germany) for plasma DNA, following standard procedures.

Target sequencing
Target sequencing was performed using the Ion Ampli-Seq™ Comprehensive Cancer Panel, which comprises all exons from 409 genes associated with different types of tumors (AmpliSeq, Ion Torrent™). This panel, based on multiplex PCR, was performed with as little as 40 ng of DNA from the tumor sample. Library was prepared based on Ion AmpliSeq™ Library Preparation protocol and sequenced at Ion Proton™ platform (Ion Torrent™), according to the manufacturer's instructions.

Whole-exome sequencing
Whole-exome sequencing of the tumor and leukocyte DNA samples from the patient and his mother were performed using the TargetSeq™ Exome Enrichment Kit (Life Technologies), followed by paired-end sequencing (75 × 50) in SOLiD 5500xl System (Life Technologies). Leukocyte DNA sample from his father was submitted for wholeexome sequencing using Ion Xpress™ Plus Fragment Library kit and Ion TargetSeq™ Exome Enrichment Kit (Life Technologies), followed by sequencing at Ion Proton™ platform (Ion Torrent™), according to the manufacturer's instructions. Single-end sequencing was performed on an Ion PI™ Chip v2 with 200 pb sequencing kit (Ion Torrent™).

Whole-genome sequencing
Whole-genome sequencing was performed by mate-paired DNA libraries prepared using the 5500 SOLiD™ Mate-Paired Library Construction Kit (Life Technologies), following the manufacturer's instructions. Briefly, genomic DNA from tumor and leukocyte of the patient was sheared with a Covaris sonicator into approximately 2 Kb fragments, circularized with mate-paired adaptors, nick-translated and digested, incorporated with sequencing adaptors and individual barcodes (distinct barcodes were used for the tumor and leukocyte DNA), and submitted to emulsion PCR. Mate-paired sequencing (60 × 60) was performed in a SOLiD 5500xl System (Life Technologies).

Sanger sequencing
For validation, primers were designed flanking the variants, in order to generate fragments of nearly 400 bp. PCR reactions were performed using GoTaq® Green Master Mix (Promega, Madison, WI, USA), using 15 ng of DNA, with 300 nM of each primer, for a final reaction volume of 20 μL. Approximately 200 ng of PCR-amplified fragments were purified with ExoSAP-IT (USB Corporation, Cleveland, OH, USA) and sequenced in both directions. All alterations were evaluated in all four samples (mother, father, and patient's leukocyte and tumor). Products were analyzed using an ABI 3130xl DNA sequencer (Applied Biosystems, Foster City, CA, USA), and sequences were aligned with the respective gene reference sequence using CLC Genomics Workbench Software (QIAGEN, Hilden, Germany).

Comparative genome hybridization based on microarrays (array CGH)
Comparative genomic hybridization based on microarrays was performed in a commercial whole-genome 180 K platform containing 180,000 oligonucleotide probes (Agilent Technologies; design 22060), using DNA from the tumor sample. Reference DNA was a commercially available human pool of samples from multiple anonymous healthy donors (Promega Corporation). Technical procedures are described in Torrezan et al. [11]. Hybridization and washing were performed as recommended by the manufacturer. Scanned images were processed using Feature Extraction 10.7.3.1 software (Agilent Technologies), and array CGH analysis was conducted with Nexus Copy Number software 7.0 (Biodiscovery). We used the FASST2 segmentation algorithm, according to the following settings: minimum of five consecutive probes (effective resolution of 70 Kb for CNA calling), significance threshold set at 10 −8 , and threshold log 2 Cy3/Cy5 of 0.33 and −0.3 for gains for loss, respectively, and 1.2 and −1.1 for high copy number gains and homozygous losses, respectively. All copy number alterations are reported in the Database of Genomic Variants [12].

Whole-exome sequencing
Sequencing reads from SOLiD 5500xl System were mapped to the reference genome (GRCh37/hg19) with Lifescope (LifeScope™ Genomic Analysis Software v2.5.1). Sequencing reads from Ion Proton™ were mapped to the reference genome (GRCh37/hg19) with TMAP (torrent mapper 4.2.18). Sequence variants (SNVs and indels) were identified following the GATK protocol vs3.2-2-gec30cee [13] for SOLiD 5500XL System and with Torrent Variant Caller 4.0-5, followed by GATK protocol vs3.2-2-gec30cee for Ion Proton reads. Variants were annotated using SnpEff version 3.5d (build 2014-03-05) [14] and an in-house developed script. Identified variants were compared to dbNSFP version 2.4 [15,16]; COSMIC v69 [17]; 1000 Genomes [18]; NHLBI GO Exome Sequencing Project version ESP6500SI-V2 [19]; HapMap [20]; and dbSNP version 138 [21,22] for further annotation. Somatic variants were defined for regions with a minimum coverage of 10× for both tumor and leukocyte sample from the patient and minimum variant frequency of 20% in the tumor only. De novo variants were defined with a minimum coverage of 10× for leukocyte samples of the patient and his parents and with a minimum variant frequency of 30%. To identify rare polymorphisms inherited in homozygosity, we selected variants with a minimum coverage of 10× for leukocyte samples of the patient and his parents, detected in heterozygosis in both parents and homozygosis in the patient presenting a minor allele frequency ≤10% in the public databases (1000 Genome [18], NHLBI GO Exome Sequencing Project [19], and HapMap [20]). To identify germline compound heterozygosis cases, we identified genes with two distinct heterozygous mutations in the patient, where each variant was exclusively present in one of his parents in heterozygosity and detected in regions with a minimum coverage of 10× for all leukocyte samples. We discarded variants that are detected with a minor allele frequency above 10% in the 1000 Genome Project [18]. Copy number alterations were detected using the bioinformatics packages Excavator (version 2.2) [23] and cn.mops (version 1.8.9) [24] by comparing exome data from the tumor to leukocyte from the patient. For visualization, we used circos 0.67-7 package [25].

Ingenuity Pathway Analysis
We applied the core analysis of Ingenuity Pathway Analysis (IPA) system (QIAGEN, Germantown, MD, USA) to identify gene interaction networks.

Digital droplet PCR
Digital droplet PCR assays were carried out using the QX200™ Droplet Digital™ (ddPCR™) System (Bio-Rad). A primer-probe assay labeled with FAM was designed for the amplification of wild-type WT1 gene and a primerprobe assay labeled with HEX was designed for the amplification of the gene fusion event (EWS-WT1). For the amplification reaction, we used 1× ddPCR Supermix, 1× primer-probe assay (FAM), 1× primer-probe assay (HEX), and 4 μl of DNA. Droplet generation, PCR amplification, and droplet counting were performed following the manufacturer's recommendations.
DNA samples from tumor tissue (60 ng-6 pg) and from leukocytes (6 ng) were used as positive and negative controls, respectively. We loaded 4 μl of 30 μl (1/8) of cfDNA samples from the patient. We also used DNA and cell-free DNA extracted from leukocytes and plasma, respectively, from healthy donors as negative controls and performed non-template control. All reactions were performed with at least two replicates.

Results
In this study, we performed a comprehensive genomic profiling of one case of desmoplastic small round cell tumor (DSRCT) by a combination of targeted sequencing, array CGH, whole genome, and whole exome applied in a family based format. The patient studied here is a 26-yearold male, who presented at A.C. Camargo Cancer Center in November, 2011, with a large abdominal mass and a small nodule on the pelvic region. Staging images showed that the disease was limited to the abdominal cavity. He underwent a CT-guided biopsy that revealed a desmoplastic small round cell tumor (DSRCT), showing positivity for EMA, desmin, and nuclear staining for WT1 (carboxy-terminus antibody). FISH analysis was positive for EWS translocation.
The patient started systemic treatment with 4 cycles of vincristine, cyclophosphamide, and doxorubicin (VAC) and alternated with ifosfamide, carboplatin, and etoposide (ICE). After 4 cycles, the best response was stable disease, with minor reduction in the tumor dimensions. He underwent complete surgical cytoreduction, with resection of the large mass and resection of peritoneal implant on the pelvic region, and hyperthermic intraperitoneal chemotherapy (HIPEC), with cisplatin and doxorubicin. The patient presented a complete recovery from these procedures. After that, the patient received four more cycles of chemotherapy and total abdominal irradiation (total of 30 Gy). After a follow-up of 48 months since surgery, the patient is asymptomatic, with no signs of disease (Additional file 1: Figure S1).
The comprehensive genomic profiling was carried out in multiple fronts. To identify actionable mutations, we carried out targeted sequencing of the most important actionable genes in the tumor sample. Further, to identify genes and mutations possibly involved with tumor onset and predisposition, we performed whole exome sequencing, using the DSRCT tumor sample and blood samples from the patient and his both parents (Additional file 2: Table S1). Array CGH was performed in tumor sample to identify structural rearrangements and copy number imbalances. Whole-genome sequencing defined one tumor marker used to precisely monitor patient after treatment.

Producing a portrait of somatically acquired variants
To identify actionable mutations or pathways related to DSRCT in this patient, we initially performed targeted sequencing using a cancer-oriented gene panel composed of 409 genes in the tumor sample (Comprehensive Cancer Panel -Thermo Scientific). Since no actionable mutation for targeted therapy was detected in the tumor, we investigated the complete landscape of somatic mutations by whole-exome sequencing (WES) the tumor and the patient's leukocyte. The analysis of the tumor revealed 15 somatic acquired mutations, 12 of which were proteinaffecting variants (validated by capillary sequencing) including one non-sense mutation in the ZNF808 gene, one nucleotide change at the 3' splice site of RIMS4, and 10 missense mutations considered disease associated by at least one pathogenicity prediction program (Table 1).
Gene ontology enrichment analysis of the 12 genes harboring protein-affecting somatic mutations revealed several biological processes such as muscle tissue/organ development (ZFPM2 and MEGF10), which is related to the mesothelium origin of this tumor, cell adhesion (DPP4, CDH9, CNTNAP4, MEGF10), response to mechanic stimulus (MEIS2, TRPA4), and response to abiotic stimulus (DPP4, MEIS2, TRPA1) (Additional file 3: Table S2). Intriguingly, biological network analysis obtained using Ingenuity Pathway Analysis (IPA) interconnected the 15 genes harboring somatic mutations in a single network associated with cell death and survival, cell damage or degeneration, and nervous system development and function (Fig. 1a).

Searching for germline variants associated with DSRCT
In an attempt to identify de novo germline mutations and inherited variants potentially associated with DSRCT, we performed whole-exome sequencing of the leukocyte DNA from the patient's parents (Additional file 2: Table S1). For selecting genetic variants, a minimum coverage of 10× in all samples with at least 30% frequency of the variant allele was considered. Based on these criteria, no de novo variants were found.
Next, we explored the possibility of finding genetic variants associated with DSRCT in an autosomal recessive model of inheritance. First, we searched for polymorphisms (MAF ≤10%) occurring in homozygosity in the patient that were inherited from both heterozygous parents. Four polymorphisms inherited in homozygosis were found in ADAMTS12, RASSF1, VEZT, and ISX genes (population frequency ranging from 2.0 to 7.8%). All polymorphisms lead to missense alterations, two of them reported as possibly disease associated by at least one pathogenicity prediction program (Table 2). Interestingly, biological network analysis showed an interconnection between the four genes in a single network from IPA, showing association with cell cycle and digestive system development and function (Fig. 1b).
Next, to identify candidate genes affected by compound heterozygosity, we looked for genes containing two distinct variants inherited independently from each parent. We could confirm 11 genes affected by germline compound heterozygous mutations, in which one variant allele was inherited from the mother and the other variant allele was inherited from the father (MAF ≤10%) ( Table 3). We speculate that these genes could be associated with risk to DSRCT development in an autosomal recessive model of inheritance.
Gene ontology analysis of these 11 genes showed enrichment of biological processes related to muscle tissue development (LAMB2, SYNE1, and TTN), morphogenesis (C2CD3 and TTN), and cell cycle (RSPH1, TTN, and SPICE1) (Additional file 4: Table S4). Functional analysis of IPA revealed that 9 of the 11 genes are interconnected in a single network related to cancer, organismal injury and abnormalities, and gastrointestinal disease.

Copy number alterations
Genomic copy number alterations (CNA) were investigated in the DSRCT tumor sample using array CGH in a 180-K platform. Few copy number alterations were detected (Fig. 2a), including aneuploidies such as gain of chromosomes 5 and 18, and 11p, 13q, and 22q deletions. Only one small focal homozygous deletion was identified in a segment of~1.3 Mb at 9p22.2 (chr9:17,106,384-18,449,088; hg19), encompassing the CNTLN and SH3GL2 genes. The detection of copy number losses affecting terminal segments of 11p and 22q suggested the presence of the chromosomal translocation t(11;22)(p13;q12) involving EWSR1 and WT1 genes.
Additionally, we used NGS data from the WES of the tumor and patient's leukocyte to search for CNAs, and 98.4% of the events detected by array CGH were validated, demonstrating the efficacy of WES for identification of a wide range of somatic events, besides point mutations and indels ( Fig. 2b; Additional file 5: Table S3).

Establishment of a personalized monitoring strategy based on detection of circulating tumor DNA (ctDNA)
The genomic translocation t(11;22)(p13;q12), which is considered the molecular hallmark of this tumor type, was initially detected by FISH analysis for diagnostic purposes. Furthermore, the same translocation was also detected in the array CGH and confirmed by whole-genome sequencing. The mate-pair approach used for whole-genome sequencing followed by validation by PCR and Sanger sequencing allowed the precise delimitation of the chromosomal breakpoints (Fig. 3). This somatic fusion event was then used as a tumor biomarker for monitoring the patient along the follow-up period using a liquid biopsy-based approach. We searched for traces of circulating tumor DNA in plasma samples collected 17, 36, 42, and 47 months after surgery. No signs of the translocation in any of the post-treatment plasma samples were detected (Fig. 3) using digital droplet PCR (ddPCR). The absence of the biomarker in any of the post-treatment plasma samples is in agreement with favorable clinical response of the patient, showing longterm disease-free survival and no sign of disease recurrence. As a positive control, we confirmed the presence of cell-free circulating DNA by detecting non-rearranged WT1 gene in all plasma samples.
Despite the unavailability of plasma samples at the moment of diagnosis and before treatment (neoadjuvant/ adjuvant chemotherapy and surgery) to be used as a baseline positive control, we carried out stringent control assays. We tested the robustness of the approach by evaluating six different amounts of tumor DNA, ranging from 60 ng to 6 pg, and searched for both the tumor-specific fusion event and for the non-rearranged WT1 gene. We achieved a linear quantification of tumor marker and wild-type DNA, and we were able to discriminate the presence of the fusion event even in extremely low quantity of DNA (6 pg) (Fig. 3d). To test the specificity of the ddPCR assay, we used leukocyte sample of the patient and also DNA from plasma and leukocyte samples of healthy donors, and the absence of detection of tumor-specific fusion confirmed the assay specificity for tumor DNA, without false positive signals (Additional file 6: Figure S2). Further, we mimicked a situation of circulating tumor DNA in body fluid by mixing 1 part of tumor DNA to 100 and 1000 parts of leukocyte DNA, then screened for the fusion event using 1 ng of this DNA mix. Typically, the detection of circulating tumor DNA has been reported as a fraction between 0.1 and 90% of plasma DNA for cancer patients, depending on tumor type and patient characteristics [27]. Here, we were able to detect the fusion event in the 1.0 and 0.1% fractions, supporting the reliability of our established approach for patient surveillance based on liquid-biopsy screening.

Discussion
DSRCT is an aggressive tumor not yet broadly investigated with the powerful genomic tools available. Mainly due to the rarity of this disease, only a few studies investigated molecular alterations in this tumor type. Shukla et al. [8] screened the occurrence of 275 COSMIC mutations in 29 oncogenes and found no alterations in any of the 24 DSRCT samples investigated. More recently, a study based on targeted exome sequencing of six adult patients with pediatric-type malignancies found AURKB and MCL1 amplifications and variants of unknown clinical significance in ARID1A and RUNX1 genes in one DSRCT [6].
Here, by performing a comprehensive screening of the genomic alterations in a family based approach, we identified somatic and germline variants possibly associated with DSRCT. In our study, the use of a commercially available gene panel did not show to be an adequate strategy. Based on this finding, one can argue that screening by commercially available gene panels is not an effective approach for most cases of rare tumors, since the targeted genes represented in these panels are usually those well characterized in common solid tumors. On the other hand, the use of WES not only revealed mutated genes but also showed robustness for detecting DNA copy number alterations. Concordance rates between WES and array CGH (the gold standard for CNA screening) were above 98% (Fig. 2). Among the CNAs not detected in the WES analysis, two of them were mapped to regions not covered by the library probes and five were low-level mosaic alterations (Additional file 3: Table S2). Thus, if we consider CNAs mapping to WES target regions and nonmosaic CNAs, concordance rates were above 99.8%.
In total, 38 somatically acquired alterations, including point mutations (15) and CNAs (23), were detected. The small number of somatic alterations identified here is in agreement with what is expected for pediatric tumors [28]. Moreover, given the occurrence of the driver EWS-WT1 fusion protein, additional oncogenic mutations for tumor onset is probably less necessary.  In blue, data from array CGH and in green data from WES. A great overlap of CNA detection can be observed using both approaches We identified 12 genes affected by somatic mutations possibly involved with the disease. These genes are involved with cellular development and morphology that are pathways in which WT1 gene plays an important role. Similarly, the genes affected by compound heterozygous mutations showed enrichment in biological processes of muscle tissue development and morphogenesis. Altogether, these data suggests that disruption of the embryonic cellular development process is involved with DSRCT onset, which is commonly seen in pediatric tumors.
Another interesting finding is that among the set of somatically mutated genes, 8 showed to be mutated in desmoplastic melanoma samples from TCGA project [29,30], in frequencies ranging from 5 to 25% of the 20 samples interrogated. This data suggests that mutation in these genes might be involved with the desmoplastic phenotype, seen in both tumor types.
Additionally, 7 out of 15 genes harboring somatic mutations (CHL1, MEGF10, MEIS2, MYH8, RIMS4, TBPL1, and ZFPM2) are regulated by the same transcription factor, LEF1 (p < 0.001 by enrichment analysis), which, in turn, is regulated by WT1 [31]. We therefore postulate that DSRCT tumors presenting increased activity of WT1 due to EWS-WT1 fusion might upregulate the expression of several genes mediated by LEF1 transcription factor. However, the accumulation of mutations in this set of genes regulated by LEF1 activation and its interrelation with the EWS-WT1 fusion protein remains to be addressed.
Finally, the definition of the precise genomic breakpoint of the t(11;22)(p13;q12) translocation by whole-genome and Sanger sequencing enabled the development of a personalized tool to precisely monitor the presence of ctDNA in plasma samples during the patient's follow-up. Detection of tumor-specific genomic rearrangements has been shown as a sensitive and specific method for monitoring of disease status of cancer patients [32][33][34][35] and has clear advantages over point mutations concerning the specificity of detection [33]. Here, we applied ctDNA screening in plasma samples collected in three clinical appointments during 3 years after surgery and, up to now, we did not detect the presence of this tumor marker. These results are in agreement with imaging exams (CT scan and PET scan) showing no signs of disease and also with good overall clinical condition of the patient and highlight the applicability of using genomic rearrangement for building personalized tool for patient surveillance. After defining the genomic breakpoint of EWS-WT1 fusion, DSRCT patients can benefit from a highly specific test that has the advantage of a rapid turnaround time and potentially higher sensitivity in detecting disease progression earlier than imaging exams or other cancer antigens measurements, as reported for other tumor types [36,37]. Monitoring through liquid biopsy is particularly attractive for solid tumors, which cannot be repeatedly sampled without more invasive procedures. Considering the rarity of this subtype of sarcoma and the lack of effective treatment, the detection of specific tumor marker and the monitoring of its persistence can improve the identification of patients with worse prognosis to tailor the treatment more properly. Thus, the perspective is to employ the approach used here for new patients and improve the outcome for those with worse prognosis.

Conclusion
To our knowledge, this is the first comprehensive genomic characterization of one DSRCT case. Continuous efforts to establish the genomic landscape of rare diseases, frequently neglected in large sequencing consortiums, are highly significant to improve the knowledge of defective pathways involved with tumor onset in general, in addition to the strong potential of revealing druggable targets for clinical use.