Protein phosphatase 1 regulatory subunit 1A in ewing sarcoma tumorigenesis and metastasis

Single-cell RNA sequencing analysis of the retina under acute high intraocular pressure

JOURNAL/nrgr/04.03/01300535-202419110-00032/figure1/v/2024-03-08T184507Z/r/image-tiff High intraocular pressure causes retinal ganglion cell injury in primary and secondary glaucoma diseases, yet the molecular landscape characteristics of retinal cells under high intraocular pressure remain unknown. Rat models of acute hypertension ocular pressure were established by injection of cross-linked hyaluronic acid hydrogel (Healaflow®). Single-cell RNA sequencing was then used to describe the cellular composition and molecular profile of the retina following high intraocular pressure. Our results identified a total of 12 cell types, namely retinal pigment epithelial cells, rod-photoreceptor cells, bipolar cells, Müller cells, microglia, cone-photoreceptor cells, retinal ganglion cells, endothelial cells, retinal progenitor cells, oligodendrocytes, pericytes, and fibroblasts. The single-cell RNA sequencing analysis of the retina under acute high intraocular pressure revealed obvious changes in the proportions of various retinal cells, with ganglion cells decreased by 23%. Hematoxylin and eosin staining and TUNEL staining confirmed the damage to retinal ganglion cells under high intraocular pressure. We extracted data from retinal ganglion cells and analyzed the retinal ganglion cell cluster with the most distinct expression. We found upregulation of the B3gat2 gene, which is associated with neuronal migration and adhesion, and downregulation of the Tsc22d gene, which participates in inhibition of inflammation. This study is the first to reveal molecular changes and intercellular interactions in the retina under high intraocular pressure. These data contribute to understanding of the molecular mechanism of retinal injury induced by high intraocular pressure and will benefit the development of novel therapies.

Combinatorial macrophage induced innate immunotherapy against Ewing sarcoma: Turning "Two Keys" simultaneously

BACKGROUND

Macrophages play important roles in phagocytosing tumor cells. However, tumors escape macrophage phagocytosis in part through the expression of anti-phagocytic signals, most commonly CD47. In Ewing sarcoma (ES), we found that tumor cells utilize dual mechanisms to evade macrophage clearance by simultaneously over-expressing CD47 and down-regulating cell surface calreticulin (csCRT), the pro-phagocytic signal. Here, we investigate the combination of a CD47 blockade (magrolimab, MAG) to inhibit the anti-phagocytic signal and a chemotherapy regimen (doxorubicin, DOX) to enhance the pro-phagocytic signal to induce macrophage phagocytosis of ES cells in vitro and inhibit tumor growth and metastasis in vivo.

METHODS

Macrophages were derived from human peripheral blood monocytes by granulocyte-macrophage colony-stimulating factor (GM-CSF) and macrophage colony-stimulating factor (M-CSF). Flow cytometry- and microscopy-based in-vitro phagocytosis assays were performed to evaluate macrophage phagocytosis of ES cells. Annexin-V assay was performed to evaluate apoptosis. CD47 was knocked out by CRISPR/Cas9 approach. ES cell-based and patient-derived-xenograft (PDX)-based mouse models were utilized to assess the effects of MAG and/or DOX on ES tumor development and animal survival. RNA-Seq combined with CIBERSORTx analysis was utilized to identify changes in tumor cell transcriptome and tumor infiltrating immune cell profiling in MAG and/or DOX treated xenograft tumors.

RESULTS

We found that MAG significantly increased macrophage phagocytosis of ES cells in vitro (p < 0.01) and had significant effect on reducing tumor burden (p < 0.01) and increasing survival in NSG mouse model (p < 0.001). The csCRT level on ES cells was significantly enhanced by DOX in a dose- and time-dependent manner (p < 0.01). Importantly, DOX combined with MAG significantly enhanced macrophage phagocytosis of ES cells in vitro (p < 0.01) and significantly decreased tumor burden (p < 0.01) and lung metastasis (p < 0.0001) and extended animal survival in vivo in two different mouse models of ES (p < 0.0001). Furthermore, we identified CD38, CD209, CD163 and CD206 as potential markers for ES-phagocytic macrophages. Moreover, we found increased M2 macrophage infiltration and decreased expression of Cd209 in the tumor microenvironment of MAG and DOX combinatorial therapy treated tumors.

CONCLUSIONS

By turning "two keys" simultaneously to reactivate macrophage phagocytic activity, our data demonstrated an effective and highly translatable alternative therapeutic approach utilizing innate (tumor associated macrophages) immunotherapy against high-risk metastatic ES.

Proteomic Characterization of Undifferentiated Small Round Cell Sarcomas With EWSR1 and CIC::DUX4 Translocations Reveals Diverging Tumor Biology and Distinct Diagnostic Markers

Undifferentiated small round cell sarcomas (USRS) of bone and soft tissue are a group of tumors with heterogenic genomic alterations sharing similar morphology. In the present study, we performed a comparative large-scale proteomic analysis of USRS (n = 42) with diverse genomic translocations including classic Ewing sarcomas with EWSR1::FLI1 fusions (n = 24) or EWSR1::ERG fusions (n = 4), sarcomas with an EWSR1 rearrangement (n = 2), CIC::DUX4 fusion (n = 8), as well as tumors classified as USRS with no genetic data available (n = 4). Proteins extracted from formalin-fixed, paraffin-embedded pretherapeutic biopsies were analyzed qualitatively and quantitatively using shotgun mass spectrometry (MS). More than 8000 protein groups could be quantified using data-independent acquisition. Unsupervised hierarchical cluster analysis based on proteomic data allowed stratification of the 42 cases into distinct groups reflecting the different molecular genotypes. Protein signatures that significantly correlated with the respective genomic translocations were identified and used to generate a heatmap of all 42 sarcomas with assignment of cases with unknown molecular genetic data to either the EWSR1- or CIC-rearranged groups. MS-based prediction of sarcoma subtypes was molecularly confirmed in 2 cases where next-generation sequencing was technically feasible. MS also detected proteins routinely used in the immunohistochemical approach for the differential diagnosis of USRS. BCL11B highly expressed in Ewing sarcomas, and BACH2 as well as ETS-1 highly expressed in CIC::DUX4-associated sarcomas, were among proteins identified by the present proteomic study, and were chosen for immunohistochemical confirmation of MS data in our study cohort. Differential expressions of these 3 markers in the 2 genetic groups were further validated in an independent cohort of n = 34 USRS. Finally, our proteomic results point toward diverging signaling pathways in the different USRS subgroups.

Role of Protein Phosphatases in Tumor Angiogenesis: Assessing PP1, PP2A, PP2B and PTPs Activity

Tumor angiogenesis, the formation of new blood vessels to support tumor growth and metastasis, is a complex process regulated by a multitude of signaling pathways. Dysregulation of signaling pathways involving protein kinases has been extensively studied, but the role of protein phosphatases in angiogenesis within the tumor microenvironment remains less explored. However, among angiogenic pathways, protein phosphatases play critical roles in modulating signaling cascades. This review provides a comprehensive overview of the involvement of protein phosphatases in tumor angiogenesis, highlighting their diverse functions and mechanisms of action. Protein phosphatases are key regulators of cellular signaling pathways by catalyzing the dephosphorylation of proteins, thereby modulating their activity and function. This review aims to assess the activity of the protein tyrosine phosphatases and serine/threonine phosphatases. These phosphatases exert their effects on angiogenic signaling pathways through various mechanisms, including direct dephosphorylation of angiogenic receptors and downstream signaling molecules. Moreover, protein phosphatases also crosstalk with other signaling pathways involved in angiogenesis, further emphasizing their significance in regulating tumor vascularization, including endothelial cell survival, sprouting, and vessel maturation. In conclusion, this review underscores the pivotal role of protein phosphatases in tumor angiogenesis and accentuate their potential as therapeutic targets for anti-angiogenic therapy in cancer.

Cordycepin Triphosphate as a Potential Modulator of Cellular Plasticity in Cancer via cAMP-Dependent Pathways: An In Silico Approach

Cordycepin, or 3'-deoxyadenosine, is an adenosine analog with a broad spectrum of biological activity. The key structural difference between cordycepin and adenosine lies in the absence of a hydroxyl group at the 3' position of the ribose ring. Upon administration, cordycepin can undergo an enzymatic transformation in specific tissues, forming cordycepin triphosphate. In this study, we conducted a comprehensive analysis of the structural features of cordycepin and its derivatives, contrasting them with endogenous purine-based metabolites using chemoinformatics and bioinformatics tools in addition to molecular dynamics simulations. We tested the hypothesis that cordycepin triphosphate could bind to the active site of the adenylate cyclase enzyme. The outcomes of our molecular dynamics simulations revealed scores that are comparable to, and superior to, those of adenosine triphosphate (ATP), the endogenous ligand. This interaction could reduce the production of cyclic adenosine monophosphate (cAMP) by acting as a pseudo-ATP that lacks a hydroxyl group at the 3' position, essential to carry out nucleotide cyclization. We discuss the implications in the context of the plasticity of cancer and other cells within the tumor microenvironment, such as cancer-associated fibroblast, endothelial, and immune cells. This interaction could awaken antitumor immunity by preventing phenotypic changes in the immune cells driven by sustained cAMP signaling. The last could be an unreported molecular mechanism that helps to explain more details about cordycepin's mechanism of action.

cAMP-PKA/EPAC signaling and cancer: the interplay in tumor microenvironment

Cancer is a complex disease resulting from abnormal cell growth that is induced by a number of genetic and environmental factors. The tumor microenvironment (TME), which involves extracellular matrix, cancer-associated fibroblasts (CAF), tumor-infiltrating immune cells and angiogenesis, plays a critical role in tumor progression. Cyclic adenosine monophosphate (cAMP) is a second messenger that has pleiotropic effects on the TME. The downstream effectors of cAMP include cAMP-dependent protein kinase (PKA), exchange protein activated by cAMP (EPAC) and ion channels. While cAMP can activate PKA or EPAC and promote cancer cell growth, it can also inhibit cell proliferation and survival in context- and cancer type-dependent manner. Tumor-associated stromal cells, such as CAF and immune cells, can release cytokines and growth factors that either stimulate or inhibit cAMP production within the TME. Recent studies have shown that targeting cAMP signaling in the TME has therapeutic benefits in cancer. Small-molecule agents that inhibit adenylate cyclase and PKA have been shown to inhibit tumor growth. In addition, cAMP-elevating agents, such as forskolin, can not only induce cancer cell death, but also directly inhibit cell proliferation in some cancer types. In this review, we summarize current understanding of cAMP signaling in cancer biology and immunology and discuss the basis for its context-dependent dual role in oncogenesis. Understanding the precise mechanisms by which cAMP and the TME interact in cancer will be critical for the development of effective therapies. Future studies aimed at investigating the cAMP-cancer axis and its regulation in the TME may provide new insights into the underlying mechanisms of tumorigenesis and lead to the development of novel therapeutic strategies.

A humanized orthotopic mouse model for preclinical evaluation of immunotherapy in Ewing sarcoma

The advent of novel cancer immunotherapy approaches is revolutionizing the treatment for cancer. Current small animal models for most cancers are syngeneic or genetically engineered mouse models or xenograft models based on immunodeficient mouse strains. These models have been limited in evaluating immunotherapy regimens due to the lack of functional human immune system. Development of animal models for bone cancer faces another challenge in the accessibility of tumor engraftment sites. Here, we describe a protocol to develop an orthotopic humanized mouse model for a bone and soft tissue sarcoma, Ewing sarcoma, by transplanting fresh human cord blood CD34+ hematopoietic stem cells into young NSG-SGM3 mice combined with subsequent Ewing sarcoma patient derived cell engraftment in the tibia of the humanized mice. We demonstrated early and robust reconstitution of human CD45+ leukocytes including T cells, B cells, natural killer cells and monocytes. Ewing sarcoma xenograft tumors successfully orthotopically engrafted in the humanized mice with minimal invasive procedures. We validated the translational utility of this orthotopic humanized model by evaluating the safety and efficacy of an immunotherapy antibody, magrolimab. Treatment with magrolimab induces CD47 blockade resulting in significantly decreased primary tumor growth, decreased lung metastasis and prolonged animal survival in the established humanized model. Furthermore, the humanized model recapitulated the dose dependent toxicity associated with the CD47 blockade as observed in patients in clinical trials. In conclusion, this orthotopic humanized mouse model of Ewing sarcoma represents an improved platform for evaluating immunotherapy in bone and soft tissue sarcoma, such as Ewing sarcoma. With careful design and optimization, this model is generalizable for other bone malignancies.

Histopathology and levels of proteins in plasma associate with survival after colorectal cancer diagnosis

BACKGROUND

The TNM system is used to assess prognosis after colorectal cancer (CRC) diagnosis. Other prognostic factors reported include histopathological assessments of the tumour, tumour mutations and proteins in the blood. As some of these factors are strongly correlated, it is important to evaluate the independent effects they may have on survival.

METHODS

Tumour samples from 2162 CRC patients were visually assessed for amount of tumour stroma, severity of lymphocytic infiltrate at the tumour margins and the presence of lymphoid follicles. Somatic mutations in the tumour were assessed for 2134 individuals. Pre-surgical levels of 4963 plasma proteins were measured in 128 individuals. The associations between these features and prognosis were inspected by a Cox Proportional Hazards Model (CPH).

RESULTS

Levels of stroma, lymphocytic infiltration and presence of lymphoid follicles all associate with prognosis, along with high tumour mutation burden, high microsatellite instability and TP53 and BRAF mutations. The somatic mutations are correlated with the histopathology and none of the somatic mutations associate with survival in a multivariate analysis. Amount of stroma and lymphocytic infiltration associate with local invasion of tumours. Elevated levels of two plasma proteins, CA-125 and PPP1R1A, associate with a worse prognosis.

CONCLUSIONS

Tumour stroma and lymphocytic infiltration variables are strongly associated with prognosis of CRC and capture the prognostic effects of tumour mutation status. CA-125 and PPP1R1A may be useful prognostic biomarkers in CRC.

Prognostic value of GLCE and infiltrating immune cells in Ewing sarcoma

Background

The prognostic value of D-glucuronyl C5-epimerase (GLCE) and mast cell infiltration in Ewing sarcoma (ES) has not been well specified and highlighted, which may facilitate survival prediction and treatment.

Methods

Several qualified datasets were downloaded from the GEO website. Common differentially expressed genes between normal subjects and ES patients in GSE17679, GSE45544, and GSE68776 were identified and screened by multiple algorithms to find hub genes with prognostic value. The prognostic value of 64 infiltrating cells was also explored. A prognostic model was established and then validated with GSE63155 and GSE63156. Finally, functional analysis was performed.

Results

GLCE and mast cell infiltration were screened as two indicators for a prognostic model. The Kaplan‒Meier analysis showed that patients in the low GLCE expression, mast cell infiltration and risk score groups had poorer outcomes than patients in the high GLCE expression, mast cell infiltration and risk score groups, both in the training and validation sets. Scatter plots and heatmaps also indicated the same results. The concordance indices and calibration analyses indicated a high prediction accuracy of the model in the training and validation sets. The time-dependent receiver operating characteristic analyses suggested high sensitivity and specificity of the model, with area under the curve values between 0.76 and 0.98. The decision curve analyses suggested a significantly higher net benefit by the model than the treat-all and treat-none strategies. Functional analyses suggested that glycosaminoglycan biosynthesis-heparan sulfate/heparin, the cell cycle and microRNAs in cancer were upregulated in ES patients.

Conclusions

GLCE and mast cell infiltration are potential prognostic indicators in ES. GLCE may affect the proliferation, angiogenesis and metastasis of ES by affecting the biosynthesis of heparan sulfate and heparin.

Targeted Therapy for EWS-FLI1 in Ewing Sarcoma

Ewing sarcoma (EwS) is a rare and predominantly pediatric malignancy of bone and soft tissue in children and adolescents. Although international collaborations have greatly improved the prognosis of most EwS, the occurrence of macrometastases or relapse remains challenging. The prototypic oncogene EWS-FLI1 acts as an aberrant transcription factor that drives the cellular transformation of EwS. In addition to its involvement in RNA splicing and the DNA damage response, this chimeric protein directly binds to GGAA repeats, thereby modifying the transcriptional profile of EwS. Direct pharmacological targeting of EWS-FLI1 is difficult because of its intrinsically disordered structure. However, targeting the EWS-FLI1 protein complex or downstream pathways provides additional therapeutic options. This review describes the EWS-FLI1 protein partners and downstream pathways, as well as the related target therapies for the treatment of EwS.

Molecular and biologic biomarkers of Ewing sarcoma: A systematic review

With an annual incidence of less than 1%, Ewing sarcoma mainly occurs in children and young adults. It is not a frequent tumor but is the second most common bone malignancy in children. It has a 5-year survival rate of 65-75%; however, it has a poor prognosis when it relapses in patients. A genomic profile of this tumor can potentially help identify poor prognosis patients earlier and guide their treatment. A systematic review of the articles concerning genetic biomarkers in Ewing sarcoma was conducted using the Google Scholar, Cochrane, and PubMed database. There were 71 articles discovered. Numerous diagnostic, prognostic, and predictive biomarkers were found. However, more research is necessary to confirm the role of some of the mentioned biomarkers. .

Identification of Methylation Signatures and Rules for Sarcoma Subtypes by Machine Learning Methods

Sarcoma, the second common type of solid tumor in children and adolescents, has a wide variety of subtypes that are often not properly diagnosed at an early stage, leading to late metastases and causing serious loss of life and property to patients and families. It exhibits a high degree of heterogeneity at the cellular, molecular, and epigenetic levels, where DNA methylation has been proposed to play a role in the diagnosis of sarcoma subtypes. Thus, this study is aimed at finding potential biomarkers at the DNA methylation level to distinguish different sarcoma subtypes. A machine learning process was designed to analyse sarcoma samples, each of which was represented by lots of methylation sites. Irrelevant sites were removed using the Boruta method, and remaining sites related to the target variables were kept for further analyses. Afterward, three feature ranking methods (LASSO, LightGBM, and MCFS) were adopted to rank these features, and six classification models were constructed by combining incremental feature selection and two classification algorithms (decision tree and random forest). Among these models, the performance of RF model was higher than that of DT model under all three ranking conditions. The specific expression of genes obtained from the annotation of highly correlated methylation site features, such as PRKAR1B, INPP5A, and GLI3, was proven to be associated with sarcoma by publications. Moreover, the quantitative rules obtained by decision tree algorithm helped us to understand the essential differences between various sarcoma types and classify sarcoma subtypes, providing a new means of clinical identification and determining new therapeutic targets.

Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

Regulation of Metastasis in Ewing Sarcoma

Ewing sarcoma (EwS) is a type of bone and soft tissue tumor in children and adolescents. Over 85% of cases are caused by the expression of fusion protein EWSR1-FLI1 generated by chromosome translocation. Acting as a potent chimeric oncoprotein, EWSR1-FLI1 binds to chromatin, changes the epigenetic states, and thus alters the expression of a large set of genes. Several studies have revealed that the expression level of EWSR1-FLI1 is variable and dynamic within and across different EwS cell lines and primary tumors, leading to tumoral heterogeneity. Cells with high EWSR1-FLI1 expression (EWSR1-FLI1-high) proliferate in an exponential manner, whereas cells with low EWSR1-FLI1 expression (EWSR1-FLI1-low) tend to have a strong propensity to migrate, invade, and metastasize. Metastasis is the leading cause of cancer-related deaths. The continuous evolution of EwS research has revealed some of the molecular underpinnings of this dissemination process. In this review, we discuss the molecular signatures that contribute to metastasis.

Human EWS-FLI protein recapitulates in Drosophila the neomorphic functions that induce Ewing sarcoma tumorigenesis

Ewing sarcoma (EwS) is a human malignant tumor typically driven by the Ewing sarcoma-Friend leukemia integration (EWS-FLI) fusion protein. A paucity of genetically modified animal models, partially owed to the high toxicity of EWS-FLI, hinders research on EwS. Here, we report a spontaneous mutant variant, EWS-FLI_1FS, that circumvents the toxicity issue in Drosophila. Through proteomic and genomic analyses, we show that human EWS-FLI_1FS interacts with the Drosophila homologues of EWS-FLI human protein partners, including core subunits of chromatin remodeling complexes, the transcription machinery, and the spliceosome; brings about a massive dysregulation of transcription that affects a significant fraction of known targets of EWS-FLI in human cells; and modulates splicing. We also show that EWS-FLI_1FS performs in Drosophila the two major neomorphic activities that it is known to have in human cells: activation of transcription from GGAA microsatellites and out competition of ETS transcription factors. We conclude that EWS-FLI_1FS reproduces in Drosophila the known oncogenic activities of EWS-FLI that drive EwS tumorigenesis in humans. These results open up an unprecedented opportunity to investigate EWS-FLI's oncogenic pathways in vivo in a genetically tractable organism.

AT-hook DNA-binding motif-containing protein one knockdown downregulates EWS-FLI1 transcriptional activity in Ewing's sarcoma cells

Ewing's sarcoma is the second most common bone malignancy in children or young adults and is caused by an oncogenic transcription factor by a chromosomal translocation between the EWSR1 gene and the ETS transcription factor family. However, the transcriptional mechanism of EWS-ETS fusion proteins is still unclear. To identify the transcriptional complexes of EWS-ETS fusion transcription factors, we applied a proximal labeling system called BioID in Ewing's sarcoma cells. We identified AHDC1 as a proximal protein of EWS-ETS fusion proteins. AHDC1 knockdown showed a reduced cell growth and transcriptional activity of EWS-FLI1. AHDC1 knockdown also reduced BRD4 and BRG1 protein levels, both known as interacting proteins of EWS-FLI1. Our results suggest that AHDC1 supports cell growth through EWS-FLI1.

Intra-tumor heterogeneity and prognostic risk signature for hepatocellular carcinoma based on single-cell analysis

Intra-tumor heterogeneity poses a serious challenge in the treatment of cancer, including hepatocellular carcinoma (HCC). Recent developments in single-cell RNA sequencing (scRNA-seq) make it possible to examine the heterogeneity of tumor cells. The Gene Expression Omnibus (GEO) database was retrieved to obtain scRNA-seq data of 13 HCC and 8 para cancer samples, and the cells were clustered using FindNeighbors and FindClusters functions. Cell subsets were defined using the "Enricher" function of the clusterProfiler package. Monocle was used to predict cell developmental trajectory. The LIMMA package included in the R program was utilized to detect differentially expressed genes (DEGs) between HCC and paracancerous tissues. Univariate Cox analysis and Least Absolute and Selection Operator (Lasso) Cox regression analysis were conducted to establish a risk assessment model. Thirteen cell subpopulations were identified from the sequencing data of 64,634 single cells. Four cell subgroups (dendritic cells, hepatocytes, liver bud hepatic cells, and liver progenitor cells) in tumor tissues were highly enriched. Between HCC and para cancer tissues, 3024 DEGs were identified, and 641 were specific markers of four cell subgroups. To develop a prognostic risk model, 9 genes among the 641 genes were selected. In the training and validation sets, the model demonstrated a higher 5-year AUC and independently served as a prognostic marker as confirmed by multivariate and univariate Cox analyses. This study revealed the characteristics of different cell subpopulations of immune cells and tumor cells from the HCC microenvironment. We established a novel nine-gene prognostic model to determine the death risk of HCC patients. The discoveries in this research improved the current knowledge of HCC heterogeneity and may inspire future research.

EWS/FLI mediated reprogramming of 3D chromatin promotes an altered transcriptional state in Ewing sarcoma

Ewing sarcoma is a prototypical fusion transcription factor-associated pediatric cancer that expresses EWS/FLI or a highly related FET/ETS chimera. EWS/FLI dysregulates transcription to induce and maintain sarcomagenesis, but the mechanisms utilized are not fully understood. We therefore sought to define the global effects of EWS/FLI on chromatin conformation and transcription in Ewing sarcoma cells using a well-validated ‘knock-down/rescue’ model of EWS/FLI function in combination with next generation sequencing assays to evaluate how the chromatin landscape changes with loss, and recovery, of EWS/FLI expression. We found that EWS/FLI (and EWS/ERG) genomic localization is largely conserved across multiple patient-derived Ewing sarcoma cell lines. This EWS/FLI binding signature is associated with establishment of topologically-associated domain (TAD) boundaries, compartment activation, enhancer-promoter looping that involve both intra- and inter-TAD interactions, and gene activation. In addition, EWS/FLI co-localizes with the loop-extrusion factor cohesin to promote chromatin loops and TAD boundaries. Importantly, local chromatin features provide the basis for transcriptional heterogeneity in regulation of direct EWS/FLI target genes across different Ewing sarcoma cell lines. These data demonstrate a key role of EWS/FLI in mediating genome-wide changes in chromatin configuration and support the notion that fusion transcription factors serve as master regulators of three-dimensional reprogramming of chromatin.

Exploring prognostic value and regulation network of PPP1R1A in hepatocellular carcinoma

Novel and accurate biomarkers are needed for early detection and progression evaluation of hepatocellular carcinoma (HCC). Protein phosphatase 1 regulatory subunit 1A (PPP1R1A) has been studied in cancer biology; however, the expression pattern and biological function of PPP1R1A in HCC are unclear. The differentially expressed genes (DEGs) in HCC were screened by The Cancer Genome Atlas (TCGA) database. Real-time PCR and immunohistochemistry (IHC) assay were used to detect the expression of PPP1R1A in BALB/c mice, human normal tissues and corresponding tumor tissues, especially HCC. Then, Kaplan-Meier analysis of patients with HCC was performed to evaluate the relationship between PPP1R1A expression and prognosis. The transcriptional regulatory network of PPP1R1A was constructed based on the differentially expressed mRNAs, microRNAs and transcription factors (TFs). To explore the downstream regulation of PPP1R1A, the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis and immune infiltration score were performed. A total of 4 DEGs were screened out. PPP1R1A was differentially distributed and expressed in BALB/c mice and human tissues. PPP1R1A expression was higher in normal tissues than that in tumor tissues, and patients with higher PPP1R1A expression had better clinical outcome in HCC. In addition, we constructed miR-21-3p/TAL1/PPP1R1A transcriptional network. Furthermore, PPP1R1A may modulate the activation of PI3K-Akt pathway, cell cycle, glycogen metabolism and the recruitment of M2 macrophage in HCC. This study may help to clarify the function and mechanism of PPP1R1A in HCC and provide a potential biomarker for tumor prevention and treatment.

Whole genome analysis reveals the genomic complexity in metastatic cutaneous squamous cell carcinoma

Metastatic cutaneous squamous cell carcinoma (CSCC) is a highly morbid disease requiring radical surgery and adjuvant therapy, which is associated with a poor prognosis. Yet, compared to other advanced malignancies, relatively little is known of the genomic landscape of metastatic CSCC. We have previously reported the mutational signatures and mutational patterns of CCCTC-binding factor (CTCF) regions in metastatic CSCC. However, many other genomic components (indel signatures, non-coding drivers, and structural variants) of metastatic CSCC have not been reported. To this end, we performed whole genome sequencing on lymph node metastases and blood DNA from 25 CSCC patients with regional metastases of the head and neck. We designed a multifaceted computational analysis at the whole genome level to provide a more comprehensive perspective of the genomic landscape of metastatic CSCC. In the non-coding genome, 3′ untranslated region (3′UTR) regions of EVC (48% of specimens), PPP1R1A (48% of specimens), and ABCA4 (20% of specimens) along with the tumor-suppressing long non-coding RNA (lncRNA) LINC01003 (64% of specimens) were significantly functionally altered (Q-value < 0.05) and represent potential non-coding biomarkers of CSCC. Recurrent copy number loss in the tumor suppressor gene PTPRD was observed. Gene amplification was much less frequent, and few genes were recurrently amplified. Single nucleotide variants driver analyses from three tools confirmed TP53 and CDKN2A as recurrently mutated genes but also identified C9 as a potential novel driver in this disease. Furthermore, indel signature analysis highlighted the dominance of ID signature 13 (ID13) followed by ID8 and ID9. ID9 has previously been shown to have no association with skin melanoma, unlike ID13 and ID8, suggesting a novel pattern of indel variation in metastatic CSCC. The enrichment analysis of various genetically altered candidates shows enrichment of “TGF-beta regulation of extracellular matrix” and “cell cycle G1 to S check points.” These enriched terms are associated with genetic instability, cell proliferation, and migration as mechanisms of genomic drivers of metastatic CSCC.

One oncogene, several vulnerabilities: EWS/FLI targeted therapies for Ewing sarcoma

EWS/FLI is the defining mutation of Ewing sarcoma. This oncogene drives malignant transformation and progression and occurs in a genetic background characterized by few other recurrent cooperating mutations. In addition, the tumor is absolutely dependent on the continued expression of EWS/FLI to maintain the malignant phenotype. However, EWS/FLI is a transcription factor and therefore a challenging drug target. The difficulty of directly targeting EWS/FLI stems from unique features of this fusion protein as well as the network of interacting proteins required to execute the transcriptional program. This network includes interacting proteins as well as upstream and downstream effectors that together reprogram the epigenome and transcriptome. While the vast number of proteins involved in this process challenge the development of a highly specific inhibitors, they also yield numerous therapeutic opportunities. In this report, we will review how this vast EWS-FLI transcriptional network has been exploited over the last two decades to identify compounds that directly target EWS/FLI and/or associated vulnerabilities.

A Seven-Gene Signature to Predict Prognosis of Patients With Hepatocellular Carcinoma

Purpose: Hepatocellular carcinoma (HCC) is one of the most prevalent malignant diseases worldwide and has a poor prognosis. Gene-based prognostic models have been reported to predict the overall survival of patients with HCC. Unfortunately, most of the genes used in earlier prognostic models lack prospective validation and, thus, cannot be used in clinical practice.

Methods: Candidate genes were selected from GEPIA (Gene Expression Profiling Interactive Analysis), and their associations with patients’ survival were confirmed by RT-PCR using cDNA tissue microarrays established from patients with HCC after radical resection. A multivariate Cox proportion model was used to calculate the coefficient of corresponding gene. The expression of seven genes of interest (MKI67, AR, PLG, DNASE1L3, PTTG1, PPP1R1A, and TTR) with two reference genes was defined to calculate a risk score which determined groups of different risks.

Results: Our risk scoring efficiently classified patients (n = 129) with HCC into a low-, intermediate-, and high-risk group. The three groups showed meaningful distinction of 3-year overall survival rate, i.e., 88.9, 74.5, and 20.6% for the low-, intermediate-, and high-risk group, respectively. The prognostic prediction model of risk scores was subsequently verified using an independent prospective cohort (n = 77) and showed high accuracy.

Conclusion: Our seven-gene signature model performed excellent long-term prediction power and provided crucially guiding therapy for patients who are not a candidate for surgery.

Modulation of serine/threonine-protein phosphatase 1 (PP1) complexes: A promising approach in cancer treatment

Cancer is the second leading cause of death worldwide. Despite the availability of numerous therapeutic options, tumor heterogeneity and chemoresistance have limited the success of these treatments, and the development of effective anticancer therapies remains a major focus in oncology research. The serine/threonine-protein phosphatase 1 (PP1) and its complexes have been recognized as potential drug targets. Research on the modulation of PP1 complexes is currently at an early stage, but has immense potential. Chemically diverse compounds have been developed to disrupt or stabilize different PP1 complexes in various cancer types, with the objective of inhibiting disease progression. Beneficial results obtained in vitro now require further pre-clinical and clinical validation. In conclusion, the modulation of PP1 complexes seems to be a promising, albeit challenging, therapeutic strategy for cancer.

Comprehensive Genomic and Transcriptomic Analysis of Three Synchronous Primary Tumours and a Recurrence from a Head and Neck Cancer Patient

Synchronous primary malignancies occur in a small proportion of head and neck squamous cell carcinoma (HNSCC) patients. Here, we analysed three synchronous primaries and a recurrence from one patient by comparing the genomic and transcriptomic profiles among the tumour samples and determining the recurrence origin. We found remarkable levels of heterogeneity among the primary tumours, and through the patterns of shared mutations, we traced the origin of the recurrence. Interestingly, the patient carried germline variants that might have predisposed him to carcinogenesis, together with a history of alcohol and tobacco consumption. The mutational signature analysis confirmed the impact of alcohol exposure, with Signature 16 present in all tumour samples. Characterisation of immune cell infiltration highlighted an immunosuppressive environment in all samples, which exceeded the potential activity of T cells. Studies such as the one described here have important clinical value and contribute to personalised treatment decisions for patients with synchronous primaries and matched recurrences.

Therapeutic Potential of NTRK3 Inhibition in Desmoplastic Small Round Cell Tumor

PURPOSE

Desmoplastic small round cell tumor (DSRCT) is a highly lethal intra-abdominal sarcoma of adolescents and young adults. DSRCT harbors a t(11;22)(p13:q12) that generates the EWSR1-WT1 chimeric transcription factor, the key oncogenic driver of DSRCT. EWSR1-WT1 rewires global gene expression networks and activates aberrant expression of targets that together mediate oncogenesis. EWSR1-WT1 also activates a neural gene expression program.

EXPERIMENTAL DESIGN

Among these neural markers, we found prominent expression of neurotrophic tyrosine kinase receptor 3 (NTRK3), a druggable receptor tyrosine kinase. We investigated the regulation of NTRK3 by EWSR1-WT1 and its potential as a therapeutic target in vitro and in vivo, the latter using novel patient-derived models of DSRCT.

RESULTS

We found that EWSR1-WT1 binds upstream of NTRK3 and activates its transcription. NTRK3 mRNA is highly expressed in DSRCT compared with other major chimeric transcription factor-driven sarcomas and most DSRCTs are strongly immunoreactive for NTRK3 protein. Remarkably, expression of NTRK3 kinase domain mRNA in DSRCT is also higher than in cancers with NTRK3 fusions. Abrogation of NTRK3 expression by RNAi silencing reduces growth of DSRCT cells and pharmacologic targeting of NTRK3 with entrectinib is effective in both in vitro and in vivo models of DSRCT.

CONCLUSIONS

Our results indicate that EWSR1-WT1 directly activates NTRK3 expression in DSRCT cells, which are dependent on its expression and activity for growth. Pharmacologic inhibition of NTRK3 by entrectinib significantly reduces growth of DSRCT cells both in vitro and in vivo, providing a rationale for clinical evaluation of NTRK3 as a therapeutic target in DSRCT.

KDM5A and PHF2 positively control expression of pro-metastatic genes repressed by EWS/Fli1, and promote growth and metastatic properties in Ewing sarcoma

Ewing sarcoma is an aggressive malignant neoplasm with high propensity for metastasis and poor clinical outcomes. The EWS/Fli1 oncofusion is the disease driver in > 90% of cases, but presents a difficult therapeutic target. Moreover, EWS/Fli1 plays a complex role in disease progression, with inhibitory effects on critical steps of metastasis. Like many other pediatric cancers, Ewing sarcoma is a disease marked by epigenetic dysregulation. Epigenetic mechanisms present alternative targeting opportunities, but their contributions to Ewing sarcoma metastasis and disease progression remain poorly understood. Here, we show that the epigenetic regulators KDM5A and PHF2 promote growth and metastatic properties in Ewing sarcoma, and, strikingly, activate expression many pro-metastatic genes repressed by EWS/Fli1. These genes include L1CAM, which is associated with adverse outcomes in Ewing sarcoma, and promotes migratory and invasive properties. KDM5A and PHF2 retain their growth promoting effects in more metastatically potent EWS/Fli1^low cells, and PHF2 promotes both invasion and L1CAM expression in this cell population. Furthermore, KDM5A and PHF2 each contribute to the increased metastatic potency of EWS/Fli1^low cells in vivo. Together, these studies identify KDM5A and PHF2 as novel disease-promoting factors, and potential new targets, in Ewing sarcoma, including the more metastatically potent EWS/Fli1^low cell population.

Controlling Ser/Thr protein phosphatase PP1 activity and function through interaction with regulatory subunits

Protein phosphatase 1 is a major Ser/Thr protein phosphatase activity in eukaryotic cells. It is composed of a catalytic polypeptide (PP1C), with little substrate specificity, that interacts with a large variety of proteins of diverse structure (regulatory subunits). The diversity of holoenzymes that can be formed explain the multiplicity of cellular functions under the control of this phosphatase. In quite a few cases, regulatory subunits have an inhibitory role, downregulating the activity of the phosphatase. In this chapter we shall introduce PP1C and review the most relevant families of PP1C regulatory subunits, with particular emphasis in describing the structural basis for their interaction.

Merkel Cell Polyomavirus and Merkel Cell Carcinoma

Viruses are the cause of approximately 15% of all human cancers. Both RNA and DNA human tumor viruses have been identified, with Merkel cell polyomavirus being the most recent one to be linked to cancer. This virus is associated with about 80% of Merkel cell carcinomas, a rare, but aggressive cutaneous malignancy. Despite its name, the cells of origin of this tumor may not be Merkel cells. This review provides an update on the structure and life cycle, cell tropism and epidemiology of the virus and its oncogenic properties. Putative strategies to prevent viral infection or treat virus-positive Merkel cell carcinoma patients are discussed.

Protein phosphatase 1 regulatory subunit 1A regulates cell cycle progression in Ewing sarcoma

Ewing sarcoma (ES) is a malignant pediatric bone and soft tissue tumor. Patients with metastatic ES have a dismal outcome which has not been improved in decades. The major challenge in the treatment of metastatic ES is the lack of specific targets and rational combinatorial therapy. We recently found that protein phosphatase 1 regulatory subunit 1A (PPP1R1A) is specifically highly expressed in ES and promotes tumor growth and metastasis in ES. In the current investigation, we show that PPP1R1A regulates ES cell cycle progression in G1/S phase by down-regulating cell cycle inhibitors p21^Cip1 and p27^Kip1, which leads to retinoblastoma (Rb) protein hyperphosphorylation. In addition, we show that PPP1R1A promotes normal transcription of histone genes during cell cycle progression. Importantly, we demonstrate a synergistic/additive effect of the combinatorial therapy of PPP1R1A and insulin-like growth factor 1 receptor (IGF-1R) inhibition on decreasing ES cell proliferation and migration in vitro and limiting xenograft tumor growth and metastasis in vivo. Taken together, our findings suggest a role of PPP1R1A as an ES specific cell cycle modulator and that simultaneous targeting of PPP1R1A and IGF-1R pathways is a promising specific and effective strategy to treat both primary and metastatic ES.

Exosomes in Bone Sarcomas: Key Players in Metastasis

Bone sarcomas are rare cancers which often present with metastatic disease and are still associated with poor survival rates. Studies in the last decade have identified that exosomes, a type of extracellular vesicle released by cells, play an important role in tumour progression and dissemination. Through the transfer of their cargo (RNAs, proteins, and lipids) across cells, they are involved in cellular cross-talk and can induce changes in cellular behaviour. Exosomes have been shown to be important in metastasis organotropism, induction of angiogenesis and vascular permeability, the education of cells towards a pro-metastatic phenotype or the interaction between stromal and tumour cells. Due to the importance exosomes have in disease progression and the high incidence of metastasis in bone sarcomas, recent studies have evaluated the implications of these extracellular vesicles in bone sarcomas. In this review, we discuss the studies that evaluate the role of exosomes in osteosarcoma, Ewing sarcoma, and preliminary data on chondrosarcoma.

A Novel lncRNA HOXC-AS3 Acts as a miR-3922-5p Sponge to Promote Breast Cancer Metastasis

Purpose: The function of long noncoding RNAs (lncRNA) in breast cancer metastasis remains largely unknown. In this work, the role of HOXC-AS3 in breast cancer progression was investigated.Methods: By using Cancer Genome Atlas (TCGA) Database, we investigated the expression of HOXC-AS3 in breast cancer and explored the association between HOXC-AS3 expression and prognosis. Then, we studied the biological function of HOXC-AS3 in cell migration and invasion both in vitro and in vivo. Furthermore, the target miRNA of HOXC-AS3, and the target mRNA of miR-3922-5p were proved.Results: HOXC-AS3 is aberrantly overexpressed in breast cancers especially the HER2⁺ type. Moreover, high expression of HOXC-AS3 has a relationship with poor clinical outcomes of breast cancer. In addition, HOXC-AS3 regulates cell Invasion and migration both in vitro and in vivo. Our results demonstrated that miR-3922-5p was a direct target of HOXC-AS3, and PPP1R1A was a target of miR-3922-5p in breast cancer.Conclusions: The novel lncRNA HOXC-AS3 acts as a miR-3922-5p sponge to upregulate PPP1R1A protein expression, and thus results in promoting breast cancer metastasis. HOXC-AS3 could be a novel therapeutic target for breast cancer therapeutics.

The Role of C-X-C Chemokine Receptor Type 4 (CXCR4) in Cell Adherence and Spheroid Formation of Human Ewing's Sarcoma Cells under Simulated Microgravity

We studied the behavior of Ewing's Sarcoma cells of the line A673 under simulated microgravity (s-µg). These cells express two prominent markers-the oncogene EWS/FLI1 and the chemokine receptor CXCR4, which is used as a target of treatment in several types of cancer. The cells were exposed to s-µg in a random-positioning machine (RPM) for 24 h in the absence and presence of the CXCR4 inhibitor AMD3100. Then, their morphology and cytoskeleton were examined. The expression of selected mutually interacting genes was measured by qRT-PCR and protein accumulation was determined by western blotting. After 24 h incubation on the RPM, a splitting of the A673 cell population in adherent and spheroid cells was observed. Compared to 1 g control cells, EWS/FLI1 was significantly upregulated in the adherent cells and in the spheroids, while CXCR4 and CD44 expression were significantly enhanced in spheroids only. Transcription of CAV-1 was upregulated and DKK2 and VEGF-A were down-regulated in both, adherent in spheroid cells, respectively. Regarding, protein accumulation EWS/FLI1 was enhanced in adherent cells only, but CD44 decreased in spheroids and adherent cells. Inhibition of CXCR4 did not change spheroid count, or structure. Under s-µg, the tumor marker EWS/FLI1 is intensified, while targeting CXCR4, which influences adhesion proteins, did not affect spheroid formation.

Inhibition of protein phosphatase 1 stimulates noncanonical ER stress eIF2α activation to enhance fisetin-induced chemosensitivity in HDAC inhibitor-resistant hepatocellular carcinoma cells

Hepatocellular carcinoma (HCC) is a common fatal type of malignant tumor that has highly metastatic and recurrent properties. Fisetin is a natural flavonoid found in various vegetables and fruits which exhibits anti-cancer and anti-inflammatory properties, as well as other effects. Thus, we hypothesized that fisetin can act as an adjuvant therapy in cancer or drug-resistant cancer cells, and further investigated the molecular mechanisms underlying the development of drug-resistance in HCC cells. We found that fisetin effectively inhibited the cell viability of not only parental cells but also histone deacetylase inhibitors-resistant (HDACis-R) cells and enhanced the chemosensitivity of HCC cells. Interestingly, fisetin did not induce cell apoptosis through the activation of the endoplasmic reticulum (ER) stress sensor of protein kinase R (PKR)-like endoplasmic reticulum kinase, but rather through the non-canonical pathway of the protein phosphatase 1 (PP1)-mediated suppression of eIF2α phosphorylation. Moreover, fisetin-induced cell apoptosis was reversed by treatment with PP1 activator or eIF2α siRNA in HCC cells. Based on these observations, we suggest that PP1-eIF2α pathways are significantly involved in the effect of fisetin on HCC apoptosis. Thus, fisetin may act as a novel anticancer drug and new chemotherapy adjuvant which can improve the efficacy of chemotherapeutic agents and diminish their side-effects.

Genome-wide analysis of canine oral malignant melanoma metastasis-associated gene expression

Oral malignant melanoma (OMM) is the most common canine melanocytic neoplasm. Overlap between the somatic mutation profiles of canine OMM and human mucosal melanomas suggest a shared UV-independent molecular aetiology. In common with human mucosal melanomas, most canine OMM metastasise. There is no reliable means of predicting canine OMM metastasis, and systemic therapies for metastatic disease are largely palliative. Herein, we employed exon microarrays for comparative expression profiling of FFPE biopsies of 18 primary canine OMM that metastasised and 10 primary OMM that did not metastasise. Genes displaying metastasis-associated expression may be targets for anti-metastasis treatments, and biomarkers of OMM metastasis. Reduced expression of CXCL12 in the metastasising OMMs implies that the CXCR4/CXCL12 axis may be involved in OMM metastasis. Increased expression of APOBEC3A in the metastasising OMMs may indicate APOBEC3A-induced double-strand DNA breaks and pro-metastatic hypermutation. DNA double strand breakage triggers the DNA damage response network and two Fanconi anaemia DNA repair pathway members showed elevated expression in the metastasising OMMs. Cross-validation was employed to test a Linear Discriminant Analysis classifier based upon the RT-qPCR-measured expression levels of CXCL12, APOBEC3A and RPL29. Classification accuracies of 94% (metastasising OMMs) and 86% (non-metastasising OMMs) were estimated.

IPP-1 controls Akt/CREB phosphorylation extension in A2a adenosine receptor signaling cascade in MIN6 pancreatic β-cell line

Signaling through A_2a adenosine receptor specifically prevent pancreatic β-cells (PBCs) loses under diabetogenic conditions. However, signaling mediators of this receptor in PBCs remained unidentified. Thus, we aimed to investigate the possible involvement of PKA/Akt/IPP-1/CREB pathway in MIN6 β-cells. In addition, we investigated IPP-1 role in A_2a receptor signaling pathway. The expression of A_2a receptor in MIN6 cell line was evaluated by RT-PCR and its functionality confirmed by quantification of cAMP in response to the CGS 21680, an A_2a receptor agonist. MTT and Brdu assays were used to evaluate cell viability and proliferation, respectively. PKA activity and insulin release were evaluated using ELISA methods. P-Akt/Akt, p-IPP-1/IPP-1, and p-CREB/CREB levels were assessed using western blotting. IPP-1 knock down assessments was performed using specific siRNA. Our result revealed that MIN6 cells express A_2a receptor which actively increased cAMP levels (with EC₅₀ = 2.41 µM) and PKA activity. Activation of this receptor increased cell viability, proliferation and insulin release. Moreover, we mentioned A_2a receptor stimulation increased p-Akt, p-IPP-1, and p-CREB levels in dose (max at 10 µM of CGS 21680) and time (max at 30 min after CGS 21680 treatment) dependent manner. Interestingly, herein, we found in IPP-1 knocked down cells, A_2a receptor failed to activate Akt and CREB. Altogether, we mentioned that in MIN6 cells A_2a receptor increase cell viability, proliferation and insulin release through PKA/Akt/IPP-1/CREB signaling pathway. In addition, we conclude A_2a receptor signaling through this pathway is dependent to activation of IPP-1.

Adaptively capturing the heterogeneity of expression for cancer biomarker identification

Background

Identifying cancer biomarkers from transcriptomics data is of importance to cancer research. However, transcriptomics data are often complex and heterogeneous, which complicates the identification of cancer biomarkers in practice. Currently, the heterogeneity still remains a challenge for detecting subtle but consistent changes of gene expression in cancer cells.

Results

In this paper, we propose to adaptively capture the heterogeneity of expression across samples in a gene regulation space instead of in a gene expression space. Specifically, we transform gene expression profiles into gene regulation profiles and mathematically formulate gene regulation probabilities (GRPs)-based statistics for characterizing differential expression of genes between tumor and normal tissues. Finally, an unbiased estimator (aGRP) of GRPs is devised that can interrogate and adaptively capture the heterogeneity of gene expression. We also derived an asymptotical significance analysis procedure for the new statistic. Since no parameter needs to be preset, aGRP is easy and friendly to use for researchers without computer programming background. We evaluated the proposed method on both simulated data and real-world data and compared with previous methods. Experimental results demonstrated the superior performance of the proposed method in exploring the heterogeneity of expression for capturing subtle but consistent alterations of gene expression in cancer.

Conclusions

Expression heterogeneity largely influences the performance of cancer biomarker identification from transcriptomics data. Models are needed that efficiently deal with the expression heterogeneity. The proposed method can be a standalone tool due to its capacity of adaptively capturing the sample heterogeneity and the simplicity in use.

Software availability

The source code of aGRP can be downloaded from https://github.com/hqwang126/aGRP.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2437-2) contains supplementary material, which is available to authorized users.