Integrative Bayesian Analysis Identifies Rhabdomyosarcoma Disease Genes

Developmental Heterogeneity of Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric embryonal solid tumor and the most common pediatric soft tissue sarcoma. The histology and transcriptome of RMS resemble skeletal muscle progenitor cells that have failed to terminally differentiate. Thus, RMS is typically thought to arise from corrupted skeletal muscle progenitor cells during development. However, RMS can occur in body regions devoid of skeletal muscle, suggesting the potential for nonmyogenic cells of origin. Here, we discuss the interplay between RMS driver mutations and cell(s) of origin with an emphasis on driving location specificity. Additionally, we discuss the mechanisms governing RMS transformation events and tumor heterogeneity through the lens of transcriptional networks and epigenetic control. Finally, we reimagine Waddington's developmental landscape to include a plane of transformation connecting distinct lineage landscapes to more accurately reflect the phenomena observed in pediatric cancers.

The MYC-YBX1 Circuit in Maintaining Stem-like Vincristine-Resistant Cells in Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a pediatric soft tissue sarcoma that causes significant devastation, with no effective therapy for relapsed disease. The mechanisms behind treatment failures are poorly understood. Our study showed that treatment of RMS cells with vincristine led to an increase in CD133-positive stem-like resistant cells. Single cell RNAseq analysis revealed that MYC and YBX1 were among the top-scoring transcription factors in CD133-high expressing cells. Targeting MYC and YBX1 using CRISPR/Cas9 reduced stem-like characteristics and viability of the vincristine-resistant cells. MYC and YBX1 showed mutual regulation, with MYC binding to the YBX1 promoter and YBX1 binding to MYC mRNA. The MYC inhibitor MYC361i synergized with vincristine to reduce tumor growth and stem-like cells in a zebrafish model of RMS. MYC and YBX expression showed a positive correlation in RMS patients, and high MYC expression correlated with poor survival. Targeting the MYC-YBX1 axis holds promise for improving survival in RMS patients.

TWIST2-mediated chromatin remodeling promotes fusion-negative rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a common soft tissue sarcoma in children that resembles developing skeletal muscle. Unlike normal muscle cells, RMS cells fail to differentiate despite expression of the myogenic determination protein MYOD. The TWIST2 transcription factor is frequently overexpressed in fusion-negative RMS (FN-RMS). TWIST2 blocks differentiation by inhibiting MYOD activity in myoblasts, but its role in FN-RMS pathogenesis is incompletely understood. Here, we show that knockdown of TWIST2 enables FN-RMS cells to exit the cell cycle and undergo terminal myogenesis. TWIST2 knockdown also substantially reduces tumor growth in a mouse xenograft model of FN-RMS. Mechanistically, TWIST2 controls H3K27 acetylation at distal enhancers by interacting with the chromatin remodelers SMARCA4 and CHD3 to activate growth-related target genes and repress myogenesis-related target genes. These findings provide insights into the role of TWIST2 in maintaining an undifferentiated and tumorigenic state of FN-RMS and highlight the potential of suppressing TWIST2-regulated pathways to treat FN-RMS.

Identify potential driver genes for PAX-FOXO1 fusion-negative rhabdomyosarcoma through frequent gene co-expression network mining

Background

Rhabdomyosarcoma (RMS) is a soft tissue sarcoma usually originated from skeletal muscle. Currently, RMS classification based on PAX-FOXO1 fusion is widely adopted. However, compared to relatively clear understanding of the tumorigenesis in the fusion-positive RMS, little is known for that in fusion-negative RMS (FN-RMS).

Methods

We explored the molecular mechanisms and the driver genes of FN-RMS through frequent gene co-expression network mining (fGCN), differential copy number (CN) and differential expression analyses on multiple RMS transcriptomic datasets.

Results

We obtained 50 fGCN modules, among which five are differentially expressed between different fusion status. A closer look showed 23% of Module 2 genes are concentrated on several cytobands of chromosome 8. Upstream regulators such as MYC, YAP1, TWIST1 were identified for the fGCN modules. Using in a separate dataset we confirmed that, comparing to FP-RMS, 59 Module 2 genes show consistent CN amplification and mRNA overexpression, among which 28 are on the identified chr8 cytobands. Such CN amplification and nearby MYC (also resides on one of the above cytobands) and other upstream regulators (YAP1, TWIST1) may work together to drive FN-RMS tumorigenesis and progression. Up to 43.1% downstream targets of Yap1 and 45.8% of the targets of Myc are differentially expressed in FN-RMS vs. normal comparisons, which also confirmed the driving force of these regulators.

Discussion

We discovered that copy number amplification of specific cytobands on chr8 and the upstream regulators MYC, YAP1 and TWIST1 work together to affect the downstream gene co-expression and promote FN-RMS tumorigenesis and progression. Our findings provide new insights for FN-RMS tumorigenesis and offer promising targets for precision therapy. Experimental investigation about the functions of identified potential drivers in FN-RMS are in progress.

Nuclear PTEN's Functions in Suppressing Tumorigenesis: Implications for Rare Cancers

Phosphatase and tensin homolog (PTEN) encodes a tumor-suppressive phosphatase with both lipid and protein phosphatase activity. The tumor-suppressive functions of PTEN are lost through a variety of mechanisms across a wide spectrum of human malignancies, including several rare cancers that affect pediatric and adult populations. Originally discovered and characterized as a negative regulator of the cytoplasmic, pro-oncogenic phosphoinositide-3-kinase (PI3K) pathway, PTEN is also localized to the nucleus where it can exert tumor-suppressive functions in a PI3K pathway-independent manner. Cancers can usurp the tumor-suppressive functions of PTEN to promote oncogenesis by disrupting homeostatic subcellular PTEN localization. The objective of this review is to describe the changes seen in PTEN subcellular localization during tumorigenesis, how PTEN enters the nucleus, and the spectrum of impacts and consequences arising from disrupted PTEN nuclear localization on tumor promotion. This review will highlight the immediate need in understanding not only the cytoplasmic but also the nuclear functions of PTEN to gain more complete insights into how important PTEN is in preventing human cancers.

Methyltransferase Inhibition Enables Tgf Driven Induction of CDKN2A and B in Cancer Cells

CDKN2A/B deletion or silencing is common across human cancer, reinforcing the general importance of bypassing its tumor suppression in cancer formation or progression. In rhabdomyosarcoma (RMS) and neuroblastoma, two common childhood cancers, the three CDKN2A/B transcripts are independently expressed to varying degrees, but one, ARF, is uniformly silenced. Although TGFβ induces certain CDKN2A/B transcripts in HeLa cells, it was unable to do so in five tested RMS lines unless the cells were pretreated with a broadly acting methyltransferase inhibitor, DZNep, or one targeting EZH2. CDKN2A/B induction by TGFβ correlated with de novo appearance of three H3K27Ac peaks within a 20 kb cis element ∼150 kb proximal to CDKN2A/B. Deleting that segment prevented their induction by TGFβ but not a basal increase driven by methyltransferase inhibition alone. Expression of two CDKN2A/B transcripts was enhanced by dCas9/CRISPR activation targeting either the relevant promoter or the 20 kb cis elements, and this "precise" manipulation diminished RMS cell propagation in vitro. Our findings show crosstalk between methyltransferase inhibition and TGFβ-dependent activation of a remote enhancer to reverse CDKN2A/B silencing. Though focused on CDKN2A/B here, such crosstalk may apply to other TGFβ-responsive genes and perhaps govern this signaling protein's complex effects promoting or blocking cancer.

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.

Non-coding RNA in rhabdomyosarcoma progression and metastasis

Rhabdomyosarcoma (RMS) is a soft tissue sarcoma of skeletal muscle differentiation, with a predominant occurrence in children and adolescents. One of the major challenges facing treatment success is the presence of metastatic disease at the time of diagnosis, commonly associated with the more aggressive fusion-positive subtype. Non-coding RNA (ncRNA) can regulate gene transcription and translation, and their dysregulation has been associated with cancer development and progression. MicroRNA (miRNA) are short non-coding nucleic acid sequences involved in the regulation of gene expression that act by targeting messenger RNA (mRNA), and their aberrant expression has been associated with both RMS initiation and progression. Other ncRNA including long non-coding RNA (lncRNA), circular RNA (circRNA) and ribosomal RNA (rRNA) have also been associated with RMS revealing important mechanistic roles in RMS biology, but these studies are still limited and require further investigation. In this review, we discuss the established roles of ncRNA in RMS differentiation, growth and progression, highlighting their potential use in RMS prognosis, as therapeutic agents or as targets of treatment.

Synthetic essentiality between PTEN and core dependency factor PAX7 dictates rhabdomyosarcoma identity

PTEN promoter hypermethylation is nearly universal and PTEN copy number loss occurs in ~25% of fusion-negative rhabdomyosarcoma (FN-RMS). Here we show Pten deletion in a mouse model of FN-RMS results in less differentiated tumors more closely resembling human embryonal RMS. PTEN loss activated the PI3K pathway but did not increase mTOR activity. In wild-type tumors, PTEN was expressed in the nucleus suggesting loss of nuclear PTEN functions could account for these phenotypes. Pten deleted tumors had increased expression of transcription factors important in neural and skeletal muscle development including Dbx1 and Pax7. Pax7 deletion completely rescued the effects of Pten loss. Strikingly, these Pten;Pax7 deleted tumors were no longer FN-RMS but displayed smooth muscle differentiation similar to leiomyosarcoma. These data highlight how Pten loss in FN-RMS is connected to a PAX7 lineage-specific transcriptional output that creates a dependency or synthetic essentiality on the transcription factor PAX7 to maintain tumor identity.

DCX and CRABP2 are candidate genes for differential diagnosis between pre-chemotherapy embryonic and alveolar rhabdomyosarcoma in pediatric patients

IMPORTANCE

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in children. More than 90% of cases are classified as embryonic RMS (ERMS) or alveolar RMS (ARMS). ERMS has a worse prognosis than ARMS. Early differential diagnosis is of paramount importance for optimization of treatment.

OBJECTIVE

To identify genes that are differentially expressed between ARMS and ERMS, which can be used for accurate rhabdomyosarcoma classification.

METHODS

Three Gene Expression Omnibus datasets composed of ARMS and ERMS samples were screened and 35 differentially expressed genes (DEGs) were identified. Receiver operating characteristic curve analysis and area under the curve analysis was performed for these 35 DEGs and seven candidate genes with the best differential expression scores between ARMS and ERMS were determined. The expression of these seven candidate genes was validated by immunohistochemical analysis of pre-chemotherapy ARMS and ERMS specimens.

RESULTS

The levels of DCX and CRABP2 were confirmed to be remarkably different between paraffin-embedded ARMS and ERMS tissues, while EGFR abundance was only marginally different between these two RMS subtypes.

INTERPRETATION

DCX and CRABP2 are potential biomarkers for distinguishing ARMS from ERMS in pre-chemotherapy pediatric patients.

TPQCI: A topology potential-based method to quantify functional influence of copy number variations

Copy number variation (CNV) is a major type of chromosomal structural variation that play important roles in many diseases including cancers. Due to genome instability, a large number of CNV events can be detected in diseases such as cancer. Therefore, it is important to identify the functionally important CNVs in diseases, which currently still poses a challenge in genomics. One of the critical steps to solve the problem is to define the influence of CNV. In this paper, we provide a topology potential based method, TPQCI, to quantify this kind of influence by integrating statistics, gene regulatory associations, and biological function information. We used this metric to detect functionally enriched genes on genomic segments with CNV in breast cancer and multiple myeloma and discovered biological functions influenced by CNV. Our results demonstrate that, by using our proposed TPQCI metric, we can detect disease-specific genes that are influenced by CNVs. Source codes of TPQCI are provided in Github (https://github.com/usos/TPQCI).

Epigenetic regulator BMI1 promotes alveolar rhabdomyosarcoma proliferation and constitutes a novel therapeutic target

Rhabdomyosarcoma (RMS) is an aggressive pediatric soft tissue sarcoma. There are two main subtypes of RMS, alveolar rhabdomyosarcoma (ARMS) and embryonal rhabdomyosarcoma. ARMS typically encompasses fusion‐positive rhabdomyosarcoma, which expresses either PAX3‐FOXO1 or PAX7‐FOXO1 fusion proteins. There are no targeted therapies for ARMS; however, recent studies have begun to illustrate the cooperation between epigenetic proteins and the PAX3‐FOXO1 fusion, indicating that epigenetic proteins may serve as targets in ARMS. Here, we investigate the contribution of BMI1, given the established role of this epigenetic regulator in sustaining aggression in cancer. We determined that BMI1 is expressed across ARMS tumors, patient‐derived xenografts, and cell lines. We depleted BMI1 using RNAi and inhibitors (PTC‐209 and PTC‐028) and found that this leads to a decrease in cell growth/increase in apoptosis in vitro, and delays tumor growth in vivo. Our data suggest that BMI1 inhibition activates the Hippo pathway via phosphorylation of LATS1/2 and subsequent reduction in YAP levels and YAP/TAZ target genes. These results identify BMI1 as a potential therapeutic vulnerability in ARMS and warrant further investigation of BMI1 in ARMS and other sarcomas.

Applications of CRISPR/Cas9 in the research of malignant musculoskeletal tumors

Background

Malignant tumors of the musculoskeletal system, especially osteosarcoma, Ewing sarcoma and rhabdomyosarcoma, pose a major threat to the lives and health of adolescents and children. Current treatments for musculoskeletal tumors mainly include surgery, chemotherapy, and radiotherapy. The problems of chemotherapy resistance, poor long-term outcome of radiotherapy, and the inherent toxicity and side effects of chemical drugs make it extremely urgent to seek new treatment strategies.

Main text

As a potent gene editing tool, the rapid development of CRISPR/Cas9 technology in recent years has prompted scientists to apply it to the study of musculoskeletal tumors. This review summarizes the application of CRISPR/Cas9 technology for the treatment of malignant musculoskeletal tumors, focusing on its essential role in the field of basic research.

Conclusion

CRISPR, has demonstrated strong efficacy in targeting tumor-related genes, and its future application in the clinical treatment of musculoskeletal tumors is promising.

HDAC6 promotes growth, migration/invasion, and self-renewal of rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a devastating pediatric sarcoma. The survival outcomes remain poor for patients with relapsed or metastatic disease. Effective targeted therapy is lacking due to our limited knowledge of the underlying cellular and molecular mechanisms leading to disease progression. In this study, we used functional assays in vitro and in vivo (zebrafish and xenograft mouse models) to demonstrate the crucial role of HDAC6, a cytoplasmic histone deacetylase, in driving RMS tumor growth, self-renewal, and migration/invasion. Treatment with HDAC6-selective inhibitors recapitulates the HDAC6 loss-of-function phenotypes. HDAC6 regulates cytoskeletal dynamics to promote tumor cell migration and invasion. RAC1, a Rho family GTPase, is an essential mediator of HDAC6 function, and is necessary and sufficient for RMS cell migration and invasion. High expression of RAC1 correlates with poor clinical prognosis in RMS patients. Targeting the HDAC6-RAC1 axis represents a promising therapeutic option for improving survival outcomes of RMS patients.

SNAIL Promotes Metastatic Behavior of Rhabdomyosarcoma by Increasing EZRIN and AKT Expression and Regulating MicroRNA Networks

Rhabdomyosarcoma (RMS) is a predominant soft tissue tumor in children and adolescents. For high-grade RMS with metastatic involvement, the 3-year overall survival rate is only 25 to 30%. Thus, understanding the regulatory mechanisms involved in promoting the metastasis of RMS is important. Here, we demonstrate for the first time that the SNAIL transcription factor regulates the metastatic behavior of RMS both in vitro and in vivo. SNAIL upregulates the protein expression of EZRIN and AKT, known to promote metastatic behavior, by direct interaction with their promoters. Our data suggest that SNAIL promotes RMS cell motility, invasion and chemotaxis towards the prometastatic factors: HGF and SDF-1 by regulating RHO, AKT and GSK3β activity. In addition, miRNA transcriptome analysis revealed that SNAIL-miRNA axis regulates processes associated with actin cytoskeleton reorganization. Our data show a novel role of SNAIL in regulating RMS cell metastasis that may also be important in other mesenchymal tumor types and clearly suggests SNAIL as a promising new target for future RMS therapies.

Bayesian Modeling Identifies PLAG1 as a Key Regulator of Proliferation and Survival in Rhabdomyosarcoma Cells

We recently developed a novel computational algorithm that incorporates Bayesian methodology to identify rhabdomyosarcoma disease genes whose expression level correlates with copy-number variations, and we identified PLAG1 as a candidate oncogenic driver. Although PLAG1 has been shown to contribute to other type of cancers, its role in rhabdomyosarcoma has not been elucidated. We observed that PLAG1 mRNA is highly expressed in rhabdomyosarcoma and is associated with PLAG1 gene copy-number gain. Knockdown of PLAG1 dramatically decreased cell accumulation and induced apoptosis in rhabdomyosarcoma cells, whereas its ectopic expression increased cell accumulation in vitro and as a xenograft and promoted G₁ to S-phase cell-cycle progression. We found that PLAG1 regulates IGF2 expression and influences AKT and MAPK pathways in rhabdomyosarcoma, and IGF2 partially rescues cell death triggered by PLAG1 knockdown. The expression level of PLAG1 correlated with the IC₅₀ of rhabdomyosarcoma cells to BMS754807, an IGF receptor inhibitor. IMPLICATIONS: Our data demonstrate that PLAG1 contributes to proliferation and survival of rhabdomyosarcoma cells at least partially by inducing IGF2, and this new understanding may have the potential for clinical translation.

Insights into pediatric rhabdomyosarcoma research: Challenges and goals

Overall survival rates for pediatric patients with high-risk or relapsed rhabdomyosarcoma (RMS) have not improved significantly since the 1980s. Recent studies have identified a number of targetable vulnerabilities in RMS, but these discoveries have infrequently translated into clinical trials. We propose streamlining the process by which agents are selected for clinical evaluation in RMS. We believe that strong consideration should be given to the development of combination therapies that add biologically targeted agents to conventional cytotoxic drugs. One example of this type of combination is the addition of the WEE1 inhibitor AZD1775 to the conventional cytotoxic chemotherapeutics, vincristine and irinotecan.

Twist2 amplification in rhabdomyosarcoma represses myogenesis and promotes oncogenesis by redirecting MyoD DNA binding

Li et al. show that TWIST2 shapes the epigenetic landscape to drive chromatin opening at oncogenic loci and chromatin closing at myogenic loci. These epigenetic changes redirect MyoD binding from myogenic genes toward oncogenic, metabolic, and growth genes.

Rhabdomyosarcoma (RMS) is an aggressive pediatric cancer composed of myoblast-like cells. Recently, we discovered a unique muscle progenitor marked by the expression of the Twist2 transcription factor. Genomic analyses of 258 RMS patient tumors uncovered prevalent copy number amplification events and increased expression of TWIST2 in fusion-negative RMS. Knockdown of TWIST2 in RMS cells results in up-regulation of MYOGENIN and a decrease in proliferation, implicating TWIST2 as an oncogene in RMS. Through an inducible Twist2 expression system, we identified Twist2 as a reversible inhibitor of myogenic differentiation with the remarkable ability to promote myotube dedifferentiation in vitro. Integrated analysis of genome-wide ChIP-seq and RNA-seq data revealed the first dynamic chromatin and transcriptional landscape of Twist2 binding during myogenic differentiation. During differentiation, Twist2 competes with MyoD at shared DNA motifs to direct global gene transcription and repression of the myogenic program. Additionally, Twist2 shapes the epigenetic landscape to drive chromatin opening at oncogenic loci and chromatin closing at myogenic loci. These epigenetic changes redirect MyoD binding from myogenic genes toward oncogenic, metabolic, and growth genes. Our study reveals the dynamic interplay between two opposing transcriptional regulators that control the fate of RMS and provides insight into the molecular etiology of this aggressive form of cancer.

Rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in children and represents a high-grade neoplasm of skeletal myoblast-like cells. Decades of clinical and basic research have gradually improved our understanding of the pathophysiology of RMS and helped to optimize clinical care. The two major subtypes of RMS, originally characterized on the basis of light microscopic features, are driven by fundamentally different molecular mechanisms and pose distinct clinical challenges. Curative therapy depends on control of the primary tumour, which can arise at many distinct anatomical sites, as well as controlling disseminated disease that is known or assumed to be present in every case. Sophisticated risk stratification for children with RMS incorporates various clinical, pathological and molecular features, and that information is used to guide the application of multifaceted therapy. Such therapy has historically included cytotoxic chemotherapy as well as surgery, ionizing radiation or both. This Primer describes our current understanding of RMS epidemiology, disease susceptibility factors, disease mechanisms and elements of clinical care, including diagnostics, risk-based care of newly diagnosed and relapsed disease and the prevention and management of late effects in survivors. We also outline potential opportunities to further translate new biological insights into improved clinical outcomes.