Integration of genomic copy number variations and chemotherapy-response biomarkers in pediatric sarcoma

Germline Genetic Mutations in Adult Patients with Sarcoma: Insight into the Middle East Genetic Landscape

Data on germline mutations in soft tissue and bone sarcomas are scarce. We sought to identify the prevalence of germline mutations in adult sarcoma patients treated at a tertiary cancer center. Newly diagnosed patients were offered germline genetic testing via an 84-gene panel. The prevalence of pathogenic germline variants (PGVs) and their association with disease-, and patient- related factors are reported. A total of 87 patients were enrolled, the median age was 48 (19-78) years, and 47 (54%) were females. Gastrointestinal stromal tumors (n = 12, 13.8%), liposarcoma (n = 10, 11.5%), and Ewing sarcoma (n = 10, 11.5%) were the main subtypes. A total of 20 PGVs were detected in 18 (20.7%) patients. Variants of uncertain significance, in the absence of PGVs, were detected in 40 (45.9%) patients. Young age (p = 0.031), presence of a second primary cancer (p = 0.019), and female gender (p = 0.042) were correlated with the presence of PGVs. All identified PGVs have potential clinical actionability and cascade testing, and eight (44.44%) suggested eligibility for a targeted therapy. Almost one in five adult patients with soft tissue and bone sarcomas harbor pathogenic or likely pathogenic variants. Many of these variants are potentially actionable, and almost all have implications on cancer screening and family counselling. In this cohort from the Middle East, younger age, presence of a second primary tumor, and female gender were significantly associated with higher PGVs rates. Larger studies able to correlate treatment outcomes with genetic variants are highly needed.

Tyrosine kinase inhibitors in osteosarcoma: Adapting treatment strategiesa

Osteosarcoma (OS) is an aggressive primary bone malignancy that metastasizes rapidly. The standard of care has changed little over the previous four decades, and survival rates have plateaued. In this context, tyrosine kinase inhibitors (TKIs) emerge as potential treatments. A literature search was conducted to collect data related to receptor tyrosine kinase genetic alterations and expression in OS specimens. Gene amplification and protein expression of these receptors were linked to prognosis and tumor behavior. Relevant TKIs were evaluated as monotherapies and as parts of combination therapies. Certain TKIs, such as apatinib, regorafenib, and cabozantinib, present a potential therapeutic avenue for OS patients, especially when combined with chemotherapy. Producing long-lasting responses and enhancing quality of life remain key goals in OS treatment. To this effect, optimizing the use of TKIs by identifying biomarkers predictive of response and assessing promising TKIs in larger-scale trials to validate the efficacy and safety outcomes relative to these drugs reported in phase II clinical trials. To this effect, it is necessary to identify biomarkers predictive of response to TKIs in larger-scale trials and to validate the efficacy and safety of these drugs reported in phase II clinical trials.

Integrative Multi-OMICs Identifies Therapeutic Response Biomarkers and Confirms Fidelity of Clinically Annotated, Serially Passaged Patient-Derived Xenografts Established from Primary and Metastatic Pediatric and AYA Solid Tumors

Establishment of clinically annotated, molecularly characterized, patient-derived xenografts (PDXs) from treatment-naïve and pretreated patients provides a platform to test precision genomics-guided therapies. An integrated multi-OMICS pipeline was developed to identify cancer-associated pathways and evaluate stability of molecular signatures in a panel of pediatric and AYA PDXs following serial passaging in mice. Original solid tumor samples and their corresponding PDXs were evaluated by whole-genome sequencing, RNA-seq, immunoblotting, pathway enrichment analyses, and the drug−gene interaction database to identify as well as cross-validate actionable targets in patients with sarcomas or Wilms tumors. While some divergence between original tumor and the respective PDX was evident, majority of alterations were not functionally impactful, and oncogenic pathway activation was maintained following serial passaging. CDK4/6 and BETs were prioritized as biomarkers of therapeutic response in osteosarcoma PDXs with pertinent molecular signatures. Inhibition of CDK4/6 or BETs decreased osteosarcoma PDX growth (two-way ANOVA, p < 0.05) confirming mechanistic involvement in growth. Linking patient treatment history with molecular and efficacy data in PDX will provide a strong rationale for targeted therapy and improve our understanding of which therapy is most beneficial in patients at diagnosis and in those already exposed to therapy.

Nefarious NTRK oncogenic fusions in pediatric sarcomas: Too many to Trk

Neurotrophic Tyrosine Receptor Kinase (NTRK) genes undergo chromosomal translocations to create novel open reading frames coding for oncogenic fusion proteins; the N-terminal portion, donated by various partner genes, becomes fused to the tyrosine kinase domain of either NTRK1, NTRK2, or NTRK3. NTRK fusion proteins have been identified as driver oncogenes in a wide variety of tumors over the past three decades, including Pediatric Gliomas, Papillary Thyroid Carcinoma, Spitzoid Neoplasms, Glioblastoma, and additional tumors. Importantly, NTRK fusions function as drivers of pediatric sarcomas, accounting for approximately 15% of childhood cancers including Infantile Fibrosarcoma (IFS), a subset of pediatric soft tissue sarcoma (STS). While tyrosine kinase inhibitors (TKIs), such as larotrectinib and entrectinib, have demonstrated profound results against NTRK fusion-positive cancers, acquired resistance to these TKIs has resulted in the formation of gatekeeper, solvent-front, and compound mutations. We present a comprehensive compilation of oncogenic fusions involving NTRKs focusing specifically on pediatric STS, examining their biological signaling pathways and mechanisms of activation. The importance of an obligatory dimerization or multimerization domain, invariably donated by the N-terminal fusion partner, is discussed using characteristic fusions that occur in pediatric sarcomas. In addition, examples are presented of oncogenic fusion proteins in which the N-terminal partners may contribute additional biological activities beyond an oligomerization domain. Lastly, therapeutic approaches to the treatment of pediatric sarcoma will be presented, using first generation and second-generation agents such as selitrectinib and repotrectinib.

Replication collisions induced by de-repressed S-phase transcription are connected with malignant transformation of adult stem cells

Transcription replication collisions (TRCs) constitute a major intrinsic source of genome instability but conclusive evidence for a causal role of TRCs in tumor initiation is missing. We discover that lack of the H4K20-dimethyltransferase KMT5B (also known as SUV4-20H1) in muscle stem cells de-represses S-phase transcription by increasing H4K20me1 levels, which induces TRCs and aberrant R-loops in oncogenic genes. The resulting replication stress and aberrant mitosis activate ATR-RPA32-P53 signaling, promoting cellular senescence, which turns into rapid rhabdomyosarcoma formation when p53 is absent. Inhibition of S-phase transcription ameliorates TRCs and formation of R-loops in Kmt5b-deficient MuSCs, validating the crucial role of H4K20me1-dependent, tightly controlled S-phase transcription for preventing collision errors. Low KMT5B expression is prevalent in human sarcomas and associated with tumor recurrence, suggesting a common function of KMT5B in sarcoma formation. The study uncovers decisive functions of KMT5B for maintaining genome stability by repressing S-phase transcription via control of H4K20me1 levels.

An international working group consensus report for the prioritization of molecular biomarkers for Ewing sarcoma

The advent of dose intensified interval compressed therapy has improved event-free survival for patients with localized Ewing sarcoma (EwS) to 78% at 5 years. However, nearly a quarter of patients with localized tumors and 60–80% of patients with metastatic tumors suffer relapse and die of disease. In addition, those who survive are often left with debilitating late effects. Clinical features aside from stage have proven inadequate to meaningfully classify patients for risk-stratified therapy. Therefore, there is a critical need to develop approaches to risk stratify patients with EwS based on molecular features. Over the past decade, new technology has enabled the study of multiple molecular biomarkers in EwS. Preliminary evidence requiring validation supports copy number changes, and loss of function mutations in tumor suppressor genes as biomarkers of outcome in EwS. Initial studies of circulating tumor DNA demonstrated that diagnostic ctDNA burden and ctDNA clearance during induction are also associated with outcome. In addition, fusion partner should be a pre-requisite for enrollment on EwS clinical trials, and the fusion type and structure require further study to determine prognostic impact. These emerging biomarkers represent a new horizon in our understanding of disease risk and will enable future efforts to develop risk-adapted treatment.

Single-cell transcriptional profiling reveals heterogeneity and developmental trajectories of Ewing sarcoma

PURPOSE

Ewing sarcoma (EwS) is an aggressive malignant neoplasm composed of small round cells. The heterogeneity and developmental trajectories of EwS are uncertain.

METHODS

Single-cell RNA sequencing was performed on 4 EwS tumor tissue samples, and 3 transcriptional atlases were generated. K-nearest neighbor algorithm was used to predict the origin of tumor cells at single-cell resolution. Monocle2 package was used to perform pseudotime trajectory analysis in tumor cells. Differentially expressed genes were compared against those in all other clusters via the FindMarkers function, and then they were subjected to GO analysis using clusterProfiler package.

RESULTS

Combined with the results of k-nearest neighbor algorithm and pseudotime trajectory analysis in tumor cells, we thought meningeal EwS originated from neural crest cells during epithelial to mesenchymal transition and simulated the process of neural crest cell lineage differentiation. But for perirenal EwS and spinal EwS, we hypothesized that after the neural crest cell lineage mutated into them, the tumor cells did not maintain the differentiation trajectory of neural crest cell lineage, and the development trajectory of tumor cells became chaotic. GO analysis results showed that interferon signaling pathway-related biological processes play an essential role in the tumorigenesis and tumor progression process of EwS, and among these biological processes genes, JAK1 gene up-regulated most significantly and highly expressed in all tumor cells. Ruxolitinib was used to explore the function of JAK1. Targeting JAK1 can promote apoptosis of EwS tumor cells, inhibit the migration and invasion of EwS tumor cells, and inhibit cell proliferation by inducing cell cycle S phase arrest.

CONCLUSION

EwS was derived from neural crest cell lineage with variable developmental timing of oncogenic conversion, and the JAK1 might be a candidate for therapeutic targets of EwS.

Establishment and characterization of NCC-UPS4-C1: a novel cell line of undifferentiated pleomorphic sarcoma from a patient with Li-Fraumeni syndrome

Li-Fraumeni syndrome (LFS) is a cancer predisposition syndrome caused by a germline mutation of the TP53. The lifetime risk of cancer in individuals with LFS is ≥ 70% for men and ≥ 90% for women. Undifferentiated pleomorphic sarcoma (UPS) is one of the core cancers associated with LFS. UPS is a subtype of undifferentiated soft tissue sarcoma that shows no identifiable line of differentiation. The standard curative treatment for UPS is complete surgical resection. However, local recurrence and distant metastasis to the lung can usually be found after resection of the UPS. Therefore, a novel treatment strategy for patients with UPS is required. Although well characterized, patient-derived tumor cell lines facilitate the high-throughput screening of a large number of drugs, and no sarcoma cell lines derived from a patient with LFS have been registered in public cell banks. Thus, this study aimed to establish a novel, well-characterized UPS cell line from a patient with LFS. From surgically resected UPS tumor tissues, we established the first UPS cell line from a patient with LFS and named it NCC-UPS4-C1. NCC-UPS4-C1 harbored copy number alterations and had the TP53 tumor suppressor gene mutation. The cells exhibited constant cell growth and invasive ability. This well-characterized NCC-UPS4-C1 cell line was then utilized for high-throughput screening of 214 anti-cancer drugs, and two effective drugs were identified. One of the two drugs, romidepsin, was commonly effective for the NCC-UPS1-C1, NCC-UPS2-C1, and NCC-UPS3-C1 cell lines that we previously reported; a potential drug for the treatment of UPS was suggested using well-characterized UPS cell lines. These data indicate that NCC-UPS4-C1, which is the first sarcoma cell line established from a patient with LFS, enables researchers to conduct vigorous preclinical research on UPS.

Case report of recurrent fibromatosis with laryngeal involvement: Treatment based on network analyses of NGS data

Aggressive fibromatosis (AF) is a rare, benign neoplasm originating from musculoaponeurotic stromal structures characterized by aggressive growth and infiltration of local tissues. To date, only six previous cases of AF involving the larynx have been reported. The present case was that of a 70-year-old female patient with a 5-year history of hoarseness and an enlarging neck mass consistent with aggressive recurrent fibromatosis. MRI displayed a large, solid mass arising from the left anterior cervical space displacing the trachea and upper airway to the right. At one year after the initial radical resection, the patient presented with recurrence. A second radical excision was performed not including a laryngectomy. Radiation therapy was considered for possible local treatment to prevent subsequent tumor recurrence. According to the literature, AF has recurrence rates of up to 40-70% within 18 months. There is only sparse literature to guide treatment. Using the mutations detected in the patient's AF tissue, an Ingenuity Pathway Analysis (IPA) was used to guide treatment of the recurrence. In the present case, the IPA analysis indicated the use of pazopanib to treat the patient's cancer. In general, surgery appears to be the treatment of choice for head and neck AF, but the management of recurrence is controversial.

BRCA1/ATF1-Mediated Transactivation is Involved in Resistance to PARP Inhibitors and Cisplatin

Homologous recombination (HR)-deficient cells are sensitive to PARP inhibitors through a synthetic lethal effect. We previously developed an HR activity assay named Assay of Site-Specific HR Activity (ASHRA). Here, we evaluated the HR activity of 30 missense variants of BRCA1 by ASHRA and found that several BRCA1 variants showed intermediate HR activity, which was not clearly discerned by our previous analyses using a conventional method. HR activity measured by ASHRA was significantly correlated with sensitivity to olaparib. However, cells expressing the severely HR-deficient BRCA1-C61G variant were resistant to olaparib, and resistance was dependent on high expression of activating transcription factor 1 (ATF1), which binds to BRCA1 and activates the transcription of target genes to regulate cell proliferation. The BRCA1-C61G variant bound to ATF1 and stimulated ATF1-mediated transactivation similar to wild-type BRCA1. High expression of ATF1 conferred resistance to olaparib and cisplatin activating BRCA1/ATF1-mediated transcription without affecting HR activity in BRCA2-knockdown or RAD51-knockdown cells, but not in BRCA1-knockdown cells. These results suggest that ASHRA is a useful method to evaluate HR activity in cells and to predict the sensitivity to PARP inhibitors. The expression level of ATF1 might be an important biomarker of the effect of PARP inhibitors and platinum agents on HR-deficient tumors with the BRCA1-C61G variant or alteration of non-BRCA1 HR factors such as BRCA2 and RAD51.

Significance

ASHRA could evaluate HR activity in cells and predict the sensitivity to PARP inhibitors. High expression level of ATF1 may predict the resistance of BRCAness tumors with alterations of non-BRCA1 HR factors to PARP inhibitors and platinum agents.

Integrative analysis of prognostic long non-coding RNAs with copy number variation in bladder cancer

Copy number variations (CNVs), which can affect the role of long non-coding RNAs (lncRNAs), are important genetic changes seen in some malignant tumors. We analyzed lncRNAs with CNV to explore the relationship between lncRNAs and prognosis in bladder cancer (BLCA). Messenger RNA (mRNA) expression levels, DNA methylation, and DNA copy number data of 408 BLCA patients were subjected to integrative bioinformatics analysis. Cluster analysis was performed to obtain different subtypes and differently expressed lncRNAs and coding genes. Weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression gene and lncRNA modules. CNV-associated lncRNA data and their influence on cancer prognosis were assessed with Kaplan-Meier survival curve. Multi-omics integration analysis revealed five prognostic lncRNAs with CNV, namely NR2F1-AS1, LINC01138, THUMPD3-AS1, LOC101928489,and TMEM147-AS1,and a risk-score signature related to overall survival in BLCA was identified. Moreover, validated results in another independent Gene Expression Omnibus (GEO) dataset, GSE31684, were consistent with these results. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the mitogen-activated protein kinase (MAPK) signaling pathway, focal adhesion pathway, and Janus kinase-signal transducers and activators of transcription (JAK-STAT) signaling pathway were enriched in a high-risk score pattern, suggesting that imbalance in these pathways is closely related to tumor development. We revealed the prognosis-related lncRNAs by analyzing the expression profiles of lncRNAs and CNVs, which can be used as prognostic biomarkers for BLCA.

PAX3/7-FOXO1 fusion-negative alveolar rhabdomyosarcoma in Schuurs-Hoeijmakers syndrome

PAX3/7-FOXO1 fusion-negative alveolar rhabdomyosarcoma (ARMS) developed in a patient presenting with intellectual disability and dysmorphic facial features. Whole exome sequencing analysis of a germline sample identified a PACS1 c.607 C>T de novo variant and the patient was diagnosed with Schuurs-Hoeijmakers syndrome (SHS). SHS is a rare disease characterized by intellectual disability and dysmorphic facial features, among various physical abnormalities, due to PACS1 c.607 C>T de novo variant. Due to the rarity of the SHS, diagnosis based on phenotypic information is difficult. To date, there have been no previous reports describing malignancy associated with SHS. Comprehensive somatic mutation analysis revealed a unique pattern of genetic alterations in the PAX3/7-FOXO1 fusion-negative ARMS tumor, including mutations in the oncogene, HRAS; MYOD1, a molecule essential for muscle differentiation; and KMT2C and TET1, genes encoding factors involved in epigenetic regulation. Although the role of PACS1 in tumorigenesis is unclear, it is reported to function in apoptosis regulation. Our case suggests that PACS1 could have a novel role in oncogenesis.

Multimodal analysis of cell-free DNA whole-genome sequencing for pediatric cancers with low mutational burden

Sequencing of cell-free DNA in the blood of cancer patients (liquid biopsy) provides attractive opportunities for early diagnosis, assessment of treatment response, and minimally invasive disease monitoring. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we introduce an integrated genetic/epigenetic analysis method and demonstrate its utility on 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. Our method achieves sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. Moreover, we benchmark different metrics for cell-free DNA fragmentation analysis, and we introduce the LIQUORICE algorithm for detecting circulating tumor DNA based on cancer-specific chromatin signatures. Finally, we combine several fragmentation-based metrics into an integrated machine learning classifier for liquid biopsy analysis that exploits widespread epigenetic deregulation and is tailored to cancers with low mutation rates. Clinical associations highlight the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent genetic aberrations, and it renders the benefits of liquid biopsy more readily accessible for childhood cancers.

Ameloblastic fibrosarcoma of the maxilla with EGFR exon 20 insertions: Relevance of whole-exome sequencing in molecular understanding and therapeutic proposals for rare cancers

Ameloblastic fibrosarcoma (AFS) is the most common odontogenic sarcoma, but the incidence is relatively low, and its molecular biology is poorly understood. We experienced a young female patient with a rapidly growing soft tissue tumor of the left maxilla, which eventually occupied the left side of the oral cavity. Histologically, the tumor mainly consisted of a proliferation of atypical spindle to polygonal cells without any specific differentiation, but a small number of benign odontogenic epithelial foci mainly in the tumor periphery were also noted; thus, a diagnosis of AFS was made. We performed whole-exome sequencing (WES) on the tumor to investigate its molecular features and identify therapeutic options. We found that the tumor harbored EGFR exon 20 insertions and MDM2 amplification; the former may be a target for newly developed tyrosine kinase inhibitors in case of recurrence. To the best of our knowledge, this is the first case of AFS for which WES was performed and with EGFR mutation. Our case provides new genetic information on AFS and suggests that comprehensive genetic analysis can clarify the molecular biology in rare cancers, potentially leading to the proposal of therapeutic strategies.

Next-Generation Sequencing Identifies Potential Actionable Targets in Paediatric Sarcomas

Background: Bone and soft-tissue sarcomas represent 13% of all paediatric malignancies. International contributions to introduce next-generation sequencing (NGS) approaches into clinical application are currently developing. We present the results from the Precision Medicine program for children with sarcomas at a reference centre. Results: Samples of 70 paediatric sarcomas were processed for histopathological analysis, reverse transcriptase polymerase chain reaction (RT-PCR) and next-generation sequencing (NGS) with a consensus gene panel. Pathogenic alterations were reported and, if existing, targeted recommendations were translated to the clinic. Seventy paediatric patients with sarcomas from 10 centres were studied. Median age was 11.5 years (range 1–18). Twenty-two (31%) had at least one pathogenic alteration by NGS. Thirty pathogenic mutations in 18 different genes were detected amongst the 22 patients. The most frequent alterations were found in TP53, followed by FGFR4 and CTNNB1. Combining all biological studies, 18 actionable variants were detected and six patients received targeted treatment observing a disease control rate of 78%. Extrapolating the results to the whole cohort, 23% of the patients would obtain clinical benefit from this approach. Conclusions: Paediatric sarcomas have a different genomic landscape when compared to adult cohorts. Incorporating NGS targets into paediatric sarcomas’ therapy is feasible and allows personalized treatments with clinical benefit in the relapse setting.

Identification of 6 gene markers for survival prediction in osteosarcoma cases based on multi-omics analysis

Multiple-omics sequencing information with high-throughput has laid a solid foundation to identify genes associated with cancer prognostic process. Multiomics information study is capable of revealing the cancer occurring and developing system according to several aspects. Currently, the prognosis of osteosarcoma is still poor, so a genetic marker is needed for predicting the clinically related overall survival result. First, Office of Cancer Genomics (OCG Target) provided RNASeq, copy amount variations information, and clinically related follow-up data. Genes associated with prognostic process and genes exhibiting copy amount difference were screened in the training group, and the mentioned genes were integrated for feature selection with least absolute shrinkage and selection operator (Lasso). Eventually, effective biomarkers received the screening process. Lastly, this study built and demonstrated one gene-associated prognosis mode according to the set of the test and gene expression omnibus validation set; 512 prognosis-related genes (P < 0.01), 336 copies of amplified genes (P < 0.05), and 36 copies of deleted genes (P < 0.05) were obtained, and those genes of the mentioned genomic variants display close associations with tumor occurring and developing mechanisms. This study generated 10 genes for candidates through the integration of genomic variant genes as well as prognosis-related genes. Six typical genes (i.e. MYC, CHIC2, CCDC152, LYL1, GPR142, and MMP27) were obtained by Lasso feature selection and stepwise multivariate regression study, many of which are reported to show a relationship to tumor progressing process. The authors conducted Cox regression study for building 6-gene sign, i.e. one single prognosis-related element, in terms of cases carrying osteosarcoma. In addition, the samples were able to be risk stratified in the training group, test set, and externally validating set. The AUC of five-year survival according to the training group and validation set reached over 0.85, with superior predictive performance as opposed to the existing researches. Here, 6-gene sign was built to be new prognosis-related marking elements for assessing osteosarcoma cases' surviving state.

Detection of Structural Variants in Circulating Cell-Free DNA from Sarcoma Patients Using Next Generation Sequencing

Circulating tumour DNA (ctDNA) analysis using next generation sequencing (NGS) is being implemented in clinical practice for treatment stratification and disease monitoring. However, using ctDNA to detect structural variants, a common occurrence in sarcoma, can be challenging. Here, we use a sarcoma-specific targeted NGS panel to identify translocations and copy number variants in a cohort of 12 tissue specimens and matched circulating cell-free DNA (cfDNA) from soft tissue sarcoma patients, including alveolar rhabdomyosarcoma (n = 2), Ewing's Sarcoma (n = 2), synovial sarcoma (n = 2), extraskeletal myxoid chondrosarcoma (n = 1), clear cell sarcoma (n = 1), undifferentiated round cell sarcoma (n = 1), myxoid liposarcoma (n = 1), alveolar soft part cell sarcoma (n = 1) and dedifferentiated liposarcoma (n = 1). Structural variants were detected in 11/12 (91.6%) and 6/12 (50%) of tissue and plasma samples, respectively. Structural variants were detected in cfDNA at variant allele frequencies >0.2% with an average sequencing depth of 1026×. The results from this cohort show clinical potential for using NGS in ctDNA to aid in the diagnosis and clinical monitoring of sarcomas and warrant additional studies in larger cohorts.

Sarcoma treatment in the era of molecular medicine

Sarcomas are heterogeneous and clinically challenging soft tissue and bone cancers. Although constituting only 1% of all human malignancies, sarcomas represent the second most common type of solid tumors in children and adolescents and comprise an important group of secondary malignancies. More than 100 histological subtypes have been characterized to date, and many more are being discovered due to molecular profiling. Owing to their mostly aggressive biological behavior, relative rarity, and occurrence at virtually every anatomical site, many sarcoma subtypes are in particular difficult-to-treat categories. Current multimodal treatment concepts combine surgery, polychemotherapy (with/without local hyperthermia), irradiation, immunotherapy, and/or targeted therapeutics. Recent scientific advancements have enabled a more precise molecular characterization of sarcoma subtypes and revealed novel therapeutic targets and prognostic/predictive biomarkers. This review aims at providing a comprehensive overview of the latest advances in the molecular biology of sarcomas and their effects on clinical oncology; it is meant for a broad readership ranging from novices to experts in the field of sarcoma.

Systems Biology Approach Identifies Prognostic Signatures of Poor Overall Survival and Guides the Prioritization of Novel BET-CHK1 Combination Therapy for Osteosarcoma

Osteosarcoma (OS) patients exhibit poor overall survival, partly due to copy number variations (CNVs) resulting in dysregulated gene expression and therapeutic resistance. To identify actionable prognostic signatures of poor overall survival, we employed a systems biology approach using public databases to integrate CNVs, gene expression, and survival outcomes in pediatric, adolescent, and young adult OS patients. Chromosome 8 was a hotspot for poor prognostic signatures. The MYC-RAD21 copy number gain (8q24) correlated with increased gene expression and poor overall survival in 90% of the patients (n = 85). MYC and RAD21 play a role in replication-stress, which is a therapeutically actionable network. We prioritized replication-stress regulators, bromodomain and extra-terminal proteins (BETs), and CHK1, in order to test the hypothesis that the inhibition of BET + CHK1 in MYC-RAD21+ pediatric OS models would be efficacious and safe. We demonstrate that MYC-RAD21+ pediatric OS cell lines were sensitive to the inhibition of BET (BETi) and CHK1 (CHK1i) at clinically achievable concentrations. While the potentiation of CHK1i-mediated effects by BETi was BET-BRD4-dependent, MYC expression was BET-BRD4-independent. In MYC-RAD21+ pediatric OS xenografts, BETi + CHK1i significantly decreased tumor growth, increased survival, and was well tolerated. Therefore, targeting replication stress is a promising strategy to pursue as a therapeutic option for this devastating disease.

Identification of potential prognostic small nucleolar RNA biomarkers for predicting overall survival in patients with sarcoma

Objective

The main purpose of the present study is to screen prognostic small nucleolar RNA (snoRNA) markers using the RNA‐sequencing (RNA‐seq) dataset of The Cancer Genome Atlas (TCGA) sarcoma cohort.

Methods

The sarcoma RNA‐seq dataset comes from the TCGA cohort. A total of 257 sarcoma patients were included into the prognostic analysis. Multiple bioinformatics analysis methods for functional annotation of snoRNAs and screening of targeted drugs, including biological network gene ontology tool, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA), and connectivity map (CMap) are used.

Results

We had identified 15 snoRNAs that were significantly related to the prognosis of sarcoma and constructed a prognostic signature based on four prognostic snoRNA (U3, SNORA73B, SNORD46, and SNORA26) expression values. Functional annotation of these four snoRNAs by their co‐expression genes suggests that some of them were closely related to cell cycle‐related biological processes and tumor‐related signaling pathways, such as Wnt, mitogen‐activated protein kinase, target of rapamycin, and nuclear factor‐kappa B signaling pathway. GSEA of the risk score suggests that high risk score phenotype was significantly enriched in cell cycle‐related biological processes, protein SUMOylation, DNA replication, p53 binding, regulation of DNA repair, and DNA methylation, as well as Myc, Wnt, RB1, E2F, and TEL pathways. Then we also used the CMap online tool to screen five targeted drugs (rilmenidine, pizotifen, amiprilose, quipazine, and cinchonidine) for this risk score model in sarcoma.

Conclusion

Our study have identified 15 snoRNAs that may be serve as novel prognostic biomarkers for sarcoma, and constructed a prognostic signature based on four prognostic snoRNA expression values.

Sequential Whole Exome Sequencing Reveals Somatic Mutations Associated with Platinum Response in NSCLC

Purpose

Resistance is one of the main limitations of successful platinum treatment in non-small-cell lung cancer (NSCLC) patients. In this study, we aimed to identify somatic mutations associated with platinum response.

Patients and Methods

A total of 57 patients who received platinum-based chemotherapy only and 13 patients who received neoadjuvant chemotherapy (NAC) were enrolled. Somatic mutations were obtained from targeted and whole exome sequencing (WES).

Results

Somatic mutations in a total of 225 genes were observed. Nonsynonymous variants in EGFR, TTN, TP53 and KRAS, and copy number variations (SCNVs) in chromosome 8q24.3 and 22q11.21 were identified to be associated with platinum response. Based on these mutations, the mutational signature associated with the failure of DNA double-strand break and calcium signaling pathways were identified to be associated with platinum response. Besides, we observed a decrease in tumor mutational burden after chemotherapy. We also evaluated the mutation spectrum consistency between cell-free DNA (cfDNA) and tissue DNA. Somatic mutations detected in cfDNA were consistent with that in tDNA, which indicated that plasma might be used for somatic mutation detection.

Conclusion

These results support that somatic mutations can affect platinum drug response and provide potential clinical biomarkers for NSCLC treatment.

The Osteosarcoma Microenvironment: A Complex But Targetable Ecosystem

Osteosarcomas are the most frequent primary bone sarcomas, affecting mainly children, adolescents, and young adults, and with a second peak of incidence in elderly individuals. The current therapeutic management, a combined regimen of poly-chemotherapy and surgery, still remains largely insufficient, as patient survival has not improved in recent decades. Osteosarcomas are very heterogeneous tumors, both at the intra- and inter-tumor level, with no identified driver mutation. Consequently, efforts to improve treatments using targeted therapies have faced this lack of specific osteosarcoma targets. Nevertheless, these tumors are inextricably linked to their local microenvironment, composed of bone, stromal, vascular and immune cells and the osteosarcoma microenvironment is now considered to be essential and supportive for growth and dissemination. This review describes the different actors of the osteosarcoma microenvironment and gives an overview of the past, current, and future strategies of therapy targeting this complex ecosystem, with a focus on the role of extracellular vesicles and on the emergence of multi-kinase inhibitors.

Genomics and Therapeutic Vulnerabilities of Primary Bone Tumors

Osteosarcoma, Ewing sarcoma and chondrosarcoma are rare diseases but the most common primary tumors of bone. The genes directly involved in the sarcomagenesis, tumor progression and treatment responsiveness are not completely defined for these tumors, and the powerful discovery of genetic analysis is highly warranted in the view of improving the therapy and cure of patients. The review summarizes recent advances concerning the molecular and genetic background of these three neoplasms and, of their most common variants, highlights the putative therapeutic targets and the clinical trials that are presently active, and notes the fundamental issues that remain unanswered. In the era of personalized medicine, the rarity of sarcomas may not be the major obstacle, provided that each patient is studied extensively according to a road map that combines emerging genomic and functional approaches toward the selection of novel therapeutic strategies.

Bioinformatics analysis of the expression of HOXC13 and its role in the prognosis of breast cancer

The homeobox (HOX) genes, a class of transcription factors, are known to promote embryonic development and induce tumor formation. To date, the HOXA and HOXB gene families have been reported to be associated with breast cancer. However, the expression and exact role of homeobox C13 (HOXC13) in breast cancer has not yet been investigated. In the present study, the HOXC13 expression in human breast cancer was evaluated using the Oncomine database and Cancer Cell Line Encyclopedia (CCLE). Next, the Gene expression-based Outcome for Breast cancer online database, cBioportal, University of California Santa Cruz Xena browser and bc-GenExMinerv were used to explore the specific expression of HOXC13 in breast cancer. The methylation and mutation status of HOXC13 in breast cancer was then validated using the CCLE and cBioportal databases. Finally, the co-expression of HOX transcript antisense RNA (HOTAIR) and HOXC13 in breast cancer were analyzed and their impact on clinical prognosis determined. It was found that the expression of HOXC13 was high in breast cancer compared with other types of cancer, such as gastric cancer and colon cancer. Following co-expression analysis, a significant positive association was identified between HOTAIR and HOXC13. An association between HOTAIR and HOXC13, and lymph node and distant metastasis recurrence was also revealed during the development of breast cancer. Of note, survival analysis showed that high expression of HOTAIR and HOXC13 predicted poor prognosis. These findings revealed that HOXC13 plays an important role in the progression of breast cancer. However, the specific mechanism needs to be confirmed by subsequent experiments.

A Shallow Convolutional Learning Network for Classification of Cancers Based on Copy Number Variations

Genomic copy number variations (CNVs) are among the most important structural variations. They are linked to several diseases and cancer types. Cancer is a leading cause of death worldwide. Several studies were conducted to investigate the causes of cancer and its association with genomic changes to enhance its management and improve the treatment opportunities. Classification of cancer types based on the CNVs falls in this category of research. We reviewed the recent, most successful methods that used machine learning algorithms to solve this problem and obtained a dataset that was tested by some of these methods for evaluation and comparison purposes. We propose three deep learning techniques to classify cancer types based on CNVs: a six-layer convolutional net (CNN6), residual six-layer convolutional net (ResCNN6), and transfer learning of pretrained VGG16 net. The results of the experiments performed on the data of six cancer types demonstrated a high accuracy of 86% for ResCNN6 followed by 85% for CNN6 and 77% for VGG16. The results revealed a lower prediction accuracy for one of the classes (uterine corpus endometrial carcinoma (UCEC)). Repeating the experiments after excluding this class reveals improvements in the accuracies: 91% for CNN6 and 92% for Res CNN6. We observed that UCEC and ovarian serous carcinoma (OV) share a considerable subset of their features, which causes a struggle for learning in the classifiers. We repeated the experiment again by balancing the six classes through oversampling of the training dataset and the result was an enhancement in both overall and UCEC classification accuracies.

An update on emerging drugs in osteosarcoma: towards tailored therapies?

Introduction: Current treatment of conventional and non-conventional high-grade osteosarcoma (HGOS) is based on the surgical removal of primary tumor and, when possible, of metastases and local reccurrence, together with systemic pre- and post-operative chemotherapy with drugs that have been used since decades. Areas covered: This review is intended to summarize the new agents and therapeutic strategies that are under clinical evaluation in HGOS, with the aim to increase the cure probability of this highly malignant bone tumor, which has not significantly improved during the last 30-40 years. The list of drugs, compounds and treatment modalities presented and discussed here has been generated by considering only those that are included in presently ongoing and recruiting clinical trials, or which have been completed in the last 2 years with reported results, on the basis of the information obtained from different and continuously updated databases. Expert opinion: Despite HGOS is a rare tumor, several clinical trials are presently evaluating different treatment strategies, which may hopefully positively impact on the outcome of patients who experience unfavorable prognosis when treated with conventional therapies.

Identifying subpathway signatures for individualized anticancer drug response by integrating multi-omics data

Background

Individualized drug response prediction is vital for achieving personalized treatment of cancer and moving precision medicine forward. Large-scale multi-omics profiles provide unprecedented opportunities for precision cancer therapy.

Methods

In this study, we propose a pipeline to identify subpathway signatures for anticancer drug response of individuals by integrating the comprehensive contributions of multiple genetic and epigenetic (gene expression, copy number variation and DNA methylation) alterations.

Results

Totally, 46 subpathway signatures associated with individual responses to different anticancer drugs were identified based on five cancer-drug response datasets. We have validated the reliability of subpathway signatures in two independent datasets. Furthermore, we also demonstrated these multi-omics subpathway signatures could significantly improve the performance of anticancer drug response prediction. In-depth analysis of these 46 subpathway signatures uncovered the essential roles of three omics types and the functional associations underlying different anticancer drug responses. Patient stratification based on subpathway signatures involved in anticancer drug response identified subtypes with different clinical outcomes, implying their potential roles as prognostic biomarkers. In addition, a landscape of subpathways associated with cellular responses to 191 anticancer drugs from CellMiner was provided and the mechanism similarity of drug action was accurately unclosed based on these subpathways. Finally, we constructed a user-friendly web interface-CancerDAP (http://bio-bigdata.hrbmu.edu.cn/CancerDAP/) available to explore 2751 subpathways relevant with 191 anticancer drugs response.

Conclusions

Taken together, our study identified and systematically characterized subpathway signatures for individualized anticancer drug response prediction, which may promote the precise treatment of cancer and the study for molecular mechanisms of drug actions.

Electronic supplementary material

The online version of this article (10.1186/s12967-019-2010-4) contains supplementary material, which is available to authorized users.

The International Conference on Intelligent Biology and Medicine (ICIBM) 2018: genomics meets medicine

During June 10–12, 2018, the International Conference on Intelligent Biology and Medicine (ICIBM 2018) was held in Los Angeles, California, USA. The conference included 11 scientific sessions, four tutorials, one poster session, four keynote talks and four eminent scholar talks that covered a wide range of topics ranging from 3D genome structure analysis and visualization, next generation sequencing analysis, computational drug discovery, medical informatics, cancer genomics to systems biology. While medical genomics has always been a main theme in ICIBM, this year we for the first time organized the BMC Medical Genomics Supplement for ICIBM. Here, we describe 15 ICIBM papers selected for publishing in BMC Medical Genomics.